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

Size: px
Start display at page:

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

Transcription

1 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 Physics, Hong Kong January, University Paul Sabatier Toulouse III, IMT. 2 INRIA Rennes - Bretagne Atlantique & ENS Cachan Bretagne, IRMAR. 3 CNRS & University of Rennes I & ENS Cachan Bretagne, IRMAR. 1 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

2 Outline 1 Problem and objectives A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

3 Numerical simulation of plasmas Different scales in plasmas, for example collisions parameterized by the Knudsen number ε different models. Kinetic models Particles represented by a distribution function f (x, v, t). Solving the Vlasov equation (with source term S (ε)) t f + v x f + q m (E + v B) vf = S (ε) coupled to Maxwell equations or Poisson equation. In 3D 7 variables: 3 in space, 3 in velocity and the time heavy computations. 3 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

4 Fluid models Moment equations on physical quantities linked to f (density, mean velocity, temperature, etc.). Smaller cost, but lost of precision. General difficulties Find a well adapted model for our system, with a good precision/cost ratio. If two scales in the same simulation: two schemes with an interface? Difficult to deal with the interface! 4 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

5 Vlasov-Poisson-BGK system Collisions between particles parameterized by the Knudsen number ε. The Vlasov-Poisson-BGK equations (on [0, L x ] R R +, L x R + ) t f + v x f + E v f = 1 ε Q (f ), x E = 1 + f dv. Periodic conditions on f and E. Zero mean condition on E and initial conditions Lx 0 E (x, t) dx = 0, t 0, f (x, v, 0) = f 0 (x, v), x [0, L x ], v R. A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

6 Q (f ) is the BGK collisions operator (Bhatnagar-Gross-Krook): Q (f ) = M (U) f, M (U) the Maxwellian having the same first three moments than f, denoted by U. We denote by N (L Q ) = Span { M, vm, v 2 M }, the kernel of the linearized operator L Q of Q. 6 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

7 Objectives Problem and objectives Construction of a stable and consistent at the ε 0 limit scheme, of constant cost with respect to ε: an asymptotic preserving (AP) scheme 4. Reduction of the numerical cost: micro-macro decomposition with particles. Especially at the ε 0 limit since collisions bring f to its Maxwellian : f M (U) 0 when ε 0. 4 Jin, SIAM JSC A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

