Large-Eddy Simulation of Supercritical Mixing Layers
|
|
- Amelia Robinson
- 5 years ago
- Views:
Transcription
1 Sonderforschungsbereich/Transregio 4 Annual Report Large-Eddy Simulation of Supercritical Mixing Layers By C.A. Niedermeier, M. Jarczyk, S. Hickel, N.A. Adams AND M. Pfitzner Institute of Aerodynamics and Fluid Mechanics, Technische Universität München Boltzmannstr. 1, 8748 Garching bei München Within the development of a reliable numerical tool for the simulation of a whole rocket combustion chamber, real gas thermodynamics have been implemented into two CFD codes, the in-house code INCA of the Institute of Aerodynamics and Fluid Mechanics at the Technische Universität München and OpenFOAM, used by the Institute for Thermodynamics at the Universität der Bundeswehr München. The present work is part of the validation process, where Large Eddy Simulations (LES) with different subgridscale (SGS) models are conducted with these codes. Predictions of the implicit LES method ALDM and two explicit one equation models dynamic and non-dynamic are compared with results for Direct Numerical Simulations (DNS) of a non-reacting, transitional, temporal mixing layer of counterflowing oxygen and hydrogen. All models are able to reproduce the DNS data very well. A future task will be the validation of the SGS models for real gas jet flows and subsequent simulations of the mixing and combustion process of coaxial injectors. 1. Introduction The technology of cryogenic rocket combustion engines has been used successfully for many years and is well known today. Demand for increasing rocket performance and reliability is challenging in particular if restricted budgets and shortened development cycles have to be considered. For this reason, the application of computational methods in the development process increases steadily and therefore also the demand for computational fluid dynamics (CFD) tools that are able to simulate the flow at rocket combustor conditions. The pressure in modern combustion chambers often exceeds MPa, causing many propellants, like oxygen and hydrogen, to become supercritical at injection. In this high pressure environment, molecular interactions significantly affect the fluid properties. Hence, a real gas equation of state and suitable relations for the transport properties have to be used for the numerical simulation. We apply the cubic Peng-Robinson (PR) [1] equation of state (EOS) combined with an empirical volume correction method to obtain the thermodynamic behaviour of the real gas [2]. The ultimate goal of our work is to analyze the propellant injection into a rocket combustion chamber through multiple coaxial injectors and the subsequent combustion process. As this process is strongly three-dimensional and unsteady, Large Eddy Simulation (LES) is the most suitable method. Institute for Thermodynamics, Universität der Bundeswehr München, Werner-Heisenberg-Weg 39, 877 Neubiberg
2 136 C.A. Niedermeier, M. Jarczyk, S. Hickel, N.A. Adams & M. Pfitzner As a first step towards this goal, a non-reacting, transitional, temporal mixing layer of hydrogen and oxygen studied by Okong o et al. [3] is used in the present work as test case for validating the two CFD codes. The Institute for Thermodynamics of the Universität der Bundeswehr München uses the multi-component solver of the Open- FOAM package with an explicit one equation model (dynamic (DOEM) and non-dynamic (OEM)) as subgrid-scale (SGS) model, whereas the simulations carried out at the Institute of Aerodynamics and Fluid Mechanics at the Technische Universität München (TUM) are conducted by means of the in-house CFD code INCA using an implicit LES (ILES) method, namely the Adaptive Local Deconvolution Method (ALDM). This method gives very good results for compressible flows of ideal gases [4, ] and its usage is now extended to real gas flows. Additionally, Direct Numerical Simulations (DNS) are conducted with this code to have an own reference database for comparison with the LES results. The paper is organized as follows: In Sec. 2, the underlying conservation equations which are solved by our codes are briefly depicted. Secs. 3 and 4 describe the SGS and thermodynamic modeling approaches, while in Secs. and 6 the transport coefficients and the boundary and initial conditions of the simulations are given. The results consisting of flow visualizations of different variables at the transitional time are shown in Sec. 7. In Sec. 8, we draw conclusions from our results and give an outlook to future work. 2. Conservation Equations The general conservations equations for the H 2 /O 2 binary mixture under consideration are given by Okong o et al. [3] as ρe t t (ρu i ) t ρy O t ρ t + (ρu j) =, (2.1) + (ρu iu j +pδ ij ) + (ρy Ou j ) + [(ρe t +p)u j ] = τ ij, (2.2) = J Oj, (2.3) = q IKj + τ iju i. (2.4) Here x j, j {1,2,3} is a cartesian coordinate, t is the time, ρ is the density, u i is the velocity component in direction i {1,2,3}, e t = e + u i u i /2 is the total energy, p is the thermodynamic pressure and and Y O is the oxygen mass fraction. q IKj is the Irwing-Kirkwood (IK) form of the heat flux vector [6], J Oj is the oxygen mass flux in direction j and τ ij is the Newtonian viscous stress tensor. A detailed description of the flux calculation can be found in [3] and [7]. This flux formulation includes the Dufour and Soret effect as well as diffusion due to pressure gradients. In the investigated test cases, this very detailed diffusion modeling is needed in the DNS for the reproduction of the flow field to be as accurate as possible. The CFD codes follow two different approaches. INCA is a density-based finite volume method that directly solves Eqs. (2.1)-(2.4), whereas the OpenFOAM solver uses a pressure based PISO (Pressure Implicit with Splitting of Operators) algorithm proposed by Issa et al. [8, 9].
3 LES of Supercritical Mixing Layers 137 The PISO algorithm for compressible flows used in OpenFOAM is an extension of the incompressible approach. As in the incompressible limit the continuity equation becomes an auxiliary condition more than a conservation equation, it is no longer reasonable to solve the system of equations above. Therefore, pressure based algorithms have been developed, where a conservation equation for the pressure (developed from the discretized momentum and continuity equation) is solved instead of Eq. (2.1). A special issue here is the nonlinearity of the equation of state needed in high pressure environments, which has to be accounted for in the solution approach. Energy conservation in OpenFOAM is realized by solving the enthalpy equation (ρh) t + (ρhu j) = q IKj + Dp Dt with h = e t + p ρ (2.) instead of the total energy equation (2.4). The viscous heating term τ ij u i / is neglected here as it is insignificant in the subsonic flows studied in this work and known to negatively affect stability in pressure based solution approaches. 3. Turbulence Modeling and Numerics As an intermediate approach between Reynolds-averaged Navier-Stokes (RANS) and DNS, LES resolves the large scales of turbulence while the small scales below the grid width have to be modeled. In this paper, we apply two distinct modeling approaches: Implicit LES (ILES) with the Adaptive Local Deconvolution Method (ALDM) and explicit LES with a one-equation model in dynamic (DOEM) and non-dynamic form (OEM). The basic idea of ILES is to directly use the truncation error of the unmodified conservation equations (2.1)-(2.4). ALDM incorporates free parameters in the discretization scheme which can be used to control the truncation error. A physically motivated implicit SGS model that is consistent with turbulence theory is obtained through parameter calibration, see Ref. [4]. ALDM is implemented in INCA for Cartesian collocated grids and used to discretize the convective terms of the Navier-Stokes equations (see Ref. [] for a detailed description). The diffusive terms are discretized by 2 nd order centered differences and a 3 rd order explicit Runge-Kutta method is used for time integration. For the DNS conducted with INCA in the present work, the convective terms are discretized by 6 th order centered differences with all other properties of the computational framework being the same as for the ILES. For explicit LES, a compressible version of the OEM and DOEM is used in OpenFOAM []. An additional transport equation for the turbulent kinetic energy is solved in this model, which is then used to calculate the eddy viscosity. With OpenFOAM, 2 nd order centered differences are applied for spatial discretization and a 2 nd order backward differencing scheme is used for time integration. 4. Thermodynamic Modeling To account for real gas effects, all thermodynamic properties are calculated as the sum of an ideal reference value and a departure function that accounts for real gas effects. For closing the system of conservation equations, enthalpy h and constantpressure specific heat c p must be provided in OpenFOAM, while internal energy e and
4 138 C.A. Niedermeier, M. Jarczyk, S. Hickel, N.A. Adams & M. Pfitzner constant-volume specific heat c v must be provided in INCA. These are defined as p ( ( ) ) Vm h(t,p) = h (T)+ V m T dp, (4.1) T p p e(t,v) = e (T)+ V V ( p T T c p (T,p) = c v (T,p) c v (T,V) = ( e T ( ) ) p dv, (4.2) T V ( ) 2 p T ( ) V p V ) V T, (4.3). (4.4) Here the subscript refers to the ideal reference state at low pressure. In the present work, the departure functions on the right hand side are determined by the Peng- Robinson (PR) EOS [1] p = RT V b a(t) V 2 +2Vb b 2. (4.) V is the molar volume and R is the universal gas constant with a value of R = J/molK. The constants a(t) and b are calculated from empirical relations. a(t) accounts for attractive forces between the molecules in the fluid and is calculated from the empirical equation ( ( a(t) =.4723 R2 T 2 c p c 1 κ )) 2 T 1, (4.6) where κ = ω.26992ω 2 is a function of the acentric factor ω. The effects of the reduction of free volume by the particular volume of the molecules are taken into account via b =.77796RT c /p c. T c and p c are the critical temperatures and pressures of the modeled species (H 2 : 33.1 K / 1.3 MPa, O 2 : 14.6 K /.4 MPa). For multi-component mixtures, the extended Corresponding States Principle is applied, where a mixture of a fixed composition is assumed to behave like a pure fluid. The mixture properties are also calculated using the PR equation of state with parameters a and b calculated from real gas mixing rules a = x i x j a ij, (4.7) i j T c b = i x i b i. (4.8) Here a i and b i are the pure component parameters, and a ij is calculated from a ij = ai a j (1 k ij ) with the binary interaction coefficient k ij, which is set to zero for the investigations presented in this work. As the Peng-Robinson equation of state is known to be not very accurate in predicting the density in transcritical regions, an empirical correction method established by Harstad et al. [2] is chosen for the final implementation.
5 LES of Supercritical Mixing Layers 139 FIGURE 1. Mixing layer configuration: fluid 1 = oxygen, fluid 2 = hydrogen.. Transport Coefficients The transport coefficients are calculated with empirical correlations given by Okong o et al. [3] as µ = µ R (T/[.(T 1 +T 2 )]).7, (.1) Sc = µ/(ρα D D) = ( Y O.186YO 2.268YO 6 ) [ 1+(88.6/T) 1.], (.2) Pr = µc p /(mλ) = 1.33/T.1. (.3) The Reynolds number Re of the mixing layer is adjusted using the reference viscosity µ R, which can be calculated from Re =.(ρ 1 +ρ 2 ) U δ ω, µ R. (.4) The investigations presented in this work are carried out at Re = 7, which is sufficient for the mixing layer s transition [3]. ρ 1 and ρ 2 are the densities of the two different layers of the shear flow, U is the initial velocity difference across the layer ( U = U 1 U 2 ) and δ ω, is the initial vorticity thickness (δ ω, = 6.89 mm). 6. Initial and boundary conditions We chose the configuration of a temporal oxygen/hydrogen mixing layer investigated by Okong o et al. [3] as a validation test case. The domain is periodic in x 1 - and - direction and of outflow type in the -direction. The layer is not symmetric in as the layer growth is expected to be significantly larger on the hydrogen side. Table 1 summarizes the the mixing layer s mean flow properties. The free stream densities (ρ 1 and ρ 2 ) are calculated at the free stream temperatures (T 1 and T 2 ) and the initial uniform pressure p using the specified EOS. The initial mean profiles for streamwise velocity, temperature and species mass fractions are prescribed using an error function profile erf( π ). To initiate the layer roll up, vorticity distur-
6 14 C.A. Niedermeier, M. Jarczyk, S. Hickel, N.A. Adams & M. Pfitzner Mean quantity = (O 2) = (H 2) ū 1, m/s ρ, kg/m p, atm T, K 4 6 Y O 1 TABLE 1. Mixing layer mean flow properties. bances are superimposed on the mean velocity profile. Their wavelength corresponds to the most unstable incompressible mode [11]. These streamwise and spanwise disturbances are specified as ω 1 (, ) = F 3D λ 1 U Γ 1 f 2 ( )f 3 ( ), (6.1) λ 3 U ω 3 (x 1, ) = F 2D f 1 (x 1 )f 2 ( ). (6.2) Γ 3 The functions f 1, f 2 and f 3 are are given by 2 ( ) ( πx1 πx1 f 1 (x 1 ) = A i sin 2 i +A 3 sin π ), (6.3) λ 1 8λ 1 2 i= f 2 ( ) = exp [ λ 3 π ( x2 δ ω, ) 2 ], (6.4) ( ) ( ) 2πx3 πx3 f 3 ( ) = B sin +B 1 sin. (6.) The relative amplitudes of the vorticity disturbances are chosen asf 2D =.1 andf 3D =.. The coefficients A i and B i have the values A = 1, A 1 =., A 2 = A 3 =.3, B = 1 and B 1 =.2. The perturbation wave length is λ 1 = 7.29δ ω, and λ 3 =.6λ 1, corresponding to the domain size which is chosen to cover four wavelengths each in streamwise and spanwise direction. The circulations Γ 1 and Γ 3 are calculated by integratingf 2 f 3 andf 1 f 2 over the respective planes and dividing by twice the number of wavelengths in x 1 and once the number of wavelengths in, respectively. The initial velocity fluctuations can then finally be determined by solving the corresponding Poisson equations for the vorticities. Table 2 summarizes the domain sizes and grid resolutions which are used for the DNS and LES. L 3 7. Results In this section we compare the results of our reference DNS with ALDM, OEM and DOEM. Of course, an LES cannot recover all the features of a much higher resolved
7 LES of Supercritical Mixing Layers 141 Cases Domain size Grid resolution L 1 L 2 L 3,m 3 N 1 N 2 N 3 DNS LES TABLE 2. Domain size and grid resolutions for the H 2/O 2 mixing layer. DNS by definition. Therefore, the DNS results must be filtered and coarsened to the LES grid resolution for comparison with the LES results. We obtain the filtered DNS (FDNS) as reference solution by top-hat filtering over cubes of 4x4x4 cells. Figure 2 (left) shows the oxygen mass fraction in the - -plane at x 1 =.176 m. All SGS models are able to reproduce the mushroom-like structures in the O 2 mass fraction, only the fact that the leftmost structure already crosses the domain boundary in the FDNS is not covered to the full extent by the different LES. The tilt of the structures is best reproduced with ALDM and OEM, while the penetration of hydrogen into the oxygen side of the layer is a bit overestimated by all SGS models. The amount of unmixed oxygen remaining between the mushroom-like structures fits the FDNS very well in all LES, especially with ALDM, where these oxygen streaks are most prominent. The streamwise vorticity (Fig. 2 (right)) gives similar results. Shape, magnitude and position of the upper vortices obtained with all models are very close to each other and agree well with the FDNS results. While the minimum values are best captured with DOEM, ALDM gives the best result regarding the shape of the areas of positive streamwise vorticity. Only the fact that the positive peak values differ from each other in the FDNS is captured to the full extent by neither of the SGS models, but besides that the overall agreement is very good. The results for the non-dimensional spanwise vorticity (Fig. 3 (right)) in the x plane at =.6 m show comparable behaviour. The LES are able to reproduce many of the geometrical features of the FDNS contour as well as the overall position. The magnitude is generally in the right range as well. For example, ALDM is able to capture the maximum values in the middle and near the right boundary of the domain, only the area of maximum vorticity near the left boundary of the domain is overestimated. For OEM and DOEM the situation is contrariwise. While the vorticity magnitude near the left boundary matches the FDNS reference better, the maxima in the middle and near the right boundary are missing. That the results obtained with OpenFOAM do not match the shape of the high vorticity region in the lower part of the flow field is probably caused by boundary interactions.
8 142 C.A. Niedermeier, M. Jarczyk, S. Hickel, N.A. Adams & M. Pfitzner (a) (b) (c) (d) (e) (f) (g) (h) FIGURE 2. Y O for (a) FDNS, (c) ALDM, (e) OEM and (g) DOEM and ω 1δ ω,/ U for (b) FDNS, (d) ALDM, (f) OEM and (h) DOEM at t = 1 in the --plane at x 1 =.176 m.
9 LES of Supercritical Mixing Layers (a) (b) (c) x /δ 1 ω, (d) x /δ 1 ω, (e) x /δ 1 ω, (f) x /δ 1 ω, (g) x /δ 1 ω, (h) x /δ 1 ω, (i) x /δ 1 ω, (j) x /δ 1 ω, (k) x /δ 1 ω, (l) x x x 1 FIGURE 3. Y O for (a) FDNS, (d) ALDM, (g) OEM and (j) DOEM, ρ in kg/m 4 for (b) FDNS, (e) ALDM, (h) OEM and (k) DOEM and ω 3δ ω,/ U for (c) FDNS, (f) ALDM, (i) OEM and (l) DOEM at t = 1 in the x 1--plane at =.6 m.
10 144 C.A. Niedermeier, M. Jarczyk, S. Hickel, N.A. Adams & M. Pfitzner Regions with so-called high density gradient magnitude (HDGM) are of special interest among the geometrical features of the mixing layer. They redistribute turbulent energy from the normal to the tangential direction [12], therefore being crucial for the correct simulation of mixing and combustion. As can be seen in Fig. 3 (middle), again all models are able to reproduce the most important geometrical features of the FDNS solution s density gradient. The pressure based nature of the OpenFOAM approach (OEM and DOEM) obviously leads to higher density fluctuations, leading to more distorted results than with INCA (FDNS and ALDM). Besides that, the LES results are quite similar to each other with the strong density gradient in the lower left part of the domain being best captured with ALDM. Also, the area of the maximum density gradient in the middle is the smallest for ALDM, leading to the best match with the FDNS reference. The results for the O 2 mass fraction in the x 1 - -plane at =.6 m (Fig. 3 (left)) are similar to the previous findings. The general flow characteristics are met well with all SGS models. The substantial amount of hydrogen being present in the core of the vortex is well captured especially by ALDM and DOEM. Again, the results obtained with OpenFOAM look more distorted, but besides that both codes give comparable, good results. 8. Conclusions and outlook We simulated a non-reacting, transitional, temporal mixing layer of counterflowing hydrogen and oxygen by means of LES with different SGS models. The ILES method ALDM was used as well as the explicit SGS models OEM and DOEM. Additionally, a DNS of the configuration was simulated as a reference for comparison with the LES results. All SGS models generally show a good agreement with the FDNS results in the instantaneous flow visualizations. They are able to capture the most important geometrical flow features with their current parameterization. For a few distinct features, an advantage of ALDM and DOEM over OEM can be seen, but for a final assessment additional fully turbulent test cases will be used. There we expect the differences between the SGS models to become more prominent than in the transitional test case presented in this work. The TUM in-house code INCA showed its capability of performing high quality realgas simulations (DNS and LES). These yield the basis for the usage of ALDM for realgas flows and thereby for the very detailed investigation of supercritical flow phenomena. Additionally, OpenFOAM qualified as a promising industry-oriented tool for the final goal of simulating reacting flows in a rocket combustion chamber. Future work will focus on further validation of the SGS modeling for real-gas flows and later on the simulation of supercritical jet injection. The next step will be the investigation of fully turbulent oxygen/hydrogen and nitrogen/nitrogen mixing layers. Additionally, test cases like coaxial jet flow will be used to support this validation and optimization process. Acknowledgments Financial support has been provided by the German Research Foundation (Deutsche Forschungsgemeinschaft DFG) in the framework of the Sonderforschungsbereich Transregio 4. Computational resources have been provided by the Leibniz-Rechenzentrum München (LRZ).
11 References LES of Supercritical Mixing Layers 14 [1] PENG, D.-Y. AND ROBINSON, D.P. (1976). A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam., 1(1), [2] HARSTAD, K., MILLER, R. AND BELLAN, J. (1997). Efficient High-Pressure State Equations. AIChE Journal, 43(6), [3] OKONG O, N., HARSTAD, K. AND BELLAN, J. (22). Direct Numerical Simulations of O 2 /H 2 Temporal Mixing Layers Under Supercritical Conditions. AIAA Journal, 4(), [4] HICKEL, S. AND LARSSON, J. (28). An adaptive local deconvolution model for compressible turbulence. Proceedings of the CTR Summer Program 28, [] HICKEL, S. (211). Implicit subgrid-scale modeling for Large Eddy Simulation of compressible flows and shock turbulence interaction. Submitted to Phys. Fluids. [6] SARMAN, S. AND EVANS, D. (1992). Heat flux and mass diffusion in binary Lennard-Jones mixtures. Phys. Review A: General Physics, 4(4), [7] MILLER, R., HARSTAD, K. AND BELLAN, J. (21). Direct numerical simulations of supercritical fluid mixing layers applied to heptane-nitrogen. J. Fluid Mech., 436, [8] ISSA, R. (198). Solution of the Implicitly Discretized Fluid Flow Equations by Operator-Splitting. J. Comp. Phys., 62, 4 6. [9] ISSA, R., AHMADI-BEFRUI, B., BESHAY, K. AND GROSMAN, A. (1991). Solution of the Implicitly Discretized Reacting Flow Equations by Operator Splitting. J. Comp. Phys., 93, [] SCHUMANN, U. (197). Linear stability of finite difference equations for threedimensional flow problems. J. Comp. Phys., 18, [11] MOSER, R. AND ROGERS, M. (1991). Mixing transition and the cascade to small scales in a plane mixing layer. Phys. Fluids A, 3, [12] TASKINOGLU, E. AND BELLAN, J. (2). A posteriori study using a DNS database describing fluid disintegration and binary-species mixing under supercritical pressure: heptane and nitrogen. J. Fluid Mech., 64,
Large Eddy Simulation of Supercritical Nitrogen Jets
Sonderforschungsbereich/Transregio 40 Annual Report 2012 145 Large Eddy Simulation of Supercritical Nitrogen Jets By M. Jarczyk, M. Pfitzner, C. A. Niedermeier, S. Hickel AND N. A. Adams Institut für Thermodynamik,
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationCFD Simulation of high pressure real gas flows
CFD Simulation of high pressure real gas flows on the progress form art to physics :o) Maria Magdalena Poschner Prof. Dr. rer. nat. Pfitzner UNIVERSITÄT DER BUNDESWEHR MÜNCHEN Fakultät für Luft und Raumfahrt
More informationSteady Laminar Flamelet Modeling for turbulent non-premixed Combustion in LES and RANS Simulations
Sonderforschungsbereich/Transregio 4 Annual Report 213 139 Steady Laminar Flamelet Modeling for turbulent non-premixed Combustion in LES and RANS Simulations By H. Müller, C. A. Niedermeier, M. Pfitzner
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationModelling of subgrid-scale phenomena in supercritical transitional mixing layers: an aprioristudy
J. Fluid Mech. (27), vol. 593, pp. 57 91. c 27 Cambridge University Press doi:1.117/s221127875 Printed in the United Kingdom 57 Modelling of subgrid-scale phenomena in supercritical transitional mixing
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 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 informationAn evaluation of LES for jet noise prediction
Center for Turbulence Research Proceedings of the Summer Program 2002 5 An evaluation of LES for jet noise prediction By B. Rembold, J. B. Freund AND M. Wang Large-eddy simulation (LES) is an attractive
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationA parametrized non-equilibrium wall-model for large-eddy simulations
Center for Turbulence Research Proceedings of the Summer Program 212 127 A parametrized non-equilibrium wall-model for large-eddy simulations By S. Hickel, E. Touber, J. Bodart AND J. Larsson Wall-models
More informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationLarge-eddy simulation of cryogenic jet injection at supercritical pressures
Sonderforschungsbereich/Transregio 4 Annual Report 16 131 Large-eddy simulation of cryogenic jet injection at supercritical pressures By J. Matheis, H. Müller, S. Hickel AND M. Pfitzner Institute for Aerodynamics
More informationVelocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step
Copyright c 2004 Tech Science Press CMC, vol.1, no.3, pp.275-288, 2004 Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step B. Wang 1, H.Q. Zhang 1, C.K. Chan 2 and X.L.
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 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 informationOn the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows
Center for Turbulence Research Annual Research Briefs 1998 267 On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows By Jeffrey S. Baggett 1. Motivation and objectives
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 informationThe Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing
The Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM 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 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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationMixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin
Mixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin At very high pressures and densities, what is different and what is similar about the processes of Injection and Atomization,
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 informationWALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION
WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION Joongnyon Kim, Kyoungyoun Kim, Hyung Jin Sung Department of Mechanical Engineering, Korea Advanced Institute of Science
More informationTURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows
TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows Presented by William A. Sirignano Mechanical and Aerospace Engineering University of California
More informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationDNS, LES, and wall-modeled LES of separating flow over periodic hills
Center for Turbulence Research Proceedings of the Summer Program 4 47 DNS, LES, and wall-modeled LES of separating flow over periodic hills By P. Balakumar, G. I. Park AND B. Pierce Separating flow in
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 informationA NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS
June 30 - July 3, 2015 Melbourne, Australia 9 7B-4 A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS C.-Y. Chang, S. Jakirlić, B. Krumbein and C. Tropea Institute of Fluid Mechanics and Aerodynamics /
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationDirect numerical simulation of a turbulent reacting jet
Center for Turbulence Research Annual Research Briefs 999 59 Direct numerical simulation of a turbulent reacting jet By B. J. Boersma. Motivation and objectives Turbulent reacting jets are important in
More informationAdvanced Numerical Methods for non-premixed Flames
0 Advanced Numerical Methods for non-premixed Flames Annarita Viggiano Department of Environmental Engineering and Physics University of Basilicata Italy 1. Introduction Engine designers are under increasing
More informationPredicting natural transition using large eddy simulation
Center for Turbulence Research Annual Research Briefs 2011 97 Predicting natural transition using large eddy simulation By T. Sayadi AND P. Moin 1. Motivation and objectives Transition has a big impact
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
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 informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationLarge Eddy Simulation of Flame Flashback by Combustion Induced Vortex Breakdown
June 30 - July 3, 2015 Melbourne, Australia 9 1C-5 Large Eddy Simulation of Flame Flashback by Combustion Induced Vortex Breakdown Eike Tangermann Institut für Mathematik und Rechneranwendung Universität
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationDNS of Reacting H 2 /Air Laminar Vortex Rings
46th AIAA Aerospace Sciences Meeting and Exhibit 7-10 January 2008, Reno, Nevada AIAA 2008-508 DNS of Reacting H 2 /Air Laminar Vortex Rings Jeff Doom and Krishnan Mahesh University of Minnesota, Minneapolis,
More informationCOMPRESSIBLE TURBULENT CHANNEL AND PIPE FLOW: SIMILARITIES AND DIFFERENCES
COMPRESSIBLE TURBULENT CHANNEL AND PIPE FLOW: SIMILARITIES AND DIFFERENCES Somnath Ghosh, Rainer Friedrich Lehrstuhl für Aerodynamik, Technische Universität München Boltzmannstrasse 15, D-85748 Garching,
More informationPeriodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.
Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.
More informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More informationUncertainty quantification for RANS simulation of flow over a wavy wall
Uncertainty quantification for RANS simulation of flow over a wavy wall Catherine Gorlé 1,2,3, Riccardo Rossi 1,4, and Gianluca Iaccarino 1 1 Center for Turbulence Research, Stanford University, Stanford,
More informationJ. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and
J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationNumerical Investigation of the Transonic Base Flow of A Generic Rocket Configuration
1 Numerical Investigation of the Transonic Base Flow of A Generic Rocket Configuration A. Henze, C. Glatzer, M. Meinke, W. Schröder Institute of Aerodynamics, RWTH Aachen University, Germany March 21,
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationLarge Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena
29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied
More informationVerified Calculation of Nonlinear Dynamics of Viscous Detonation
Verified Calculation of Nonlinear Dynamics of Viscous Detonation Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and Joseph
More informationPhysics of turbulent flow
ECL-MOD 3A & MSc. Physics of turbulent flow Christophe Bailly Université de Lyon, Ecole Centrale de Lyon & LMFA - UMR CNRS 5509 http://acoustique.ec-lyon.fr Outline of the course A short introduction to
More informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationNumerical Investigation of Flow and Combustion in a Single-Element GCH4/GOX Rocket Combustor: Aspects of Turbulence Modeling
Numerical Investigation of Flow and Combustion in a Single-Element GCH4/GOX Rocket Combustor: Aspects of Turbulence Modeling A. Chemnitz, T. Sattelmayer, C. Roth, O. Haidn Technische Universität München,
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 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 informationA NEW MODELING OF SHEET CAVITATION CONSIDERING THE THERMODYNAMIC EFFECTS. Yumiko Sekino Ebara Cooporation, Ltd., Ota-ku, Tokyo JAPAN
Cav3-GS-6-3 Fifth International Symposium on Cavitation (CAV3) Osaka, Japan, November -4, 3 A NEW MODELING OF SHEET CAVITATION CONSIDERING THE THERMODYNAMIC EFFECTS Takashi Tokumasu Institute of Fluid
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
More informationL.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133, ORSAY CEDEX, FRANCE fax number:
Large Eddy Simulations of a spatially developing incompressible 3D mixing layer using the v-ω formulation. C. TENAUD, S. PELLERIN, A. DULIEU and L. TA PHUOC L.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133,
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationTurbulent Dynamics and Heat Transfer in Transcritical Channel Flow
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 Turbulent Dynamics and Heat Transfer in Transcritical Channel Flow Kukjin Kim School of Mechanical
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationCOMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL
COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL Kengo Asada 1 and Kozo Fujii 2 ABSTRACT The effects of body force distribution on the
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Ecole Centrale Paris] On: 27 January 29 Access details: Access Details: [subscription number 79383386] Publisher Taylor & Francis Informa Ltd Registered in England and
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 informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
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 informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
More information+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2):
International Conference on Methods of Aerophysical Research, ICMAR 008 LARGE EDDY SIMULATION OF TURBULENT ROUND IMPINGING JET B.B. Ilyushin, D.V. Krasinsky Kutateladze Institute of Thermophysics SB RAS
More informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationCOMMUTATION ERRORS IN PITM SIMULATION
COMMUTATION ERRORS IN PITM SIMULATION B. Chaouat ONERA, 93 Châtillon, France Bruno.Chaouat@onera.fr Introduction Large eddy simulation is a promising route. This approach has been largely developed in
More informationSurvey of Turbulent Combustion Models for Large Eddy Simulations of Propulsive Flowfields
AIAA SciTech 5-9 January 2015, Kissimmee, Florida 53rd AIAA Aerospace Sciences Meeting AIAA 2015-1379 Survey of Turbulent Combustion Models for Large Eddy Simulations of Propulsive Flowfields Downloaded
More informationHigh Speed Aerodynamics. Copyright 2009 Narayanan Komerath
Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
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 informationENHANCED HEAT TRANSFER IN LAMINAR PULSATING FLOW PAST A FLAT PLATE
HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta ENHANCED HEAT TRANSFER IN LAMINAR PULSATING FLOW PAST A FLAT PLATE Cárdenas Miranda, A.*
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 informationComputers and Mathematics with Applications. Investigation of the LES WALE turbulence model within the lattice Boltzmann framework
Computers and Mathematics with Applications 59 (2010) 2200 2214 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Investigation
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationarxiv: v1 [physics.flu-dyn] 13 Jul 2018
Thermodynamic analysis and large-eddy simulations of LOx-CH 4 and LOx-H 2 flames at high pressure Christoph Traxinger a,, Julian Zips a, Michael Pfitzner a a Institute for Thermodynamics, Bundeswehr University
More informationOn the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows
PHYSICS OF FLUIDS VOLUME 11, NUMBER 5 MAY 1999 On the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows L. Shao a) Laboratoire de Mécanique des Fluides
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationENGINEERING MECHANICS 2012 pp Svratka, Czech Republic, May 14 17, 2012 Paper #195
. 18 m 2012 th International Conference ENGINEERING MECHANICS 2012 pp. 309 315 Svratka, Czech Republic, May 14 17, 2012 Paper #195 NUMERICAL SIMULATION OF TRANSITIONAL FLOWS WITH LAMINAR KINETIC ENERGY
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay
More informationAtmospheric Boundary Layer Studies with Unified RANS-LES and Dynamic LES Methods
5st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 7 - January 23, Grapevine (Dallas/Ft. Worth Region), Texas AIAA 23-747 Atmospheric Boundary Layer Studies with
More informationThe Johns Hopkins Turbulence Databases (JHTDB)
The Johns Hopkins Turbulence Databases (JHTDB) HOMOGENEOUS BUOYANCY DRIVEN TURBULENCE DATA SET Data provenance: D. Livescu 1 Database Ingest and Web Services: C. Canada 1, K. Kalin 2, R. Burns 2 & IDIES
More informationNeeds work : define boundary conditions and fluxes before, change slides Useful definitions and conservation equations
Needs work : define boundary conditions and fluxes before, change slides 1-2-3 Useful definitions and conservation equations Turbulent Kinetic energy The fluxes are crucial to define our boundary conditions,
More informationIMPLEMENTATION OF POD AND DMD METHODS IN APACHE SPARK FRAMEWORK FOR SIMULATION OF UNSTEADY TURBULENT FLOW IN THE MODEL COMBUSTOR
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More information