On a New Diagram Notation for the Derivation of Hyperbolic Moment Models

Size: px
Start display at page:

Download "On a New Diagram Notation for the Derivation of Hyperbolic Moment Models"

Transcription

1 On a New Diagram Notation for the Derivation of Hyperbolic Moment Models Julian Koellermeier, Manuel Torrilhon, Yuwei Fan March 17th, 2017 Stanford University J. Koellermeier 1 / 57 of Hyperbolic Moment Models

2 Germany s largest technical university 45k students, 10k staff J. Koellermeier 2 / 57 of Hyperbolic Moment Models

3 Outline Conclusion J. Koellermeier 3 / 57 of Hyperbolic Moment Models

4 Introduction Motivation Rarefied Gas Dynamics Boltzmann Equation Aim Derive hyperbolic PDE systems for rarefied gas flows Extension of standard fluid dynamic equations Reentry flows Micro channel flows J. Koellermeier 4 / 57 of Hyperbolic Moment Models

5 Introduction Motivation Rarefied Gas Dynamics Boltzmann Equation Aim Derive hyperbolic PDE systems for rarefied gas flows Extension of standard fluid dynamic equations Reentry flows Micro channel flows Importance of Hyperbolicity Physical solutions with bounded propagation speeds Well-posedness and stability of the solution J. Koellermeier 4 / 57 of Hyperbolic Moment Models

6 Rarefied Gas Dynamics Motivation Rarefied Gas Dynamics Boltzmann Equation Goal solve and simulate flow problems involving rarefied gases Knudsen number distinguish flow regimes by orders of the Knudsen number Kn = mean free path length reference length = λ L J. Koellermeier 5 / 57 of Hyperbolic Moment Models

7 Boltzmann Transport Equation Motivation Rarefied Gas Dynamics Boltzmann Equation t f (t, x, c) + c i f (t, x, c) = S(f ) x i PDE for particles probability density function f (t, x, c) Describes change of f due to transport and collisions Collision operator S Usually a 7-dimensional phase space Free motion Collision Free path J. Koellermeier 6 / 57 of Hyperbolic Moment Models

8 Boltzmann Transport Equation 1D Motivation Rarefied Gas Dynamics Boltzmann Equation f (t, x, c) + c f (t, x, c) = S(f ) t x PDE for particles probability density function f (t, x, c) Describes change of f due to transport and collisions Collision operator S 3-dimensional phase space Free motion Collision Free path J. Koellermeier 7 / 57 of Hyperbolic Moment Models

9 Macroscopic Quantities Motivation Rarefied Gas Dynamics Boltzmann Equation f (t, x, c) is related to the macroscopic quantities density ρ(t, x), velocity v(t, x), temperature θ(t, x) ρ(t, x) = ρ(t, x)v(t, x) = 1 2 ρ(t, x)θ(t, x) ρ(t, x)v(t, x)2 = R R R f (t, x, c) dc cf (t, x, c) dc 1 2 c2 f (t, x, c) dc J. Koellermeier 8 / 57 of Hyperbolic Moment Models

10 Discretization in Velocity Space Grad HME QBME Review of Hyperbolic Moment Models J. Koellermeier 9 / 57 of Hyperbolic Moment Models