8 Micro-macro model Problem and objectives PIC method Projection step Macro part Micro-macro decomposition 56 f = M (U) + g: t M (U)+v x M (U)+E v M (U)+ t g +v x g +E v g = 1 ε g. Transport operator T = v x +E v : t M (U) + T M (U) + t g + T g = 1 ε g. Hypothesis: first 3 moments of g are zeros: mg := m (v) g (x, v, t) dv = 0, with m (v) := ) T (1, v, v 2. 2 Verified at the numerical level? If not, we have to impose it. 5 Lemou, Mieussens, SIAM JSC Liu, Yu, CMP A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

9 Micro-macro equations PIC method Projection step Macro part Let Π M the orthogonal projection in L 2 ( M 1 dv ) on N (L Q ): Π M (ϕ) = 1 ρ + [ (v u) (v u) ϕ ϕ + T ( v u 2 1 ) ( v u 2 2T 2 T ) ] 1 ϕ M. Properties: (I Π M ) ( t M) = Π M (g) = Π M ( t g) = 0. By applying (I Π M ) to Vlasov-BGK micro equation on g t g + (I Π M ) T (M + g) = 1 ε g. 9 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

10 PIC method Projection step Macro part By applying Π M to Vlasov-BGK macro equation on M (U) t M + Π M T (M + g) = 0. And by taking its first 3 moments t U + x F (U) + x vm (v) g = S (U), with F (U) the Euler flux and S (U) a source term. Micro-macro system 7 { t g + (I Π M ) T (M + g) = 1 ε g t U + x F (U) + x vm (v) g = S (U) equivalent to Vlasov-BGK. 7 Bennoune, Lemou, Mieussens, JCP A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

11 Algorithm Problem and objectives PIC method Projection step Macro part 1. Solving the micro part by a Particle-In-Cell (PIC) method. 2. Projection step to numerically force to zero the first three moments of g (matching procedure 8 ). 3. Solving the macro part by a finite volume scheme (mesh on x), with a source term dependent on g. 1-3 coupling: similarities with the δf method 9 but here: AP scheme. 8 Degond, Dimarco, Pareschi, IJNMF, Brunner, Valeo, Krommes, Phys. of Plasmas A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

12 1. PIC method on g Problem and objectives PIC method Projection step Macro part Solving by a PIC method the equation t g + (I Π M ) T (M + g) = 1 ε g t g + T g = (I Π M ) T M + Π M T g g ε =: S g, coupled to the Poisson equation for the electric field. Model: having N p numerical particles, with position x k, velocity v k and weight ω k, g is approximated by N p g Np (x, v, t) = ω k (t) δ (x x k (t)) δ (v v k (t)). k=1 2 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

13 PIC algorithm Problem and objectives PIC method Projection step Macro part Mesh generation Electric field computing on the mesh Initialization of positions and velocities of particles Computation of charge density on the mesh (deposition) Interpolation of electric field on the particles Evolution of the weights (if source term) Movement of particles 13 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

14 PIC method Projection step Macro part Equation with source term splitting. 1. Solving the transport part t g + T g = 0. Displacement of particles thanks to the equations of motion dx k dt (t) = v k (t), dv k dt (t) = E (x k (t), t). Second order Verlet scheme (for example): v n+ 1 2 k = v n k + t 2 E n (x n k ) x n+1 k = xk n + tv n+ 1 2 k v n+1 k = v n+ 1 2 k + t 2 E n+1 ( x n+1 ) k. 14 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

15 PIC method Projection step Macro part 2. Solving the source part t g = (I Π M ) (T M) + Π M (T g) 1 ε g. Equation on weights of the form dω k dt (t) = α k (t) ω k (t) ε where α k is the weight associated to (I Π M ) (T M) + Π M (T g). To have an AP scheme, we make the stiff term implicit: ω n+1 k = ω n k + t αn k ωn+1 t kε. 15 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

16 2. Projection step At time n + 1, we have N p k=1 Problem and objectives PIC method Projection step Macro part g n+1 (x, v) = ω n+1 k δ ( x x n+1 ) ( ) ( k δ v v n+1 k g x, v, t n+1 ). Micro-macro structure: mg ( x, v, t n+1) = 0. A priori not true for g n+1 (x, v) correction We compute U g := mg n+1 0 on each cell x i : U g (x i ) = mg n+1 xi = ω k m (v k ). k/x k [x i,x i+1 ] 2. We seek h N (L Q )= Span { M, vm, v 2 M } h (x, v) = λ (x) m (v) M (x, v) s.t. U g (x i ) = mh (x i, v). 10 C., Crouseilles, Lemou, KRM A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

17 PIC method Projection step Macro part Solving N x linear systems (3 3) U g (x i ) = A i λ i with A i a matrix containing moments of M and λ i λ (x i ). Computation of weights γ k associated to h on each particle. 3. Correction of the weights in order to preserve the micro-macro structure ω new k ω k γ k. By construction : mg n+1,new xi = U g (x i ) mh (x i, v) = 0. Remark: correction of order 1 (by using linear shape functions) (3N x 3N x ) system. 17 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

18 3. Macro part Problem and objectives PIC method Projection step Macro part Solving the equation t U + x F (U) = S (U) x vm (v) g =: S (U, g). Finite volume method U n+1 i = Ui n t x ( ) Fi+1/2 n F i 1/2 n t S i n, with Rusanov flux Fi+1/2 n = 1 ( ( ) F U n 2 i+1 + F (U n i ) a i+1/2 (U i+1 U i ) ), where a i+1/2 = max j=i,i+1 ( abs R ( JF ( xj ))), R (JF ) being the eigenvalues of the Jacobian of F. Computation of the moments of g to evaluate S n i = S (U n i ) ( vmg n+1 xi+1/2 vmg n+1 xi 1/2 x ). 18 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

19 Landau damping Problem and objectives Landau damping Computational cost Initial distribution function: f (x, v, 0) = 1 2π exp ( v 2 2 ) (1 + α cos (kx)). Micro-macro initializations: 1 + α cos (kx) U (x) = α cos (kx) and g (x, v, t = 0) = 0. Parameters: α = 0.01, k = 0.5, L x = 2π/k. Observation: electrical energy E (t) = E (t, x) 2 dx. 19 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

20 AP property Problem and objectives Landau damping Computational cost Convergence to Euler-Poisson at the ε 0 limit? E (t), comparison with Euler-Poisson. N x = 128, N p = , different values of ε log E(t) L Time t Good convergence AP scheme. Euler MiMa, eps= MiMa, eps=0.1 MiMa, eps=0.3 MiMa, eps=1 0 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

21 Importance of the projection step Landau damping Computational cost ρ (x), with and without correction. ε = 1, N p = at t = Without correction Correction of order 0 Correction of order Charge density 1-rho Position x Correction the instabilities disappear. 1 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

22 Numerical-noise reduction Landau damping Computational cost ρ (x), comparison with a PIC method on the whole f. ε = 1, N x = 128, N p = , at t = PIC BGK MiMa Charge density 1-rho Position x Micro-macro decomposition noise due to the PIC method reduced. 2 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

23 Computational cost Problem and objectives Landau damping Computational cost E (t), comparison with a PIC method on the whole f. ε = 10 ε = 10 4 N p Time MiMa s. PIC-BGK s. N p Time MiMa s. PIC-BGK s. 23 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

24 Conclusion Problem and objectives Landau damping Computational cost Micro-macro decomposition for Vlasov-Poisson-BGK using a PIC method. Projection step to numerically force the first 3 moments of g to zero. AP scheme. Reduction of the numerical noise due to the PIC method because only on g. Reduction of the cost compared to grid methods at the ε 0 limit because few particles are sufficient. Also true when g (x, v, 0) A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

25 Landau damping Computational cost - M. Benoune, M. Lemou, L. Mieussens: Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier-Stokes asymptotics, J. Comput. Phys. 227, pp (2008). - S. Brunner, E. Valeo, J.A. Krommes: Collisional delta-f scheme with evolving background for transport time scale simulations, Phys. of Plasmas 12, (1999). - A. C., N. Crouseilles, M. Lemou: Micro-macro decomposition for Vlasov-BGK equation using particles, Kinetic and Related Models 5, pp (2012). - N. Crouseilles, M. Lemou: An asymptotic preserving scheme based on a micro-macro decomposition for collisional Vlasov equations: diffusion and high-field scaling limits, Kinetic and Related Models 4, pp (2011). - P. Degond, G. Dimarco, L. Pareschi: The moment guided Monte Carlo method, International Journal for Numerical Methods in Fluids 67, pp (2011). - S. Jin: Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations, SIAM J. Sci. Comput. 21, pp (1999). - M. Lemou, L. Mieussens: A new asymptotic preserving scheme based on micro-macro formulation for linear kinetic equations in the diffusion limit, SIAM J. Sci. Comput. 31, pp (2008). - T.-P. Liu, S.-H. Yu: Boltzmann Equation: Micro-Macro Decompositions and Positivity of Shock Profiles, Comm. in Math. Phys. 246, pp (2004). - B. Yan, S. Jin: A successive penalty-based asymptotic-preserving scheme for kinetic equations, SIAM J. Sci. Comput., to appear. Thank you for your attention! 25 A. Crestetto, N. Crouseilles, M. Lemou Micro-macro scheme for Vlasov-Poisson-BGK

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. Workshop Asymptotic-Preserving schemes, Porquerolles.

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

Modèles hybrides préservant l asymptotique pour les plasmas

Modèles hybrides préservant l asymptotique pour les plasmas Modèles hybrides préservant l asymptotique pour les plasmas Anaïs Crestetto 1, Nicolas Crouseilles 2, Fabrice Deluzet 1, Mohammed Lemou 3, Jacek Narski 1 et Claudia Negulescu 1. Groupe de Travail Méthodes

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

Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling

Micro-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 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 Rennes, 14ème Journée de l équipe

More information

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

Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles

Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto, Nicolas Crouseilles, Mohammed Lemou To cite this version: Anaïs Crestetto, Nicolas Crouseilles,

More information

Hybrid and Moment Guided Monte Carlo Methods for Kinetic Equations

Hybrid 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 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

Uncertainty 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 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 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

Uncertainty 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. 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 information

Dissertation. presented to obtain the. Université Paul Sabatier Toulouse 3. Mention: Applied Mathematics. Luc MIEUSSENS

Dissertation. 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 information

Uncertainty Quantification for multiscale kinetic equations with high dimensional random inputs with sparse grids

Uncertainty Quantification for multiscale kinetic equations with high dimensional random inputs with sparse grids Uncertainty Quantification for multiscale kinetic equations with high dimensional random inputs with sparse grids Shi Jin University of Wisconsin-Madison, USA Kinetic equations Different Q Boltmann Landau

More information

High Order Semi-Lagrangian WENO scheme for Vlasov Equations

High Order Semi-Lagrangian WENO scheme for Vlasov Equations High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder

More information

A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas

A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University

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

Boundary Value Problems and Multiscale Coupling Methods for Kinetic Equations SCHEDULE

Boundary Value Problems and Multiscale Coupling Methods for Kinetic Equations SCHEDULE Boundary Value Problems and Multiscale Coupling Methods for Kinetic Equations April 21-24, 2016 Department of Mathematics University of Wisconsin-Madison SCHEDULE Thursday, April 21 Friday, April 22 Saturday,

More information

arxiv: v1 [math.na] 7 Nov 2018

arxiv: v1 [math.na] 7 Nov 2018 A NUMERICAL METHOD FOR COUPLING THE BGK MODEL AND EULER EQUATION THROUGH THE LINEARIZED KNUDSEN LAYER HONGXU CHEN, QIN LI, AND JIANFENG LU arxiv:8.34v [math.na] 7 Nov 8 Abstract. The Bhatnagar-Gross-Krook

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

Adaptive semi-lagrangian schemes for transport

Adaptive semi-lagrangian schemes for transport for transport (how to predict accurate grids?) Martin Campos Pinto CNRS & University of Strasbourg, France joint work Albert Cohen (Paris 6), Michel Mehrenberger and Eric Sonnendrücker (Strasbourg) MAMCDP

More information

Controlling numerical dissipation and time stepping in some multi-scale kinetic/fluid simulations

Controlling numerical dissipation and time stepping in some multi-scale kinetic/fluid simulations 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

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

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

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

Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equations

Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equations Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equations Irene M. Gamba, Shi Jin, and Liu Liu Abstract In this

More information

Semi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations

Semi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.

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

An Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation

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

STABLE STEADY STATES AND SELF-SIMILAR BLOW UP SOLUTIONS

STABLE STEADY STATES AND SELF-SIMILAR BLOW UP SOLUTIONS STABLE STEADY STATES AND SELF-SIMILAR BLOW UP SOLUTIONS FOR THE RELATIVISTIC GRAVITATIONAL VLASOV- POISSON SYSTEM Mohammed Lemou CNRS and IRMAR, Rennes Florian Méhats University of Rennes 1 and IRMAR Pierre

More information

Anomalous transport of particles in Plasma physics

Anomalous transport of particles in Plasma physics Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University

More information

Exponential methods for kinetic equations

Exponential methods for kinetic equations Exponential methods for kinetic equations Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara, Italy http://utenti.unife.it/lorenzo.pareschi/ lorenzo.pareschi@unife.it Joint research

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

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

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

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

Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck

Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN

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

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

Finite 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 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 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

Landau Damping Simulation Models

Landau Damping Simulation Models Landau Damping Simulation Models Hua-sheng XIE (u) huashengxie@gmail.com) Department of Physics, Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027, P.R.China Oct. 9, 2013

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

analysis for transport equations and applications

analysis for transport equations and applications Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg

More information

The Moment Guided Monte Carlo Method

The Moment Guided Monte Carlo Method The Moment Guided Monte Carlo Method Pierre Degond,, Giacomo Dimarco,,3 and Lorenzo Pareschi Université de Toulouse; UPS, INSA, UT, UTM ; Institut de Mathématiques de Toulouse ; F-3 Toulouse, France. CNRS;

More information

Lattice Bhatnagar Gross Krook model for the Lorenz attractor

Lattice Bhatnagar Gross Krook model for the Lorenz attractor Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,

More information

Exponential Runge-Kutta for inhomogeneous Boltzmann equations with high order of accuracy

Exponential Runge-Kutta for inhomogeneous Boltzmann equations with high order of accuracy Exponential Runge-Kutta for inhomogeneous Boltzmann equations with high order of accuracy Qin Li, Lorenzo Pareschi Abstract We consider the development of exponential methods for the robust time discretization

More information

Kinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence

Kinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence 2013 US-Japan JIFT workshop on New Aspects of Plasmas Kinetic Simulation NIFS, November 22-23, 2013 Kinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence cf. Revisit for Landau

More information

Hypocoercivity for kinetic equations with linear relaxation terms

Hypocoercivity for kinetic equations with linear relaxation terms Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT

More information

ASYMPTOTIC PRESERVING (AP) SCHEMES FOR MULTISCALE KINETIC AND HYPERBOLIC EQUATIONS: A REVIEW

ASYMPTOTIC PRESERVING (AP) SCHEMES FOR MULTISCALE KINETIC AND HYPERBOLIC EQUATIONS: A REVIEW ASYMPTOTIC PRESERVING (AP) SCHEMES FOR MULTISCALE KINETIC AND HYPERBOLIC EQUATIONS: A REVIEW SHI JIN Contents 1. Introduction 1 2. Hyperbolic systems with stiff relaxations 3 3. Kinetic equations: the

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

Accurate representation of velocity space using truncated Hermite expansions.

Accurate representation of velocity space using truncated Hermite expansions. Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli

More information

Moments conservation in adaptive Vlasov solver

Moments conservation in adaptive Vlasov solver 12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,

More information

arxiv: v1 [math.na] 25 Oct 2018

arxiv: v1 [math.na] 25 Oct 2018 Multi-scale control variate methods for uncertainty quantification in kinetic equations arxiv:80.0844v [math.na] 25 Oct 208 Giacomo Dimarco and Lorenzo Pareschi October 26, 208 Abstract Kinetic equations

More information

A successive penalty based Asymptotic-Preserving scheme for kinetic equations

A successive penalty based Asymptotic-Preserving scheme for kinetic equations A successive penalty based Asymptotic-Preserving scheme for kinetic equations Bokai Yan Shi Jin September 30, 202 Abstract We propose an asymptotic-preserving AP) scheme for kinetic equations that is efficient

