CFD modeling of single-phase flow in a packed bed with MRI validation

Computational fluid dynamics (CFD) simulations are becoming increasingly widespread with the advent of more powerful computers and more sophisticated software. The aim of these developments is to facilitate more accurate reactor design and optimization methods compared to traditional lumped-parameter models. However, in order for CFD to be a trusted method, it must be validated using experimental data acquired at sufficiently high spatial resolution. This article validates an in-house CFD code by comparison with flow-field data obtained using magnetic resonance imaging (MRI) for a packed bed with a particle-to-column diameter ratio of 2. Flows characterized by inlet Reynolds numbers, based on particle diameter, of 27, 55, 111, and 216 are considered. The code used employs preconditioning to directly solve for pressure in low-velocity flow regimes. Excellent agreement was found between the MRI and CFD data with relative error between the experimentally determined and numerically predicted flow-fields being in the range of 3-9%. © 2012 American Institute of Chemical Engineers (AIChE).

The design of multi-element airfoils through multi-objective optimization techniques

This paper presents the development and the application of a multi-objective optimization framework for the design of two-dimensional multi-element high-lift airfoils. An innovative and efficient optimization algorithm, namely Multi-Objective Tabu Search (MOTS), has been selected as core of the framework. The flow-field around the multi-element configuration is simulated using the commercial computational fluid dynamics (cfd) suite Ansys cfx. Elements shape and deployment settings have been considered as design variables in the optimization of the Garteur A310 airfoil, as presented here. A validation and verification process of the cfd simulation for the Garteur airfoil is performed using available wind tunnel data. Two design examples are presented in this study: a single-point optimization aiming at concurrently increasing the lift and drag performance of the test case at a fixed angle of attack and a multi-point optimization. The latter aims at introducing operational robustness and off-design performance into the design process. Finally, the performance of the MOTS algorithm is assessed by comparison with the leading NSGA-II (Non-dominated Sorting Genetic Algorithm) optimization strategy. An equivalent framework developed by the authors within the industrial sponsor environment is used for the comparison. To eliminate cfd solver dependencies three optimum solutions from the Pareto optimal set have been cross-validated. As a result of this study MOTS has been demonstrated to be an efficient and effective algorithm for aerodynamic optimizations. Copyright © 2012 Tech Science Press.

The Effect of Inlet Guide Vanes on Inlet Flow Distortion Transfer and Transonic Fan Stability

An investigation was carried out into the effects of variable inlet guide vanes (VIGVs) on the performance and stability margin of a transonic fan in the presence of inlet flow distortion. The study was carried out using computational fluid dynamics (CFD) and validated with experimental data. The capability of CFD to predict the changes in performance with or without VIGVs in the presence of an inlet flow distortion is assessed. Results show that the VIGVs improve the performance and stability margin and do so by reducing the amount of swirl at inlet to the rotor component of the fan.

Initial phase modelling in numerical explosion applied to process safety

The utilisation of computational fluid dynamics (CFD) in process safety has increased significantly in recent years. The modelling of accidental explosion via CFD has in many cases replaced the classical Multi Energy and Brake Strehlow methods. The benefits obtained with CFD modelling can be diminished if proper modelling of the initial phase of explosion is neglected. In the early stages of an explosion, the flame propagates in a quasi-laminar regime. Proper modelling of the initial laminar phase is a key aspect in order to predict the peak pressure and the time to peak pressure. The present work suggests a modelling approach for the initial laminar phase in explosion scenarios. Findings are compared with experimental data for two classical explosion test cases which resemble the common features in chemical process areas (confinement and congestion). A detailed analysis of the threshold for the transition from laminar to turbulent regime is also carried out. The modelling is implemented in a fully 3D Navier-Stokes compressible formulation. Combustion is treated using a laminar flamelet approach based on the Bray, Moss and Libby (BML) formulation. A novel modified porosity approach developed for the unstructured solver is also considered. Results agree satisfactorily with experiments and the modelling is found to be robust. © 2013 The Institution of Chemical Engineers.

Analysis of high pressure premixed flames using Equivalent Reactor Networks for predicting NOx emissions