11 A Transformed Velocity Variable Discretization in Velocity Space Grad HME QBME f c ξ(t, x, c) := c v(t, x) θ(t, x) = f (ξ) = ρ ) exp ( ξ2 2πθ 2 10 f (c) = ρ ) (c v)2 exp ( 2πθ 2θ c J. Koellermeier 10 / 57 of Hyperbolic Moment Models

12 A Transformed Velocity Variable Discretization in Velocity Space Grad HME QBME f c ξ(t, x, c) := c v(t, x) θ(t, x) = f (ξ) = ρ ) exp ( ξ2 2πθ 2 10 f (c) = ρ ) (c v)2 exp ( 2πθ 2θ c Transformed velocity space reduces numerical complexity J. Koellermeier 10 / 57 of Hyperbolic Moment Models

13 A Transformed Velocity Variable Discretization in Velocity Space Grad HME QBME f c ξ(t, x, c) := c v(t, x) θ(t, x) = f (c) = ρ ) (c v)2 exp ( 2πθ 2θ c f (ξ) = ρ ) exp ( ξ2 2πθ 2 Transformed velocity space reduces numerical complexity J. Koellermeier 10 / 57 of Hyperbolic Moment Models

14 Discretization in Velocity Space Grad HME QBME A Transformed Velocity Variable f HcL c v (t, x) ξ(t, x, c) := p θ(t, x) = f HxL ρ (c v )2 f (c) = exp 2θ 2πθ c ρ ξ f (ξ) = exp 2 2πθ Transformed velocity space reduces numerical complexity J. Koellermeier 11 / 57 of Hyperbolic Moment Models x

15 Model Order Reduction Discretization in Velocity Space Grad HME QBME Ansatz: Expansion f (t, x, c) = M α=0 f α (t, x)φ v,θ α (ξ) Basis function: is weighted Hermite polynomial φ v,θ α (ξ) = 1 ) θ α+1 2 exp ( ξ2 H α (ξ) 2π 2 J. Koellermeier 12 / 57 of Hyperbolic Moment Models

16 Model Order Reduction Discretization in Velocity Space Grad HME QBME Ansatz: Expansion f (t, x, c) = M α=0 f α (t, x)φ v,θ α (ξ) Basis function: is weighted Hermite polynomial φ v,θ α (ξ) = 1 ) θ α+1 2 exp ( ξ2 H α (ξ) 2π 2 Reduction of Complexity One PDE for f (t, x, c) that is high-dimensional System of PDEs for ρ(t, x), v(t, x), θ(t, x), f α (t, x) that is low-dimensional J. Koellermeier 12 / 57 of Hyperbolic Moment Models

17 Grad s Method [Grad, 1949] Discretization in Velocity Space Grad HME QBME Galerkin Approach Standard method Multiplication with test function and integration Grad result t u M + A Grad x u M = 0, A Grad = u 4 = (ρ, v, θ, f 3, f 4 ) T v ρ θ ρ v θ v ρ 0 ρθ 0 4f 3 2 v 4 f 3θ 3f ρ 5f θ v Loss of hyperbolicity J. Koellermeier 13 / 57 of Hyperbolic Moment Models

18 Discretization in Velocity Space Grad HME QBME Hyperbolic Moment Equations (HME) [Cai et al., 2012] Modification of equations Based on Grad s method Modification of last equation to achieve hyperbolicity HME result t u M + A HME x u M = 0, A HME = u 4 = (ρ, v, θ, f 3, f 4 ) T v ρ θ ρ v θ v ρ 0 ρθ 0 4f 3 2 v 4 ρ 0 f 3 θ v f 3θ Globally hyperbolic for every state vector u M J. Koellermeier 14 / 57 of Hyperbolic Moment Models

19 Discretization in Velocity Space Grad HME QBME Quadrature-Based Moment Equations (QBME) [JK, 2013] Quadrature-Based Projection Approach Based on Grad s method Substitution of integrals by Gaussian quadrature QBME result t u M + A QBME x u M = 0, A QBME = u 4 = (ρ, v, θ, f 3, f 4 ) T v ρ θ ρ v θ v ρ 0 ρθ 0 4f f 4 θ v 4 f 3θ ρ 5f 4 f 3 θ+ 15f 4 ρθ v Globally hyperbolic for every state vector u M J. Koellermeier 15 / 57 of Hyperbolic Moment Models

20 Additions and Developments Discretization in Velocity Space Grad HME QBME Hyperbolic Moment Models Extension to n-d Framework for hyperbolic models... Applications Relativistic Boltzmann equation Quantum Gas... J. Koellermeier 16 / 57 of Hyperbolic Moment Models

21 Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations J. Koellermeier 17 / 57 of Hyperbolic Moment Models

22 Derivation of Moment System Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations f (t, x, c) = ( fα t v f α 1 f α 2 2 M α=0 f α (t, x)φ v,θ α (ξ) f (t, x, c) + c f (t, x, c) = 0 t x + v f α x + θ f α 1 + (α + 1) f α+1 x x + t + (vf α 1 + θf α 2 + (α + 1)f α ) v x + θ t +1 2 (vf α 2 + θf α 3 + (α + 1)f α 1 ) θ x ) φ v,θ α (ξ) = 0 J. Koellermeier 18 / 57 of Hyperbolic Moment Models

23 Derivation of Moment System (2) Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations J. Koellermeier 19 / 57 of Hyperbolic Moment Models

24 Towards a Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations Previous derivations Concise and mathematically elegant way to derive models Very theoretical and no direct insight into equations Diagram-Based Notation Better understanding of hyperbolic regularization Derivation of new models J. Koellermeier 20 / 57 of Hyperbolic Moment Models

25 Derivation 1 Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations Equation and ansatz f (t, x, c) = t f (t, x, c) + c x f (t, x, c) = 0 M α=0 f α (t, x)φ v,θ c v(t, x) α (ξ), ξ = θ(t, x) s (f α φ α (ξ)) = s f α φ α (ξ) + f α s φ α (ξ) = s f α φ α (ξ) + f α ( θ φ α (ξ) s θ + ξ φ α (ξ) s ξ) ( = s f α φ α (ξ) + f α ( θ φ α (ξ) s θ + ξ φ α (ξ) 1 s v ξ ) θ 2θ sθ J. Koellermeier 21 / 57 of Hyperbolic Moment Models

26 Derivation 2 Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations Equation and ansatz f (t, x, c) = t f (t, x, c) + c x f (t, x, c) = 0 M α=0 f α (t, x)φ v,θ c v(t, x) α (ξ), ξ = θ(t, x) ( s (f α φ α (ξ)) = s f α φ α (ξ) + f α ( θ φ α (ξ) s θ + ξ φ α (ξ) ) c f α φ α (ξ) = ( v + θξ f α φ α (ξ) 1 θ s v ξ 2θ sθ ) J. Koellermeier 22 / 57 of Hyperbolic Moment Models

27 Choice of Basis Function Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations weighted Hermite polynomial φ α (ξ) = 1 ) θ α+1 2 H α (ξ) exp ( ξ2 2π 2 1. ξ derivative: ξ φ α(ξ) = θφ α+1 (ξ) 2. θ derivative: θ φ α(ξ) = α+1 2θ φ α(ξ) 3. multiplication with ξ: ξφ α (ξ) = θφ α+1 (ξ) + α θ φ α 1 (ξ) J. Koellermeier 23 / 57 of Hyperbolic Moment Models

28 ξ derivative diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ξ φ α(ξ) = θφ α+1 (ξ) J. Koellermeier 24 / 57 of Hyperbolic Moment Models

29 θ derivative diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations θ φ α(ξ) = α + 1 2θ φ α(ξ) J. Koellermeier 25 / 57 of Hyperbolic Moment Models

30 ξ multiplication diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ξφ α (ξ) = θφ α+1 (ξ) + α θ φ α 1 (ξ) J. Koellermeier 26 / 57 of Hyperbolic Moment Models

31 Diagram Notation Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations Aim 1 Extend diagram notation from basis function to whole equation 2 Split terms of equation into steps 3 Use single diagram for each step t f (t, x, c) + c x f (t, x, c) = 0 s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( c f α φ α (ξ) = v + ) θξ f α φ α (ξ) 3 ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 27 / 57 of Hyperbolic Moment Models

32 1 diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ( f α ξ φ α (ξ) = f α ) θ φ α+1 (ξ) J. Koellermeier 28 / 57 of Hyperbolic Moment Models

33 2 diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ( f α φ α 1 s v ξ ) ( ) θ 2θ sθ = f α α sθ 2 θ 3 φ α 1 sv φ α sθ θ 2 θ φ α+1 J. Koellermeier 29 / 57 of Hyperbolic Moment Models

34 2a diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ( s f α ) φ α (ξ) = s f α φ α (ξ) J. Koellermeier 30 / 57 of Hyperbolic Moment Models

35 2b diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ( f α θ φ α (ξ) s θ = f α α + 1 ) 2θ sθ s φ α (ξ) J. Koellermeier 31 / 57 of Hyperbolic Moment Models

36 3 diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ( v + ) θξ f α φ α (ξ) = f α (θφ α 1 (ξ) + vφ α (ξ) + αφ α+1 (ξ)) J. Koellermeier 32 / 57 of Hyperbolic Moment Models

37 Grad diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 33 / 57 of Hyperbolic Moment Models

38 Grad diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 33 / 57 of Hyperbolic Moment Models

39 Grad diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 33 / 57 of Hyperbolic Moment Models

40 Grad diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 33 / 57 of Hyperbolic Moment Models

41 Grad diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 33 / 57 of Hyperbolic Moment Models

42 Grad diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 33 / 57 of Hyperbolic Moment Models

43 Single Basis Function Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations ( (1) : f α φ α ( (2) : f α φ α θ α + 1 ) 2θ xθ u f α φ α, ) ( ) α θ 3 xθ u f α φ α J. Koellermeier 34 / 57 of Hyperbolic Moment Models

44 Grad diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations J. Koellermeier 35 / 57 of Hyperbolic Moment Models

45 HME diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations J. Koellermeier 36 / 57 of Hyperbolic Moment Models

46 QBME diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations J. Koellermeier 37 / 57 of Hyperbolic Moment Models

47 Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations Simplified Hyperbolic Moment Equations (SHME) QBME and HME cut off terms during derivation Corresponds to new formulas, different adaptivity How much adaptivity is necessary? ( s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 1 s v ξ ) θ 2θ sθ s f α φ α (ξ) ( c f α φ α (ξ) = v + ) θξ f α φ α (ξ) Neglect adaptivity of s ( ) Keep adaptivity of c ( ) J. Koellermeier 38 / 57 of Hyperbolic Moment Models

48 SHME diagram Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations s (f α φ α (ξ)) = s f α φ α (ξ) + f α θ φ α (ξ) s θ + f α ξ φ α (ξ) 2a + 2b ( 1 s v ξ ) θ 2θ sθ J. Koellermeier 39 / 57 of Hyperbolic Moment Models

49 Preliminaries Diagram Notation Simplified Hyperbolic Moment Equations Simplified Hyperbolic Moment Equations (SHME) Simplified derivation Neglecting parts of adaptivity SHME result v ρ t u M + A SHME x u M = 0, θ ρ v A SHME = 0 2θ v u 4 = (ρ, v, θ, f 3, f 4 ) T ρ 0 ρθ v θ v Globally hyperbolic (equilibrium) J. Koellermeier 40 / 57 of Hyperbolic Moment Models

50 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock J. Koellermeier 41 / 57 of Hyperbolic Moment Models

51 Conservative PDE systems Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock Standard conservative PDE system t u + x F(u) = 0 J. Koellermeier 42 / 57 of Hyperbolic Moment Models

52 Conservative PDE systems Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock Standard conservative PDE system Basic Finite Volume scheme u n+1 i t u + x F(u) = 0 = u n i t x Numerical flux F needed i+ 1 2 Conservation property by design ( ) F F i+ 1 i J. Koellermeier 42 / 57 of Hyperbolic Moment Models

53 Non-conservative PDE systems Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock Non-conservative PDE system t u + A(u) x u = 0 J. Koellermeier 43 / 57 of Hyperbolic Moment Models

54 Non-conservative PDE systems Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock Non-conservative PDE system t u + A(u) x u = 0 Standard conservative PDE system t u + x F(u) = 0 Can be written in conservative form iff A(u) = F(u) u In general no flux function available What about partially conservative systems? Special numerical methods are needed J. Koellermeier 43 / 57 of Hyperbolic Moment Models

55 Numerical Methods Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock Wave Propagation scheme [LeVeque, 1997] 2 nd order upwind type scheme Castro scheme [Castro, Pares, 2004] 1 st order upwind type scheme PRICE-C scheme [Canestrelli, 2009] 1 st order centered scheme J. Koellermeier 44 / 57 of Hyperbolic Moment Models

56 1D Shock Tube Test Case Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock Riemann problem with BGK collision operator t w + A(w) x w = 1 Pw, x [ 2, 2] τ ρ L = 7, ρ R = 1 Variable vector w = (ρ, u, θ, f 3, f 4 ) Relaxation time τ = Kn ρ non-linear J. Koellermeier 45 / 57 of Hyperbolic Moment Models

57 Model Equations Grad model HME model A Grad = A HME = Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock v ρ θ ρ v θ v ρ 0 ρθ 0 4f 3 2 v 4 f 3θ 3f ρ 5f θ v v ρ θ ρ v θ v ρ 0 ρθ 0 4f 3 2 v 4 ρ 0 f 3 θ v f 3θ J. Koellermeier 46 / 57 of Hyperbolic Moment Models

58 Model Equations 2 Grad model QBME model A Grad = A QBME = Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock v ρ θ ρ v θ v ρ 0 ρθ 0 4f 3 2 v 4 f 3θ 3f ρ 5f θ v v ρ θ ρ v θ v ρ 0 ρθ 0 4f f 4 θ v 4 f 3θ ρ 5f 4 f 3 θ+ 15f 4 ρθ v J. Koellermeier 47 / 57 of Hyperbolic Moment Models

59 Model Equations 3 Grad model SHME model A Grad = Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock v ρ θ ρ v θ v ρ 0 ρθ 0 4f 3 2 v 4 f 3θ 3f ρ 5f θ v v ρ θ ρ v A SHME = 0 2θ v ρ 0 ρθ v θ v J. Koellermeier 48 / 57 of Hyperbolic Moment Models

60 Method Comparison, Kn = 0.05 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock 7 PRICE UPRICE Castro WP J. Koellermeier 49 / 57 of Hyperbolic Moment Models

61 Method Comparison, Kn = 0.5 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock 7 PRICE UPRICE Castro WP J. Koellermeier 50 / 57 of Hyperbolic Moment Models

62 Model Comparison, Kn = 0.05 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock 7 HME SHME DVM QBME J. Koellermeier 51 / 57 of Hyperbolic Moment Models

63 Model Comparison, Kn = 0.5 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock 7 HME SHME DVM QBME J. Koellermeier 52 / 57 of Hyperbolic Moment Models

64 Averaged Solution, Kn = 0.5 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock 7 QBME8 QBME9 QBME8+9 DVM J. Koellermeier 53 / 57 of Hyperbolic Moment Models

65 Model Convergence, Kn = 0.5 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock 10 1 rho u p theta Figure: HME and QBME Figure: SHME J. Koellermeier 54 / 57 of Hyperbolic Moment Models

66 2D bow shock,ma = 3, Kn = Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock ρ Figure: 1st order Figure: 2nd order J. Koellermeier 55 / 57 of Hyperbolic Moment Models

67 2D bow shock,ma = 3, Kn = 0.5 Non-conservative Numerical Methods 1D Shock Tube Test Case 2D Bow Shock ρ Figure: 1st order Figure: 2nd order J. Koellermeier 56 / 57 of Hyperbolic Moment Models

68 Summary and Further Work Summary Hyperbolic Moment Equations Hierarchical models for rarefied gases Different models are available Diagram notation Reveal details of derivation Allow derivation of new models Further work M-adaptivity J. Koellermeier 57 / 57 of Hyperbolic Moment Models

69 Summary and Further Work Summary Hyperbolic Moment Equations Hierarchical models for rarefied gases Different models are available Diagram notation Reveal details of derivation Allow derivation of new models Further work M-adaptivity Thank you for your attention J. Koellermeier 57 / 57 of Hyperbolic Moment Models

70 References Summary J. Koellermeier, R.P. Schaerer and M. Torrilhon A Framework for Hyperbolic Approximation of Kinetic Equations Using Quadrature-Based Projection Methods, Kinet. Relat. Mod. 7(3) (2014), J. Koellermeier, M. Torrilhon Hyperbolic Moment Equations Using Quadrature-Based Projection Methods, Proceedings of the 29th Symposium on Rarefied Gas Dynamics, AIP Conf. Proc (2014), , Y. Fan, J. Koellermeier, J. Li, R. Li, M. Torrilhon A framework on the globally hyperbolic moment method for kinetic equations using operator projection method J. Stat. Phys., 162(2) (2016), Z. Cai, Y. Fan and R. Li Globally hyperbolic regularization of Grad s moment system, Comm. Pure Appl. Math., 67(3) (2014), H. Grad On the kinetic theory of rarefied gases, Comm. Pure Appl. Math., 2(4) (1949), J. Koellermeier 58 / 57 of Hyperbolic Moment Models

Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form

Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative

More information

Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models

Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Julian Koellermeier, Manuel Torrilhon October 12th, 2015 Chinese Academy of Sciences, Beijing Julian Koellermeier,

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Using Quadrature-Based Projection Methods Julian Köllermeier, Manuel Torrilhon 29th International Symposium on Rarefied Gas Dynamics, Xi an July 14th, 2014 Julian Köllermeier, Manuel Torrilhon 1 / 12 Outline

More information

Simplified Hyperbolic Moment Equations

Simplified Hyperbolic Moment Equations Simplified Hyperbolic Moment Equations Julian Koellermeier and Manuel Torrilhon Abstract Hyperbolicity is a necessary property of model equations for the solution of the BGK equation to achieve stable

More information

REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS

REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS 1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,

More information

Hyperbolic Approximation of Kinetic Equations Using Quadrature-Based Projection Methods

Hyperbolic Approximation of Kinetic Equations Using Quadrature-Based Projection Methods Master Thesis Hyperbolic Approximation of Kinetic Equations Using Quadrature-Based Projection Methods Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational

More information

Physical Modeling of Multiphase flow. Boltzmann method

Physical Modeling of Multiphase flow. Boltzmann method with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory

More information

A Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids

A Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids Diese Arbeit wurde vorgelegt am Lehrstuhl für Mathematik (MathCCES) A Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids CES Seminar Computational Engineering

More information

Comparison of maximum entropy and quadrature-based moment closures for shock transitions prediction in one-dimensional gaskinetic theory

Comparison of maximum entropy and quadrature-based moment closures for shock transitions prediction in one-dimensional gaskinetic theory Comparison of maximum entropy and quadrature-based moment closures for shock transitions prediction in one-dimensional gaskinetic theory Jérémie Laplante and Clinton P. T. Groth Citation: AIP Conference

More information

Stochastic Particle Methods for Rarefied Gases

Stochastic 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 information

Atmospheric Entry. Technology, Mathematical Model and Simulation

Atmospheric Entry. Technology, Mathematical Model and Simulation Atmospheric Entry Technology, Mathematical Model and Simulation Julian Köllermeier RWTH Aachen, August 26th 2016 Outline 1. Introduction to Atmospheric Reentry 2. Rarefied Gases: From Science Fiction to

More information

Fluid Equations for Rarefied Gases

Fluid Equations for Rarefied Gases 1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel

More information

Fluid Dynamics from Kinetic Equations

Fluid Dynamics from Kinetic Equations Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the

More information

Fluid Equations for Rarefied Gases

Fluid Equations for Rarefied Gases 1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel

More information

ON HYPERBOLICITY OF 13-MOMENT SYSTEM. Zhenning Cai. Yuwei Fan. Ruo Li. (Communicated by Yan Guo)

ON HYPERBOLICITY OF 13-MOMENT SYSTEM. Zhenning Cai. Yuwei Fan. Ruo Li. (Communicated by Yan Guo) Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 20XX doi:10.3934/xx.xx.xx.xx pp. X XX ON HYPERBOLICITY OF 13-MOMENT SYSTEM Zhenning Cai CAPT & School of Mathematical Sciences Peking University

More information

arxiv: v2 [math-ph] 29 Nov 2010

arxiv: v2 [math-ph] 29 Nov 2010 An Efficient NRxx Method for Boltzmann-BGK Equation Zhenning Cai, Ruo Li and Yanli Wang arxiv:1011.5789v2 [math-ph] 29 Nov 2010 November 6, 2018 Abstract In [3], we proposed a numerical regularized moment

More information

Kinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods

Kinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods Kinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering

More information

A hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models

A hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models A hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models Alessandro Alaia, Gabriella Puppo May 31, 2011 Abstract In this work we present a non stationary domain

More information

Monte Carlo methods for kinetic equations

Monte Carlo methods for kinetic equations Monte Carlo methods for kinetic equations Lecture 4: Hybrid methods and variance reduction Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/

More information

Connection Between the Lattice Boltzmann Equation and the Beam Scheme

Connection Between the Lattice Boltzmann Equation and the Beam Scheme NASA/CR-1999-209001 ICASE Report No. 99-10 Connection Between the Lattice Boltzmann Equation and the Beam Scheme Kun Xu Hong Kong University of Science and Technology, Kowloon, Hong Kong Li-Shi Luo ICASE,

More information

SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS

SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen

More information

The Regularized 13 Moment Equations for Rarefied and Vacuum Gas Flows

The Regularized 13 Moment Equations for Rarefied and Vacuum Gas Flows xxxx The Regularized 13 Moment Equations for Rarefied and Vacuum Gas Flows Henning Struchtrup Peyman Taheri Anirudh Rana University of Victoria, Canada Manuel Torrilhon RWTH Aachen Goal: Accurate and efficient

More information

Divergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF)

Divergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF) Divergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF) arxiv:1608.07869, arxiv:1711.01657, arxiv:1709.06644 Frankfurt University 1.February.2018 Preview