More information

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

Overview of Accelerated Simulation Methods for Plasma Kinetics

Overview of Accelerated Simulation Methods for Plasma Kinetics Overview of Accelerated Simulation Methods for Plasma Kinetics R.E. Caflisch 1 In collaboration with: J.L. Cambier 2, B.I. Cohen 3, A.M. Dimits 3, L.F. Ricketson 1,4, M.S. Rosin 1,5, B. Yann 1 1 UCLA Math

More information

Lecture 4: The particle equations (1)

Lecture 4: The particle equations (1) Lecture 4: The particle equations (1) Presenter: Mark Eric Dieckmann Department of Science and Technology (ITN), Linköping University, Sweden July 17, 2014 Overview We have previously discussed the leapfrog

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

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

Computer simulation of multiscale problems

Computer simulation of multiscale problems Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20

More information

CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire

CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite

More information

A Discontinuous Galerkin Method for Vlasov Systems

A Discontinuous Galerkin Method for Vlasov Systems A Discontinuous Galerkin Method for Vlasov Systems P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu http://www.ph.utexas.edu/

More information

Stability of Mach Configuration

Stability of Mach Configuration Stability of Mach Configuration Suxing CHEN Fudan University sxchen@public8.sta.net.cn We prove the stability of Mach configuration, which occurs in moving shock reflection by obstacle or shock interaction

