Controlling numerical dissipation and time stepping in some multi-scale kinetic/fluid simulations
|
|
- Melinda Heath
- 5 years ago
- Views:
Transcription
1 Controlling numerical dissipation and time stepping in some multi-scale kinetic/fluid simulations Jian-Guo Liu Department of Physics and Department of Mathematics, Duke University Collaborators: Pierre Degond, Jeff Haack, Shi Jin, Luc Mieussens
2 Outline of talk Two numerical difficulties in multi-scale fluids/kinetic simulations Stiffness in time step due to fast waves In the Eulerian description, numerical dissipation is amplified due to the dynamics in small scale motions such as fast waves and small mean free path. I will focus on three examples All speed scheme for compressible flows (with Jaff Haack and Shi Jin) Effective local kinetic upscaling (with Pierre Degond and Luc Mieussens) Linear kinetic equation in diffusion region (with Luc Mieussens) In connection with neutron transport equation, I will also show some PDE estimates on short path statics (with Pierre Degond)
3 Asymptotic Preserving (AP) schemes definition: a scheme uniformly stable and accurate w.r.t ε ( t and x independent of ε) consequence: in the fluid regime (ε 1): scheme consistent with the fluid equation with only slight implicit time discretization to ensure stability and gain efficiency (remove fast variables) numerical dissipation is controlled in the leading variables (remove 1/ɛ factor in numerical dissipation, or appear only in the space orthogonal to the leading variables) preserving important conservation laws key: decomposition equations
4 Euler equation for compressible flows, 1757 Leonhard Euler ρ t + (ρu) = 0 conservation of mass (ρu) t + (ρu u) + p = f conservation of momentum ρ density, u velocity, f long range force Pre-thermodynamics era (1st and 2nd law of thermal dynamics not known yet), equation of state not well defined When Mach number is not small, conservation of energy becomes significant (Ma 2 = ρ 0U 2 p 0 = 2 kinetic energy γ 1 internal energy ). Scaling invariant under t Tt, x Lx, u Uu, L = UT Universal equation: from system as small as cold Fermi gas ( 100 microns) to galaxy.
5 Euler equation for compressible flows, 1757 Leonhard Euler ρ t + (ρu) = 0 conservation of mass (ρu) t + (ρu u) + p = f conservation of momentum ρ density, u velocity, f long range force Pre-thermodynamics era (1st and 2nd law of thermal dynamics not known yet), equation of state not well defined When Mach number is not small, conservation of energy becomes significant (Ma 2 = ρ 0U 2 p 0 = 2 kinetic energy γ 1 internal energy ). Scaling invariant under t Tt, x Lx, u Uu, L = UT Universal equation: from system as small as cold Fermi gas ( 100 microns) to galaxy.
6 Euler eq for incompressible flows, 1752, 1761 Leonhard Euler Dimensionless rescaling of compressible Euler eq: ρ t + (ρu) = 0 (ρu) t + (ρu u) + 1 Ma 2 p = f Mach number: Ma = U/c 0, c 0 = p0 ρ 0, ρ 0 and p 0 : dimensional scaling factors for pressure and density. As Ma 0, we eliminate the unimportant acoustic waves asymptotic expansion ρ = ρ 0 + Ma 2 ρ +, p = p 0 + Ma 2 p + ρ 0 ( u t u = 0 + u u) + p = f incompressibility momentum eq, Newton s 2nd law
7 1827, Navier-Stokes equations, 1845 C.-L. Navier Sir G.G. Stokes The dimensionless isentropic compressible Navier-Stokes equations are t ρ + m = 0 t m + (ρu u + 1ɛ ) 2 p(ρ) = (µσ), σ = u + u T 2 d u The incompressible Navier-Stokes equations are ρ( u t + u u) + p = µ u + f in Ω momentum u = 0 in Ω incompressibility The only dimensionless quantity: Reynolds number, Re = ρul/µ.
8 Hyperbolic decomposition t U + F (U) + QU = 0: Sub-systems U t + F (U) = 0, U t + QU = 0 are hyperbolic. i.e., they are dynamically stable sub-systems. F (U) is nonlinear, but it is non-stiff, any stable hyperbolic solver shall work, there is no excessive numerical dissipation. QU is stiff, but it is a linear system with constant coefficients! High order discretization with implicit non-dissipative scheme is used to stiffness in time stepping and introduce NO numerical dissipation. Automatically under-resolve small-scale details when it has no effect on bulk behavior. Asymptotic preserving: physical asymptotic behavior realized in numerical discretization, enable to do all speed computation.
9 Semi-implicit Scheme: Un+1 U n t + AU n + BU n+1 = 0 If both schemes U n+1 U n + AU n = 0, t U n+1 U n + BU n+1 = 0 t are stable, then the original scheme is also stable. Scheme U n+1 = (I ta)u n is stable iff I ta 1 + c t Scheme U n+1 = (I + tb) 1 U n is stable iff (I + tb) c t Combined scheme U n+1 = (I + tb) 1 (I ta)u n is stable iff (I + tb) 1 (I ta) 1 + C t
10 Hyperbolic decomposition t U + x F (U) + Q x U = 0: U = ( ρ m t ρ + α x m + (1 α) x m = 0 (1) ( t m + x ρu 2 + p ɛ (ρ) ) + a(t) ɛ 2 xρ = 0 (2) p(ρ) a(t)ρ p ɛ (ρ) = ɛ 2, a(t) = min p (ρ(x, t)) x ) ( ) ( 0 1 α, F (U) =, Q = Jacobian and Eigenvaues ( F 0 α (U) = p ɛ(ρ) u 2 2u αm ρu 2 + p ɛ (ρ) a(t) ɛ 2 0 ) x ), λ ± = u ± (1 α)u 2 + α p ɛ(ρ) p ɛ(ρ) 0 implies hyperbolicity, p ɛ(ρ) is bounded implies no stiffness.
11 A stable scheme Denote D u = (u j u j 1 )/h, D + u = (u j+1 u j )/h and ( ) 0 (1 α)d Q = a n D ɛ (3) F ± = 1 2 ( F (U) ± λ(u) U), { } λ(u) = max u ± (1 α)u 2 + α p(ρ) U n+1 U n + D F + (U n ) + D + F (U n ) + QU n+1 = 0 t This has a natural extension to a second order MUSCL scheme.
12 A stable scheme U n+1 U n + D F + (U n ) + D + F (U n ) + QU n+1 = 0 t ( ) 0 (1 α)d Q = a n D ɛ We can reformulate the density equation as (4) ρ n+1 2ρ n + ρ n 1 t 2 (1 α)an ɛ 2 D D + ρ n+1 = Ex(U n, U n 1 ) We can use a Fast Helmholtz Solver! After finding ρ n+1, update m n+1 explicitly.
13 Asymptotic preserving: discrete incompressible limit The incompressible limit scheme becomes u n+1 u n t α h u n + (1 α) h un+1 = 0 + h (u n u n ) + 1 ρ c + h pn+1 2 = µ h u n Note: The incompressible pressure 1 ρ c p 2 is obtained by p 2 = ρ c p (ρ c ) ρ.
14 1D Riemann Problem, Mach = 1, α =.25 %,-./012131$"!!%)(!-!!# $"' $ #"' #!!"#!"$!"%!"&!"'!"(!")!"*!"+ # #"' $"!!%)(!-!!# #!"'!!!"#!"$!"%!"&!"'!"(!")!"*!"+ #
15 2 D Riemann Problem, Mach = 1, α =.75! =
16 Driven Cavity flow, Re=4000, Mach =.1 and eps =.1, T=20, Re=4000, alpha =.01, beta = 1 #,-./0/#!!& 1/2/0/$!1/3,/0/&!!!1/45-64/0/,-. $ 1/7,84/0/# 0.9!"+ 0.8!"* 0.7!") 0.6!"( 0.5!"' 0.4!"& 0.3!"% 0.2!"$ 0.1!"# Top velocity U = 16x 2 (1 x 2 )!!!"#!"$!"%!"&!"'!"(!")!"*!"+ #
17 Backward Facing Step Flow, Re=100, Mach = AP Scheme captures correct reattachment point for low Reynolds number
18 Efficient localized kinetic upscaling [Degond-Liu-Mieussens, MMS 06] the micro-macro model (for Boltzmann): t U + x F (U) + x vmg = 0, t g + (I Π M(U) )(v x g) = Lg + Q(g, g) the localized micro-macro model t U + x F (U) + x vmg K = 0, (I Π M(U) )(v x M(U)) t g K + h(i Π M(U) )(v x g K ) = hlg K + hq(g K, g K ) h(i Π M(U) )(v x M(U)) simple idea: localization of g use the cutoff function h: determine Fluid zones : set h = 0 determine Kinetic zones : set h = 1 define buffer zones between (F) and (K) zones: set 0 h 1
19 Numerical tests: hydrodynamical scale [Degond-Liu-Mieussens, MMS 06] Sod test (density and h function) for the BGK equation with dynamic localization. 6 x 10 6 Density Sod Test for t=0.018s 1.2 Interface position, Local Knudsen number and equilibrium parameter for t=0.018 s transition function Knudsen x 10 Q density(kg/m 3 ) x h, Kn x 10, Q kinetic model macroscopic model micmac x(m) x(m)
20 Numerical tests: diffusion scale [Degond-Liu-Mieussens, MMS 06] Temperature for the radiative-heat transfer model, static localization. 1 Kinetic!=1 Fluid!= buffer kinetic model fluid model (diff micro-macro mo temperature x
21 Neutron transport equation in diffusion time scale ε t f + v f = 1 ε Lf εσ af + εs, Lf (v) = s(v, v )(f (v ) f (v))dv S 2 where v S 2 angular direction, σ total cross section, σ a absorption cross section, S internal source
22 Micro-macro decomposition. slab geometry t f + 1 ε v xf = σ ε 2 Lf σ af + S v [ 1, 1] is the x cosine direction. micro-macro decomposition: f = ρ + εg, ρ = [f ] = f (v) dv, [g] = 0 decomposition of the kinetic equation: t ρ + x [vg] = σ a ρ + S, t g + 1 ε (I [.])(v xg) = σ ε 2 Lg 1 ε 2 v xρ this formulation is well suited to derive the diffusion limit
23 Discretization and diffusion limit. the scheme: ρ n+1 i ρ n i t g n+1 i+ 2 1 g n i+ 1 2 t 2 + 4v σ i+ 1 = 2 g n+1 i+ 1 2 g n+1 i 1 2 x + 1 (I [.]) ε ε 2 Lg n+1 i = σ ai ρ n+1 i + S i, 0 g + i+ 1 2 g n i 1 2 x 1 ε 2 v ρn i+1 ρn i x g n + v i+ 2 3 g n 1 i+ 1 2 A x limit ε = 0: the scheme gives ρ n+1 i ρ n i t κ ρn i+1 2ρn i + ρ n i 1 x 2 = σ ai ρ n+1 i + S i, time explicit scheme with 3-point stencil
24 Analysis (Mieussens + L) the scheme is uniformly stable: under the CFL constraint: ρ n + ε g n C ( ρ 0 + ε g 0 ) t 2 3 ( 1 x 2 ) 3 2κ + ε x small ε: t 2 x 2 9 2κ (CFL for diffusion) large ε: t 2 3ε x (CFL for convection) the scheme is uniformly accurate: ρ(t n ) ρ n + ε g(t n ) g n C( t + x 2 + ε x)
25 Mean chord length for convex domain l = [ ] S dv 2 Γ,v n 0 l(s, v)v nds [ ] S dv Ω = 2 Γ ds S 2,v n 0 v ndv Γ dsπ = 4 Ω Γ Uniform and isotropic flux of particles entering the domain Ω across the surface Γ (neutron transport, Dirac, Fuchs, 1943) Integral geometry, geometric probability, Cauchy 1841, Czuber 1884)
26 Related references on mean chord length A.-L. Cauchy Paul Dirac Fuchs , Cauchy, mean projected area S = Γ /4 1884, Czuber, mean chord length l = Ω / S = 4 Ω / Γ 1943, Dirac, Fuchs etc, Declassified British Report MS-D-5 Atomic bomb: fissile material goes super-critical only depending on the mean chord of the piece of uranium (not mass as normally thought). Probability of reaction = l σ, σ: cross section 1958, A.M. Weinberg, E.P. Wigner, The Physical Theory of Neutron Chain Reactors. 1981, J N Bardsley and A Dubi, SIAM Appl. Math. 2003, S. Blanco and R. Fournier, Europhys. Lett.
27 orresponding to the first way of reasoning. that corresponding to the first way of r viously, these two images are both meaningful, but the corresponding processesobviously, are these two images are bo aneously under way and compensate each other exactly. The first image concerns simultaneously long under way and compens Heuristic arguments, Blanco-Fournier, 2003 A! << R R(m) R! 1 ~ R For small mean free path: B A Fig. 3 B2! 2!1 Monte Carlo simulations. Simulated average trajectory length inside a circle of radius Fig. R 2 for Monte Carlo simulations. Simulat sects with distinct random-walk characteristics: mean free path λ = m (measured two insects with distinct random-walk cha ant species Messor sancta L.); mean free path λ = m (roughly estimated for for the the ant species Messor sancta L.); m frequent direction changes ach species Blattella germanica L.). In both simulations, the asymmetry parameter cockroach of the species Blattella germanica L.) scattering phase function is g = 0.5. single scattering phase function is g = 0.5. Trajectory homothety. Illustration of the quasi-homothety of short trajectories. The Fig. second 3 Trajectory homothety. Illustration ory corresponds to a mean free path half as long as that of the first trajectory (λ 2 = trajectory λ corresponds to a mean free pat 1/2), for same statistical choices, and neglecting circle curvature, leads to a trajectory length which half as for same statistical choices, and negl For any isotropic uniform incidence, the average 2 L long (L 1/2). 2 length L 1/2). of B1 <L> R(m) Fig. 2 Once ant entered domain, it travels a longer distance due to On the other hand, most ants have an exit point very close to their entry point trajectories through the domain is independent of the characteristics of the diffusion process
28 Heuristic arguments Blanco-Fournier 2003, Bardsley-Dubi 1981 Steady neutron transport equation v f = 1 s(v, v )(f (v ) f (v))dv ε S 2 with BC: f = φ 0 for x Γ, v n 0. Constant solution (statistical equilibrium): f (x, v) φ 0. Rate at which reactions occur R = σ f (x, v)dxdv = σφ 0 4π Ω Ω S 2 R = σ l in-flux = σ l φ 0 Γ S 2,v n 0 v ndsdv = σ l φ 0 π Γ l = 4 Ω / Γ for convex domain, any mean free path
29 A. Mazzolo: Properties of diffusive random walks in bounded dom Heuristic arguments Mazzolo 2004, Blanco-Fournier l n α n τ n+1 for small τ n l r n 1 = l n, n path lengths l of first exit trajectories starting uniformly and isotropically at the boundary Fig. 1 Example of two-dimensional random walks generated with a mean free path a non-convex body (non-concentric hollow disk): solid lines represent random walk domain, while dotted lines represent random walks starting on the surface. path lengths r of first exit trajectories starting uniformly and isotropically within domain Random-walks properties. In the following, we consider constant-speed whose probabilistic law for the independent random jumps has an exponentia diffusion approximation t ρ = τ 3 ρ with BC ρ + 2τ p(l) =1/λ exp[ l/λ]. 3 nρ = 0 and uniform initial data ρ t=0 = ρ 0. where λ is the mean free path. In reactor physics, this corresponds to the monocinetic neutrons in a purely diffusing medium. First, we select a point tributed inside K, a bounded body in R 3 of volume V and surface S. Then allowed to move according to eq. (1) in a direction that is uniformly distrib
30 Some estimates (Degond + L) Theorem u τ = 1, in Ω, u τ + τ n u τ = 0, on Γ Ω 2 Γ τ + u0 2 u τ Ω 2 Ω Γ τ + u 2 where u 0 and v 0 are solutions of u 0 = 1, in Ω, u 0 = 0, on Γ u = 1, in Ω, n u = Ω / Γ, on Γ, Proposition inf u H 1, R H(u) inf H(u), Γ u=0 u H0 1 H(u) The equality holds if and only if Ω is a sphere. Γ u = 0 ( ) 1 2 u 2 u Ω
31 Expansion for τ K 0 u τ = Ω Γ τ + u + τ k v k k=1 u = 1, in Ω, n u = Ω / Γ, on Γ, v k = 0, in Ω, n v k = v k 1, on Γ, Γ Γ u = 0 v k = 0 Convergence analysis: v k Γ C (I n n) v k Γ K 0 n v k Γ = K 0 v k 1 Γ v k 2 Ω = Γ (u τ u ) v k n v k K 0 v k 1 2 Γ K0 τ K 0 u Γ
32 A technique lemma Lemma d dτ Ω u τ = Ω 2 / Γ + n u τ + Ω / Γ (u τ u 0 ) 2 + τ n u τ 2 2 Γ τ n u Γ 0 u n u τ + Ω / Γ 2 = u 2 Γ 0 2 Γ
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 informationHypocoercivity 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 informationMicro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling
Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2, Giacomo Dimarco 3 et Mohammed Lemou 4 Saint-Malo, 14 décembre 2017 1 Université de
More informationUne 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 informationAn 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 informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationKinetic/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. Workshop Asymptotic-Preserving schemes, Porquerolles.
More informationUne 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 Rennes, 14ème Journée de l équipe
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationA 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 informationComputational Astrophysics 7 Hydrodynamics with source terms
Computational Astrophysics 7 Hydrodynamics with source terms Oscar Agertz Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term
More informationModels 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 informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationEquivalence 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 informationThe 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 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 informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationFluid 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 informationPhysical 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 informationFluid 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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationChapter 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 informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationNumerical 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 informationAsymptotic-Preserving scheme based on a Finite Volume/Particle-In-Cell coupling for Boltzmann- BGK-like equations in the diffusion scaling
Asymptotic-Preserving scheme based on a Finite Volume/Particle-In-Cell coupling for Boltzmann- BGK-like equations in the diffusion scaling Anaïs Crestetto, Nicolas Crouseilles, Mohammed Lemou To cite this
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 informationDissertation. presented to obtain the. Université Paul Sabatier Toulouse 3. Mention: Applied Mathematics. Luc MIEUSSENS
Dissertation presented to obtain the HABILITATION À DIRIGER DES RECHERCHES Université Paul Sabatier Toulouse 3 Mention: Applied Mathematics by Luc MIEUSSENS Contributions to the numerical simulation in
More informationin 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 informationUncertainty Quantification for multiscale kinetic equations with random inputs. Shi Jin. University of Wisconsin-Madison, USA
Uncertainty Quantification for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Where do kinetic equations sit in physics Kinetic equations with applications
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationFluid 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 informationMicroscopically 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ρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationSECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS. Min Tang. (Communicated by the associate editor name)
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: http://aimsciences.org pp. X XX SECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS Min Tang Department of Mathematics
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationA 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 informationKinetic 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 informationA CLASS OF ASYMPTOTIC PRESERVING SCHEMES FOR KINETIC EQUATIONS AND RELATED PROBLEMS WITH STIFF SOURCES
A CLASS OF ASYMPTOTIC PRESERVING SCHEMES FOR KINETIC EQUATIONS AND RELATED PROBLEMS WITH STIFF SOURCES FRANCIS FILBET AND SHI JIN Abstract. In this paper, we propose a general framework to design asymptotic
More informationFinite volumes schemes preserving the low Mach number limit for the Euler system
Finite volumes schemes preserving the low Mach number limit for the Euler system M.-H. Vignal Low Velocity Flows, Paris, Nov. 205 Giacomo Dimarco, Univ. de Ferrara, Italie Raphael Loubere, IMT, CNRS, France
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationAn 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 informationFluid-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 informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationMagnetized Target Fusion: Insights from Mathematical Modelling
: Insights from Mathematical Modelling, UBC, PhD Candidate - Applied Math Brian Wetton, UBC, Supervisor Michael Ward, UBC, Committee Member Sandra Barsky, General Fusion, Committee Member November 28,
More informationDirect 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 informationApplications 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 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 informationUncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs.
Uncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Shanghai Jiao Tong University,
More informationComparison 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 informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationFourier 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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationHydrodynamic 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 informationHybrid and Moment Guided Monte Carlo Methods for Kinetic Equations
Hybrid and Moment Guided Monte Carlo Methods for Kinetic Equations Giacomo Dimarco Institut des Mathématiques de Toulouse Université de Toulouse France http://perso.math.univ-toulouse.fr/dimarco giacomo.dimarco@math.univ-toulouse.fr
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationOn the Asymptotic Preserving property of the Unified Gas Kinetic Scheme for the diffusion limit of linear kinetic models
On the Asymptotic Preserving property of the Unified Gas Kinetic Scheme for the diffusion limit of linear kinetic models Luc Mieussens To cite this version: Luc Mieussens. On the Asymptotic Preserving
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 informationYAN 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 informationAn asymptotic preserving scheme in the drift limit for the Euler-Lorentz system. Stéphane Brull, Pierre Degond, Fabrice Deluzet, Marie-Hélène Vignal
1 An asymptotic preserving scheme in the drift limit for the Euler-Lorentz system. Stéphane Brull, Pierre Degond, Fabrice Deluzet, Marie-Hélène Vignal IMT: Institut de Mathématiques de Toulouse 1. Introduction.
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationAn Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation
An Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation Wei Ren a, Hong Liu a,, Shi Jin b,c a J C Wu Center for Aerodynamics, School of Aeronautics and Aerospace, Shanghai Jiao Tong University,
More informationLOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 3, pp. 1 8. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S
More informationA model of alignment interaction for oriented particles with phase transition
A model of alignment interaction for oriented particles with phase transition Amic Frouvelle ACMAC Joint work with Jian-Guo Liu (Duke University, USA) and Pierre Degond (Institut de Mathématiques de Toulouse,
More informationarxiv: v1 [math.ap] 28 Apr 2009
ACOUSTIC LIMIT OF THE BOLTZMANN EQUATION: CLASSICAL SOLUTIONS JUHI JANG AND NING JIANG arxiv:0904.4459v [math.ap] 28 Apr 2009 Abstract. We study the acoustic limit from the Boltzmann equation in the framework
More informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description of many particle
More informationEffective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body
Effective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body I.G. Brykina a, B.V. Rogov b, I.L. Semenov c, and G.A. Tirskiy a a Institute of Mechanics, Moscow
More informationAsymptotic-Preserving Particle-In-Cell method for the Vlasov-Poisson system near quasineutrality
Asymptotic-Preserving Particle-In-Cell method for the Vlasov-Poisson system near quasineutrality Pierre Degond 1, Fabrice Deluzet 2, Laurent Navoret 3, An-Bang Sun 4, Marie-Hélène Vignal 5 1 Université
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationDSMC 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 informationNon-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 informationMonte 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 informationA model for a network of conveyor belts with discontinuous speed and capacity
A model for a network of conveyor belts with discontinuous speed and capacity Adriano FESTA Seminario di Modellistica differenziale Numerica - 6.03.2018 work in collaboration with M. Pfirsching, S. Goettlich
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationLecture 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 informationFluid 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 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 informationAnisotropic 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 informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows Kun Xu a,, Juan-Chen Huang b a Mathematics Department Hong Kong University of Science and Technology Clear Water Bay, Kowloon, Hong Kong b
More informationMACROSCOPIC FLUID MODELS WITH LOCALIZED KINETIC UPSCALING EFFECTS
MACROSCOPIC FLUID MODELS WITH LOCALIZED KINETIC UPSCALING EFFECTS Pierre Degond, Jian-Guo Liu 2, Luc Mieussens Abstract. This paper presents a general methodology to design macroscopic fluid models that
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationFundamentals of Acoustics
Fundamentals of Acoustics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
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 informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationA Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More informationEulerian interface-sharpening methods for hyperbolic problems
Eulerian interface-sharpening methods for hyperbolic problems Application to compressible multiphase flow Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan 11:55-12:25, March 05,
More informationImplicit kinetic relaxation schemes. Application to the plasma physic
Implicit kinetic relaxation schemes. Application to the plasma physic D. Coulette 5, E. Franck 12, P. Helluy 12, C. Courtes 2, L. Navoret 2, L. Mendoza 2, F. Drui 2 ABPDE II, Lille, August 2018 1 Inria
More information