More information

Bulk equations and Knudsen layers for the regularized 13 moment equations

Bulk equations and Knudsen layers for the regularized 13 moment equations Continuum Mech. Thermon. 7 9: 77 89 DOI.7/s6-7-5- ORIGINAL ARTICLE Henning Struchtrup Toby Thatcher Bulk equations and Knudsen layers for the regularized 3 moment equations Received: 9 December 6 / Accepted:

More information

A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows

A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National

More information

The Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016

The Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016 The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical

More information

A Comparative Study Between Cubic and Ellipsoidal Fokker-Planck Kinetic Models

A Comparative Study Between Cubic and Ellipsoidal Fokker-Planck Kinetic Models A Comparative Study Between Cubic and Ellipsoidal Fokker-Planck Kinetic Models Eunji Jun German Aerospace Center (DLR), 3773 Göttingen, Germany M. Hossein Gorji Computational Mathematics and Simulation

More information

Generalized Gas Dynamic Equations

Generalized Gas Dynamic Equations 47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 2009 Orlando Florida AIAA 2009-672 Generalized Gas Dynamic Equations Kun Xu Mathematics Department

More information

Microscopically Implicit-Macroscopically Explicit schemes for the BGK equation

Microscopically Implicit-Macroscopically Explicit schemes for the BGK equation Microscopically Implicit-Macroscopically Explicit schemes for the BGK equation Sandra Pieraccini, Gabriella Puppo July 25, 2 Abstract In this work a new class of numerical methods for the BGK model of