More information

MACROSCOPIC FLUID MODELS WITH LOCALIZED KINETIC UPSCALING EFFECTS

MACROSCOPIC 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 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

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

What place for mathematicians in plasma physics

What place for mathematicians in plasma physics What place for mathematicians in plasma physics Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 15-19 September 2008 Eric Sonnendrücker (U. Strasbourg) Math

More information

Implicit kinetic relaxation schemes. Application to the plasma physic

Implicit 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

Summer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations

Summer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations 1856-30 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The particle-in-cell M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The particle-in-cell

More information

Multi-water-bag model and method of moments for Vlasov

Multi-water-bag model and method of moments for Vlasov and method of moments for the Vlasov equation 1, Philippe Helluy 2, Nicolas Besse 3 FVCA 2011, Praha, June 6. 1 INRIA Nancy - Grand Est & University of Strasbourg - IRMA crestetto@math.unistra.fr (PhD

More information

Scalable Methods for Kinetic Equations

Scalable Methods for Kinetic Equations Scalable Methods for Kinetic Equations Oak Ridge National Laboratory October 19-23, 2015 Poster Titles and Abstracts 1 Multiscale and Multiphysics Simulations with GenASiS Reuben D. Budiardja University

More information

A class of asymptotic-preserving schemes for the Fokker-Planck-Landau equation

A class of asymptotic-preserving schemes for the Fokker-Planck-Landau equation A class of asymptotic-preserving schemes for the Fokker-Planck-Landau equation Shi Jin Bokai Yan October 2, 200 Abstract We present a class of asymptotic-preserving AP) schemes for the nonhomogeneous Fokker-Planck-Landau

