Multi-water-bag model and method of moments for Vlasov
|
|
- Patricia Deborah Miller
- 5 years ago
- Views:
Transcription
1 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 thesis supported by DGA) 2 University of Strasbourg - IRMA 3 University Henri Poincaré Nancy I - IECN & IJL and method of moments for Vlasov
2 Summary: 1 Introduction and method of moments for Vlasov
3 Vlasov-Poisson 1D system: t f (x, v, t) + v x f (x, v, t) + q m E (x, t) v f (x, v, t) = 0 x E (x, t) = q m where n 0 is the constant density of ions. ( f (x, v, t) dv n0 (x) ), Difficulties: - in 3D = 7 variables: 3 positions, 3 velocities and time = very heavy computations, - no dissipation = the solution may become turbulent (filaments, vortices, etc.). Obectives: - reduce the dimension of the problem, - keep enough information about the solution. and method of moments for Vlasov
4 We consider f piecewise constant in v: f (x, v, t) = f (x, v, t) = k A, for v [ v k (x, t), v + k (x, t)] =1 N ( A (H =1 v + (x, t) v ) ( H where H is the Heaviside function: { 0 if u < 0, H (u) = 1 if u > 0, v )) (x, t) v, N is an integer, v ± (x, t) are velocities and A are constants, for = 1,..., N. and method of moments for Vlasov
5 3 pairs of bags model in phase space (left) and corresponding distribution function (right). Idea of the model: evolve the v ± (x, t) instead of f. Resolution of transport equations for the v ± by a finite volume method (Godunov for example): t v ± + v ± x v ± q E = 0, = 1,..., N. m and method of moments for Vlasov
6 Definition: moment of order k, in v, of a function f integrable in v: M k (x, t) = + v k f (x, v, t) dv. Moments system (by taking the 2N first moments of the Vlasov equation): t M 0 + x M 1 = 0 t M k + x M k+1 k q m EM k 1 = 0, for k = 1,..., 2N 1 t E (x, t) = q m (M 1 (0, t) M 1 (x, t)). and method of moments for Vlasov
7 Remarks: - v no longer involved = one dimension less, - N can be chosen according to the expected precision and the computation capabilities, - unknown: M 0,..., M 2N 1, and M 2N = need of a closure relation. Closure relation: we have to choose a representation of f, for example by the multi-water-bag model. Moments are then written: M k (x, t) = N =1 v +k+1 A = expression of M 2N as a function of v ±. (x, t) v k+1 (x, t), k = 0,..., 2N k + 1 and method of moments for Vlasov
8 Algorithm Introduction Moments, E and v ± (x, t) known at time n. Computation of moments and E at time n + 1 with a kinetic scheme (natural expression of fluxes): - upwind discretization for the Vlasov equation: f n+1 i fi n + 1 ( v + ( fi n f n ) i 1 + v ( fi+1 n f n )) i E n t x i v fi n = 0, - multiplication by v k and integration in v: 1 v k ( f n+1 i f n ) i dv t + 1 ( ( v +k+1 f n i f n ) i 1 + v k+1 ( fi+1 n f n ) ) i dv x E n i v k v f n i dv = 0, and method of moments for Vlasov
9 - finite volume scheme: M n+1 k,i M n k,i t - expression of the flux: + 1 ( F n x k,i Fk,i 1) n ke n i Mk 1,i n = 0, Fk,i n = F ( Mk,i n, ) N ( Mn k,i+1 = A + =1 ( N A =1 v ++ v + where u + = max (u, 0) and u = min (u, 0), v +k+1 v k+1 k + 2 v + k + 2 v +k+1 v k+1 k + 2 v k + 2 ) n i ) n i+1, and method of moments for Vlasov
10 - scheme for the electric field: E n+1 i t E n i = M n 1,i M n 1,0. Computation of v ± (x, t) at time n + 1 by: - the Newton-Raphson method, iterative and hard to initialize, - the algorithm of Gosse and Runborg when A = ±1, unstable when two bags are close to each other. and method of moments for Vlasov
11 Validation Introduction Validation of the method on Landau damping and two-stream instability, by comparing it to the Particle In Cell (PIC) method. Initial condition of the Landau damping: 1 f (x, v, t = 0) = (1 + ɛ cos (kx)) e v2 2 2π E (x, t = 0) = q m ɛ sin (kx). k Initial condition of the two-stream instability: ) 1 f (x, v, t = 0) = (1 + ɛ cos (kx)) (e 2 (v v 0) e (v+v 0) 2 2 2π E (x, t = 0) = q ɛ sin (kx). m k and method of moments for Vlasov
12 Landau damping Landau damping: k=0.5, ɛ = Electrical energy (logarithmic scale) as a function of time. and method of moments for Vlasov
13 Studied test cases Initial conditions such that: - f (x, v, t = 0) is exactly described by one or two pairs of bags, - the solution destabilizes: creation of filaments, vortices, etc. Examples at initial time: and method of moments for Vlasov
14 First test case: 1 pair of bags Phase space (v as a function of x) at times 0 and 2: particles (PIC method), v + 1 and v 1 (water-bag). and method of moments for Vlasov
15 Phase space (v as a function of x) at times 5 and 10: particles (PIC method), v + 1 and v 1 (water-bag). and method of moments for Vlasov
16 Phase space (v as a function of x) at times 20 and 50: particles (PIC method), v + 1 and v 1 (water-bag). and method of moments for Vlasov
17 Second test case: 2 pairs of bags Phase space (v as a function of x) at times 0 and 2: particles (PIC method), v + 1, v 1, v + 2 and v 2 (water-bag). and method of moments for Vlasov
18 Phase space (v as a function of x) at times 3 and 4: particles (PIC method), v + 1, v 1, v + 2 and v 2 (water-bag). and method of moments for Vlasov
19 Phase space (v as a function of x) at times 5 and 10: particles (PIC method), v + 1, v 1, v + 2 and v 2 (water-bag). and method of moments for Vlasov
20 Conclusions: - monovalued solution correctly described, - v ± (x, t) cannot be multivalued = loss of information when filaments appear, - used algorithm unstable when at least two bags are close to each other (ill-conditioned matrix). Work in progress in order to improve the representation of the filaments: - stabilize the algorithm of Gosse and Runborg when several bags are close to each other, - first idea: merge these bags. and method of moments for Vlasov
21 Phase space (v as a function of x) at times 0 and 0.1: particles (PIC method), v + 1, v 1, v + 2, v 2, v + 3 and v 3 (water-bag). and method of moments for Vlasov
22 Phase space (v as a function of x) at time 0.4: particles (PIC method), v + 1, v 1, v + 2, v 2, v + 3 and v 3 (water-bag). and method of moments for Vlasov
23 Phase space (v as a function of x) at time 0.4: - left: f obtained by the PIC method (particles), - right: f reconstructed by the method of moments/multi-water-bag. and method of moments for Vlasov
24 Perspectives: - generalize the algorithm of Gosse and Runborg when the A are different of ±1, - couple our method to a PIC method: create particles when filaments appear. Thank you for your attention! and method of moments for Vlasov
25 References Introduction - N. Besse, F. Berthelin, Y. Brenier, P. Bertrand: The multi-water-bag equations for collisionless kinetic modeling, Kinetic and Related Models 2, (2009) - R.O. Fox, F. Laurent, M. Massot: Numerical simulation of spray coalescence in an Eulerian framework: Direct quadrature method of moments and multi-fluid method, Journal of Computational Physics 227, (2008) Algorithm of Gosse and Runborg - Y. Brenier, L. Corrias: A kinetic formulation for multi-branch entropy solutions of scalar conservation laws, Annales de l Institut Henri Poincaré. Analyse Non Linéaire 15, (1998) - L. Gosse, O. Runborg: Resolution of the finite Markov moment problem, Comptes Rendus Mathématique. Académie des Sciences. Paris 12, (2005) and method of moments for Vlasov
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 informationHigh 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 informationNhung 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 informationBenchmarks 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 informationBerk-Breizman and diocotron instability testcases
Berk-Breizman and diocotron instability testcases M. Mehrenberger, N. Crouseilles, V. Grandgirard, S. Hirstoaga, E. Madaule, J. Petri, E. Sonnendrücker Université de Strasbourg, IRMA (France); INRIA Grand
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. The 8th International Conference on Computational
More informationModélisation et simulation d écoulements diphasiques avec contrainte granulaire
Modélisation et simulation d écoulements diphasiques avec contrainte granulaire T. Gallouët, P. Helluy, J.-M. Hérard, J. Nussbaum LATP Marseille IRMA Strasbourg EDF Chatou ISL Saint-Louis Context Flow
More informationAdaptive simulation of Vlasov equations in arbitrary dimension using interpolatory hierarchical bases
Adaptive simulation of Vlasov equations in arbitrary dimension using interpolatory hierarchical bases Erwan Deriaz Fusion Plasma Team Institut Jean Lamour (Nancy), CNRS / Université de Lorraine VLASIX
More informationWhat 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 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 informationReduced Vlasov-Maxwell modeling
Reduced Vlasov-Maxwell modeling Philippe Helluy, Michel Massaro, Laurent Navoret, Nhung Pham, Thomas Strub To cite this version: Philippe Helluy, Michel Massaro, Laurent Navoret, Nhung Pham, Thomas Strub.
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 informationMoments 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 informationMichel 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 informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
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 informationAccurate 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 informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationComputational Methods in Plasma Physics
Computational Methods in Plasma Physics Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Purpose of Talk Describe use of numerical methods to solve simple problem in plasma
More informationHybrid Simulation Method ISSS-10 Banff 2011
Hybrid Simulation Method ISSS-10 Banff 2011 David Burgess Astronomy Unit Queen Mary University of London With thanks to Dietmar Krauss-Varban Space Plasmas: From Sun to Earth Space Plasma Plasma is (mostly)
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 1 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationA new fully two-dimensional conservative semi-lagrangian method: applications on polar grids, from diocotron instability to ITG turbulence
A new fully two-dimensional conservative semi-lagrangian method: applications on polar grids, from diocotron instability to ITG turbulence Nicolas Crouseilles, Pierre Glanc, Sever Adrian Hirstoaga, Eric
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 informationLandau 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 informationEulerian calculations of wave breaking and multi-valued solutions in a traveling wave tube
Eulerian calculations of wave breaking and multi-valued solutions in a traveling wave tube J.G. Wöhlbier, S. Jin, 2 and S. Sengele 3 Los Alamos National Laboratory, ISR-6, Radio Frequency and Accelerator
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationWater-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes
Water-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes PhD work directed by N. Besse with invaluable input from P. Bertrand, E. Gravier, P. Morel,R.
More informationKinetic 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 informationPh.D. Qualifying Exam Friday Saturday, January 3 4, 2014
Ph.D. Qualifying Exam Friday Saturday, January 3 4, 2014 Put your solution to each problem on a separate sheet of paper. Problem 1. (5166) Assume that two random samples {x i } and {y i } are independently
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 informationHamiltonian fluid closures of the Vlasov-Ampère equations: from water-bags to N moment models
Hamiltonian fluid closures of the Vlasov-Ampère equations: from water-bags to N moment models M Perin Cristel Chandre P.J. Morrison E Tassi To cite this version: M Perin Cristel Chandre P.J. Morrison E
More informationFinite Volume Schemes for Vlasov
Finite Volume Schemes for Vlasov Nicolas Crouseilles, Pierre Glanc, Michel Mehrenberger, Christophe Steiner To cite this version: Nicolas Crouseilles, Pierre Glanc, Michel Mehrenberger, Christophe Steiner.
More informationParametric Studies of Kinetically-Enhanced Raman Backscatter and Electron Acoustic Thomson Scattering
Parametric Studies of Kinetically-Enhanced Raman Backscatter and Electron Acoustic Thomson Scattering D. J. Strozzi 1, E. A. Williams 1, A. B. Langdon 1, A. Bers 2 1 Lawrence Livermore National Laboratory,
More information13. REDUCED MHD. Since the magnetic field is almost uniform and uni-directional, the field has one almost uniform component ( B z
13. REDUCED MHD One often encounters situations in which the magnetic field is strong and almost unidirectional. Since a constant field does not produce a current density, these fields are sometimes said
More informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationEXASCALE COMPUTING. Implementation of a 2D Electrostatic Particle in Cell algorithm in UniÞed Parallel C with dynamic load-balancing
ExaScience Lab Intel Labs Europe EXASCALE COMPUTING Implementation of a 2D Electrostatic Particle in Cell algorithm in UniÞed Parallel C with dynamic load-balancing B. Verleye P. Henry R. Wuyts G. Lapenta
More informationVlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L22109, doi:10.1029/2008gl035608, 2008 Vlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field
More informationPalindromic Discontinuous Galerkin Method
Palindromic Discontinuous Galerkin Method David Coulette, Emmanuel Franck, Philippe Helluy, Michel Mehrenberger, Laurent Navoret To cite this version: David Coulette, Emmanuel Franck, Philippe Helluy,
More informationGyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas
Frank Jenko Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas Max-Planck-Institut für Plasmaphysik, Garching Workshop on Vlasov-Maxwell Kinetics WPI, Vienna,
More informationA 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 informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More informationContraction Mappings Consider the equation
Contraction Mappings Consider the equation x = cos x. If we plot the graphs of y = cos x and y = x, we see that they intersect at a unique point for x 0.7. This point is called a fixed point of the function
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 informationAdaptive 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 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 informationSummer 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 informationConservative semi-lagrangian schemes for Vlasov equations
Conservative semi-lagrangian schemes for Vlasov equations Nicolas Crouseilles Michel Mehrenberger Eric Sonnendrücker October 3, 9 Abstract Conservative methods for the numerical solution of the Vlasov
More informationWeibel instability and filamentary structures of a relativistic electron beam in plasma
Mitglied der Helmholtz-Gemeinschaft 7 th Direct Drive and Fast Ignition Workshop Prague, 3-6 May, 009 Weibel instability and filamentary structures of a relativistic electron beam in plasma Anupam Karmakar
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationComputational Methods of Physical Problems with Mixed Scales
Computational Methods of Physical Problems with Mixed Scales Shi Jin University of Wisconsin-Madison IPAM Tutorial, March 10-12, 2008 Scales and Physical Laws Figure from E & Engquist, AMS Notice Connections
More informationA mechanism for magnetic field stochastization and energy release during an edge pedestal collapse
A mechanism for magnetic field stochastization and energy release during an edge pedestal collapse S. S. Kim, Hogun Jhang, T. Rhee, G. Y. Park, R. Singh National Fusion Research Institute, Korea Acknowledgements:
More informationOn stably stratified homogeneous shear flows subjected to rotation
Center for Turbulence Research Proceedings of the Summer Program 2000 241 On stably stratified homogeneous shear flows subjected to rotation By B. A. Pettersson Reif, A.Ooi AND P. A. Durbin Theoretical
More informationHybrid Simulations: Numerical Details and Current Applications
Hybrid Simulations: Numerical Details and Current Applications Dietmar Krauss-Varban and numerous collaborators Space Sciences Laboratory, UC Berkeley, USA Boulder, 07/25/2008 Content 1. Heliospheric/Space
More informationDecay instability of electron acoustic waves
Communications in Nonlinear Science and Numerical Simulation 13 (2008) 215 220 www.elsevier.com/locate/cnsns ecay instability of electron acoustic waves F. Valentini a,b, *, T.M. O Neil b,.h.e. ubin b
More informationFloquet Theory for Internal Gravity Waves in a Density-Stratified Fluid. Yuanxun Bill Bao Senior Supervisor: Professor David J. Muraki August 3, 2012
Floquet Theory for Internal Gravity Waves in a Density-Stratified Fluid Yuanxun Bill Bao Senior Supervisor: Professor David J. Muraki August 3, 212 Density-Stratified Fluid Dynamics Density-Stratified
More informationModels for Global Plasma Dynamics
Models for Global Plasma Dynamics F.L. Waelbroeck Institute for Fusion Studies, The University of Texas at Austin International ITER Summer School June 2010 Outline 1 Models for Long-Wavelength Plasma
More informationMinimization formulation of a bi-kinetic sheath
Minimization formulation of a bi-kinetic sheath Mehdi Badsi (PHD), B. Després and Martin Campos-Pinto Laboratoire Jacques-Louis Lions UPMC-Paris VI, CNRS UMR 7598, Paris, France Sponsors: ANR Chrome, FRFCM
More informationTheory and simulations of hydrodynamic instabilities in inertial fusion
Theory and simulations of hydrodynamic instabilities in inertial fusion R. Betti Fusion Science Center, Laboratory for Laser Energetics, University of Rochester IPAM/UCLA Long Program PL2012 - March 12
More informationWave Turbulence and Condensation in an Optical Experiment
Wave Turbulence and Condensation in an Optical Experiment S. Residori, U. Bortolozzo Institut Non Linéaire de Nice, CNRS, France S. Nazarenko, J. Laurie Mathematics Institute, University of Warwick, UK
More informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
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 informationNatalia Tronko S.V.Nazarenko S. Galtier
IPP Garching, ESF Exploratory Workshop Natalia Tronko University of York, York Plasma Institute In collaboration with S.V.Nazarenko University of Warwick S. Galtier University of Paris XI Outline Motivations:
More informationLong time behavior of solutions of Vlasov-like Equations
Long time behavior of solutions of Vlasov-like Equations Emanuele Caglioti "Sapienza" Università di Roma (Mathematics) Journée Claude Bardos - Laboratoire Jacques Louis Lions - 2011 Equations Vlasov Poisson
More informationAn Eulerian method for computing multi-valued solutions of the Euler-Poisson equations and applications to wave breaking in klystrons
Submitted to Phys. Rev. E, March 3 An Eulerian method for computing multi-valued solutions of the Euler-Poisson equations and applications to wave breaking in klystrons Xiantao Li, John G. Wöhlbier, Shi
More informationSummary of the Equations of Fluid Dynamics
Reference: Summary of the Equations of Fluid Dynamics Fluid Mechanics, L.D. Landau & E.M. Lifshitz 1 Introduction Emission processes give us diagnostics with which to estimate important parameters, such
More informationThermodynamic form of the equation of motion for perfect fluids of grade n
Thermodynamic form of the equation of motion for perfect fluids of grade n Henri Gouin To cite this version: Henri Gouin. Thermodynamic form of the equation of motion for perfect fluids of grade n. Comptes
More informationComparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry
Comparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry D. D. Schnack University of Wisconsin Madison Jianhua Cheng, S. E. Parker University of Colorado
More informationIntegration of Vlasov-type equations
Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma
More informationStatistical Mechanics and Evolution Equations Rome, January 28-29, 2016
Statistical Mechanics and Evolution Equations Rome, January 28-29, 2016 A Workshop in Honor of Carlo Marchioro on the Occasion of his 70th Birthday TITLES AND ABSTRACTS Kazuo Aoki (Kyoto University, Japan)
More informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationSimulation of Plasma Wakefields and Weibel Instability of Electron Beams in Plasma Using Codes LSP and OSIRIS
Simulation of Plasma Wakefields and Weibel Instability of Electron Beams in Plasma Using Codes and OSIRIS 15 Electron density Beam density Plasma density 4 6 4 A. Solodov, C. Ren, J. Myatt, and R. Betti
More informationA fully implicit, exactly conserving algorithm for multidimensional particle-in-cell kinetic simulations
A fully implicit, exactly conserving algorithm for multidimensional particle-in-cell kinetic simulations L. Chacón Applied Mathematics and Plasma Physics Group Theoretical Division Los Alamos National
More informationComputation of the Semiclassical Limit of the Schrödinger Equation with Phase Shift by a Level Set Method
Computation of the Semiclassical Limit of the Schrödinger Equation with Phase Shift by a Level Set Method Shi Jin and Xu Yang Abstract In this paper, we show how the level set method, developed in [8,
More informationAdvection, Conservation, Conserved Physical Quantities, Wave Equations
EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary
More informationAccurate Multiple-Precision Gauss-Legendre Quadrature
Accurate Multiple-Precision Gauss-Legendre Quadrature Laurent Fousse 1,2 1 Université Henri-Poincaré Nancy 1 2 INRIA Rocquencourt 18th IEEE International Symposium on Computer Arithmetic Outline Motivation
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationHigh Resolution Wave Propagation Schemes for Two-Fluid Plasma Simulations
High Resolution Wave Propagation Schemes for Two-Fluid Plasma Simulations Ammar H. Hakim A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy University
More informationKinetics of the Raman Instability in a Laser Plasma
WDS'05 Proceedings of Contributed Papers, Part II, 383 390, 2005. ISBN 80-86732-59-2 MATFYZPRESS Kinetics of the Raman Instability in a Laser Plasma M. Mašek and K. Rohlena Institute of Physics, Academy
More informationAMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code
AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code George Wilkie (gwilkie@umd.edu) Supervisor: William Dorland (bdorland@umd.edu) October 11, 2011 Abstract Simulations of turbulent plasma in
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationModern Challenges in Nonlinear Plasma Physics A Conference Honouring the Career of Dennis Papadopoulos
Modern Challenges in Nonlinear Plasma Physics A Conference Honouring the Career of Dennis Papadopoulos Progress in Plasma Physics by Numerical Simulation: Collisionless Shocks Manfred Scholer Max-Planck-Institut
More informationBernstein-Greene-Kruskal (BGK) Modes in a Three Dimensional Unmagnetized Plasma
Bernstein-Greene-Kruskal (BGK) Modes in a Three Dimensional Unmagnetized Plasma C. S. Ng & A. Bhattacharjee Space Science Center University of New Hampshire Space Science Center Seminar, March 9, 2005
More informationLattice Boltzmann approach to liquid - vapour separation
Lattice Boltzmann approach to liquid - vapour separation T.Biciușcă 1,2, A.Cristea 1, A.Lamura 3, G.Gonnella 4, V.Sofonea 1 1 Centre for Fundamental and Advanced Technical Research, Romanian Academy Bd.
More informationStability of a plasma confined in a dipole field
PHYSICS OF PLASMAS VOLUME 5, NUMBER 10 OCTOBER 1998 Stability of a plasma confined in a dipole field Plasma Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 Received
More informationNewton s Second Law of Motion. Isaac Newton (in 1689) The Universe is Regular and Predictable. The Foundation of Science
Key Concepts and Fundamental Theorems Isaac Newton (in 689) of Atmospheric Science This is the earliest portrait of Newton to survive. The artist was Godfrey Kneller, perhaps the greatest portrait painter
More informationReview exercise
Review eercise y cos sin When : 8 y and 8 gradient of normal is 8 y When : 9 y and 8 Equation of normal is y 8 8 y8 8 8 8y 8 8 8 8y 8 8 8 8y 8 8 8 y e ln( ) e ln e When : y e ln and e Equation of tangent
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationOptimization of Particle-In-Cell simulations for Vlasov-Poisson system with strong magnetic field
Optimization of Particle-In-Cell simulations for Vlasov-Poisson system with strong magnetic field Edwin Chacon-Golcher Sever A. Hirstoaga Mathieu Lutz Abstract We study the dynamics of charged particles
More informationQuadrature-based moment methods for polydisperse multiphase flow modeling
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2013 Quadrature-based moment methods for polydisperse multiphase flow modeling Cansheng Yuan Iowa State University
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 informationShock Waves. 1 Steepening of sound waves. We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: kˆ.
Shock Waves Steepening of sound waves We have the result that the velocity of a sound wave in an arbitrary reference frame is given by: v u kˆ c s kˆ where u is the velocity of the fluid and k is the wave
More informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationRecent computational methods for high frequency waves in heterogeneous media ABSTRACT. Shi Jin. 1. Introduction
Recent computational methods for high frequency waves in heterogeneous media Shi Jin ABSTRACT In this note, we review our recent results on the Eulerian computation of high frequency waves in heterogeneous
More informationNUMERICAL RESOLUTION OF CONSERVATION LAWS WITH OPENCL
ESAIM: PROCEEDINGS, Juillet 2013, Vol. 40, p. 51-62 C. Bourdarias, S. Gerbi, Editors NUMERICAL RESOLUTION OF CONSERVATION LAWS WITH OPENCL A. Crestetto 1, P. Helluy 2 and J. Jung 3 Abstract. We present
More informationEntropy evolution and dissipation in collisionless particle-in-cell gyrokinetic simulations
Max-Planck-Insititut für Plasmaphysik Entropy evolution and dissipation in collisionless particle-in-cell gyrokinetic simulations A. Bottino Objectives Develop a numerical tool able to reproduce and predict
More informationA FIELD SPACE-BASED LEVEL SET METHOD FOR COMPUTING MULTI-VALUED SOLUTIONS TO 1D EULER-POISSON EQUATIONS
A FIELD SPACE-BASED LEVEL SET METHOD FOR COMPUTING MULTI-VALUED SOLUTIONS TO 1D EULER-POISSON EQUATIONS HAILIANG LIU AND ZHONGMING WANG Abstract. We present a field space based level set method for computing
More informationNUMERICAL SIMULATION OF DROPLET DISPERSION AND EVAPORATION WITH A MOMENTS-BASED CFD MODEL
Proceedings of COBEM 29 Copyright c 29 by ABCM 2th International Congress of Mechanical Engineering November 15-2, 29, Gramado, RS, Brazil NUMERICAL SIMULATION OF DROPLET DISPERSION AND EVAPORATION WITH
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
More information