More information

Numerical methods for kinetic equations

Numerical methods for kinetic equations Numerical methods for kinetic equations Lecture 6: fluid-kinetic coupling and hybrid methods Lorenzo Pareschi Department of Mathematics and Computer Science University of Ferrara, Italy http://www.lorenzopareschi.com

More information

arxiv: v1 [physics.comp-ph] 7 Mar 2018

arxiv: v1 [physics.comp-ph] 7 Mar 2018 A BGK model for high temperature rarefied gas flows C. Baranger 1, G. Marois 1, J. Mathé 1, J. Mathiaud 1, L. Mieussens arxiv:1803.0617v1 [physics.comp-ph] 7 Mar 018 1 CEA-CESTA 15 avenue des Sablières

More information

Kibble-Zurek dynamics and off-equilibrium scaling of critical cumulants in the QCD phase diagram

Kibble-Zurek dynamics and off-equilibrium scaling of critical cumulants in the QCD phase diagram Kibble-Zurek dynamics and off-equilibrium scaling of critical cumulants in the QCD phase diagram Raju Venugopalan BNL/Heidelberg arxiv:1310.1600 [cond-mat.stat-mech] Heidelberg Seminar, June 8 th, 2016

More information

Hydrodynamic Limits for the Boltzmann Equation

Hydrodynamic Limits for the Boltzmann Equation Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC

More information

Hyperbolic Moment Equations in Kinetic Gas Theory Based on Multi-Variate Pearson-IV-Distributions

Hyperbolic Moment Equations in Kinetic Gas Theory Based on Multi-Variate Pearson-IV-Distributions Commun. Comput. Phys. doi: 0.408/cicp.009.09.049 Vol. 7, No. 4, pp. 639-673 April 00 Hyperbolic Moment Equations in Kinetic Gas Theory Based on Multi-Variate Pearson-IV-Distributions Manuel Torrilhon Seminar

More information

Lattice Boltzmann approach to liquid - vapour separation

Lattice Boltzmann approach to liquid - vapour separation Lattice Boltzmann approach to liquid - vapour separation T.Biciușcă 1,2, A.Cristea 1, A.Lamura 3, G.Gonnella 4, V.Sofonea 1 1 Centre for Fundamental and Advanced Technical Research, Romanian Academy Bd.

More information

Stability and consistency of kinetic upwinding for advection diffusion equations

Stability and consistency of kinetic upwinding for advection diffusion equations IMA Journal of Numerical Analysis (006 6, 686 7 doi:0.093/imanum/drl005 Advance Access publication on March, 006 Stability and consistency of kinetic upwinding for advection diffusion equations MANUEL

More information

Non-equilibrium mixtures of gases: modeling and computation

Non-equilibrium mixtures of gases: modeling and computation Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present

More information

Modelling and numerical methods for the diffusion of impurities in a gas

Modelling 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 information

1 Continuum Models Some History and Background Added Value of Continuum Equations... 6

1 Continuum Models Some History and Background Added Value of Continuum Equations... 6 Re g u la riz a tio n o f Gra d s 1 3 -Mo m e nt-equ a tio n s in K in e tic Ga s Th e o ry Manuel Torrilhon Department of Mathematics & Center for Computational Engineering Science RWTH Aachen University

More information

Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases

Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases Heinrich Freistühler and Blake Temple Proceedings of the Royal Society-A May 2017 Culmination of a 15 year project: In this we propose:

More information

Numerical Simulation of the Rarefied Gas Flow through a Short Channel into a Vacuum

Numerical Simulation of the Rarefied Gas Flow through a Short Channel into a Vacuum Numerical Simulation of the Rarefied Gas Flow through a Short Channel into a Vacuum Oleg Sazhin Ural State University, Lenin av.5, 6283 Ekaterinburg, Russia E-mail: oleg.sazhin@uralmail.com Abstract. The

More information

Analytical and numerical solutions of the regularized 13 moment equations for rarefied gases

Analytical and numerical solutions of the regularized 13 moment equations for rarefied gases xxxx Analytical and numerical solutions of the regularized 13 moment equations for rarefied gases Henning Struchtrup University of Victoria, Canada Manuel Torrilhon ETH Zürich Peyman Taheri University

More information

Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion

Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 La Tremblade, Congrès SMAI 2017 5

More information

Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th

Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,

More information

Title of communication, titles not fitting in one line will break automatically

Title of communication, titles not fitting in one line will break automatically Title of communication titles not fitting in one line will break automatically First Author Second Author 2 Department University City Country 2 Other Institute City Country Abstract If you want to add

More information

An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling

An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3 Saint-Malo 13 December 2016 1 Université

More information

Chapter 1 Direct Modeling for Computational Fluid Dynamics

Chapter 1 Direct Modeling for Computational Fluid Dynamics Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the

More information

Anisotropic fluid dynamics. Thomas Schaefer, North Carolina State University

Anisotropic fluid dynamics. Thomas Schaefer, North Carolina State University Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding

More information

Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I

Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA

More information

Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model

Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model G. A. Radtke, N. G. Hadjiconstantinou and W. Wagner Massachusetts Institute of Technology,

More information

Multi species Lattice Boltzmann Models and Applications to Sustainable Energy Systems

Multi species Lattice Boltzmann Models and Applications to Sustainable Energy Systems Multi species Lattice Boltzmann Models and Applications to Sustainable Energy Systems Pietro Asinari, PhD Dipartimento di Energetica (DENER), Politecnico di Torino (POLITO), Torino 10129, Italy e-mail:

More information

Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles

Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational

More information

Towards realizable hyperbolic moment closures for viscous heat-conducting gas flows based on a maximum-entropy distribution

Towards realizable hyperbolic moment closures for viscous heat-conducting gas flows based on a maximum-entropy distribution Continuum Mech. Thermodyn. (23) 25:573 63 DOI.7/s6-2-252-y ORIGINAL ARTICLE James G. McDonald Clinton P. T. Groth Towards realizable hyperbolic moment closures for viscous heat-conducting gas flows based

More information

Lecture 5: Kinetic theory of fluids

Lecture 5: Kinetic theory of fluids Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities

More information

Nonrenormalizability of nonequilibrium quantum field theory in the classical approximation

Nonrenormalizability of nonequilibrium quantum field theory in the classical approximation Nonrenormalizability of nonequilibrium quantum field theory in the classical approximation Bin Wu IPhT, CEA/Saclay RPP 2015, Institut Henri Poincare Jan 16, 2015 T. Epelbaum, F. Gelis and B. Wu, Phys.

More information

Fluid dynamics for a vapor-gas mixture derived from kinetic theory

Fluid dynamics for a vapor-gas mixture derived from kinetic theory IPAM Workshop II The Boltzmann Equation: DiPerna-Lions Plus 20 Years (IPAM-UCLA, April 15-17, 2009) Fluid dynamics for a vapor-gas mixture derived from kinetic theory Kazuo Aoki Department of Mechanical

More information

Lattice Boltzmann Method for Moving Boundaries

Lattice Boltzmann Method for Moving Boundaries Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving

More information

Li Ruo. Adaptive Mesh Method, Numerical Method for Fluid Dynamics, Model Reduction of Kinetic Equation

Li Ruo. Adaptive Mesh Method, Numerical Method for Fluid Dynamics, Model Reduction of Kinetic Equation Li Ruo Contact Information Mailing Address: School of Mathematical Sciences, Peking University, Beijing, 100871, P. R. China Tel: 86-10-6276-7345(O), 86-150-1019-4872(M) E-mail: rli@math.pku.edu.cn Homepage:

More information

Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme

Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,

More information

An improved unified gas-kinetic scheme and the study of shock structures

An improved unified gas-kinetic scheme and the study of shock structures IMA Journal of Applied Mathematics (2011) 76, 698 711 doi:10.1093/imamat/hxr002 Advance Access publication on March 16, 2011 An improved unified gas-kinetic scheme and the study of shock structures KUN

More information

Statistical mechanics of random billiard systems

Statistical mechanics of random billiard systems Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore

More information

A hybrid kinetic-fluid model for solving the gas dynamics Boltzmann-BGK equation.

A hybrid kinetic-fluid model for solving the gas dynamics Boltzmann-BGK equation. A hybrid kinetic-fluid model for solving the gas dynamics Boltzmann-BGK equation. N. Crouseilles a, b, P. Degond a, and M. Lemou a a MIP, UMR CNRS 564, UFR MIG, Université Paul Sabatier 8, route de Narbonne,

More information

Exponential moments for the homogeneous Kac equation

Exponential moments for the homogeneous Kac equation Exponential moments for the homogeneous Kac equation Maja Tasković University of Pennsylvania Young Researchers Workshop: Stochastic and deterministic methods in kinetic theory, Duke University, November

More information

On the Boltzmann equation: global solutions in one spatial dimension

On the Boltzmann equation: global solutions in one spatial dimension On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators

More information

Traffic flow on networks: conservation laws models. Benedetto Piccoli I.A.C. C.N.R. Rome

Traffic flow on networks: conservation laws models. Benedetto Piccoli I.A.C. C.N.R. Rome Traffic flow on networks: conservation laws models Benedetto Piccoli I.A.C. C.N.R. Rome Conservation laws on networks u t + f(u) x=0 Dynamics at nodes? 1. The only conservation at nodes does not determine

More information

Symmetry of the linearized Boltzmann equation: Entropy production and Onsager-Casimir relation

Symmetry of the linearized Boltzmann equation: Entropy production and Onsager-Casimir relation Symmetry of the linearized Boltzmann equation: Entropy production and Onsager-Casimir relation Shigeru TAKATA ( 髙田滋 ) Department of Mechanical Engineering and Science, (also Advanced Research Institute

More information

Direct Modeling for Computational Fluid Dynamics

Direct Modeling for Computational Fluid Dynamics Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.

More information

Lattice Boltzmann Method

Lattice Boltzmann Method Chapter 5 Lattice Boltzmann Method 5.1 Pseudo-kinetic computational methods Understanding the flow properties of complex fluids in complex topologies is of importance to technological applications and

More information

Fourier Law and Non-Isothermal Boundary in the Boltzmann Theory

Fourier Law and Non-Isothermal Boundary in the Boltzmann Theory in the Boltzmann Theory Joint work with Raffaele Esposito, Yan Guo, Rossana Marra DPMMS, University of Cambridge ICERM November 8, 2011 Steady Boltzmann Equation Steady Boltzmann Equation v x F = Q(F,

More information

Hilbert Sixth Problem

Hilbert Sixth Problem Academia Sinica, Taiwan Stanford University Newton Institute, September 28, 2010 : Mathematical Treatment of the Axioms of Physics. The investigations on the foundations of geometry suggest the problem:

More information

DSMC simulations of thermal escape

DSMC simulations of thermal escape DSMC simulations of thermal escape Alexey N. Volkov, R.E. Johnson, O.J. Tucker, J.T. Erwin Materials Science & Engineering University of Virginia, USA Financial support is provided by NASA through Planetary

More information

Adiabatic piston - Draft for MEMPHYS

Adiabatic piston - Draft for MEMPHYS diabatic piston - Draft for MEMPHYS 4th ovember 2006 Christian GRUBER (a) and nnick LESE (b) (a) Institut de héorie des Phénomènes Physiques École Polytechnique Fédérale de Lausanne, CH-05 Lausanne, Switzerland

More information

The Boltzmann Equation and Its Applications

The Boltzmann Equation and Its Applications Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC

More information

Heat transfer in micro devices packaged in partial vacuum

Heat transfer in micro devices packaged in partial vacuum Heat transfer in micro devices packaged in partial vacuum AnirudhRana 1, Manuel Torrilhon 2, HenningStruchtrup 1 1 DepartmentofMechanicalEngineering, University of Victoria, Victoria BC, Canada 2 DepartmentofMathematics,

More information

Kinetic relaxation models for reacting gas mixtures

Kinetic 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 information

in Bounded Domains Ariane Trescases CMLA, ENS Cachan

in Bounded Domains Ariane Trescases CMLA, ENS Cachan CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline

More information

Predicting Breakdown of the Continuum Equations Under Rarefied Flow Conditions

Predicting Breakdown of the Continuum Equations Under Rarefied Flow Conditions Predicting Breakdown of the Continuum Equations Under Rarefied Flow Conditions Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109 Abstract. The breakdown of the

More information

Analysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium

Analysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium Analysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium Adriano W. da Silva May 30, 2013 Abstract A binary gaseous mixture with reversible reaction of type A + A = B

More information

arxiv: v1 [physics.flu-dyn] 14 Feb 2017

arxiv: v1 [physics.flu-dyn] 14 Feb 2017 This draft was prepared using the LaTeX style file belonging to the Journal of Fluid Mechanics 1 arxiv:170.0410v1 [physics.flu-dyn] 14 Feb 017 A comment on An improved macroscale model for gas slip flow

More information

Scaling Parameters in Rarefied Flow and the Breakdown of the Navier-Stokes Equations Mechanical Engineering Research Report No: 2004/09

Scaling Parameters in Rarefied Flow and the Breakdown of the Navier-Stokes Equations Mechanical Engineering Research Report No: 2004/09 Scaling Parameters in Rarefied Flow and the Breakdown of the Navier-Stokes Equations Mechanical Engineering Research Report No: 2004/09 Michael Macrossan, Centre for Hypersonics, University of Queensland

More information

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso

More information

Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture

Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture E.V. Kustova, E.A. Nagnibeda, Yu.D. Shevelev and N.G. Syzranova Department of Mathematics and Mechanics,

More information

arxiv: v1 [math.na] 14 Aug 2018

arxiv: v1 [math.na] 14 Aug 2018 An Efficient Steady-State Solver for icroflows with High-Order oment odel Zhicheng Hu, Guanghui Hu August 17, 218 arxiv:188.5472v1 [math.na] 14 Aug 218 Abstract In [Z. Hu, R. Li, and Z. Qiao. Acceleration

More information

Estimations of Rotational Relaxation Parameters in Diatomic Gases

Estimations of Rotational Relaxation Parameters in Diatomic Gases Estimations of Rotational Relaxation Parameters in Diatomic Gases Vladimir V. Riabov Department of Mathematics and Computer Science, Rivier College, 420 S. Main St., Nashua, NH 03060, USA Abstract: The

More information

Applications of the compensated compactness method on hyperbolic conservation systems

Applications of the compensated compactness method on hyperbolic conservation systems Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,

More information

Discrete-velocity models and numerical schemes for the Boltzmann-BGK equation in plane and axisymmetric geometries

Discrete-velocity models and numerical schemes for the Boltzmann-BGK equation in plane and axisymmetric geometries to appear in the Journal of Computational Physics Discrete-velocity models and numerical schemes for the Boltzmann-BGK equation in plane and axisymmetric geometries Luc Mieussens Mathématiques Appliquées

More information

The propagation of chaos for a rarefied gas of hard spheres

The propagation of chaos for a rarefied gas of hard spheres The propagation of chaos for a rarefied gas of hard spheres Ryan Denlinger 1 1 University of Texas at Austin 35th Annual Western States Mathematical Physics Meeting Caltech February 13, 2017 Ryan Denlinger

More information

LATTICE BOLTZMANN MODELLING OF PULSATILE FLOW USING MOMENT BOUNDARY CONDITIONS

LATTICE BOLTZMANN MODELLING OF PULSATILE FLOW USING MOMENT BOUNDARY CONDITIONS 6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 5 June 28, Glasgow, UK LATTICE BOLTZMANN MODELLING OF PULSATILE FLOW USING MOMENT

More information

12. MHD Approximation.

12. MHD Approximation. Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal

More information

TAU Extensions for High Enthalpy Flows. Sebastian Karl AS-RF

TAU Extensions for High Enthalpy Flows. Sebastian Karl AS-RF TAU Extensions for High Enthalpy Flows Sebastian Karl AS-RF Contents Motivation Extensions available in the current release: Numerical schemes for super- and hypersonic flow fields Models for gas mixtures,

More information

Why Should We Be Interested in Hydrodynamics?

Why Should We Be Interested in Hydrodynamics? Why Should We Be Interested in Hydrodynamics? Li-Shi Luo Department of Mathematics and Statistics Center for Computational Sciences Old Dominion University, Norfolk, Virginia 23529, USA Email: lluo@odu.edu

More information

A Very Brief Introduction to Conservation Laws

A Very Brief Introduction to Conservation Laws A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation

More information

Lecture 6 Gas Kinetic Theory and Boltzmann Equation

Lecture 6 Gas Kinetic Theory and Boltzmann Equation GIAN Course on Rarefied & Microscale Gases and Viscoelastic Fluids: a Unified Framework Lecture 6 Gas Kinetic Theory and Boltzmann Equation Feb. 23 rd ~ March 2 nd, 2017 R. S. Myong Gyeongsang National

More information

Computational Methods of Physical Problems with Mixed Scales

Computational Methods of Physical Problems with Mixed Scales Computational Methods of Physical Problems with Mixed Scales Shi Jin University of Wisconsin-Madison IPAM Tutorial, March 10-12, 2008 Scales and Physical Laws Figure from E & Engquist, AMS Notice Connections

More information

Application of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows

Application of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows Application of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows Timothy R. Deschenes and Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor,

More information

YAN GUO, JUHI JANG, AND NING JIANG

YAN GUO, JUHI JANG, AND NING JIANG LOCAL HILBERT EXPANSION FOR THE BOLTZMANN EQUATION YAN GUO, JUHI JANG, AND NING JIANG Abstract. We revisit the classical ork of Caflisch [C] for compressible Euler limit of the Boltzmann equation. By using

More information

Models of collective displacements: from microscopic to macroscopic description

Models of collective displacements: from microscopic to macroscopic description Models of collective displacements: from microscopic to macroscopic description Sébastien Motsch CSCAMM, University of Maryland joint work with : P. Degond, L. Navoret (IMT, Toulouse) SIAM Analysis of

More information

Rigid Body Motion in a Special Lorentz Gas

Rigid Body Motion in a Special Lorentz Gas Rigid Body Motion in a Special Lorentz Gas Kai Koike 1) Graduate School of Science and Technology, Keio University 2) RIKEN Center for Advanced Intelligence Project BU-Keio Workshop 2018 @Boston University,

More information

Unified Gas-kinetic Scheme for the Study of Non-equilibrium Flows

Unified Gas-kinetic Scheme for the Study of Non-equilibrium Flows Unified Gas-kinetic Scheme for the Study of Non-equilibrium Flows by WANG, Ruijie A Thesis Submitted to The Hong Kong University of Science and Technology in Partial Fulfillment of the Requirements for

More information