More information

Quantum Hydrodynamics models derived from the entropy principle

Quantum Hydrodynamics models derived from the entropy principle 1 Quantum Hydrodynamics models derived from the entropy principle P. Degond (1), Ch. Ringhofer (2) (1) MIP, CNRS and Université Paul Sabatier, 118 route de Narbonne, 31062 Toulouse cedex, France degond@mip.ups-tlse.fr

More information

A HIERARCHY OF HYBRID NUMERICAL METHODS FOR MULTI-SCALE KINETIC EQUATIONS

A HIERARCHY OF HYBRID NUMERICAL METHODS FOR MULTI-SCALE KINETIC EQUATIONS A HIERARCHY OF HYBRID NUMERICAL METHODS FOR MULTI-SCALE KINETIC EQUATIONS FRANCIS FILBET AND THOMAS REY Abstract. In this paper, we construct a hierarchy of hybrid numerical methods for multi-scale kinetic

More information

On a class of implicit-explicit Runge-Kutta schemes for stiff kinetic equations preserving the Navier-Stokes limit

On a class of implicit-explicit Runge-Kutta schemes for stiff kinetic equations preserving the Navier-Stokes limit On a class of implicit-eplicit Runge-Kutta schemes for stiff kinetic equations preserving the Navier-Stokes limit Jingwei Hu Xiangiong Zhang June 8, 17 Abstract Implicit-eplicit (IMEX) Runge-Kutta (RK)

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

