Minimization formulation of a bi-kinetic sheath
|
|
- Ami James
- 5 years ago
- Views:
Transcription
1 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 29 septembre 2014 Minimization formulation of a bi-kinetic sheath p. 1 / 25
2 Section 1 Minimization formulation of a bi-kinetic sheath p. 2 / 25
3 What is a sheath? Sheath = boundary layer for a quasi-neutral plasma (i +, e ) facing a ( isolated, metallic) wall. Important parameter ε λ D << 1. Physics : due to the high mobility ratio (m i /me 2000 >> 1), the local quasi-neutrality in the core is highly violated at the wall, and the plasma continuously flows towards the wall. This is kind of a paradox since there is no exterior forces, and the trivial solution n i n e = φ = 0 is not the physical one. Terminology : core, Debye length, pre-sheath, sheath, Boltzmanian e (n e = n 0 e φ ), Bohm criterion, floating potential, ambipolar hypothesis,... From Stangeby : The Plasma Boundary of Magnetic Fusion Devices (2011). Fundamental for applications : Tokamaks (ITER), ionic engine (satellites),... Minimization formulation of a bi-kinetic sheath p. 2 / 25
4 From Chen The physical solution is n i = n i2 > n e. Here E = φ 0 Condition on the Mach number : M > 1. Kinetic Bohm criterion : 0 f given i v 2 dv < 1. Minimization formulation of a bi-kinetic sheath p. 3 / 25
5 Some references Physics : Chen ( to Plasma Physics 74 ) ; Review by Riemann (The Bohm criterion and sheath formation 91 ), Baalrud-Hegna (Kinetic theory of the presheath and the Bohm criterion, Manfredi-Devaux, Magnetized plasma-wall transition. Consequences for wall sputtering and erosion, Stangeby The Plasma Boundary of Magnetic Fusion Devices, Math : Raviart- Greengard, (A Boundary-Value problem for the stationary Vlasov-Poisson equations : The Plane Diode, 1990) P-H. Maire, Fluid models with a plasma and b.c., PHD thesis 96 (with Sentis and Golse) H. Guillard, M 1 in fluid models at the SOL (2012). Additional references : the Milne problem Marschak 47 : variational equation for the mean field. Modern treatment in Golse and al,... Minimization formulation of a bi-kinetic sheath p. 4 / 25
6 Our angle of attack Do not use Maxwellian e which ends up in a compatibility issue with the ambipolar hypothesis, and usually with different models in different parts of the domain. = Following Stangeby (see below), use truncated Maxwellians for e Do not use fluid equations for ions : = start from bi-kinetic (Vlasov) + Poisson equation Study exact stationary solutions with convenient B.C., eliminate everything f i, f e, n i, n e in function of the electric potential φ Solve the equation for φ Minimization formulation of a bi-kinetic sheath p. 5 / 25
7 Section 2 Minimization formulation of a bi-kinetic sheath p. 6 / 25
8 Equations Domain is (x, v) [0, 1] R v x f i φ (x) v f i = 0, Vlasov i + v x f e + m i φ m (x) v f e = 0, e Vlasov e ε 2 φ (x) f i dv + f edv = 0, Poisson, (ε λ D ) vfi dv vf edv = 0, Ambipolarity, (J = te = 0) f i (0, v) = f given i (v), v > 0, f i (1, v) = 0, v < 0, f e(0, v) = n 0 me m i e me v 2 2m i, v > 0, Scaled Maxwellian at entrance f e(1, v) = αf e(1, v), v < 0, Material dependent reemission of e with parameter 0 α = α given < 1, φ(0) = 0, φ(1) = φ w, Just a convention, Floating potential, unknown at this stage. Additional condition is E = φ 0 Warning : n 0 and φ w are unknowns at this level. Minimization formulation of a bi-kinetic sheath p. 6 / 25
9 Trajectories Assume the physical condition E = φ > 0 e i + Char. : 1 2 v 2 m i φ(x) = cst m e Char. : 1 2 v 2 + φ(x) = cst Sep. : 1 2 v 2 m i φ(x) = m i φ m e m w e Sep. : 1 2 v 2 + φ(x) = 0 Minimization formulation of a bi-kinetic sheath p. 7 / 25
10 ( v Notice for example that me 2 2 m ) i me φ m i R e me 2m i dv = 2πe φ. One gets after more computations n e(x) = R fedv = n 0 ( 2πe φ(x) (1 α) γ e = R vfedv = (1 α)n 0 mi m e e φw, 2φ w e v 2 2 Formulas ) v v 2 +2φ(x) dv, n i (x) = R f i dv = 0 f given v i (v) dv, v 2 2φ(x) γ i = R vf i dv = 0 f given i (v)vdv. One still has two degrees of freedom. In order to simplify the mathematical analysis, we consider two equations { γi = γ e Ambipolarity, n i (0) = n e(0) Local neutrality at entrance. Elimination of n 0 yields the resulting non linear equation with φ w only γ e n = γ i e(0) n i (0) Minimization formulation of a bi-kinetic sheath p. 8 / 25
11 Equation for φ w The equation writes W(φ w ) = b, < φ w 0, with ( mi ) W(ψ) = f given i (v)dv e ψ + m e 0 0 f given i (v)vdv e v 2 2ψ 2 dv and b = 2π 1 α 0 f given i (v)vdv. A first result : The function W being monotone (increasing) there exists a unique φ w if and only 0 f given i (v)vdv 0 f given i (v)dv 2m i 1 α πm e 1 + α. Corollary : No solutions for α = 1 (that would imply Boltzmanian e ). Corollary : If the is a =, then φ w = 0 and the solution is trivial. Corollary : One has the bound [ ln ( 2π ) ln ( 2 + mi m e 0 0 f given (v)dv i f given (v)vdv i )] ln(1 α) φ w 0. Minimization formulation of a bi-kinetic sheath p. 9 / 25
12 Non linear Poisson equation Now that φ w and n 0 are uniquely determined, one can solves { ε 2 φ (x) = n i (φ) n e(φ) := q (φ), φ(0) = 0 and φ(1) = φ w, where the potential is q(φ) = f given i (v)v v 2 2φdv 0 ( 2πe +n φ 2 0 (1 α) 2 v ) v 2 + 2φdv 0. 2φw e v Set with the space 1 J ε(φ) = ε 2 φ (x) q(φ(x))dx, φ V V = { v H 1 (0, L), 0 = v(0) v v(l) = φ w }. Proposition (evident) : There exists a minimum of the problem J ε(φ ε) J ε(φ), φ ε V, φ V. Minimization formulation of a bi-kinetic sheath p. 10 / 25
13 Preliminary properties of the minimizer Since 0 q by construction, the inequality J ε(φ ε) J(xφ w ) yields a bound on the H 1 norm : φ ε H 1 (0,1) = O ( ε 1). Minimizers are non increasing 0 x y 1 = 0 φ ε(x) φ ε(y) φ w. This is proved by comparison methods. If 0 > x > 1 = 0 > φ ε(x) > φ w, the minimizer is C 2 solution of the non linear Poisson equation. Since now ε 2 φ ε (x) = n i n e, it yields n i n e H 1 (0,1) = O(ε) which shows quasi-neutrality with a boundary layer (equal to the physical sheath) of size ε in a weak norm. Minimization formulation of a bi-kinetic sheath p. 11 / 25
14 Bohm criterion The interesting physical situation is when n i n e. At entrance (x = 0) (n i n e) (φ) = q (φ) = q (0)φ + O(φ 2 ) Since φ 0, a necessary condition is q (0) 0 (Bohm criterion) Complete Bohm crit. : 0 f given i (v) dv 2π + (1 α) 2φw v 2 0 f given i (v)dv e v2 2 dv v 2. 2π (1 α) 2φw e v2 2 dv Proposition : If this inequality is true, then n i n e everywhere in [0, 1]. Notice the simplified Bohm criterion 0 f given i (v) dv v 2 0 f given i (v)dv < 1. Minimization formulation of a bi-kinetic sheath p. 12 / 25
15 Compatiblity issue We need at the same time two inequalities. First one says : low velocities 0 f given i (v)vdv φ w : 0 f given i (v)dv < 2m i 1 α πm e 1 + α. Second one says : high velocities Simplified Bohm criterion : Lemma : There exists non trivial f given i conditions if 1 < 2mi 1 α πm e 1+α. 0 f given i (v) dv v 2 0 f given i (v)dv < 1. L 1 (R + ) L (R + ) that satisfy all Proof of the necessary condition : based on double Cauchy-Schwarz inequality ( 4 ( f given i (v)dv) f given i (v) dv ) 2 ( ) 2 f given i (v)vdv 0 ( 0 ) ( f given i (v)dv 0 0 v 0 f given i (v) dv ) ( 2 v 2 f given i (v)vdv). 0 Minimization formulation of a bi-kinetic sheath p. 13 / 25
16 Mobility ratio Quite fortunately the high mass ratio 1 << m i 3000 (for Z = 1.5) yields m e 2m i 43.7 πm e Proposition : There exists a important interval in which sheath solutions exist 0 α α c 0.93 These are also exact stationary solutions of bi-kinetic+maxwell equations. Minimization formulation of a bi-kinetic sheath p. 14 / 25
17 Summary of the claim Claim : Assume 0 α α c, ρ 0 = n i (0) n 0 (0) = 0, f given i I ad (ρ 0, α) and ε > 0. Assume the Bohm criterion. Then the Vlasov-Poisson-Ampere formulation of the sheath is well-posed in the following sense : There exists a unique strictly concave and decreasing φ ε C 2 [0, 1] solution of the non linear Poisson equation with b.c. φ(0) = 0 and φ(1) = φ w < 0, where the equation on φ w does not depend on ε. The densities f i, f e ( L 1 L ) ( (0, 1) R +) are weak solutions : note that f e is not a strong solution since it is a truncated Maxwellian. There is a sheath (i.e. a boundary layer) : φ ε H 1 (0,1) = O ( ε 1) and n i n e H 1 (0,1) = O(ε). The local charge is positive : ni(x) > n e(x) for x (0, 1]. Minimization formulation of a bi-kinetic sheath p. 15 / 25
18 Section 3 Minimization formulation of a bi-kinetic sheath p. 16 / 25
19 Algorithm Numerical algorithm developed by Mehdi. Choose a f given i that satisfies all conditions and inequalities. Use high order quadratures to compute q. Solve the equation on φ w Use a gradient algorithm to compute φ. Converge until φ n+1 φ n ɛ. Minimization formulation of a bi-kinetic sheath p. 16 / 25
20 The potential q for the chosen data Here φ w 2.4. Seems q is convex in this case : characterization still an open problem. Minimization formulation of a bi-kinetic sheath p. 17 / 25
21 Ion distribution : ε = 1, 0.1, 0.01 and zoom Note ε is a scaling factor : y = x ε Minimization formulation of a bi-kinetic sheath p. 18 / 25
22 Electron distribution : ε = 1, 0.1, 0.01 and zoom Minimization formulation of a bi-kinetic sheath p. 19 / 25
23 n i and n e ε = 1 ε = 0.1 ε = 0.01 ε = 0.01 (zoom) Minimization formulation of a bi-kinetic sheath p. 20 / 25
24 φ and u = fi vdv fi dv φ u ε = 1 ε = 0.1 ε = 0.01 Minimization formulation of a bi-kinetic sheath p. 21 / 25
25 Influence of the reemission factor Seems important in view of engineering and applications α (n i, n e) is increasing function. Minimization formulation of a bi-kinetic sheath p. 22 / 25
26 Section 4 Minimization formulation of a bi-kinetic sheath p. 23 / 25
27 What is a RF sheath? It is an new problem discussed in plasma physics community : Heuraux (U. Nancy), Colas (CEA-IRFM), Myra-D ippolito,.... Heating of a fusion plasma can be described by a radio-frequency wave φ rf cos ωt. PDE model is linear Maxwell s equations by with the cold plasma dielectric tensor, but......sheath modify strongly the boundary conditions. The boundary condition might be non linear due to the sheath. Minimization formulation of a bi-kinetic sheath p. 23 / 25
28 A fundamental simple idea Boltzmannian electrons are n e = n 0 e φ. A RF wave modifies the electric potential φ tot = φ + φ rf cos(ωt) so that averaging in time modify the apparent electronic density since 1 T e φ+φ rf cos(ωt) dt = e φ 1 T e φ rf cos(ωt) dt = e φ I 0 (φ rf ), T = 2π T 0 T 0 ω where the rectification factor I 0 (φ rf ) 1 is the modified Bessel function of the first kind. Finally n e = n apparent 0 e φ, n apparent 0 = n 0 I 0 (φ rf ). = It is possible to add the multiplicative rectification factor in the previous equations. = At first order of approximation : φ w and φ are the same. I do not know if this conclusion is physically correct for RF sheath. Minimization formulation of a bi-kinetic sheath p. 24 / 25
29 A sheath solution has been constructed, starting from stationary solutions of two species Vlasov + Poisson. A well-posed non linear minimization equation shows up after reduction for the mean field (Poisson). It is reminiscent of the seminal Marshak s proposition of a variational formulation of the Milne s problem. is very robust. More to be done to compare with physical solutions. This can be a basis for further mathematical developments on kinetic sheath solutions. Preprint soon online. Minimization formulation of a bi-kinetic sheath p. 25 / 25
Numerical stability of plasma sheath
Numerical stability of plasma sheath Mehdi Badsi, Michel Mehrenberger, Laurent Navoret To cite this version: Mehdi Badsi, Michel Mehrenberger, Laurent Navoret. Numerical stability of plasma sheath. 2018.
More informationScattering structure of Vlasov equations around inhomogeneous Boltzmannian states
B. Després (LJLL/UPMC and IUF Scattering structure of Vlasov s around inhomogeneous Boltzmannian states B. Després (LJLL/UPMC and IUF Collisionless Boltzmann (Vlasov Equation 207 p. / 2 Goal Model problem
More informationSheath stability conditions for PIC simulation
Plasma_conference (2011/11/24,Kanazawa,Japan) PIC Sheath stability conditions for PIC simulation Osaka Prefecture. University H. Matsuura, and N. Inagaki, Background Langmuir probe is the most old but
More informationQuasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths
in the presence of plasma sheaths Department of Mathematical Sciences Ulsan National Institute of Science and Technology (UNIST) joint work with Masahiro Suzuki (Nagoya) and Chang-Yeol Jung (Ulsan) The
More informationVlasov simulations of plasma-wall interactions in a magnetized and weakly collisional plasma
PHYSICS OF PLASMAS 13, 083504 2006 Vlasov simulations of plasma-wall interactions in a magnetized and weakly collisional plasma S. Devaux a Laboratoire de Physique des Milieux Ionisés et Applications,
More informationKinetic Theory of the Presheath and the Bohm Criterion
UW-CPTC 10-3 Kinetic Theory of the Presheath and the Bohm Criterion S D Baalrud and C C Hegna Department of Engineering Physics, University of Wisconsin-Madison, 1500 Engineering Drive, Madison, WI 53706,
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 informationEquilibrium and transport in Tokamaks
Equilibrium and transport in Tokamaks Jacques Blum Laboratoire J.-A. Dieudonné, Université de Nice Sophia-Antipolis Parc Valrose 06108 Nice Cedex 02, France jblum@unice.fr 08 septembre 2008 Jacques Blum
More informationMultiscale modelling of sheath physics in edge transport codes
EUROFUSION WPPFC-CP(6) 6599 N Mellet et al. Multiscale modelling of sheath physics in edge transport codes Preprint of Paper to be submitted for publication in Proceedings of 26th IAEA Fusion Energy Conference
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 informationOne dimensional hybrid Maxwell-Boltzmann model of shearth evolution
Technical collection One dimensional hybrid Maxwell-Boltzmann model of shearth evolution 27 - Conferences publications P. Sarrailh L. Garrigues G. J. M. Hagelaar J. P. Boeuf G. Sandolache S. Rowe B. Jusselin
More informationModélisation de sources plasma froid magnétisé
Modélisation de sources plasma froid magnétisé Gerjan Hagelaar Groupe de Recherche Energétique, Plasma & Hors Equilibre (GREPHE) Laboratoire Plasma et Conversion d Énergie (LAPLACE) Université Paul Sabatier,
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 informationDevelopment of a stable coupling of the Yee scheme with linear current
Development of a stable coupling of the Yee scheme with linear current Martin Campos Pinto (LJLL), Bruno Després (LJLL) Stéphane Heuraux (IJL), Filipe Da Silva (IPFN+IST) 15 octobre 2013 Munich: Pereverzev
More informationTWO-SCALE SIMULATION OF MAXWELL S EQUATIONS
ESAIM: PROCEEDINGS, February 27, Vol.6, 2-223 Eric Cancès & Jean-Frédéric Gerbeau, Editors DOI:.5/proc:27 TWO-SCALE SIMULATION OF MAWELL S EQUATIONS Hyam Abboud,Sébastien Jund 2,Stéphanie Salmon 2, Eric
More informationParallel transport and profile of boundary plasma with a low recycling wall
1 TH/P4-16 Parallel transport and profile of boundary plasma with a low recycling wall Xian-Zhu Tang 1 and Zehua Guo 1 1 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, U.S.A.
More informationKinetic theory of ions in the magnetic presheath
Kinetic theory of ions in the magnetic presheath Alessandro Geraldini 1,2, Felix I. Parra 1,2, Fulvio Militello 2 1. Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, Oxford
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
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 informationElectromagnetic waves in magnetized plasma The dispersion relation
Electromagnetic waves in magnetized plasma The dispersion relation Bruno Després (LJLL-UPMC) Electromagnetic waves in magnetized plasma The dispersion relation p. 1 / 32 Vlasov-Maxwell Vectors are in bold.
More informationProgress on Quantitative Modeling of rf Sheaths
Progress on Quantitative Modeling of rf Sheaths D. A. D Ippolito, J. R. Myra, H. Kohno and J. C. Wright Lodestar Research Corporation, Boulder, Colorado, 80301 May, 2011 Prepared for the 19th Topical Conference
More informationEnergy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas
Energy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas Yingda Cheng Andrew J. Christlieb Xinghui Zhong March 18, 2014 Abstract In this paper, we apply our recently developed energy-conserving
More informationNonlinear instability of periodic BGK waves for Vlasov-Poisson system
Nonlinear instability of periodic BGK waves for Vlasov-Poisson system Zhiwu Lin Courant Institute Abstract We investigate the nonlinear instability of periodic Bernstein-Greene-Kruskal(BGK waves. Starting
More informationA Kinetic Theory of Planar Plasma Sheaths Surrounding Electron Emitting Surfaces
A Kinetic Theory of Planar Plasma Sheaths Surrounding Electron Emitting Surfaces J. P. Sheehan1, I. Kaganovich2, E. Barnat3, B. Weatherford3, H. Wang2, 4 1 2 D. Sydorenko, N. Hershkowitz, and Y. Raitses
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
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 informationRobustness for a Liouville type theorem in exterior domains
Robustness for a Liouville type theorem in exterior domains Juliette Bouhours 1 arxiv:1207.0329v3 [math.ap] 24 Oct 2014 1 UPMC Univ Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005, Paris,
More informationBeams and magnetized plasmas
Beams and magnetized plasmas 1 Jean-Pierre BOEUF LAboratoire PLAsma et Conversion d Energie LAPLACE/ CNRS, Université Paul SABATIER, TOULOUSE Beams and magnetized plasmas 2 Outline Ion acceleration and
More informationEDP with strong anisotropy : transport, heat, waves equations
EDP with strong anisotropy : transport, heat, waves equations Mihaï BOSTAN University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Nachos team INRIA Sophia Antipolis, 3/07/2017 Main goals Effective
More informationKinetic Theory of Instability-Enhanced Collisions and Its Application to Langmuir s Paradox and the Multi-Species Bohm Criterion
Kinetic Theory of Instability-Enhanced Collisions and Its Application to Langmuir s Paradox and the Multi-Species Bohm Criterion Scott D. Baalrud in collaboration with Chris C. Hegna and James D. Callen
More information1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),
Real Analysis 2, Math 651, Spring 2005 April 26, 2005 1 Real Analysis 2, Math 651, Spring 2005 Krzysztof Chris Ciesielski 1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer
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 informationA note on the plasma sheath and the Bohm Criterion
A note on the plasma sheath and the Bohm Criterion G.D. Severn Dept. of Physics, University of San Diego, San Diego CA 92110 (Dated: April 6, 2006) PACS numbers: 52.27.Aj, 52.27.Cm The word sheath in connection
More informationIntroduction to Plasma Physics
Mitglied der Helmholtz-Gemeinschaft Introduction to Plasma Physics CERN School on Plasma Wave Acceleration 24-29 November 2014 Paul Gibbon Outline Lecture 1: Introduction Definitions and Concepts Lecture
More informationGyrokinetic simulations with GYSELA: Main current issues in physics & numerics
Gyrokinetic simulations with GYSELA: Main current issues in physics & numerics Y. Sarazin, Y. Asahi 2, N. Bouzat, G. Dif-Pradalier, P. Donnel, C. Ehrlacher, C. Emeriau 3, X. Garbet, Ph. Ghendrih, V. Grandgirard,
More informationModelling of JT-60U Detached Divertor Plasma using SONIC code
J. Plasma Fusion Res. SERIES, Vol. 9 (2010) Modelling of JT-60U Detached Divertor Plasma using SONIC code Kazuo HOSHINO, Katsuhiro SHIMIZU, Tomonori TAKIZUKA, Nobuyuki ASAKURA and Tomohide NAKANO Japan
More informationAn Improved RF-Sheath Boundary Condition and Implications for ICRF Modeling
An Improved RF-Sheath Boundary Condition and Implications for ICRF Modeling J. R. Myra Lodestar Research Corporation, Boulder, CO, USA H. Kohno Kyushu Institute of Technology, Iizuka, Fukuoka, Japan September
More informationPhysique des plasmas radiofréquence Pascal Chabert
Physique des plasmas radiofréquence Pascal Chabert LPP, Ecole Polytechnique pascal.chabert@lpp.polytechnique.fr Planning trois cours : Lundi 30 Janvier: Rappels de physique des plasmas froids Lundi 6 Février:
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 informationOn the locality of parallel transport of heat carrying electrons in the SOL
P1-068 On the locality of parallel transport of heat carrying electrons in the SOL A.V. Chankin* and D.P. Coster Max-Planck-Institut für Pasmaphysik, 85748 Garching, Germany Abstract A continuum Vlasov-Fokker-Planck
More informationanalysis 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 informationModel adaptation in hierarchies of hyperbolic systems
Model adaptation in hierarchies of hyperbolic systems Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France February 15th, 2012 DFG-CNRS Workshop Nicolas Seguin (LJLL, UPMC) 1 / 29 Outline of the
More informationMagnetized ion collection by oblique surfaces including self-consistent drifts: Mach-probes of arbitrary shape.
1 Magnetized ion collection by oblique surfaces including self-consistent drifts: Mach-probes of arbitrary shape I H Hutchinson Plasma Science and Fusion Center and and Engineering Department MIT APS DPP
More informationCharacteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges
Plasma Science and Technology, Vol.6, No.6, Jun. 204 Characteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges M. M. HATAMI, A. R. NIKNAM 2 Physics Department
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 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 information4 Hilbert spaces. The proof of the Hilbert basis theorem is not mathematics, it is theology. Camille Jordan
The proof of the Hilbert basis theorem is not mathematics, it is theology. Camille Jordan Wir müssen wissen, wir werden wissen. David Hilbert We now continue to study a special class of Banach spaces,
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationL-H transitions driven by ion heating in scrape-off layer turbulence (SOLT) model simulations
L-H transitions driven by ion heating in scrape-off layer turbulence (SOLT) model simulations D.A. Russell, D.A. D Ippolito and J.R. Myra Research Corporation, Boulder, CO, USA Presented at the 015 Joint
More informationThe Invalidity of a Mach Probe Model
The Invalidity of a Mach Probe Model I. H. Hutchinson Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA, USA Abstract Despite its recent application to the interpretation
More informationMultiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions
Multiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions March 6, 2013 Contents 1 Wea second variation 2 1.1 Formulas for variation........................
More informationPRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING
PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING Second Edition MICHAEL A. LIEBERMAN ALLAN J, LICHTENBERG WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC PUBLICATION CONTENTS PREFACE xrrii PREFACE
More informationKineMa. Numerical Modeling of Kinetic Magnetized Plasmas. April 3rd-7th 2017 IESC, Cargèse
KineMa Numerical Modeling of Kinetic Magnetized Plasmas April 3rd-7th 2017 IESC, Cargèse Program 2 Abstracts of lectures Josh Burby (Courant Institute, New York) Discrete collisionless kinetics: how far
More informationL(q, q) = m 2 q2 V (q) 2 m + V (q)
Lecture 7 Phase Space, Part 1 MATH-GA 71.1 Mechanics 1 Phase portraits 1.1 One dimensional system Consider the generic one dimensional case of a point mass m described by a generalized coordinate q and
More informationModeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range
Modeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range Dr.-Ing. Frank H. Scharf CST of America What is a plasma? What is a plasma? Often referred to as The fourth
More informationOn a Class of Multidimensional Optimal Transportation Problems
Journal of Convex Analysis Volume 10 (2003), No. 2, 517 529 On a Class of Multidimensional Optimal Transportation Problems G. Carlier Université Bordeaux 1, MAB, UMR CNRS 5466, France and Université Bordeaux
More informationAPPLICATIONS OF DIFFERENTIABILITY IN R n.
APPLICATIONS OF DIFFERENTIABILITY IN R n. MATANIA BEN-ARTZI April 2015 Functions here are defined on a subset T R n and take values in R m, where m can be smaller, equal or greater than n. The (open) ball
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 informationContents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments
Contents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments IEC:! IEC at UIUC modified into a space thruster.! IEC has
More informationplasmas Lise-Marie Imbert-Gérard, Bruno Després. July 27th 2011 Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris.
Lise-Marie Imbert-Gérard, Bruno Després. Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris. July 27th 2011 1 Physical and mathematical motivations 2 approximation of the solutions 3 4 Plan
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
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 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 informationOn the harmonic Boltzmannian waves in laser-plasma interaction
On the harmonic Boltzmannian waves in laser-plasma interaction Mihai Bostan 1 and Simon Labrunie 1 Laboratoire de Mathématiques de Besançon, UM CNS 663, Université de Franche-Comté, 16 route de Gray, 53
More informationThe Physics of Field-Aligned Currents
The Physics of Field-Aligned Currents Andrew N. Wright UNIVERSITY OF ST ANDREWS Magnetospheric Current Circuit There is a rich structure of currents flowing parallel (j ) and perpendicular (j ) to the
More informationEXISTENCE OF GLOBAL WEAK SOLUTIONS TO A HYBRID VLASOV-MHD MODEL FOR MAGNETIZED PLASMAS
EXISTENCE OF GLOBAL WEAK SOLUTIONS TO A HYBRID VLASOV-MHD MODEL FOR MAGNETIZED PLASMAS BIN CHENG, ENDRE SÜLI, AND CESARE TRONCI Abstract. We prove the global-in-time existence of large-data finite-energy
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 informationLectures on basic plasma physics: Introduction
Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters
More informationHigh-order energy conserving, (discontinuous) finite-element algorithms for (gyro) kinetic simulations of plasmas
High-order energy conserving, (discontinuous) finite-element algorithms for (gyro) kinetic simulations of plasmas Ammar H. Hakim 1 G. Hammett 1 E. Shi 1,2 1 Princeton Plasma Physics Laboratory, Princeton,
More informationLecture 8 - Carrier Drift and Diffusion (cont.) September 21, 2001
6.720J/3.43J - Integrated Microelectronic Devices - Fall 2001 Lecture 8-1 Lecture 8 - Carrier Drift and Diffusion (cont.) September 21, 2001 Contents: 1. Non-uniformly doped semiconductor in thermal equilibrium
More informationMATHS 730 FC Lecture Notes March 5, Introduction
1 INTRODUCTION MATHS 730 FC Lecture Notes March 5, 2014 1 Introduction Definition. If A, B are sets and there exists a bijection A B, they have the same cardinality, which we write as A, #A. If there exists
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationIf g is also continuous and strictly increasing on J, we may apply the strictly increasing inverse function g 1 to this inequality to get
18:2 1/24/2 TOPIC. Inequalities; measures of spread. This lecture explores the implications of Jensen s inequality for g-means in general, and for harmonic, geometric, arithmetic, and related means in
More informationDevelopment of LH wave fullwave simulation based on FEM
Development of LH wave fullwave simulation based on FEM S. Shiraiwa and O. Meneghini on behalf of LHCD group of Alacator C-Mod PSFC, MIT 2010/03/10 San-Diego, CA Special acknowledgements : R. Parker, P.
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory CENTRAL PEAKING OF MAGNETIZED GAS DISCHARGES Francis F. Chen and Davide Curreli LTP-1210 Oct. 2012 Electrical Engineering Department Los Angeles, California
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationEP2Plus: a hybrid plasma. plume/spacecraft. interaction code. F. Cichocki, M. Merino, E. Ahedo
EP2Plus: a hybrid plasma plume/spacecraft interaction code F. Cichocki, M. Merino, E. Ahedo 24 th SPINE meeting ESTEC, Noordwijk, October 23 rd, 2017 Contents Introduction to EP2PLUS Overall structure
More informationOn Dust Particle Dynamics in Tokamak Edge Plasma
On Dust Particle Dynamics in Tokamak Edge Plasma Sergei Krasheninnikov University of California, San Diego, USA With contributions from Y. Tomita 1, R. D. Smirnov, and R. K. Janev 1 1 National Institute
More informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More informationSensors Plasma Diagnostics
Sensors Plasma Diagnostics Ken Gentle Physics Department Kenneth Gentle RLM 12.330 k.gentle@mail.utexas.edu NRL Formulary MIT Formulary www.psfc.mit.edu/library1/catalog/ reports/2010/11rr/11rr013/11rr013_full.pdf
More informationContinuum Electrostatics for Ionic Solutions with Nonuniform Ionic Sizes
Continuum Electrostatics for Ionic Solutions with Nonuniform Ionic Sizes Bo Li Department of Mathematics and Center for Theoretical Biological Physics (CTBP) University of California, San Diego, USA Supported
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 informationAnomalous 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 informationLecture 7 - Carrier Drift and Diffusion (cont.) February 20, Non-uniformly doped semiconductor in thermal equilibrium
6.720J/3.43J - Integrated Microelectronic Devices - Spring 2007 Lecture 7-1 Lecture 7 - Carrier Drift and Diffusion (cont.) February 20, 2007 Contents: 1. Non-uniformly doped semiconductor in thermal equilibrium
More informationElectromagneic Waves in a non- Maxwellian Dusty Plasma
Electromagneic Waves in a non- Maxwellian Dusty Plasma Nazish Rubab PhD student, KF University Graz IWF-OEAW Graz 26 January, 2011 Layout Dusty Plasma Non-Maxwellian Plasma Kinetic Alfven Waves Instability
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory Equilibrium theory for plasma discharges of finite length Francis F. Chen and Davide Curreli LTP-6 June, Electrical Engineering Department Los Angeles, California
More informationThe Stationary Boltzmann equation for a two component gas in the slab.
The Stationary oltzmann equation for a two component gas in the slab. Stéphane rull* Abstract The stationary oltzmann equation for hard forces in the context of a two component gas is considered in the
More informationFull Wave Propagation Modelling in View to Integrated ICRH Wave Coupling/RF Sheaths Modelling
EUROFUSION CP(15)03/04 J. Jacquot et al. Full Wave Propagation Modelling in View to Integrated ICRH Wave Coupling/RF Sheaths Modelling (27th April 2015 29th April 2015) Lake Arrowhead, California, USA
More informationAngular momentum preserving CFD on general grids
B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration
More informationMODERN PHYSICS OF PLASMAS (19 lectures)
UNIT OF STUDY OUTLINE (PHYS 3021, 3921, 3024, 3924, 3025, 3925) MODERN PHYSICS OF PLASMAS (19 lectures) Course coordinator and principal lecturer: Dr Kostya (Ken) Ostrikov Lecturer (normal student stream,
More informationPhysics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection
1 ITR/P1-37 Physics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection J.P. Boeuf a, G. Fubiani a, G. Hagelaar a, N. Kohen a, L. Pitchford a, P. Sarrailh a, and A. Simonin
More informationTHEOREMS, ETC., FOR MATH 515
THEOREMS, ETC., FOR MATH 515 Proposition 1 (=comment on page 17). If A is an algebra, then any finite union or finite intersection of sets in A is also in A. Proposition 2 (=Proposition 1.1). For every
More informationON THE VLASOV-MAXWELL SYSTEM WITH A STRONG EXTERNAL MAGNETIC FIELD
ON THE VLASOV-MAXWELL SYSTEM WITH A STONG EXTENAL MAGNETIC FIELD Francis Filbet, Tao Xiong, Eric Sonnendr ucker To cite this version: Francis Filbet, Tao Xiong, Eric Sonnendr ucker. ON THE VLASOV-MAXWELL
More informationPlasma waves in the fluid picture I
Plasma waves in the fluid picture I Langmuir oscillations and waves Ion-acoustic waves Debye length Ordinary electromagnetic waves General wave equation General dispersion equation Dielectric response
More informationExistence result for the coupling problem of two scalar conservation laws with Riemann initial data
Existence result for the coupling problem of two scalar conservation laws with Riemann initial data Benjamin Boutin, Christophe Chalons, Pierre-Arnaud Raviart To cite this version: Benjamin Boutin, Christophe
More informationProjection Theorem 1
Projection Theorem 1 Cauchy-Schwarz Inequality Lemma. (Cauchy-Schwarz Inequality) For all x, y in an inner product space, [ xy, ] x y. Equality holds if and only if x y or y θ. Proof. If y θ, the inequality
More informationParallel Kelvin-Helmholtz instability in edge plasma
Journal of Physics: Conference Series OPEN ACCESS Parallel Kelvin-Helmholtz instability in edge plasma To cite this article: H Guillard et al 2014 J. Phys.: Conf. Ser. 561 012009 View the article online
More informationPositive ion flux from a low-pressure electronegative discharge
Plasma Sources Sci. Technol. 8 (1999) 457 462. Printed in the UK PII: S0963-0252(99)03994-8 Positive ion flux from a low-pressure electronegative discharge T E Sheridan, P Chabert and R W Boswell Space
More information