This paper describes a computational study of lean premixed high pressure methane-air flames, using Computational Fluid Dynamics (CFD) together with a reactor network approach. A detailed chemical reaction mechanism is employed to predict pollutant concentrations, placing emphasis on nitrogen oxide emissions. The reacting flow field is divided into separate zones in which homogeneity of the physical and chemical conditions prevails. The defined zones are interconnected forming an Equivalent Reactor Network (ERN). Three flames are examined for which experimental data is available. Flame A is characterised by an equivalence ratio of 0.43 while Flames B and C are richer with equivalence ratios of 0.5 and 0.56 respectively. Computations are performed for a range of operating conditions, quantifying the effect in the emitted NOx levels. Model predictions are compared against the available experimental data. Sensitivity analysis is performed to investigate the effect of the network size, in order to define the optimum number of reactors for accurate predictions of the species mass fractions. © 2012 Elsevier Ltd. All rights reserved.

Off-wall boundary conditions for turbulent flows obtained from buffer-layer minimal flow units

There is strong evidence that the transport processes in the buffer region of wall-bounded turbulence are common across various flow configurations, even in the embryonic turbulence in transition (Park et al., Phys. Fl. 24). We use this premise to develop off-wall boundary conditions for turbulent simulations. Boundary conditions are constructed from DNS databases using periodic minimal flow units and reduced order modeling. The DNS data was taken from a channel at Reτ=400 and a zero-pressure gradient transitional boundary layer (Sayadi et al., submitted to J. Fluid Mech.). Both types of boundary conditions were first tested on a DNS of the core of the channel flow with the aim of extending their application to LES and to spatially evolving flows.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

Application of scalar dissipation rate model to siemens DLE combustors

The standard design process for the Siemens Industrial Turbomachinery, Lincoln, Dry Low Emissions combustion systems has adopted the Eddy Dissipation Model with Finite Rate Chemistry for reacting computational fluid dynamics simulations. The major drawbacks of this model have been the over-prediction of temperature and lack of species data limiting the applicability of the model. A novel combustion model referred to as the Scalar Dissipation Rate Model has been developed recently based on a flamelet type assumption. Previous attempts to adopt the flamelet philosophy with alternative closure models have failed, with the prediction of unphysical phenomenon. The Scalar Dissipation Rate Model (SDRM) was developed from a physical understanding of scalar dissipation rate, signifying the rate of mixing of hot and cold fluids at scales relevant to sustain combustion, in flames and was validated using direct numerical simulations data and experimental measurements. This paper reports on the first industrial application of the SDRM to SITL DLE combustion system. Previous applications have considered ideally premixed laboratory scale flames. The industrial application differs significantly in the complexity of the geometry, unmixedness and operating pressures. The model was implemented into ANSYS-CFX using their inbuilt command language. Simulations were run transiently using Scale Adaptive Simulation turbulence model, which switches between Large Eddy Simulation and Unsteady Reynolds Averaged Navier Stokes using a blending function. The model was validated in a research SITL DLE combustion system prior to being applied to the actual industrial geometry at real operating conditions. This system consists of the SGT-100 burner with a glass square-sectioned combustor allowing for detailed diagnostics. This paper shows the successful validation of the SDRM against time averaged temperature and velocity within measurement errors. The successful validation allowed application of the SDRM to the SGT-100 twin shaft at the relevant full load conditions. Limited validation data was available due to the complexity of measurement in the real geometry. Comparison of surface temperatures and combustor exit temperature profiles showed an improvement compared to EDM/FRC model. Furthermore, no unphysical phenomena were predicted. This paper presents the successful application of the SDRM to the industrial combustion system. The model shows a marked improvement in the prediction of temperature over the EDM/FRC model previously used. This is of significant importance in the future applications of combustion CFD for understanding of hardware mechanical integrity, combustion emissions and dynamics of the flame. Copyright © 2012 by ASME.

An approach to multi-fidelity optimization of aeroengine compression systems