Michel Mehrenberger 1 & Eric Sonnendrücker 2 ECCOMAS 2016

Michel Mehrenberger 1 & Eric Sonnendrücker 2 ECCOMAS 2016 for Vlasov type for Vlasov type 1 2 ECCOMAS 2016 In collaboration with : Bedros Afeyan, Aurore Back, Fernando Casas, Nicolas Crouseilles, Adila Dodhy, Erwan Faou, Yaman Güclü, Adnane Hamiaz, Guillaume

More information

Numerical methods for plasma physics in collisional regimes

Numerical methods for plasma physics in collisional regimes J. Plasma Physics (15), vol. 81, 358116 c Cambridge University Press 1 doi:1.117/s3778176 1 Numerical methods for plasma physics in collisional regimes G. Dimarco 1,Q.Li,L.Pareschi 1 and B. Yan 3 1 Department

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

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

Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas

Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Michael Barnes Plasma Science & Fusion Center Massachusetts Institute of Technology Collaborators:

More information

Nhung Pham 1, Philippe Helluy 2 and Laurent Navoret 3

Nhung Pham 1, Philippe Helluy 2 and Laurent Navoret 3 ESAIM: PROCEEDINGS AND SURVEYS, September 014, Vol. 45, p. 379-389 J.-S. Dhersin, Editor HYPERBOLIC APPROXIMATION OF THE FOURIER TRANSFORMED VLASOV EQUATION Nhung Pham 1, Philippe Helluy and Laurent Navoret

