EVALUATION OF EULER-EULER AND EULER-LAGRANGE STRATEGIES FOR LARGE-EDDY SIMULATIONS OF TURBULENT REACTING FLOWS
|
|
- Kerrie Phillips
- 5 years ago
- Views:
Transcription
1 EVALUATION OF EULER-EULER AND EULER-LAGRANGE STRATEGIES FOR LARGE-EDDY SIMULATIONS OF TURBULENT REACTING FLOWS M. García 1, Y. Sommerer 1, T. Schönfeld 1 and T. Poinsot 1 CERFACS - Toulouse, France IMFT - Toulouse, France Abstract Keywords This paper describes a study of Euler-Lagrange two-phase flow solver for Large- Eddy Simulations (LES). In modern LES for reacting two-phase flows, the coupling between the dispersed and the gaseous phase is not only a physical problem where all submodels for drop evaporation, combustion, wall interaction must be specified. It is also a numerical problem that, especially for strong coupling between the two phases, must be performed at each iteration of an unsteady computation and run efficiently on massively parallel machines. These issues are investigated here by implementing a droplet solver into a LES compressible code and assessing the results in terms of physics but also of numerical efficiency. The tests are performed for dispersed particles in homogeneous turbulence: results show that the code scales very well when the particles are homogeneously distributed but leads to large numerical losses if particles are located in one part only of the computational domain (as it is often the case for reacting flows). Results demonstrate the need for dynamic load balancing in such cases. Parallel computing, unstructured grids, two-phase flows, Lagrangian method, solid particles. INTRODUCTION Today, the RANS (Reynolds Averaged Navier Stokes) equations are routinely solved to design combustion chambers, for both gaseous and liquid fuels. Recently, in order to provide better accuracy for the prediction of mean flows but also to give access to unsteady phenomena occurring in combustion devices (such as instabilities, flashback or quenching), Large-Eddy Simulation (LES) has been extended to reacting flows. The success of these approaches for gaseous flames in the last years [1,, 3, 4, 5, 6, 7, 8, 9, 10, 11] is a clear illustration of their potential. LES gives access to the large scale structures of the flow, which reduces the importance of modelling and naturally captures a significant part of the physics controlling these flames. Even though LES has already demonstrated its potential for gaseous flames, its extension to two-phase flames is still largely to be done. There are many reasons which slow down the application of LES to two-phase reacting flows:
2 the physical submodels required to describe the atomization of a liquid fuel jet, the dispersion of droplets, their interaction with walls, evaporation and combustion are as difficult to build in LES as in RANS because they are essentially subgrid scale phenomena for which even the basic mechanisms are often not well understood. the numerical implementation of two-phase flow LES remains a challenge. The equations for both the gaseous and the dispersed phases must be solved together at each time step in a strongly coupled manner. This differs from classical RANS where the resolution of the two phases can be done in a weak procedure, bringing first the gas flow to convergence, then the droplets and finally iterating until convergence of both phases. This paper mainly focuses on the second issue: the efficiency of the coupled solvers for the two phases, especially in the context of parallel computers. Recent tests performed at CERFACS show that a modern LES code can run gaseous reacting flows with speedups of the order of 4900 on 5000 processors (see Events section). Maintaining a similar parallel efficiency for a two-phase flow solver raises questions which have never been addressed before. The first question is related, for such architectures, to the paradigm used to describe the two-phase flow: most RANS codes use Euler- Lagrange (EL) methods in which the flow is solved using an Eulerian method and the particles are tracked using a Lagrangian approach. An alternative technique is to use two-fluid models in which both the gas and the dispersed phases are solved using an Eulerian method (Euler-Euler or EE). The drawbacks and advantages of each method are discussed below: The EL approach is more intuitive because it highlights trajectories of particles (or droplets) which appeal to the CFD user. However, these particles are often not real particles but parcels which can contain thousands of droplets because present computers can not handle the millions of droplets which are created by standard fuel injectors every second. Therefore models are required for these parcels which sometimes makes modelling more difficult [1]. The EE approach requires an initial basic modelling effort which is larger than for the EL method [13]. When a two-fluid model is used, the EE model faces difficulties in handling droplet clouds with extended size distributions because it basically assumes that all droplets in a given cell move at the same speed.
3 In very dense zones where the topology of the flow might differ from a cloud of droplets, the EE approach is often preferred. For RANS computations, EE techniques are commonly used for fluidized beds [14, 15] or for chemical reactors [16, 17, 18]. In terms of computer implementation the EL approach is not well-suited to parallel computers: since two different solvers must be coupled, the complexity of the implementation on a parallel computer increases drastically compared to a single-phase code. Two methods may be used for LES: (1) task parallelization in which certain processors compute the gas flow and others compute the droplets flow and () domain partitioning in which droplets are computed together with the gas flow on geometrical subdomains mapped on parallel processors. Droplets must then be exchanged between processors when leaving a subdomain to enter an adjacent domain. For LES, it is easy to show that only domain partitioning is efficient on large grids because task parallelization would require the communication of very large three-dimensional data sets at each iteration between all processors. However, codes based on domain partitioning are difficult to optimize on massively parallel architectures when droplets are clustered in one part of the domain (typically, near the fuel injectors). Moreover, the distribution of droplets may change during the computation: for a gas turbine reignition sequence, for example, the chamber is filled with droplets when the ignition begins thus ensuring an almost uniform droplet distribution; these droplets then evaporate rapidly during the computation, leaving droplets only in the near injector regions. This leads to a poor speedup on a parallel machine if the domain is decomposed in the same way for the entire computation. As a result, dynamic load balancing strategies are required to redecompose the domain during the computation itself to preserve a high parallel efficiency [19]. On the other hand, EE techniques are naturally parallel because the flow and the droplets are solved using the same solver [0]. The history of RANS development has shown that both EE and EL are useful and either is found today in most commercial codes. Moreover, coupling strategies between EE and EL methods within the same application are considered for certain cases. In the present work, we focus on the EL technique for LES and investigate the key issue of parallel efficiency. This question has already been identified by other groups as a critical issue to run LES for reacting two-phase flows on massively parallel computers [19].
4 In the present paper, a Lagrangian solver is first coupled with a compressible LES code for reacting flows on hybrid grids (the AVBP code of CER- FACS: see Parallelization is performed using domain partitioning. The partitioning algorithm is the same as the one used for gaseous flows and does not take into account the additional computing cost induced by the droplets. The implementation is first verified by computing the dispersion of non evaporating droplets distributed homogeneously in homogeneous turbulence which is a classical test case [0, 1, ]. For this homogeneous case, the parallel efficiency is tested. More realistic cases are then tested where the droplets are distributed non-uniformly in the domain: this situation is more representative of real combustion chambers where droplets will be present mainly near the fuel injector. For this type of loading, adequate speedups are more difficult to obtain, which confirms the need for dynamic load balancing. NUMERICAL AND PHYSICAL MODELLING This section describes briefly the governing equations used in the AVBP solver for the gaseous and dispersed phases for the two-phase flow Euler-Euler [3] and Euler-Lagrangian model for dilute flows with one-way coupling between the two phases. The particle-tracking scheme and the interpolation algorithm are also included at the end of the section. Two-phase flow EE model Gaseous phase The continuity, momentum and energy equations for the gaseous phase are: (ρ g E g ) t (ρ g u g,i ) t ρ g t + (ρ gu g,j ) = 0 (1) + (ρ gu g,i u g,j ) = p g x i + τ g,ij () + (u g,j(ρ g E g + p g )) = (τ g,iju g,i ) q g,j. (3) Pressure is obtained from the equation of state p g = ρ g rt g where r is computed from the mass fractions r = R/W with W the air molecular weight. The Newtonian viscous tensor is defined as: ( ug,i τ g,ij = ρ g ν g + u g,j ) x i 3 δ u g,i ij (4) x i
5 and the heat diffusion is calculated from the classical Fourier law: q g,j = λ g T g. (5) Dispersed phase The continuity, momentum and energy equations for the dispersed phase are: (α p ρ p ) t (α p ρ p u p,i ) t (α p ρ p E QB ) t + (α pρ p u p,j ) = 0 (6) + (α pρ p u p,i u p,j ) = p QB x i + (α pτ QB,ij ) + I p,i (7) + (α pρ p u p,j E QB ) = α pρ p τ p E QB ( p QB δ ij + α p τ QB,ij ) u p,i + ( ) α p ρ p κ QB E QB (8) where α p = nπd 3 /6 with n the number of particles by unit volume, τ p is the particles relaxation time (to be defined in the following section) and κ QB = 5τ p E QB /3 where E QB δqp is the quasi-brownian energy which accounts for the uncorrelated particles velocity as described in [13]. The quasi-brownian pressure and stress tensor p QB and τ QB,ij are defined by analogy with the kinetic theory of gases. p QB is calculated through an equivalent equation of state: p QB = 3 α pρ p E QB (9) and τ QB,ij has a definition similar to the gaseous viscous tensor: ( up,i τ QB,ij = ρ p ν QB + u p,j ) x i 3 δ u p,i ij x i (10) where ν QB = τ p E QB /3 is the quasi-brownian viscosity and I p,i is the standard drag force (see [3] for details).
6 Two-phase flow EL model Gaseous phase The compressible Navier-Stokes equations are solved in their conservative form for the gaseous phase. The continuity, momentum and energy equations are the same as for the gaseous phase of the EE model (see Eqs. 1, and 3). Dispersed phase The dispersed phase consists of particles which are assumed to be rigid spheres. If the density of particles is much larger than the fluid density, the forces acting on particles reduce to drag and gravity [1] (though the latter is not considered in this study). Particle size is assumed small compared to the smallest turbulence length scale, the Kolmogorov scale, while collisions among particles are neglected. The analysis presented in this work is limited to one-way coupling from gas-to-particle phases with the assumption of small particle loadings. The particle equations of motion can then be written for a single particle as: x p,i t u p,i t = u p,i (11) = u p,i ũ g,i τ p (1) where ũ g,i is the fluid velocity at the position of the particle assuming that the flow field is locally undisturbed by the presence of this particle [4, 5]. Due to the very small droplet Reynolds number measured in the simulation, the particle relaxation time τ p is defined as the characteristic time for Stokes drag: τ p = d pρ p 18µ g. (13) Locating particles in elements of arbitrary shape In order to interpolate the discreted flow fields from the grid point onto the particle position (inside a cell) it is necessary to know the cell the particle is located in as well as the particle exact location within the cell. This is straightforward for regular, uniform Cartesian grids where the physical coordinates can be transformed into a uniform computational space. However,
7 this is no longer valid for unstructured grids. In order to decide if the particle is located inside the element or not, the scalar product is taken between the vector starting from the vertex of the element to the particle position and the inward normal vector of the corresponding edge. Then, the particle is inside the element if all the scalar products of each edge are positive. If any of them is negative, the particle falls outside the element in that direction. Search algorithm for particles on unstructured grids During each time step of the simulation, the particles change their position. But before they can contribute information to the grid or sample field information from it, their new host element (or cell) must be identified. Therefore, they must be traced to the grid and the idea is to exploit as an initial, most likely, guess that host element before the particle was moved. Under the reasonable condition that a particle can only travel to a neighboring cell during one single time step, it is sufficient to limit the search to the cell the particle was located in at the end of the previous timestep and to its adjacent cells. More precisely, the closest point of the mesh to the particle location is evaluated and only the elements surrounding to that point are considered. In case of search failure, the elements surrounding all the close points are considered. Such an approach is usually referred to as modified brute-force algorithm [6] and has been retained for this study. Interpolation of gaseous-phase properties For this particular case, hexahedral elements are used, so that, flow variables such as the velocity vector or the fluid viscosity, are evaluated at particle position via trilinear Lagrange interpolation from the values at the cell vertices. The AVBP solver The AVBP solver is a finite volume code based on a cell-vertex formulation. It solves the laminar and turbulent compressible Navier-Stokes equations in two and three space dimensions for hybrid and unstructured grids. Steady state or unsteady flows may be simulated, and the variations of molecular weights and heat capacities with temperature and mixture composition are accounted for. In this work, a Lax-Wendroff scheme was used for the numerical discretization in space and time.
8 APPLICATION TO HOMOGENEOUS ISOTROPIC TURBU- LENCE WITH PARTICLES The bulk of the simulations for this study was performed on a 64 3 uniform grid, the length of the computational domain is π 10 3 m in the three directions and has a cell size of x = m. Parallel simulations were performed on a Compaq AlphaServer SC with, 4, 8, and 1 processors and.48 million particles trajectories were computed at each time step corresponding approx. to 10 particles per cell. This number is a minimum value if proper sampling (specially for combustion) must be achieved in each cell. Initialization of the gaseous phase The gaseous phase is initialized with a divergence free velocity obeying a Passot-Pouquet spectrum [7] for the kinetic energy: ( ) k 4 E(k) = C e (k/ke). (14) ke Here k is the wave length and k e corresponds to the most energetic wave length. The gaseous Navier-Stokes equations are run up to t 0 = s which corresponds to half an eddy turnover time. This is necessary to ensure that the generated homogeneous isotropic turbulence (HIT) flow is indeed a solution of the Navier-Stokes equations before inserting the dispersed phase. Table 1. Flow field parameters at time t 0. N x/η u g [m/s] qg [m /s ] ε g [m /s 3 ] τ ε [s] Table 1 summarizes the main characteristics of the gaseous field obtained at time t 0. In this table, η is the Kolmogorov scale, u g is the fluctuating turbulent velocity, and qg and ε g are the gaseous kinetic energy and dissipation rate defined as: q g = 1 u g,iu g,i, ε g = ν g ug,i u g,i (15) respectively, and ν g is the kinematic viscosity defined by ν g = m /s. The last two parameters are used to form the Eulerian time macro-scale or turnover time, τ ε = q g/ε g. Averages of turbulence quantities obtained over the computational domain are denoted.
9 Figure 1. Snapshot (top) and D slice (bottom) of the initial particle distribution for the homogeneous (left) and non-homogeneous (right) cases. Initialization of the dispersed phase The dynamics of the particles in HIT depends on the Stokes number which is the ratio between the particle relaxation time and a characteristic time scale of turbulence. In the present study we define a Stokes number based on the turnover time, St ε = τp τ. For all simulations the particles diameter, d ε p, was set to m which gives τ p = s and St ε = This small value for St ε indicates that the particles initial velocity is close to the gaseous phase velocity so that it is a reasonable assumption to initialize the dispersed phase velocity field with the gaseous one. For EE simulations, this implies that the quasi-brownian energy is initialized with a value close to zero.
10 In both approaches, particles are randomly placed within the computational domain and periodic boundary conditions are used in all directions. Figure 1 shows the initial particle distribution for the two cases studied in this work: the homogeneous distribution, where all processors share the same number of particles, the non-homogeneous distribution, where only one processor contains all particles. These two cases correspond to limit situations: for the homogeneous loading case, high parallel efficiency should be obtained for the EL approach since all processors will have the same number of particles. On the other hand, the non-homogeneous case mimics a real combustor in which a few processors only located in the fuel injectors will contain most particles while others will have no particles at all. Comparison of dispersed two-phase flows Figure reports the temporal evolution of the kinetic energy of Lagrangian particles: q p Lagrangian = 1 N p u p,i u p,i (16) N p i=1 where N p is the total number of particles. This is compared with the particle kinetic energy obtained in the Eulerian simulation (note that prime in Eq. 15 is omitted here as in HIT u p,i = 0): q p Eulerian = 1 u p,i u p,i. (17) Note that qp Eulerian is obtained by adding the correlated part (mesoscopic energy) q p, and the uncorrelated part (or quasi-brownian energy) δqp [13]. The results show an overall good agreement between the two approaches and confirm that both methods provide the same result for this academic test case. As expected for this range of Stokes numbers [0], the particle kinetic energy qp obtained with both methods is larger than the gaseous phase kinetic energy qg.
11 Kinetic energy [m / s ] q p EE q p EL q g EE or EL δq p 0 3 t [s] 4 5 6x10-5 Figure. Time evolution of the dispersed phase kinetic energy q p, computed with the Euler- Euler model (squares), and with the Euler-Lagrange model (solid line). Gaseous phase kinetic energy q g (circles). Quasi-Brownian kinetic energy δq p (crosses). Analysis of code scalability In this section, the scalability of the Euler-Lagrange model is analyzed by means of two basic parameters used to measure the efficiency of parallel implementation: the speedup, S run (N procs ), and the reference single-phase CPU time ratio, T run (N procs ). The former is defined as the ratio between the CPU time of a simulation with processors and the CPU time of a simulation with a given number of processors N procs : S run (N procs ) = T run() T run (N procs ). (18) The latter is defined as the ratio between the CPU time of a simulation with a given number of procs and the CPU time of the reference single-phase simulation with processors: T run (N procs ) = T run(n procs ) T single phase (). (19)
12 The values of these two parameters are summarized in Tables and 3 for the homogeneous and non-homogeneous cases (note that in the EE model, the CPU time is the same for both cases since the computational cost of this approach does not depend on the number of particles). Table. Summary of the speedup of EE and EL models. N procs Ideal scaling Single-phase Two-phase EE Two-phase EL Homog Two-phase EL Non-Homog Table 3. Summary of the CPU time ratios of EE and EL models. N procs Single-phase Two-phase EE Two-phase EL Homog Two-phase EL Non-Homog Homogeneous particle loading For the homogeneous case Figure 3 exhibits an excelent speedup for both single and two-phase models. More precisely, the good parallel implementation of the two-phase EL model can be observed from the super-linear values obtained with more than 4 processors. Figure 3 (bottom) indicates that taking into account the dispersed phase results in about twice and a half the CPU time the single-phase simulation for the EE model and about twice the CPU time for the EL approach. The differences between both models increase with the number of processors, thus, the higher the number of processors, the faster the EL model, leading to differences from 18% for a two-processors simulation up to 55% for the simulation performed with 1 processors. This is due to a
13 lower ratio of particles by processor which implies less computational operations, less memory requirements and, therefore, less cache memory defaults. These effects are stronger than the increase of the communication cost related to a higher number of particles crossing processors domain. Although the current tests have been performed over only a few processors, the results reported in Figure 3 show encouraging scaling of the EL homogeneous case. The super-linear speedup obtained for the EL homogeneous case confirms the improved behaviour of the EL approach when increasing the number of partitions. Non-homogeneous particle loading In the non-homogeneous case, the dispersed phase in the EL model is only computed by a single processor (the master node in the present study). This leads to a significant increase in CPU time compared to the homogeneous case where all parallel processors execute the same number of instructions for both gaseous and dispersed phases. The main reasons contributing to this increase are: 1) a higher number of numerical operations per processor related to the solution of Eqs. 11 and 1, such as search and interpolation procedures; ) higher memory requirements and 3) communication cost related to the flux of particles between processors. On the other hand, the EE approach does not suffer from this load balancing problem. This is confirmed in Figure 4 which shows that the speedup of the non-homogeneous two-phase EL approach drifts rapidly far from the ideal, linear curve. The drop of performance in this case is not related to large communications costs between processors but merely to the parallel load imbalance. Thus, an efficient parallel implementation of an EL model can only be achieved by applying the same principles used in parallel computing so as to reduce the computational time required to complete a task: namely, decomposition of a data set into smaller subsets executing the same instructions on separate processors. This points out the need of dynamic load balancing for real two-phase flow LES in complex geometry (see also [19]), e.g., by using multi-constraint partitioning algorithms which take into account particles loading on each processor.
14 Ideal scaling Single phase Two-phase EE Two-phase EL Speedup, S run (N procs ) N procs 3 CPU time ratio, T run (N procs ) Single-phase Two-phase EE Two-phase EL N procs Figure 3. Homogeneous particle distribution. Top: Speedup of the single and two-phase flow models, Eq. 18. Bottom: CPU time ratio with respect to the single-phase simulation, Eq. 19 (all scalings are reported relative to the corresponding two-processors simulation).
15 Ideal scaling Single phase Two-phase EE Two-phase EL Speedup, S run (N procs ) N procs 3 CPU time ratio, T run (N procs ) Single-phase Two-phase EE Two-phase EL N procs Figure 4. Non-homogeneous particle distribution. Top: Speedup of the single and twophase flow models, Eq. 18. Bottom: CPU time ratio with respect to the single-phase simulation, Eq. 19 (all scalings are reported relative to the corresponding two-processors simulation).
16 CONCLUSIONS AND PERSPECTIVES In the present work Euler-Euler and Euler-Lagrange formulation for twophase reacting flows were compared in terms of parallel efficiency. Tests were performed for dispersed particles in homogeneous isotropic turbulence: results show that the code scales very well when the particles are homogeneously distributed but leads to large numerical losses if particles are located in only one partition of the computational domain. As a consequence, a strategy to control the load of particles before and during the simulations must be considered in order to avoid big differences in the number of particles between processors, e.g., with multi-constraint partitioning algorithms to take into account particles information. Future work includes the integration of particle/mesh load balancing capabilities and the development of additional modules for evaporation and two-way coupling. REFERENCES [1] D. Caraeni, C. Bergström, and L. Fuchs. Modeling of liquid fuel injection, evaporation and mixing in a gas turbine burner using large eddy simulation. Flow, Turb. and Combustion, 65:3 44, 000. [] V.K. Chakravarthy and S. Menon. Subgrid modeling of turbulent premixed flames in the flamelet regime. Flow, Turb. and Combustion, 65: , 000. [3] O. Colin, F. Ducros, D. Veynante, and T. Poinsot. A thickened flame model for large eddy simulations of turbulent premixed combustion. Phys. Fluids, 1: , 000. [4] H. Forkel and J. Janicka. Large-eddy simulation of a turbulent hydrogen diffusion flame. Flow, Turb. and Combustion, 65: , 000. [5] H. Pitsch and L. Duchamp de la Geneste. Large eddy simulation of premixed turbulent combustion using a level-set approach. Proc. of the Comb. Institute, 9, 00. [6] K. Mahesh, G. Constantinescu, and P. Moin. A numerical method for large-eddy simulation in complex geometries. J. Comput. Phys., 197:15 40, 004. [7] L. Selle, G. Lartigue, T. Poinsot, A. Kaufmann, W. Krebs, and D. Veynante. Compressible large-eddy simulation of turbulent combustion in complex geometry on unstructured meshes. Combust. Flame, 137: , 004. [8] Y. Sommerer, D. Galley, T. Poinsot, S. Ducruix, F. Lacas, and D. Veynante. Large eddy simulation and experimental study of flashback and blow-off in a lean partially premixed swirled burner. J. of Turbulence, 5, 004. [9] V. Moureau, G. Lartigue, Y. Sommerer, C. Angelberger, O. Colin, and T. Poinsot. Highorder methods for dns and les of compressible multi-component reacting flows on fixed and moving grids. J. Comput. Phys., 0: , 005.
17 [10] S. Roux, G. Lartigue, T. Poinsot, U. Meier, and C. Bérat. Studies of mean and unsteady flow in a swirled combustor using experiments, acoustic analysis and large eddy simulations. Combust. Flame, 141:40 54, 005. [11] T. Poinsot and D. Veynante. Theoretical and numerical combustion (Second Edition). R.T. Edwards, 005. [1] P. Moin. Large eddy simulation of multi-phase turbulent flows in realistic combustors. Prog. in Computational Fluid Dynamics, 4:37 40, 004. [13] P. Février, O. Simonin, and K.D. Squires. Partitioning of particle velocities in gas-solid turbulent flows into a continuous field and a spatially-uncorrelated random distribution: theoretical formalism and numerical study. J. Fluid Mech., submitted, 005. [14] H. Enwald, E. Peirano, and A.-E. Almstedt. Eulerian two-phase flow theory applied to fluidization. Int. J. Multiphase Flow, :1 66, [15] E. Peirano and B. Leckner. Fundamentals of turbulent gas-solid flows applied to circulating fluidized bed combustion. Prog. Energy Combust. Sci., 4:59 96, [16] J. Gao, C. Xu, S. Lin, and G. Yang. Simulations of gas-liquid-solid 3-phase flow and reaction in FCC riser reactors. AIChE Journal, 47(3):677 69, 001. [17] A. K. Das. Computational fluid dynamics simulation of gas-solid risers : reactive flow modelling. PhD thesis, Universiteit Gent, 00. [18] A. Gobin, H. Neau, O. Simonin, J.-R. Linas, V. Reiling, and J.-L. Sélo. Fluid dynamic numerical simulation of a gas phase polymerization reactor. International Journal for Numerical Methods in Fluids, 43: , 003. [19] F. Ham, A. Sourabh, G. Iaccarino, X. Wu, Herrmann. M., G. Constantinescu, K. Mahesh, and P. Moin. Unstructured les of reacting multiphase flows in realistic gas turbine combustors. In NASA AMES/Stanford University Center for Turbulence Research, editor, Annual Research Briefs, pages , 003. [0] A. Kaufmann, O. Simonin, and T. Poinsot. Direct Numerical Simulation of Particle-Laden Homogeneous Isotropic Turbulent Flows Using a Two-Fluid Model Formulation. In 5th International Conference on Multiphase Flow, 004. [1] M. Boivin, O. Simonin, and K.D. Squires. Direct numerical simulation of turbulence modulation by particles in isotropic turbulence. J. Fluid Mech., 375:35 63, [] P. Février. Etude numérique des effets de concentration préférentielle et de corrélation spatiale entre vitesses de particules solides en turbulence homogène isotrope stationnaire. PhD thesis, Institut National Polytechnique de Toulouse, 000. [3] A. Kaufmann. Towards Eulerian-Eulerian large eddy simulation of reactive two phase flow - TH/CFD/04/11. PhD thesis, Institut National Polytechnique de Toulouse, 004. [4] R. Gatignol. The faxén formulae for a rigid particle in an unsteady non-uniform stokes flow. J. Méch. Theor. Appl, 1:143, 1983.
18 [5] M. Maxey and J. Riley. Equation of motion for a small rigid sphere in a turbulent fluid flow. Phys. Fluids, 6:883, [6] S.V. Apte, K. Mahesh, P. Moin, and J.C. Oefelein. Subgrid modeling of turbulent premixed flames in the flamelet regime. Int. J. of Multiphase Flow, 9: , 003. [7] T. Passot and A. Pouquet. Numerical simulation of compressible homogeneous flow in the turbulent regime. J. Fluid Mech., 181: , 1987.
Evaluation of numerical strategies for LES of two-phase reacting flows
Center for Turbulence Research Proceedings of the Summer Program 2006 197 Evaluation of numerical strategies for LES of two-phase reacting flows By E. Riber,, M. García, V. Moureau, H. Pitsch, O. Simonin
More informationStudies of mean and unsteady flow in a swirled combustor using experiments, acoustic analysis and Large Eddy Simulations
Studies of mean and unsteady flow in a swirled combustor using experiments, acoustic analysis and Large Eddy Simulations S. Roux a,, G. Lartigue a, T. Poinsot a,b, U. Meier c and C. Bérat d a CERFACS,
More informationA Ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries
Center for Turbulence Research Annual Research Briefs 2005 269 A Ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries By V. Moureau, P. Minot, C. Bérat AND H. Pitsch
More informationDynamics of Lean Premixed Systems: Measurements for Large Eddy Simulation
Dynamics of Lean Premixed Systems: Measurements for Large Eddy Simulation D. Galley 1,2, A. Pubill Melsió 2, S. Ducruix 2, F. Lacas 2 and D. Veynante 2 Y. Sommerer 3 and T. Poinsot 3 1 SNECMA Moteurs,
More informationComparison between Lagrangian and Eulerian. mesoscopic modelling approaches for inertial. particles suspended in decaying isotropic.
Comparison between Lagrangian and Eulerian mesoscopic modelling approaches for inertial particles suspended in decaying isotropic turbulence A. Kaufmann a,1 M. Moreau b O. Simonin b J. Helie b,2 a CERFACS,
More informationLarge-eddy simulation of an industrial furnace with a cross-flow-jet combustion system
Center for Turbulence Research Annual Research Briefs 2007 231 Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system By L. Wang AND H. Pitsch 1. Motivation and objectives
More informationOverview of Turbulent Reacting Flows
Overview of Turbulent Reacting Flows Outline Various Applications Overview of available reacting flow models LES Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Ever-Expanding
More informationXXXVIII Meeting of the Italian Section of the Combustion Institute
Coupling a Helmholtz solver with a Distributed Flame Transfer Function (DFTF) to study combustion instability of a longitudinal combustor equipped with a full-scale burner D. Laera*, S.M. Camporeale* davide.laera@poliba.it
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal PI: Sarma L. Rani Department of Mechanical and Aerospace Engineering The
More informationProcess Chemistry Toolbox - Mixing
Process Chemistry Toolbox - Mixing Industrial diffusion flames are turbulent Laminar Turbulent 3 T s of combustion Time Temperature Turbulence Visualization of Laminar and Turbulent flow http://www.youtube.com/watch?v=kqqtob30jws
More informationExercises in Combustion Technology
Exercises in Combustion Technology Exercise 4: Turbulent Premixed Flames Turbulent Flow: Task 1: Estimation of Turbulence Quantities Borghi-Peters diagram for premixed combustion Task 2: Derivation of
More informationWall-Functions and Boundary Layer Response to Pulsating and Oscillating Turbulent Channel Flows
K. Hanjalić, Y. Nagano and S. Jakirlić (Editors) Wall-Functions and Boundary Layer Response to Pulsating and Oscillating Turbulent Channel Flows D. Panara 1, M. Porta 2,R. Dannecker 1, and B. Noll 1 1
More informationBest Practice Guidelines for Combustion Modeling. Raphael David A. Bacchi, ESSS
Best Practice Guidelines for Combustion Modeling Raphael David A. Bacchi, ESSS PRESENTATION TOPICS Introduction; Combustion Phenomenology; Combustion Modeling; Reaction Mechanism; Radiation; Case Studies;
More informationLES AND ACOUSTIC ANALYSIS OF COMBUSTION INSTABILITIES IN GASTURBINES
LES AND ACOUSTIC ANALYSIS OF COMBUSTION INSTABILITIES IN GASTURBINES T. Poinsot and L. Selle IMF and CERFACS, Toulouse, France Abstract Keywords RANS (Reynolds Averaged) techniques for combustion constitute
More informationCompressible Large Eddy Simulation of turbulent combustion in complex geometry on unstructured meshes
Compressible Large Eddy Simulation of turbulent combustion in complex geometry on unstructured meshes L. Selle a, G. Lartigue a, T. Poinsot b, R. Koch c, K.-U. Schildmacher c, W. Krebs d, P. Kaufmann d
More informationNew Developments in Large Eddy Simulation of Complex Flows
14th Australasian Fluid Mechanics Conference Adelaide University, Adelaide, Australia 1-14 December 21 New Developments in Large Eddy Simulation of Complex Flows P. Moin Center for Turbulence Research
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationA G-equation formulation for large-eddy simulation of premixed turbulent combustion
Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent
More informationFLAME WRINKLING FACTOR DYNAMIC MODELING FOR LARGE EDDY SIMULATIONS OF TURBULENT PREMIXED COMBUSTION
August 8 -, Poitiers, France FLAME WRINKLING FACTOR DYNAMIC MODELING FOR LARGE EDDY SIMULATIONS OF TURBULENT PREMIXED COMBUSTION Thomas Schmitt, Matthieu Boileau, Denis Veynante Laboratoire EMC CNRS -
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationFlame / wall interaction and maximum wall heat fluxes in diffusion burners
Flame / wall interaction and maximum wall heat fluxes in diffusion burners de Lataillade A. 1, Dabireau F. 1, Cuenot B. 1 and Poinsot T. 1 2 June 5, 2002 1 CERFACS 42 Avenue Coriolis 31057 TOULOUSE CEDEX
More informationA comparison between two different Flamelet reduced order manifolds for non-premixed turbulent flames
8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013 A comparison between two different Flamelet
More informationPairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray
Pairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray Georges Akiki, Kai Liu and S. Balachandar * Department of Mechanical &
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal and Vijaya Rani PI: Sarma L. Rani Department of Mechanical and Aerospace
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationSpectral analysis of energy transfer in variable density, radiatively heated particle-laden flows
Center for Turbulence Research Proceedings of the Summer Program 24 27 Spectral analysis of energy transfer in variable density, radiatively heated particle-laden flows By H. Pouransari, H. Kolla, J. H.
More informationMULTIPHASE FLOW MODELLING
MULTIPHASE FLOW MODELLING 1 Introduction 2 Outline Multiphase Flow Modeling Discrete phase model Eulerian model Mixture model Volume-of-fluid model Reacting Flow Modeling Eddy dissipation model Non-premixed,
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationChapter 7. Basic Turbulence
Chapter 7 Basic Turbulence The universe is a highly turbulent place, and we must understand turbulence if we want to understand a lot of what s going on. Interstellar turbulence causes the twinkling of
More informationDEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION
TMS (The Minerals, Metals & Materials Society), DEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION G.F. Yao, C. W. Hirt, and
More informationD. VEYNANTE. Introduction à la Combustion Turbulente. Dimanche 30 Mai 2010, 09h00 10h30
D. VEYNANTE Introduction à la Combustion Turbulente Dimanche 30 Mai 2010, 09h00 10h30 Introduction to turbulent combustion D. Veynante Laboratoire E.M2.C. CNRS - Ecole Centrale Paris Châtenay-Malabry France
More informationSensitivity analysis to the normal grid-resolution in a turbulent channel flow using large-eddy simulations
Center for Turbulence Research Annual Research Briefs 01 Sensitivity analysis to the normal grid-resolution in a turbulent channel flow using large-eddy simulations By J. Dombard AND G. Iaccarino 1. Motivation
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationLecture 10 Turbulent Combustion: The State of the Art
Lecture 10 Turbulent Combustion: The State of the Art 10.-1 Engineering applications are typically turbulent turbulence models These models use systematic mathematical derivations based on the Navier-
More informationNUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER
NUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER Naveen Kumar D 1*, Pradeep R 2 and Bhaktavatsala H R 3 1 Assistant Professor Department of Mechanical Engineering, M S Engineering College,
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationEvaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows
Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows by Ali Afshar A thesis submitted in conformity with the requirements for the degree of Masters
More informationTowards Large Eddy Simulation of Non-Homogeneous Particle Laden Turbulent Gas Flows Using Euler-Euler Approach
Towards Large Eddy Simulation of Non-Homogeneous Particle Laden Turbulent Gas Flows Using Euler-Euler Approach Eleonore Riber CERFACS Av. G. Coriolis 317 Toulouse riber@cerfacs.fr Mathieu Moreau IMFT Av.
More informationLarge-Eddy Simulation of Reacting Turbulent Flows in Complex Geometries
K. Mahesh Aerospace Engineering and Mechanics, University of Minnesota, 107 Akerman Hall, Minneapolis, MN 55455 G. Constantinescu Civil and Environmental Engineering, University of Iowa, Iowa City, IA
More informationThere are no simple turbulent flows
Turbulence 1 There are no simple turbulent flows Turbulent boundary layer: Instantaneous velocity field (snapshot) Ref: Prof. M. Gad-el-Hak, University of Notre Dame Prediction of turbulent flows standard
More informationA DEDICATED LES EXPERIMENTAL DATABASE FOR THE ASSESSMENT OF LES SGS MODELS: THE PULSATILE JET IMPINGEMENT IN TURBULENT CROSS FLOW
A DEDICATED LES EXPERIMENTAL DATABASE FOR THE ASSESSMENT OF LES SGS MODELS: THE PULSATILE JET IMPINGEMENT IN TURBULENT CROSS FLOW Hubert Baya Toda Energy Applications Techniques IFP Energie Nouvelles Rueil
More informationImpact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions
Journal of Physics: Conference Series PAPER OPEN ACCESS Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions To cite this article: A Rosiak and
More informationAnalysis of dynamic models for turbulent premixed combustion
Center for Turbulence Research Proceedings of the Summer Program 2012 387 Analysis of dynamic models for turbulent premixed combustion By D. Veynante, T. Schmitt, M. Boileau AND V. Moureau Very few attempts
More informationReacting Flow Modeling in STAR-CCM+ Rajesh Rawat
Reacting Flow Modeling in STAR-CCM+ Rajesh Rawat Latest Additions (v 7.02/v 7.04) Eulerian Multi-phase Reaction Model Soot Model Moment Methods PPDF Flamelet Multi-stream model Complex Chemistry Model
More informationANSYS Advanced Solutions for Gas Turbine Combustion. Gilles Eggenspieler 2011 ANSYS, Inc.
ANSYS Advanced Solutions for Gas Turbine Combustion Gilles Eggenspieler ANSYS, Inc. 1 Agenda Steady State: New and Existing Capabilities Reduced Order Combustion Models Finite-Rate Chemistry Models Chemistry
More informationSimulation of a lean direct injection combustor for the next high speed civil transport (HSCT) vehicle combustion systems
Center for Turbulence Research Annual Research Briefs 27 241 Simulation of a lean direct injection combustor for the next high speed civil transport (HSCT) vehicle combustion systems By H. El-Asrag, F.
More informationMSc. Thesis Project. Simulation of a Rotary Kiln. MSc. Cand.: Miguel A. Romero Advisor: Dr. Domenico Lahaye. Challenge the future
MSc. Thesis Project Simulation of a Rotary Kiln MSc. Cand.: Miguel A. Romero Advisor: Dr. Domenico Lahaye 1 Problem Description What is a Rotary Kiln? A Rotary Kiln is a pyroprocessing device used to raise
More informationNUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP REGIME
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00394 NUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP
More informationS. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives
Center for Turbulence Research Annual Research Briefs 2005 325 The dynamics of premixed flames propagating in non-uniform velocity fields: Assessment of the significance of intrinsic instabilities in turbulent
More informationTowards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows
Center for Turbulence Research Annual Research Briefs 009 199 Towards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows By M. Kostka, E.
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationThe effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets
ILASS Americas, 26 th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR, May 204 The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization
More informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More informationLES Approaches to Combustion
LES Approaches to combustion LES Approaches to combustion LES Approaches to Combustion W P Jones Department of Mechanical Engineering Imperial College London Exhibition Road London SW7 2AZ SIG on Combustion
More informationTransported PDF Calculations of Combustion in Compression- Ignition Engines
International Multidimensional Engine Modeling User s Group Meeting at the SAE Congress Detroit, MI 15 April 2013 Transported PDF Calculations of Combustion in Compression- Ignition Engines V. Raj Mohan
More informationA Priori Model for the Effective Lewis Numbers in Premixed Turbulent Flames
Paper # 070LT-0267 Topic: Turbulent Flames 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013.
More informationLES of the Sandia Flame D Using an FPV Combustion Model
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 82 (2015 ) 402 409 ATI 2015-70th Conference of the ATI Engineering Association LES of the Sandia Flame D Using an FPV Combustion
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationMulti-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah
Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary
More informationModeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet Model with Heat-Loss Effects
9 th U. S. National Combustion Meeting Organized by the Central States Section of the Combustion Institute May 17-20, 2015 Cincinnati, Ohio Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet
More informationNumerical Simulations of Hydrogen Auto-ignition in a Turbulent Co-flow of Heated Air with a Conditional Moment Closure
Numerical Simulations of Hydrogen Auto-ignition in a Turbulent Co-flow of Heated Air with a Conditional Moment Closure I. Stanković*, 1, A. Triantafyllidis, E. Mastorakos, C. Lacor 3 and B. Merci 1, 4
More informationModeling flame brush thickness in premixed turbulent combustion
Center for Turbulence Research Proceedings of the Summer Program 2006 299 Modeling flame brush thickness in premixed turbulent combustion By E. Knudsen, O. Kurenkov, S. Kim, M. Oberlack AND H. Pitsch Turbulent
More informationCoupling LES, radiation and structure in gas turbine simulations
Center for Turbulence Research Proceedings of the Summer Program 2010 239 Coupling LES, radiation and structure in gas turbine simulations By J. Amaya, E. Collado, B. Cuenot AND T. Poinsot Multi-physics
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationA validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved
Center for Turbulence Research Annual Research Briefs 2009 185 A validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved By E. Knudsen AND
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationComputational model for particle deposition in turbulent gas flows for CFD codes
Advanced Computational Methods and Experiments in Heat Transfer XI 135 Computational model for particle deposition in turbulent gas flows for CFD codes M. C. Paz, J. Porteiro, A. Eirís & E. Suárez CFD
More informationTurbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media
Copyright 011 Tech Science Press CMES, vol.71, no., pp.149-155, 011 Turbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media Zhi Yuan Cui 1, Jiu Hui Wu 1 and Di Chen Li 1 Abstract:
More informationLarge-eddy simulation of an evaporating and reacting spray
Center for Turbulence Research Annual Research Briefs 2008 479 Large-eddy simulation of an evaporating and reacting spray By T. Lederlin AND H. Pitsch 1. Motivation and objectives 1.1. Evaporative spray
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER
ADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER Ing. Vojtech Betak Ph.D. Aerospace Research and Test Establishment Department of Engines Prague, Czech Republic Abstract
More informationSimulation of Turbulent Lifted Flames and their Transient Propagation
25 th ICDERS August 2-7th, 2015 Leeds, UK Simulation of Turbulent Lifted Flames and their Transient Propagation S. Ruan, Z. Chen, N. Swaminathan University of Cambridge Cambridge, UK 1 Introduction Turbulent
More informationDispersed Multiphase Flow Modeling using Lagrange Particle Tracking Methods Dr. Markus Braun Ansys Germany GmbH
Dispersed Multiphase Flow Modeling using Lagrange Particle Tracking Methods Dr. Markus Braun Ansys Germany GmbH 2011 ANSYS, Inc., Markus Braun 1 Overview The Euler/Lagrange concept Breaking the barrier
More informationExperimental analysis and large eddy simulation to determine the response of non premixed flame submitted to acoustic forcing
Experimental analysis and large eddy simulation to determine the response of non premixed flame submitted to acoustic forcing B. Varoquié, J.P. Légier, F. Lacas, D. Veynante and T. Poinsot Laboratoire
More informationNumerical and physical instabilities in massively parallel LES of reacting flows
Numerical and physical instabilities in massively parallel LES of reacting flows Dr Thierry Poinsot 1 Institut de Mecanique des Fluides de Toulouse CNRS UMR 5502, FRANCE poinsot@imft.fr Summary. LES of
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationProspects for High-Speed Flow Simulations
Prospects for High-Speed Flow Simulations Graham V. Candler Aerospace Engineering & Mechanics University of Minnesota Support from AFOSR and ASDR&E Future Directions in CFD Research: A Modeling & Simulation
More informationDEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
DRAFT Proceedings of ASME IMECE: International Mechanical Engineering Conference & Exposition Chicago, Illinois Nov. 5-10, 2006 IMECE2006-14867 DEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
More informationDroplet Behavior within an LPP Ambiance
Copyright 2010 Tech Science Press FDMP, vol.6, no.4, pp.399-408, 2010 Droplet Behavior within an LPP Ambiance M. Chrigui 1,2, L. Schneider 1, A. Zghal 2, A. Sadiki 1 and J. Janicka 1 Abstract: This paper
More informationTable of Contents. Preface... xiii
Preface... xiii PART I. ELEMENTS IN FLUID MECHANICS... 1 Chapter 1. Local Equations of Fluid Mechanics... 3 1.1. Forces, stress tensor, and pressure... 4 1.2. Navier Stokes equations in Cartesian coordinates...
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationCoupling tabulated chemistry with large-eddy simulation of turbulent reactive flows
Center for Turbulence Research Proceedings of the Summer Program 2008 237 Coupling tabulated chemistry with large-eddy simulation of turbulent reactive flows By R. Vicquelin, B. Fiorina, N. Darabiha, D.
More informationTURBULENCE MODULATION IN LARGE EDDY SIMULATION OF BACKWARD-FACING STEP FLOW LADEN WITH PARTICLES
Engineering MECHANICS, Vol. 20, 2013, No. 3/4, p. 299 307 299 TURBULENCE MODULATION IN LARGE EDDY SIMULATION OF BACKWARD-FACING STEP FLOW LADEN WITH PARTICLES Jaroslav Volavý*, Miroslav Jícha* This work
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationSimulation of atomization : from DNS to industrial applications
Simulation of atomization : from DNS to industrial applications MUSAF III - 29/09/2016 D. Zuzio, J.-L. Estivalèzes, O. Rouzaud, P. Gajan, P. Villedieu PhD/postdoc : G. Blanchard, I. Marter, A. Orazzo,
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationTopics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films
Lecture Topics Transient Droplet Vaporization Convective Vaporization Liquid Circulation Transcritical Thermodynamics Droplet Drag and Motion Spray Computations Turbulence Effects Topics in Other Lectures
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationEcoulements Diphasiques
Ecoulements Diphasiques Marc Massot Professor at Ecole Centrale Paris Laboratoire EM2C - UPR CNRS 288 - Ecole Centrale Paris Visiting Professor, Department of Mechanical Engineering, Center for Turbulence
More informationApplication of Chimera Grids in Rotational Flow
CES Seminar Write-up Application of Chimera Grids in Rotational Flow Marc Schwalbach 292414 marc.schwalbach@rwth-aachen.de Supervisors: Dr. Anil Nemili & Emre Özkaya, M.Sc. MATHCCES RWTH Aachen University
More informationEffects of the preferential segregation of droplets on evaporation and turbulent mixing
Effects of the preferential segregation of droplets on evaporation and turbulent mixing Julien Reveillon, François-Xavier Demoulin To cite this version: Julien Reveillon, François-Xavier Demoulin. Effects
More informationJoint use of Compressible Large-Eddy. Simulation and Helmholtz solvers for the. analysis of rotating modes in an industrial.
Joint use of Compressible Large-Eddy Simulation and Helmholtz solvers for the analysis of rotating modes in an industrial swirled burner L. Selle a, L. Benoit a, T. Poinsot b, F. Nicoud c and W. Krebs
More informationThermoacoustic Instabilities Research
Chapter 3 Thermoacoustic Instabilities Research In this chapter, relevant literature survey of thermoacoustic instabilities research is included. An introduction to the phenomena of thermoacoustic instability
More informationBest Practice Guidelines for Computational Turbulent Dispersed Multiphase Flows. René V.A. Oliemans
Best Practice Guidelines for Computational Turbulent Dispersed Multiphase Flows René V.A. Oliemans ERCOFTAC Seminar, Innventia, Stockholm, June 7-8, 2011 1 Vermelding onderdeel organisatie Department of
More informationReynolds number scaling of inertial particle statistics in turbulent channel flows
Reynolds number scaling of inertial particle statistics in turbulent channel flows Matteo Bernardini Dipartimento di Ingegneria Meccanica e Aerospaziale Università di Roma La Sapienza Paolo Orlandi s 70th
More information