The aerodynamic design of turbomachinery presents the design optimisation community with a number of exquisite challenges. Chief among these are the size of the design space and the extent of discontinuity therein. This discontinuity can serve to limit the full exploitation of high-fidelity computational fluid dynamics (CFD): such codes require detailed geometric information often available only sometime after the basic configuration of the machine has been set by other means. The premise of this paper is that it should be possible to produce higher performing designs in less time by exploiting multi-fidelity techniques to effectively harness CFD earlier in the design process, specifically by facilitating its participation in configuration selection. The adopted strategy of local multi-fidelity correction, generated on demand, combined with a global search algorithm via an adaptive trust region is first tested on a modest, smooth external aerodynamic problem. Speed-up of an order of magnitude is demonstrated, comparable to established techniques applied to smooth problems. A number of enhancements aimed principally at effectively evaluating a wide range of configurations quickly is then applied to the basic strategy, and the emerging technique is tested on a generic aeroengine core compression system. A similar order of magnitude speed-up is achieved on this relatively large and highly discontinuous problem. A five-fold increase in the number of configurations assessed with CFD is observed. As the technique places constraints neither on the underlying physical modelling of the constituent analysis codes nor on first-order agreement between those codes, it has potential applicability to a range of multidisciplinary design challenges. © 2012 by Jerome Jarrett and Tiziano Ghisu.

Computational engineering design for micro-scale combustors

At present, optimisation is an enabling technology in innovation. Multi-objective and multi-disciplinary design tools are essential in the engineering design process, and have been applied successfully in aerospace and turbomachinery applications extensively. These approaches give insight into the design space and identify the trade-offs between the competing performance measures satisfying a number of constraints at the same time. It is anticipated here that the same benefits can be obtained for the design of micro-scale combustors. In this paper, a multi-disciplinary automated design optimisation system was developed for this purpose, which comprises a commercial computational fluid dynamics package and a multi-objective variant of the Tabu Search optimisation algorithm. The main objectives that are considered in this study are to optimise the main micro-scale combustor design characteristics and to satisfy manufacturability considerations from the very beginning of the whole design operation. Hydrogen-air combustion as well as 14 geometrical and 2 operational parameters are used to describe and model the design problem. Two illustrative test cases will be presented, in which the most important device operational requirements are optimised, and the efficiency of the developed optimisation system is demonstrated. The identification, assessment and suitability of the optimum design configurations are discussed in detail. Copyright © 2012 by ASME.

Scaling of transitionally-rough flow structures in DNSes of channels with riblets up to Reτ=550

The linear, drag-reducing effect of vanishingly small riblets breaks down once their size is in the transitionally-rough regime. We have previously reported that this breakdown is caused by the additional Reynolds stresses produced by the appearance of elongated spanwise rollers just above the riblet surface. These rollers are related with the Kelvin--Helmholtz instability of free shear layers, and to similar structures appearing over other rough and porous surfaces. However, because of the limited Reτ=180 in our previous DNSes, it could not be determined whether those structures scaled in inner or outer units. Furthermore, it is questionable if results in the transitionally-rough regime at Reτ=180 can be extrapolated to configurations of practical interest. At such small Reynolds numbers, roughness of transitional size can perturb a large portion of the boundary layer, which is not the case in most industrial and atmospheric applications. To clarify these issues we have conducted a set of DNSes at Reτ=550. Our results indicate that the spanwise rollers scale in wall units, and support the validity of the extrapolation to configurations of practical interest.

Kelvin-Helmholtz-like instability of turbulent flows over riblets

The turbulent drag reduction due to riblets is a function of their size and, for different configurations, collapses well with a length scale l+g=(A+g)1/2, based in the groove cross-section Ag. The initially linear drag reduction breaks down for l+g≈11, which agrees in our DNS with the previously reported appearance of quasi-two-dimensional spanwise rollers immediately above the riblets. They are similar to those found over porous surfaces and plant canopies, and can be traced to a Kelvin-Helmholtz-like instability associated with the relaxation of the impermeability condition for the wall-normal velocity. The extra Reynolds stress associated with them accounts quantitatively for the drag degradation. An inviscid model for the instability confirms its nature, agreeing well with the observed perturbation wavelengths and shapes. The onset of the instability is determined by a length scale L+w that, for conventional riblet geometries, is proportional to l+g. The instability onset, L+w≥4, corresponds to the empirical breakdown point l+g≈11.