More information

Postprint.

Postprint. http://www.diva-portal.org Postprint This is the accepted version of a chapter published in Domain Decomposition Methods in Science and Engineering XXI. Citation for the original published chapter: Gander,

More information

From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray

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

Research of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method

Research of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,

More information

Recent Innovative Numerical Methods

Recent Innovative Numerical Methods Recent Innovative Numerical Methods Russel Caflisch IPAM Mathematics Department, UCLA 1 Overview Multi Level Monte Carlo (MLMC) Hybrid Direct Simulation Monte Carlo (DSMC) Accelerated Molecular Dynamics

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

Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics

Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics Luigi Barletti (Università di Firenze) Carlo Cintolesi (Università di Trieste) 6th MMKT Porto Ercole, june 9th

More information

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

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction

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

Week 6 Notes, Math 865, Tanveer

Week 6 Notes, Math 865, Tanveer Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,

More information

Benchmarks in Computational Plasma Physics

Benchmarks in Computational Plasma Physics Benchmarks in Computational Plasma Physics P. Londrillo INAF, Bologna, Italie S. Landi Università di Firenze, Italie What you compute when you do computations of the Vlasov equation? Overview A short review

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

arxiv: v2 [math.na] 5 Jul 2016

arxiv: v2 [math.na] 5 Jul 2016 AN ASYMPTOTIC PRESERVING SCHEME FOR THE RELATIVISTIC VLASOV MAXWELL EQUATIONS IN THE CLASSICAL LIMIT by arxiv:1602.09062v2 [math.na] 5 Jul 2016 Nicolas Crouseilles, Lukas Einkemmer & Erwan Faou Abstract.

More information

A cut cell method for the 3D simulation of Crookes radiometer

A cut cell method for the 3D simulation of Crookes radiometer A cut cell method for the D simulation of Crookes radiometer Guillaume Dechristé and Luc Mieussens Univ. Bordeaux, IMB, UMR 5251, F-400 Talence, France. CNRS, IMB, UMR 5251, F-400 Talence, France. Abstract.

More information

Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas

Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas Vladimir Kolobov, a,b,1 Robert Arslanbekov a and Dmitry Levko a a CFD Research Corporation, Huntsville, AL 35806, USA b The

More information

Uniformly accurate averaging numerical schemes for oscillatory evolution equations

Uniformly accurate averaging numerical schemes for oscillatory evolution equations Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University

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