Eulerian Monte Carlo Method for solving joint velocity-scalar PDF: numerical aspects and validation
|
|
- Erick Nicholson
- 5 years ago
- Views:
Transcription
1 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 1/63 Eulerian Monte Carlo Method for solving joint velocity-scalar PDF: numerical aspects and validation Olivier Soulard and Vladimir Sabel nikov ONERA, Palaiseau, France CEA, Bruyères-le-Châtel, France Workshop and Minicourse Conceptual Aspects of Turbulence, Wolfgang Pauli Institute, Vienna, February 2008
2 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 2/63 Contents PART I 1 Introduction 2 Velocity-scalar PDF Modeled velocity-scalar PDF Lagrangian and Eulerian Monte Carlo methods 3 Derivation of SPDEs for solving velocity-scalar PDF 4 Implications for scalar PDFs 5 Simulation of a turbulent premixed methane flame over a backward facing step 6 Conclusions for part I
3 Velocity-scalar PDF Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 3/63
4 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 4/63 Reactive Navier-Stokes equations Navier-Stokes equations for the velocity U, the density ρ, and a turbulent reactive scalar c : ρ t + U i t + U j c t + U j (ρu j ) = 0 U i = 1 ρ c = 1 ρ P x i + ν 2 U i J j + S with P : pressure, ν: cinematic diffusion coefficient, J j : scalar diffusion term and S: chemical source terms. 1 ρ Low Mach number assumption: ρ = ρ(c) and S = S(c, ρ(c)) = S(c)
5 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 5/63 Density weighted statistics For variable density flows, it is usual to work with density-weighted statistics: p Uc (U, c) = δ(c(x, t) c)δ(u(x, t) U) f Uc (U, c) = ρ U, c ρ Favre averages are noted Q : p Uc (U, c) = ρ(c) ρ p Uc(U, c) (Low Mach assumption) Q = ρq ρ
6 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 6/63 Modeled velocity-composition PDF equation Modeled transport equation for the Favre PDF f Uc : ( ρ t f ) Uc + ( ) ρ U j fuc = ( [ ρ 1 ] ) P G jk (U k Ũ k ) f Uc + 1 U j ρ 2 ρ C 0 ɛ 2 fuc U j U j + ( ρ ω c (c c) c f ) Uc ( ρ S(c) c f ) Uc with: First line: transport in physical space, treated exactly Second line: pressure fluctuations and molecular diffusion modelled by the Generalized Langevin model. ɛ is the turbulent dissipation Third line: Scalar mixing is modeled by the IEM model. Chemistry is treated exactly Equation (1) is a Fokker-Planck equation (1)
7 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 7/63 Solving PDF equation (1) PDF equation (1) has N d = N + 7 dimensions where N is the number of scalars = large in practical applications Finite volume/element/difference methods: CPU cost increases exponentially with N d Not suitable Monte Carlo methods: CPU cost increases linearly with N d OK Two options: Lagrangian (Particle) Monte Carlo (LMC) methods Eulerian Monte Carlo (EMC) methods
8 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 8/63 Lagrangian Monte Carlo methods (1/2) PDF is represented by a set of N p stochastic particles Each particle is a sample of the physical properties of the system F Np = N p X n=1 w (n) δ(c c (n) )δ(u U (n) )δ(x x (n) ) ; ρ e f U c = F Np w (n), c (n), U (n), x (n) : mass, concentration, velocity, position of particle (n) Particles evolve according to Ito stochastics ODEs (SODEs) (Pope, 1985): dc (n) = ω c c (n) ec dt + S(c (n) )dt (2) du (n) j = 1 ρ P dt G jk (U (n) k e U k )dt + p C 0 ɛ dw (n) j (t) (3) dx (n) j = U (n) j dt (4)
9 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 9/63 Lagrangian Monte Carlo methods (2/2) W (n) j are independent Gaussian white (Brownian) noises : = δ nm δ ij dt dw (n) j dw (m) k The absence of symbol in the stochastic product (multiplicative noise) in eq. (3) denotes the Ito interpretation (no-correlation property) Numerous publications document the convergence and accuracy of LMC methods LMC used in many complex calculations (including LES) LMC implemented in commercial CFD codes
10 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 10/63 Eulerian Monte Carlo methods Principle PDF represented by stochastic Eulerian fields Each field is a sample of the physical properties of the system Each field evolves according to a stochastic PDE (SPDE) Development of EMC methods is useful and stimulating LMC/EMC competition could push both approaches forward EMC methods already designed for joint scalar PDFs Theoretical works : Valiño (98), Sabel nikov & Soulard (05) Calculations : Naud et al (04), Sabel nikov & Soulard (06), Mustata et al. (06)
11 Derivation of SPDEs for solving velocity-scalar PDFs Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 11/63
12 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 12/63 Objectives and Method The question is : How can we derive SPDEs for a stochastic velocity field U a stochastic scalar field θ so that their statistics yield the Favre PDF? To answer this question : we transpose existing LMC models to a Eulerian framework 3 steps 1. We introduce the stochastic density (necessary to convert Lagrangian statistics into Eulerian ones) 2. We use the notion of stochastic characteristics to convert Lagrangian equations (SODEs) into Eulerian ones (SPDEs) 3. We show that Eulerian PDF (of U and θ) weighted by the stochastic density is identical to Favre PDF f Uc :
13 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 13/63 Stochastic density (1/2) The following Ito SODEs are used in LMC methods (Pope (1985)) ( dθ = ω c θ θ ) dt + S(θ)dt (5) Important: du j = 1 ρ dx j = U j dt P dt G jk (U k Ũk)dt + C 0 ɛ dw j (t) (6) SODEs (5)-(7) not only define a stochastic velocity U and a stochastic scalar θ, but also a stochastic density, noted r In a Lagrangian framework, the continuity equation is given by : (7) r = r 0 j ; j = Det [j ik ], with j ik = x i x 0k j is the Jacobian and r 0 is the initial stochastic density j defines the transformation from intial position x 0 to current one x
14 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 14/63 Stochastic density (2/2) The stochastic density r is different from the physical density ρ The evolution of r is given by: dr = r div (U) dt (8) For instance, in the case of a constant density ρ = const, but r const because div(u) 0: d dt (div(u)) = U i U j x i + x i (G ij (U j U j )) P ρ x i x i C 0 ɛ ɛ x i Ẇ i (9) This property is a consequence of the closure assumption used in the Langevin model, which only imposes a continuity constraint on the mean velocity field, but not on its instantaneous value
15 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 15/63 Lagrangian - Eulerian transposition SODEs (5) - (8) are the stochastic characteristics of the following hyperbolic Ito SPDEs : r t + U r j = r (U j ) (10) θ t + U j U i t + U j θ = ω c (θ θ) + S(θ) (11) U i = 1 ρ P x i G ij (U j Ũj) + C 0 ɛ Ẇi (12) These hyperbolic equations can be rewritten also in conservative form: r t + t (rθ) + t (ru i) + (ru j ) = 0 (ru j θ) = r ω c (θ θ) + rs(θ) (14) (ru j U i ) = r ρ (13) P x i + rg ij (U j U j ) + r C 0 ɛ Ẇ j (15)
16 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 16/63 Eulerian PDF The Eulerian PDF of U and θ is defined by : f Uθ (U, Θ; x, t) = δ(u U(x, t))δ(θ θ(x, t)) (16) r varies we introduce statistics weighted by the stochastic density r f Uθ:r (U, Θ; x, t) = r U, Θ f UΘ (U, Θ; x, t) (17) r By standard techniques, one can show that weighted PDF is identical to Favre PDF f Uc : f Uθ:r = f Uc provided that r = ρ
17 Implications for scalar PDF methods Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 17/63
18 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 18/63 Implications for scalar PDF methods (1/3) EMC method has been devised for the full velocity-composition PDF How do the previous results apply to the computation of the following scalar PDF equation? ( ( ρ t f ) c + ( ρ x Ũj f ) c = ρ Γ f ) c T j + ( ρ ω c (c c) c f ) c ( ρ S(c) c f ) c scalar turbulent advection is now directly modelled by a gradient diffusion assumption Γ T = C µ k2 / ɛ is a turbulent diffusivity (18)
19 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 19/63 Implications for scalar PDF methods (2/3) Previous work by Sabel nikov & Soulard (PRE 72, , 2005): Velocity U is modeled by a gaussian decorrelated velocity An equation is sought for the stochastic scalar θ so that its unweighted statistics give immediately the Favre PDF: p θ = f c (in particular θ = c), and thus the introduction of stochastic density is not needed Sabel nikov & Soulard obtained the following SPDE for the scalar θ: θ t + U j θ = ω c (θ θ ) + S(θ) U i = Ũ i Γ T 1 ρ ρ x i 1 2 Γ T x i + V i, V i = 2Γ T Ẇ i (19)
20 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 20/63 Implications for scalar PDF methods (3/3) The symbol denotes the Stratonovitch interpretation of the multiplicative noise - stochastic part of the convection term in (19). The mean value of V j θ is not zero and is equal to: ( ) V j θ x = θ j x W j Γ T + 1 Γ T θ 2 We remind that in Ito calculus this correlation is equal zero: θ V j = 0 W Stochastic convection and SPDE (19) can be rewritten in Ito form: V j θ = V j θ ρ θ t + ρ Ũ j θ ( ) θ Γ T + ρ V j θ = Eq (19) is in non conservative form + 1 Γ T θ 2 ( ) θ ρ Γ T ω c (θ θ ) + S(θ)
21 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 21/63 Zero-correlation time limit in (13)-(15) The same reasonings can be applied to system (13)-(15): But equivalence is between weighted statistics f θ = r θ r p θ = f c To facilitate our derivation, we consider the Simplified Langevin model. The tensor G ij reduces to : G ij = 1 τ rel δ ij, where 1 τ rel = We let τ rel 0 in (13)-(15), but keep ɛτ 2 rel finite : A solution to (20) is ( ) ɛ 4 C 0 k r U j r U j r + r C 0 ɛτrelẇj 2 = 0 (20) U i = Ũi + Γ T 1 ρ ρ x i Γ T x i + 2Γ T Ẇ i
22 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 22/63 Final system of SPDEs The final system of hyperbolic SPDEs is : U i = Ũi + Γ T 1 ρ r t + t (rθ) + ρ x i Γ T x i + 2Γ T Ẇ i (21) (r U j ) = 0 (22) ((rθ) U j ) = r ω c (θ θ r ) + rs (23) Eq. (21)-(23) are statistically equivalent to PDF equation (18) This remains true even in the case when τ rel is not small These SPDEs are different from, but stochastically equivalent to those of Sabel nikov & Soulard (PR 72, , 2005) They are in conservative form may be easier to implement to CFD codes The fluctuating velocity has the property: Ũ = 0
23 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 23/63 Simulation of a backward facing step with a RANS/EMC solver Work done in collaboration with M. Ourliac (PhD student)
24 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 24/63 Configuration Stoichiometric premixed air/methane flame over a backward facing step: Inlet velocity U e 55 m/s Inlet temperature T e 525 K Pressure P 10 5 P a Equivalence ratio 1 Inlet t.k.e k e 40 m 2 /s 2 Inlet turb. dissip. ɛ e 800 m 2 /s 3
25 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 25/63 RANS & EMC solver (1/2) RANS solver: Continuity Momentum Turbulent kinetic energy Turbulent dissipation ) ρ t + ( ρ Ũj ρ Ũ i t + = 0 ) ( ρ ŨjŨi ) ρ k t + ( ρ Ũ j k ) ρ ɛ t + ( ρ Ũjɛ = P + σ ij = ( ρ ν t P r k = ( ρ ν t P r ɛ k ) + P k d k ɛ ) + P ɛ d ɛ EMC solver: Stoch. Vel. U i = Ũi + Γ T 1 Stoch. Density Mass fraction Total enthalpy r t + ρ ρ x i (r U j ) = 0 t (ry k) + t (rh t) + Γ T x i + 2Γ T Ẇ i ( ) ((ry k ) U j ) = r ω k Y k Ỹ k + rs(y, h t ) ( ((rh t ) U j ) = r ω h h t h ) t
26 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 26/63 RANS & EMC solver (2/2) Coupling between the solvers: ρ, Γ T RANS solver ρ, Ũ, k, ɛ Y k, h t EMC solver P = ρ r T EMC solver based on Stratonovitch calculus
27 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 27/63 Parameters Chemistry: CH 4 + 2O 2 + βn 2 CO 2 + 2H 2 O + βn 2 Adiabatic walls Mesh (default value) : 100 x 40 cells Stochastic fields (default value): 50 Evolution of the CPU time per iteration and per number of cells against the number of stochastic fields CPU / Iteration / Pts N c
28 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 28/63 Analysis of the stochastic error (1/2) N=5 N=10 N=25 N=50 N= <T> 1400 <T> N=5 N=10 N=25 N=50 N= Z (a) Z (b) Figure 1: Influence of the number of stochastic fields N on the mean temperature vertical profiles, at x = 0.25 m (a) and x = 0.46 m (b)
29 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 29/63 Analysis of the stochastic error (2/2) T rms 400 T rms N=5 N=10 N=25 N=50 N= N=5 N=10 N=25 N=50 N= Z (a) Z (b) Figure 2: Influence of the number of stochastic fields N on the RMS temperature vertical profiles, at x = 0.25 m (a) and x = 0.46 m (b)
30 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 30/63 Results Mean and RMS temperature fields: T: y x 0.2 Trms: y x
31 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 31/63 Experiments (Magre & Moreau 1988) Mean and RMS temperature fields:
32 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 32/63 Profiles (1/3) Comparison of mean and RMS temperature with experiment: EMC solver Experiment <T> T rms Z Mean Temperature at x=210 mm 200 EMC solver Experiment Z RMS Temperature at x=210 mm
33 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 33/63 Profiles (2/3) Comparison of mean and RMS temperature with experiment: 2200 EMC solver Experiment <T> T rms Z EMC solver Experiment Z Mean Temperature at x=250 mm RMS Temperature at x=250 mm
34 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 34/63 Profiles (3/3) Comparison of mean and RMS temperature with experiment: EMC solver Experiment 1600 <T> T rms EMC solver Experiment Z Mean Temperature at x=460 mm Z RMS Temperature at x=460 mm
35 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 35/63 PDF (1/2) Comparison between computed and measured PDFs at x=210 mm EMC solver Experiment EMC solver Experiment PDF PDF T T z = 6mm z = 8mm
36 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 36/63 PDF (2/2) Comparison between computed and measured PDFs at x=250 mm EMC solver Experiment EMC solver Experiment PDF PDF T T z = 6mm z = 8mm
37 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 37/63 Conclusions for Part I Extension of EMC method for solving velocity-scalar PDF has been proposed Appropriate numerical scheme must be designed The derived SPDEs need to be tested on simple 1D tests The overall method needs to be compared against LMC method Corresponding EMC method for solving scalar PDF has also been proposed The EMC method for scalar PDF has been applied to the computation of a turbulent premixed methane flame over a backward facing step
38 Part II Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 38/63
39 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 39/63 Contents 1 Velocity-scalar PDF Modeled velocity-scalar PDF Lagrangian and Eulerian Monte Carlo methods 2 Scalar-velocity EMC method 3 Numerical scheme 4 Validation tests Riemann problem Return to Gaussianity Auto-ignition of a Methane-air mixture 5 Conclusions and perspectives
40 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 40/63 Turbulent reacting flow We consider a turbulent reacting flow with density ρ, velocity U, total enthalpy h t, and mass fractions Y α (α = 1,, N s ) For variable density flows, it is usual to work with Favre statistics: Reynolds PDF f = δ(y α (x, t) Y α )δ(h t (x, t) h t )δ(u(x, t) U) Favre PDF f = (ρ U, Y α, h t ) ρ f Reynolds and Favre averages are noted Q and Q : Q = ρq ρ Low Mach number assumption: pressure fluctuations are neglected in the thermodynamical equations ρ and the chemical source terms are functions of the reactive scalars
41 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 41/63 Modeled scalar-velocity PDF equation Modeled transport equation in 1D for the Favre PDF f: t ( ρ f ) + x ( ρ U f ) ( ρ U = (24) [ 1 ρ ( [ 1 ρ ρ h t ( ρ C φ ω Y α P ( x C 1ω U Ũ) ] ) f ( h t h ) ] ) t f P t C hω ( Y α Ỹ α ) f + ρ Sα f First line: transport in physical space, treated exactly C 0 ρ ω k 2 f U 2 Second line: pressure fluctuations and molecular diffusion modelled by the simplified Langevin model Third and fourth lines: mixing modelled by the IEM model P : pressure, S α : source term, k: turbulent kinetic energy, ω: turbulent frequency, C # : constants )
42 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 42/63 Solving PDF equation (24) PDF equation (24) has N d = N s + 7 dimensions where N is the number of scalars = large in practical applications Finite volume/element/difference methods: CPU cost increases exponentially with N d Not suitable Monte Carlo methods: CPU cost increases linearly with N d OK Two options: Lagrangian Monte Carlo (LMC) methods Eulerian Monte Carlo (EMC) methods
43 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 43/63 Lagrangian Monte Carlo methods (1/2) PDF represented by a set of N p stochastic particles Each particle is a sample of the physical properties of the system F Np = N p X n=1 w (n) δ(u U (n) )δ(y α Y α (n) )δ(h t h (n) t )δ(x x (n) ) ; ρ f e = F Np w (n), U (n), Y α (n), h (n) t, x (n) : mass, velocity, concentration, enthalpy, position of particle (n) Particles evolve according to SODEs (Pope, 1985): du (n) j dy (n) α dh (n) t dx (n) j = 1 P q dt C 1 ω U eu dt + C 0 ω e kdw (n) j (t) (25) ρ = C φ ω Y α (n) Y f α dt + S α (n) dt (26) = 1 ρ P t dt C hω h (n) t h e t dt (27) = U (n) j dt (28)
44 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 44/63 Lagrangian Monte Carlo methods (2/2) W (n) j are independent Brownian noises : = δ nm δ ij dt dw (n) j dw (m) k The absence of symbol in the stochastic product in eq. (26) denotes the Ito interpretation Numerous publications document the convergence and accuracy of LMC methods LMC used in many complex calculations (including LES) LMC implemented in commercial CFD codes
45 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 45/63 Eulerian Monte Carlo methods Principle PDF represented by stochastic Eulerian fields Each field is a sample of the physical properties of the system Each field evolves according to an SPDE Development of EMC methods is useful and stimulating LMC/EMC competition could push both approaches forward EMC methods already designed for joint scalar PDFs Theoretical works : Valiño (98), Sabel nikov & Soulard (05) Calculations : Naud et al (04), Sabel nikov & Soulard (06)
46 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 46/63 Eulerian Monte Carlo methods EMC methods also designed for joint scalar-velocity PDFs Theoretical derivation of SPDEs : Soulard & Sabel nikov (05) But no validation nor application Purpose of this work to propose a numerical scheme for solving the SPDEs to assess its performances
47 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 47/63 Scalar-velocity SPDEs SPDEs stochastically equivalent to PDF equation (24) can be derived from SODEs (25)-(28) by the notion of stochastic characteristic: r t + (r U) = 0 x ( (r U) + r U 2 + P ) = t x t t (r Y α) + x (r UY α) = r ω ( 1 r ρ (r H) + x (r UH) = r ω ( H h t ) ) P ( ) x r C 1ω U Ũ + r ρ ( ) Y α Ỹα P t + r (29) C 0 ω kẇ (30) (31) + r S α (32) U is the stochastic velocity field, H is the stochastic total enthalpy and Y α is the stochastic mass fraction. W is a standard Brownian noise and Ẇ is its time derivative r is the stochastic density. r is different from the physical density ρ
48 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 48/63 Scalar-velocity SPDEs We note Q s the average over the stochastic fields, and Q r the average weighted by the density r: Q r = rq s r s The corresponding PDFs are respectively denoted by f s and f r. The following equivalences exist: f = f r = r U, Y α, H s ρ f s and f = ρ ρ f r = r U, Y α, H s ρ with consistency conditions: r s = ρ and Favre and Reynolds averages are given by: Q = Q r = rq s ρ and Q = r ρ Q s r ρ s = 1 f s Mean pressure is given by: P = ρr g T = rr g T s
49 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 49/63 Numerical scheme: spatial discretization We use a finite volume method we must define a monotone numerical flux Characteristic waves of SPDEs (29)-(32) are complex due to the presence of both mean and instantaneous quantities upwinding monotone fluxes (like Godunov method) are difficult to use We decide to use instead a centered monotone flux Only requires limited information on the characteristic wave system Many different centered scheme exist We retain the GFORCE flux (Toro and Titarev, 2006) To achieve 2nd order, we use the GFORCE flux in conjunction with a 2nd/3rd order WENO interpolation
50 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 50/63 Numerical scheme: temporal discretization 3 requirements for the time discretization it should be strong stability preserving (SSP): in order to maintain monotonicity it should allow for the implicit treatment of stiff source terms it should be of weak order 2 In the deterministic case, implicit-explicit (IMEX) Runge-Kutta (RK) schemes with the SSP property exist (Pareschi and Russo, 2005) IMEX schemes are an interesting alternative to splitting methods the stiff and non-stiff parts of the system are discretized together The stiff part is treated implicitly, and the non-stiff part explicitly We extend these deterministic IMEX-RK schemes to the stochastic case We procede as in Tocino et al (2002) We obtain a new class of Stochastic IMEX schemes (S-IMEX)
51 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 51/63 First test: Riemann problem (1/4) Calculations are performed on a L = 1 m domain At initial time, the domain is divided into a left and a right state: for x < 0.5m Ũ L = 0 m/s ũ 2 L = 50 m 2 /s 2 ρ L = kg/m 3 = 10 5 P a P L for x > 0.5m Ũ R = 0 m/s ũ 2 R = 0 m 2 /s 2 ρ R = kg/m 3 = P a P R The turbulent frequency is taken equal to ω = 200 s 1 N = 400 stochastic fields and N x = 320 cells are used
52 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 52/63 First test: Riemann problem (2/4) 10 Stochastic velocity 120 Stochastic velocity Solution to the Riemann problem Velocity (m/s) 0 Velocity (m/s) x (m) (a) x (m) First test: examples of stochastic velocity fields (a) at t=0 s and (b) at t = s (b)
53 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 53/63 First test: Riemann problem (3/4) Computed mean density Solution to the Riemann problem Computed mean pressure Solution to the Riemann problem Density (kg/m^3) Pressure (Pa) x (m) x (m) (a) First test: profiles of mean density (a), mean pressure(b) at t = s (b)
54 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 54/63 First test: Riemann problem (4/4) 100 Computed mean velocity Solution to the Riemann problem Velocity (m/s) Velocity variance (m^2/s^2) x (m) (c) x (m) (d) First test: profiles of mean velocity (c) and velocity variance (d) at t = s. Mean profiles are compared against the exact solution to the Riemann problem 0
55 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 55/63 Second test: Return to Gaussiannity (1/3) For this test only, we replace the Brownian coefficient by 2ωσ 2 instead of C 0 ω k σ 2 is a constant and is taken equal to 1 the velocity field PDF should tend to a Gaussian with variance σ 2 At the left boundary: Dirac distributions are imposed, with means: Ũ = 10.4m/s, ρ = 0.146kg/m 3, P = 10 5 P a At initial time, the stochastic fields are initialized with the left boundary conditions The domain has a length L = 0.04m and is discretized with N x = 40 cells The turbulent frequency is taken equal to ω = 2500s 1
56 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 56/63 Second test: Return to Gaussiannity (2/3) 15 Variance Skewness Flatness Hyper flatness Moment value 5 3 velocity x (a) Second test: (a): evolution of the velocity variance,skewness,flatness and hyper-flatness ; (b) : stochastic velocity fields from which the moments are computed x (b)
57 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 57/63 Second test: Return to Gaussiannity (3/3) Stochastic error Best fit: 3.8 N^-.72 Stochastic error Best fit: 15.8 N^-.55 Stochastic error for the velocity flatness 0.1 Stochastic error for the velocity hyper-flatness Number of stochastic fields N (a) Number of stochastic fields N (b) Second test : convergence rate for the velocity flatness (a) and hyper-flatness (b)
58 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 58/63 Third test: Auto-ignition (1/4) At inlet, a stoichiometric methane-air mixture is injected: temperature T in = 1500K pressure P in = P a Gaussian random velocity with Ũ in = 10.4ms 1 and variance ũ 2 in = 1 m 2 s 2 At outlet, zero gradients, except for the pressure which value is fixed at P in The length of the 1D domain is L = m The methane-air chemistry is dealt with a simple one step global reaction (Westbrook, 1984)
59 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 59/63 Third test: Auto-ignition (2/4) 3500 t=1 ms t=1.5 ms t=2 ms t=3 ms t=1 ms t=1.5 ms t=2 ms t=3 ms Mean temperature (K) Temperature variance (K^2) x (m) x (m) (a) (b) Third test: Time evolution of the profiles of the mean (a) and variance (b) of the temperature
60 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 60/63 Third test: Auto-ignition (3/4) Mean velocity (m/s) t=1 ms t=1.5 ms t=2 ms t=3 ms Velocity variance (m^2/s^2) t=1 ms t=1.5 ms t=2 ms t=3 ms x (m) x (m) (a) (b) Third test: Time evolution of the profiles of the mean (a) and variance (b) of the velocity
61 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 61/63 Third test: Auto-ignition (4/4) mean velocity convergence rate velocity variance convergence rate mean temperature convergence rate temperature variance convergence rate Theoretical 10 mean velocity convergence coeff. velocity variance convergence coeff. mean temperature convergence coeff. temperature variance convergence coeff. Convergence rate Convergence coefficient t (s) t (s) (a) (b) Third test: Time evolution of the convergence rates (a) and coefficients (b) for the means and variances of velocity and temperature
62 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 62/63 Third test: Auto-ignition Spatial error Best fit: ^4 N_x^-2.1 Spatial error Best fit: ^2 N_x^-1.65 Spatial error for the mean velocity Spatial error for the velocity variance Number of cells N_x Number of cells N_x (a) (b) Third test: Spatial convergence of the mean (a) and variance (b) of the velocity at t = 2.5ms 100
63 Workshop and Minicourse Conceptual Aspects of Turbulence, February 2008 p. 63/63 Conclusions for Part II A numerical scheme for solving the SPDEs obtained in Soulard and Sabel nikov has been proposed finite volume scheme based on a monotone centered second order numerical flux new weak second order Runge-Kutta scheme, whith the SSP property, and an implicit treatment of chemical source terms (S-IMEX) 1D validation tests have been performed Monotonicity was verified on a Riemann problem Return to Gaussiannity was also checked An auto-ignition problem was studied Further developments of this work : extension of the numerical method to 2D problems calculation of a practical cases
LES 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 informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
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 informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationThe Hybrid Method for the PDF Equations of Turbulent Reactive Flows: Consistency Conditions and Correction Algorithms
Journal of Computational Physics 172, 841 878 (2001) doi:10.1006/jcph.2001.6861, available online at http://www.idealibrary.com on The Hybrid Method for the PDF Equations of Turbulent Reactive Flows: Consistency
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 information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable
More informationInsights into Model Assumptions and Road to Model Validation for Turbulent Combustion
Insights into Model Assumptions and Road to Model Validation for Turbulent Combustion Venke Sankaran AFRL/RQR 2015 AFRL/RQR Basic Research Review UCLA Jan 20, 2015 AFTC PA Release# 15011, 16 Jan 2015 Goals
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationMass Transfer in Turbulent Flow
Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationThe Truth about diffusion (in liquids)
The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion
More informationOn two fractional step finite volume and finite element schemes for reactive low Mach number flows
Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes
More informationPDF modeling and simulation of premixed turbulent combustion
Monte Carlo Methods Appl. Vol. 4 No. 4 (8), pp. 343 377 DOI. / MCMA.8.6 c de Gruyter 8 PDF modeling and simulation of premixed turbulent combustion Michael Stöllinger and Stefan Heinz Abstract. The use
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationLOW TEMPERATURE MODEL FOR PREMIXED METHANE FLAME COMBUSTION
ISTP-16, 2005, PRAGUE 16TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA LOW TEMPERATURE MODEL FOR PREMIXED METHANE FLAME MBUSTION M. Forman, J.B.W.Kok,M. Jicha Department of Thermodynamics and Environmental
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
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 informationReactive Transport in Porous Media
Reactive Transport in Porous Media Daniel M. Tartakovsky Department of Mechanical & Aerospace Engineering University of California, San Diego Alberto Guadagnini, Politecnico di Milano, Italy Peter C. Lichtner,
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 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 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 informationPDF Modeling and Simulation of Premixed Turbulent Combustion
Monte Carlo Methods Appl. Vol. No. (), pp. 43 DOI 5 / MCMA.7. c de Gruyter PDF Modeling and Simulation of Premixed Turbulent Combustion Michael Stöllinger and Stefan Heinz Abstract. The use of probability
More informationLES-PDF SIMULATION OF A HIGHLY SHEARED TURBULENT PILOTED PREMIXED FLAME
LES-PDF SIMULATION OF A HIGHLY SHEARED TURBULENT PILOTED PREMIXED FLAME Abstract V. N. Prasad, K. H. Luo and W. P. Jones k.h.luo@soton.ac.uk School of Engineering Sciences, University of Southampton, Southampton
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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
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 informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
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 informationarxiv: v2 [physics.class-ph] 4 Apr 2009
arxiv:0903.4949v2 [physics.class-ph] 4 Apr 2009 Geometric evolution of the Reynolds stress tensor in three-dimensional turbulence Sergey Gavrilyuk and Henri Gouin Abstract The dynamics of the Reynolds
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
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 informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
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 informationarxiv: v1 [physics.chem-ph] 6 Oct 2011
Calculation of the Minimum Ignition Energy based on the ignition delay time arxiv:1110.1163v1 [physics.chem-ph] 6 Oct 2011 Jens Tarjei Jensen a, Nils Erland L. Haugen b, Natalia Babkovskaia c a Department
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationM U S C L E S. Modelling of UnSteady Combustion in Low Emission Systems. Project No. GRD Contract No. G4RD-CT
M U S C L E S Modelling of UnSteady Combustion in Low Emission Systems Project No. GRD-200-4098 Work Package 2 Prediction of unsteady reacting flow and validation Task 2. Prediction of lean blow out in
More informationContinuum Limit of Forward Kolmogorov Equation Friday, March 06, :04 PM
Continuum Limit of Forward Kolmogorov Equation Friday, March 06, 2015 2:04 PM Please note that one of the equations (for ordinary Brownian motion) in Problem 1 was corrected on Wednesday night. And actually
More informationA stochastic particle system for the Burgers equation.
A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u
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 informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationBoundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity
Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with
More informationMarching on the BL equations
Marching on the BL equations Harvey S. H. Lam February 10, 2004 Abstract The assumption is that you are competent in Matlab or Mathematica. White s 4-7 starting on page 275 shows us that all generic boundary
More informationAccelerating incompressible fluid flow simulations on hybrid CPU/GPU systems
Accelerating incompressible fluid flow simulations on hybrid CPU/GPU systems Yushan Wang 1, Marc Baboulin 1,2, Karl Rupp 3,4, Yann Fraigneau 1,5, Olivier Le Maître 1,5 1 Université Paris-Sud, France 2
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 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 informationInvestigation of Scalar Dissipation Rate Fluctuations in Non-Premixed Turbulent Combustion Using a Stochastic Approach
Investigation of Scalar Dissipation Rate Fluctuations in Non-Premixed Turbulent Combustion Using a Stochastic Approach Heinz Pitsch and Sergei Fedotov Flow Physics and Computation Division, Stanford University,
More informationOutline: Phenomenology of Wall Bounded Turbulence: Minimal Models First Story: Large Re Neutrally-Stratified (NS) Channel Flow.
Victor S. L vov: Israeli-Italian March 2007 Meeting Phenomenology of Wall Bounded Turbulence: Minimal Models In collaboration with : Anna Pomyalov, Itamar Procaccia, Oleksii Rudenko (WIS), Said Elgobashi
More informationConsistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames
Consistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames F. Ferraro, Y. Ge and M. Pfitzner Institut für Thermodynamik, Fakultät für Luft- und Raumfahrrttechnik Universität
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 informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
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 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 informationThis paper might be a pre-copy-editing or a post-print author-produced.pdf of an article accepted for publication. For the definitive
This paper might be a pre-copy-editing or a post-print author-produced.pdf of an article accepted for publication. For the definitive publisher-authenticated version, please refer directly to publishing
More informationSuperparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows
Superparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows June 2013 Outline 1 Filtering Filtering: obtaining the best statistical estimation of a nature
More informationFootprints: outline Üllar Rannik University of Helsinki
Footprints: outline Üllar Rannik University of Helsinki -Concept of footprint and definitions -Analytical footprint models -Model by Korman and Meixner -Footprints for fluxes vs. concentrations -Footprints
More informationKinetic relaxation models for reacting gas mixtures
Kinetic relaxation models for reacting gas mixtures M. Groppi Department of Mathematics and Computer Science University of Parma - ITALY Main collaborators: Giampiero Spiga, Giuseppe Stracquadanio, Univ.
More informationAerosol nucleation and growth in a turbulent jet using the Stochastic Fields method
Loughborough University Institutional Repository Aerosol nucleation and growth in a turbulent jet using the Stochastic Fields method This item was submitted to Loughborough University's Institutional Repository
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 informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationEulerian models. 2.1 Basic equations
2 Eulerian models In this chapter we give a short overview of the Eulerian techniques for modelling turbulent flows, transport and chemical reactions. We first present the basic Eulerian equations describing
More informationComputational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor
Computational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor SD Ravi 1, NKS Rajan 2 and PS Kulkarni 3 1 Dept. of Aerospace Engg., IISc, Bangalore, India. ravi@cgpl.iisc.ernet.in
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 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 informationProbability density function modeling of scalar mixing from concentrated sources in turbulent channel flow
PHYSICS OF FLUIDS 19, 115106 2007 Probability density function modeling of scalar mixing from concentrated sources in turbulent channel flow J. Bakosi, a P. Franzese, and Z. Boybeyi College of Science,
More informationComputation of hypersonic shock boundary layer interaction on a double wedge using a differential Reynolds Stress Model
Computation of hypersonic shock boundary layer interaction on a double wedge using a differential Reynolds Stress Model A. Bosco AICES, Templergraben 55, Aachen, 52056, Germany B. Reinartz CATS, Theaterplatz
More informationExtension to moving grids
Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical
More informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationSistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas
Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Pedro González-Casanova Henríquez Unidad de Investigación en Cómputo Aplicado DGSCA, UNAM Elvira Torondel y Ricard Garrido
More informationhydrogen auto ignition in a turbulent co flow of heated air with LES and CMC approach.
biblio.ugent.be The UGent Institutional Repository is the electronic archiving and dissemination platform for all UGent research publications. Ghent University has implemented a mandate stipulating that
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 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 informationKrylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations
Accepted Manuscript Krylov single-step implicit integration factor WENO methods for advection diffusion reaction equations Tian Jiang, Yong-Tao Zhang PII: S0021-9991(16)00029-2 DOI: http://dx.doi.org/10.1016/j.jcp.2016.01.021
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 informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationChapter 8 STATISTICAL ASSESSMENT OF NUMERICAL DIFFUSION
84 hapter 8 STATISTIAL ASSESSMENT OF NUMERIAL DIFFUSION The assessment of the performance of high-order numerical methods regarding the numerical diffusion introduced in the solution of the scalar transport
More informationFlow Structure Investigations in a "Tornado" Combustor
Flow Structure Investigations in a "Tornado" Combustor Igor Matveev Applied Plasma Technologies, Falls Church, Virginia, 46 Serhiy Serbin National University of Shipbuilding, Mikolayiv, Ukraine, 545 Thomas
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 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 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 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 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 informationChapter 13. Eddy Diffusivity
Chapter 13 Eddy Diffusivity Glenn introduced the mean field approximation of turbulence in two-layer quasigesotrophic turbulence. In that approximation one must solve the zonally averaged equations for
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 informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 05-05 Efficient Methods for Solving Advection Diffusion Reaction Equations S. van Veldhuizen ISSN 1389-6520 Reports of the Department of Applied Mathematical Analysis
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationIMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE FLUCTUATIONS
Proceedings of the Combustion Institute, Volume 9, 00/pp. 1971 1978 IMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE
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 informationFaster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs
Faster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs Christopher P. Stone, Ph.D. Computational Science and Engineering, LLC Kyle Niemeyer, Ph.D. Oregon State University 2 Outline
More information