The breakdown of the viscous regime in riblets

We investigate the mechanisms involved in the breakdown of the viscous regime in riblets, with a view to determining the point of optimum performance, where drag reduction ceases to be proportional to the riblet size. This occurs empirically for a groove cross-section $A_g^+ \approx 120^+$. To study the interaction of the riblets with the overlaying turbulent flow, we systematically conduct DNSes in a ribbed turbulent channel with increasing riblet size. The conditionally averaged crossflow above and within the grooves reveals a mean recirculation bubble that exists up to the point of viscous breakdown, isolating the groove floor from the overlying crossflow, and preventing the high momentum fluid from entering the grooves. We do not find evidence of outside vortices lodging within the grooves until $A_g^+ \approx 400$, which is well past the drag minimum, and already into the drag increasing regime. Interestingly, as the bubble breaks down, we observe that quasi-two-dimensional spanwise structures form just above the riblets, similar to those observed above porous surfaces and plant canopies, which appear to be involved in the performance degradation.

Effects of the biodiesel fuel physical properties on the swirl stabilised spray combustion characteristics

An increasin g interest in biofuel applications in modern engines requires a better understanding of biodiesel combustion behaviour. Many numerical studies have been carried out on unsteady combustion of biodiesel in situations similar to diesel engines, but very few studies have been done on the steady combustion of biodiesel in situations similar to a gas turbine combustor environment. The study of biodiesel spray combustion in gas turbine applications is of special interest due to the possible use of biodiesel in the power generation and aviation industries. In modelling spray combustion, an accurate representation of the physical properties of the fuel is a first important step, since spray formation is largely influenced by fuel properties such as viscosity, density, surface tension and vapour pressure. In the present work, a calculated biodiesel properties database based on the measured composition of Fatty Acid Methyl Esters (FAME) has been implemented in a multi-dimensional Computational Fluid Dynamics (CFD) spray simulation code. Simulations of non-reacting and reacting atmospheric-pressure sprays of both diesel and biodiesel have been carried out using a spray burner configuration for which experimental data is available. A pre-defined droplet size probability density function (pdf) has been implemented together with droplet dynamics based on phase Doppler anemometry (PDA) measurements in the near-nozzle region. The gas phase boundary condition for the reacting spray cases is similar to that of the experiment which employs a plain air-blast atomiser and a straight-vane axial swirler for flame stabilisation. A reaction mechanism for heptane has been used to represent the chemistry for both diesel and biodiesel. Simulated flame heights, spray characteristics and gas phase velocities have been found to compare well with the experimental results. In the reacting spray cases, biodiesel shows a smaller mean droplet size compared to that of diesel at a constant fuel mass flow rate. A lack of sensitivity towards different fuel properties has been observed based on the non-reacting spray simulations, which indicates a need for improved models of secondary breakup. By comparing the results of the non-reacting and reacting spray simulations, an improvement in the complexity of the physical modelling is achieved which is necessary in the understanding of the complex physical processes involved in spray combustion simulation. Copyright © 2012 SAE International.

On the two classes of global primary modal instability in laminar separation bubbles

The self-excited global instability mechanisms existing in flat-plate laminar separation bubbles are studied here, in order to shed light on the causes of unsteadiness and three- dimensionality of unforced, nominally two-dimensional separated flows. The presence of two known linear global mechanisms, namely an oscillator behavior driven by local regions of absolute inflectional instability and a centrifugal instability giving rise to a steady three- dimensionalization of the bubble, is studied in a series of model separation bubbles. Present results indicate that absolute instability, and consequently a global oscillator behavior, does not exist for two-dimensional bubbles with a peak reversed-flow velocity below 12% of the free-stream velocity. However, the three-dimensional instability becomes active for recirculation levels as low as urev ≈ 7%. These findings suggest a route to the three-dimensionality and unsteadiness observed in experiments and simulations substantially different from that usually found in the literature, in which two-dimensional vortex shedding is followed by three-dimensionalization.