Quasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths
|
|
- Carmella Floyd
- 6 years ago
- Views:
Transcription
1 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 1st Meeting of Young Researchers in PDEs Yonsei University, August 18-19, 2016
2 Euler-Poisson equations The motion of the positive ions in a plasma is governed by ρ t + (ρu) = 0, (ρu) t + (ρu u) + K ρ = ρ φ, ε φ = ρ e φ, ρ(x, t), u(x, t) and φ(x, t) represent the density, velocity of the positive ions and the electrostatic potential, respectively. K = temperature of positive ion ε is the ratio of the Debye length λd to the radius of an inner sphere R, that is, ε = λ 2 D/R 2. (Lab setting λ D 10 4 )
3 Plasma physics Plasma (in Greek anything formed ) is one of the four fundamental states of matter. Like gas, plasma does not have a definite shape or a definite volume unless enclosed in a container. Unlike gas, under the influence of a magnetic field, it may form structures such as filaments, beams and double layers. The presence of a significant number of charge carriers makes plasma electrically conductive so that it responds strongly to electromagnetic fields
4 Plasma physics Quasi-neutrality (quasi, from the Latin, as if, resembling ) describes the apparent charge neutrality of a plasma overall, while at smaller scales, the positive and negative charges making up the plasma, may give rise to charged regions and electric fields. Since electrons are very mobile, plasmas are excellent conductors of electricity, and any charges that develop are readily neutralized, and in many cases, plasmas can be treated as being electrically neutral. The distance over which quasi-neutrality may break down, is often described by the Debye length (or Debye sphere), and varies according to the physical characteristics of the plasma. The Debye length is typically less than a millimetre (ie. charged regions will not exceed a millimetre), in plasmas found in fluorescent light tubes, tokamaks (used in fusion research), and the ionosphere.
5 From the two-fluid model to our model Full two-fluid model reads as tn ± + (n ±v ±) = 0, n ±m ±( tv ± + (v ± )v ±) + T ± n ± = ±en ± φ, ε φ = 4πe(n + n ) After non-dimensionalization together with Boltzmann-Maxwell relation (m m +), one obtains EP for ion dynamics: ρ t + (ρu) = 0, (ρu) t + (ρu u) + K ρ = ρ φ, ε φ = ρ e φ.
6 Euler equations of compressible fluids: Euler equations of compressible fluids is given by ρ t + (ρu) = 0, ρ(u t + (u )u) + p = 0, S t + (u )S = 0, where P(ρ, S) = Aρ γ e S is the pressure and γ > 1 is the adiabatic index. σ := P ρ( ρ, S) = Aγ ρ γ 1 e S is the speed of sound. Note that the maximum speed of propagation of the wave of a smooth disturbance is governed by σ. Heuristically, the linear approximation of the Euler equations near the constant state (ρ, u, S) = ( ρ, 0, S) is the wave equation: ρ tt σ 2 ρ = 0. (T. Sideris 1985) C 1 solutions may blow up in finite-time. Hyperbolic + Nonlinear blow up
7 Dispersion relation EP for ion dynamics: ρ t + (ρu) = 0, (ρu) t + (ρu u) + K ρ γ = ρ φ, ε φ = ρ e φ. Linearizing about the constant state (ρ, u, φ) = (1, 0, 0), we have ρ t + u = 0, u t + K ρ = φ, ε φ = ρ + φ. Furthermore we obtain the wave equation with nonlocal term: ρ tt K ρ = (( 1 + ) 1 ρ). Dispersion relation: K K ξ 2 λ = i ξ 1 + ξ 2
8 Plasma sheaths Physical situation: a negatively charged material is immersed in a plasma and the plasma contacts with its surface. The plasma sheath is the transition from a plasma to a solid surface, where the positive ions outnumber the electrons (Langmuir, 1929) Plasma sheath is emerged. The interaction of electron and positive ion space charges in cathode sheaths, Phys. Rev. 33 (1929), pp (Bohm, 1949) Bohm criterion Minimum ionic kinetic energy for a stable sheath, in The characteristics of electrical discharges in magnetic fields, A. Guthrie and R.K.Wakerling eds., McGraw-Hill, New York, (1949),
9 Plasma sheaths The plasma sheath is a layer in a plasma which has a greater density of positive ions, and hence an overall excess positive charge, that balances an opposite negative charge on the surface of a material with which it is in contact.
10 Program: hydrodynamic description of the formation of plasma sheaths via the Euler-Poisson system Plasma sheaths (sharp transition layers near the boundary, stationary boundary layer solutions) Bohm criterion (physical criterion) Existence of the stationary solution (plasma sheath) Euler-Poisson system its quasi-neutral limiting equations (Euler system) Stability (time-asymptotic) vs blow-up Local existence (hyperbolic initial-boundary value problem) Global existence (uniform estimates) or Spectral analysis towards blow-up Quasi-neutral limit Stationary and time-dependent problem Uniform decay estimate (in terms of ε = Debye length) Correctors (information on boundary layers)
11 Euler-Poisson system and its quasi-neutral limiting eqns E-P system: ρ ε t + (ρ ε u ε ) = 0, (ρ ε u ε ) t + (ρ ε u ε u ε ) + K ρ ε = ρ ε φ ε, ε φ ε = ρ ε e φε. The associated quasi-neutral limiting equations (ε = 0) ρ 0 t + (ρ 0 u 0 ) = 0, u 0 t + u 0 u 0 + (K + 1) log ρ 0 = 0.
12 Euler-Poisson system in an annular domain Annulus domain Ω: Ω := {(x 1, x 2, x 3) R 3 1 < where L > 0 is the width of the annulus domain. We consider E-P in Ω: wehre x Ω and t 0. Initial conditions: Boundary conditions: where 1, u + and φ b are constants. ρ t + (ρu) = 0, 3 j=1 x 2 i < 1 + L}, (ρu) t + (ρu u) + K ρ = ρ φ, ε φ = ρ e φ, (ρ, u)(0, x) = (ρ 0, u 0)(x), inf x Ω ( u0 2 Kρ 0 )(x) > 0, inf ρ0(x) > 0, x Ω sup u 0 ν(x) < 0, x Ω (ρ, u)(t, 1 + L, ω) = (1, u +ω), ω := x x, φ(t, 1, ω) = φ b, φ(t, 1 + L, ω) = 0.
13 Existence of stationary radial solutions Condition (Bohm s Criterion for the annulus domain) u 2 + > a (1 + K), where a > 1 is a larger root of log(1 + L) x + log x = 0. Theorem (Existence for the stationary Q-N-L) Let u 2 + > a (1 + K) hold. Then the SLP has a unique solution ( ρ 0, ũ 0 ) C (I ) C(I ). Moreover, ρ 0 is a monotone decreasing function with the bound 1 = ρ 0 (1 + L) ρ 0 (r) ρ 0 (1). Theorem (Existence for the stationary E-P eqns.) There are constants ε 2 > 0 and δ 1 > 0 such that if ε ε 2 and φ b log ρ 0 (1) + δ 1, the SEP has a unique solution φ C 2 ([1, 1 + L]) satisfying φ(r) G(r, ρ 0 (r)) r [1, 1 + L] φ(r) log ρ 0 (r) + δ 1 r [1, 1 + L]
14 Quasi-neutral limit in the presence of plasma sheaths Theorem (Quasi-neutral limit) Under the same assumptions, the difference ( ρ ρ 0 θ 0 φ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) satisfies the decay estimates ( ρ ρ 0 θ 0 ρ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) L 2 Cε3/4, ( ρ ρ 0 θ 0 ρ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) H 1 Cε1/4, ( ρ ρ 0 θ 0 ρ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) L Cε1/2, where C > 0 is a constant independent of ε. Proposition (pointwise estimates for the correctors) Under the same assumptions, for any r I, there hold φ b φ ( 0 (1) e C 0 /ε(r 1) e ) C 0 /εl θ 0 φ (r) φ b φ 0 (1) e c 0 /ε(r 1), where c, C, c 0 and C 0 are some positive constants independent of ε. This implies that the thickness of the boundary layer is of order O( ε) = O(λ D ). By a direct calculation, it seems that the best L 2 and H 1 estimates for φ φ 0 are φ φ 0 L 2 Cε 1/4 and φ φ 0 H 1 Cε 1/4, respectively.
15 Numerical experiments in radial symmetry setting (Jung-K-Suzuki 2016)
16 Existence for the stationary Euler-Poisson system Condition Recall the physical assumptions: (Bohm criterion) u+ 2 > a (1 + K) φ b G(1, ρ 0 (1)) Theorem There are constants ε 2 > 0 and δ 1 > 0 such that if ε ε 2 and φ b log ρ 0 (1) + δ 1, the SEP has a unique solution φ C 2 ([1, 1 + L]) satisfying φ(r) G(r, ρ 0 (r)) r [1, 1 + L] φ(r) log ρ 0 (r) + δ 1 r [1, 1 + L]
17 Quasi-neutral limit problem Study the quasi-neutral limit problem (ε 0): For (ρ = ρ ε, u = u ε, φ = φ ε ), (P ε ) (ρu) r = 2 r ρu, (ρu 2 + Kρ) r = ρφ r 2 r ρu2, εr 2 (r 2 φ r ) r = (ρ e φ ), φ = φ b at r = 1, (ρ, u, φ) = (1, u +, 0) at r = 1 + L. (0.1) For (ρ = ρ 0, u = u 0, φ = φ 0 ), (P 0 ) (ρu) r = 2 r ρu, (ρu 2 + Kρ) r = ρφ r 2 r ρu2, 0 = ρ e φ, (ρ, u) = (1, u +) at r = 1 + L. (0.2) discrepancy bet n the boundary data φ b and log ρ(1) causes a sharp transition layer
18 Quasi-neutral limit in the presence of plasma sheaths Theorem (Quasi-neutral limit) Under the same assumptions, the difference ( ρ ρ 0 θ 0 φ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) satisfies the decay estimates ( ρ ρ 0 θ 0 ρ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) L 2 Cε3/4, ( ρ ρ 0 θ 0 ρ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) H 1 Cε1/4, ( ρ ρ 0 θ 0 ρ, ũ ũ0 θ 0 u, φ φ 0 θ 0 φ ) L Cε1/2, where C > 0 is a constant independent of ε, and θ 0 is a corrector which is a solution to the dominating equations: εθ 0 rr = g(r, φ0 + θ 0 ) e φ0 θ 0, θ 0 (1) = φ b + log ρ 0 (1), θ 0 (1 + L) = 0. RMK. By a direct calculation, it seems that the best L 2 and H 1 estimates for φ φ 0 are φ φ 0 L 2 Cε 1/4 and φ φ 0 H 1 Cε 1/4, respectively.
19 Pointwise estimates for the correctors Proposition (pointwise estimates) Under the same assumptions, for any r I, there hold ( φ b φ 0 (1) e C 0 /ε(r 1) e ) C 0 /εl θ 0 φ (r) φ b φ 0 (1) e c 0 /ε(r 1), where c, C, c 0 and C 0 are some positive constants independent of ε. This implies that the thickness of the boundary layer is of order O( ε) = O(λ D ).
20 Numerical Experiments for Quasi-neutral limit problem
21 Sketch of proof To find the boundary layer structures, we decompose φ = φ 0 + θ. ε d 2 dr ( φ 0 + θ) r ε d dr ( φ 0 + θ) = g(r, φ 0 + θ) e φ 0 θ. Since φ 0 is independent of ε, the terms εd 2 φ0 /dr 2 and εd φ 0 /dr are small and we drop them. In addition, to zoom the behavior of θ near r = 1, using a stretched variable r = (r 1)/ ε, we obtain d 2 d r θ + 2 d ε 2 ε r + 1 d r θ = g(r, φ 0 + θ) e φ 0 θ. (0.3) Then we look for the leading asymptotic term θ 0 φ with θ θ0 φ + εθ 1 φ. Substituting this expansion for θ in (0.3) and using the Taylor expansion for g(r, φ) and e φ about φ = φ 0 + θ 0 φ, we find that d 2 d r 2 θ0 φ = g(r, φ 0 + θ 0 φ ) φ 0 θ 0 e φ + O( ε). At the leading order, i.e. O(1), we obtain the equation for the boundary layers.
22 Sketch of proof Define a corrector function by θ(r) := (φ b φ 0 (1)) exp ( r 1 ε ) χ(r). Then θ 0 φ θ satisfies ε(θ 0 φ θ) = g(r, φ 0 + θ 0 φ ) φ 0 e θ φ 0 θ R 1, (0.4a) (θ 0 φ θ)(1) = (θ 0 φ θ)(1 + L) = 0, (0.4b) where R 1 := ( 2 εχ + εχ ) (φ b φ 0 (1))e (r 1)/ ε. On the other hand, by a straightforward calculation, one can estimate θ and R 1 as where C > 0 is a constant independent of ε. θ L 2 + R 1 L 2 Cε 1/4, (0.5)
23 By L 2 estimate in a straightforward fashion, we have 1+L 1+L ( ) ε (θ 0 φ θ) 2 dr + θ 0 φ g(r, φ 0 + θ 0 φ ) φ 0 e θ φ 0 dr L ( ) = θ g(r, φ 0 + θ 0 φ ) φ 0 e θ φ 0 g(r, φ 0 ) + e φ 0 1+L dr + ( θ + R 1)(θ 0 φ θ)dr. 1 1 This gives ε (θ 0 φ θ) 2 L 2 + θ 0 φ 2 L 2 Cε 1/4. It is easy to check that θ L 2 Cε 1/4 Combining these, we have ε (θ 0 φ ) 2 L 2 + θ 0 φ 2 L 2 Cε 1/4.
24 Moreover, define Then it satisfies w := φ φ 0 θ 0 φ. ε(r 2 w ) = r 2 {K(r, φ, φ 0 + θ 0 φ ) + J(r, φ, φ 0 + θ 0 φ )}w + R2, where R 2 = ε{2r( φ 0 + θ 0 φ ) + r 2 ( φ 0 ) }. w(1) = w(1 + L) = 0, Note that R 2 L 2 Cε 3/4. Again, by energy estimate, we arrive at ε w 2 L 2 + w 2 L 2 C R 2 2 L 2 Cε 3/2.
25 Thank you for your attention
Lecture Note 1. 99% of the matter in the universe is in the plasma state. Solid -> liquid -> Gas -> Plasma (The fourth state of matter)
Lecture Note 1 1.1 Plasma 99% of the matter in the universe is in the plasma state. Solid -> liquid -> Gas -> Plasma (The fourth state of matter) Recall: Concept of Temperature A gas in thermal equilibrium
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 informationMagnetohydrodynamic waves in a plasma
Department of Physics Seminar 1b Magnetohydrodynamic waves in a plasma Author: Janez Kokalj Advisor: prof. dr. Tomaž Gyergyek Petelinje, April 2016 Abstract Plasma can sustain different wave phenomena.
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 informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationFigure 1.1: Ionization and Recombination
Chapter 1 Introduction 1.1 What is a Plasma? 1.1.1 An ionized gas A plasma is a gas in which an important fraction of the atoms is ionized, so that the electrons and ions are separately free. When does
More informationGlobal existence for the ion dynamics in the Euler-Poisson equations
Global existence for the ion dynamics in the Euler-Poisson equations Yan Guo (Brown U), Benoît Pausader (Brown U). FRG Meeting May 2010 Abstract We prove global existence for solutions of the Euler-Poisson/Ion
More informationFundamentals of Plasma Physics
Fundamentals of Plasma Physics Definition of Plasma: A gas with an ionized fraction (n i + + e ). Depending on density, E and B fields, there can be many regimes. Collisions and the Mean Free Path (mfp)
More informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationFundamentals of Plasma Physics Basic concepts
Fundamentals of Plasma Physics Basic concepts APPLAuSE Instituto Superior Técnico Instituto de Plasmas e Fusão Nuclear Vasco Guerra Since the dawn of Mankind men has tried to understand plasma physics...
More informationA Variational Analysis of a Gauged Nonlinear Schrödinger Equation
A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological
More informationLecture 2. Introduction to plasma physics. Dr. Ashutosh Sharma
Preparation of the concerned sectors for educational and R&D activities related to the Hungarian ELI project Ion acceleration in plasmas Lecture 2. Introduction to plasma physics Dr. Ashutosh Sharma Zoltán
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 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 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 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 informationBielefeld Course on Nonlinear Waves - June 29, Department of Mathematics University of North Carolina, Chapel Hill. Solitons on Manifolds
Joint work (on various projects) with Pierre Albin (UIUC), Hans Christianson (UNC), Jason Metcalfe (UNC), Michael Taylor (UNC), Laurent Thomann (Nantes) Department of Mathematics University of North Carolina,
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 informationSpace Plasma Physics Thomas Wiegelmann, 2012
Space Plasma Physics Thomas Wiegelmann, 2012 1. Basic Plasma Physics concepts 2. Overview about solar system plasmas Plasma Models 3. Single particle motion, Test particle model 4. Statistic description
More informationThe ideal Maxwellian plasma
The ideal Maxwellian plasma Dr. L. Conde Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Plasmas are,... The plasma state of matter may be defined as a
More informationTwo-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation
Mathematics and Computers in Simulation 76 (2007) 3 7 Two-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation S. Baboolal a,, R. Bharuthram
More informationxkcd.com It IS about physics. It ALL is.
xkcd.com It IS about physics. It ALL is. Introduction to Space Plasmas! The Plasma State What is a plasma? Basic plasma properties: Qualitative & Quantitative Examples of plasmas! Single particle motion
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 informationON THE GLOBAL REGULARITY OF SUB-CRITICAL EULER-POISSON EQUATIONS WITH PRESSURE
ON THE GLOBAL REGULARITY OF SUB-CRITICAL EULER-POISSON EQUATIONS WITH PRESSURE EITAN TADMOR AND DONGMING WEI Abstract We prove that the one-dimensional Euler-Poisson system driven by the Poisson forcing
More informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
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 informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationxkcd.com It IS about physics. It ALL is.
xkcd.com It IS about physics. It ALL is. Introduction to Space Plasmas The Plasma State What is a plasma? Basic plasma properties: Qualitative & Quantitative Examples of plasmas Single particle motion
More informationBoundary Conditions for the Child Langmuir Sheath Model
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 6, DECEMBER 2000 2207 Boundary Conditions for the Child Langmuir Sheath Model Mikhail S. Benilov Abstract A collision-free space-charge sheath formed by
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationAST 553. Plasma Waves and Instabilities. Course Outline. (Dated: December 4, 2018)
AST 553. Plasma Waves and Instabilities Course Outline (Dated: December 4, 2018) I. INTRODUCTION Basic concepts Waves in plasmas as EM field oscillations Maxwell s equations, Gauss s laws as initial conditions
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationGlobal existence and asymptotic behavior of the solutions to the 3D bipolar non-isentropic Euler Poisson equation
Nonlinear Analysis: Modelling and Control, Vol., No. 3, 35 33 ISSN 139-5113 http://dx.doi.org/1.15388/na.15.3.1 Global existence and asymptotic behavior of the solutions to the 3D bipolar non-isentropic
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More information20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:
Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationNonlinear electrostatic structures in unmagnetized pair-ion (fullerene) plasmas
Nonlinear electrostatic structures in unmagnetized pair-ion (fullerene) plasmas S. Mahmood Theoretical Plasma Physics Division, PINSTECH Islamabad Collaborators: H. Saleem National Center for Physics,
More informationAnalytic Expression For the Electric Potential in the Plasma Sheath
884 IEEE TRANSACTIONS ON PLASMA SCIENCE. VOL. 17. NO. 6. DECEMBER 1989 Analytic Expression For the Electric Potential in the Plasma Sheath TERRENCE E. SHERIDAN, JR. AND JOHN A GOREE Abstract-An expression
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
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 information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
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 informationAustralian Journal of Physics
CSIRO PUBLISHING Australian Journal of Physics Volume 50, 1997 CSIRO Australia 1997 A journal for the publication of original research in all branches of physics www.publish.csiro.au/journals/ajp All enquiries
More information. (70.1) r r. / r. Substituting, we have the following equation for f:
7 Spherical waves Let us consider a sound wave in which the distribution of densit velocit etc, depends only on the distance from some point, ie, is spherically symmetrical Such a wave is called a spherical
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 informationOn Chern-Simons-Schrödinger equations including a vortex point
On Chern-Simons-Schrödinger equations including a vortex point Alessio Pomponio Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Workshop in Nonlinear PDEs Brussels, September 7
More informationGeneral relativistic computation of shocks in accretion disc and bipolar jets
General relativistic computation of shocks in accretion disc and bipolar jets Rajiv Kumar Co-author Dr. Indranil Chattopadhyay Aryabhatta Research Institute of observational sciences (ARIES) rajiv.k@aries.res.in
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 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 informationElectromagnetic Field Theory (EMT)
Electromagnetic Field Theory (EMT) Lecture # 9 1) Coulomb s Law and Field Intensity 2) Electric Fields Due to Continuous Charge Distributions Line Charge Surface Charge Volume Charge Coulomb's Law Coulomb's
More informationCluster fusion in a high magnetic field
Santa Fe July 28, 2009 Cluster fusion in a high magnetic field Roger Bengtson, Boris Breizman Institute for Fusion Studies, Fusion Research Center The University of Texas at Austin In collaboration with:
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
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 informationA quantum heat equation 5th Spring School on Evolution Equations, TU Berlin
A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationNONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE
NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE M.A. Lieberman Department of Electrical Engineering and Computer Sciences University of California Berkeley, CA 94720 Collaborators:
More informationNonlinear Wave Propagation in 1D Random Media
Nonlinear Wave Propagation in 1D Random Media J. B. Thoo, Yuba College UCI Computational and Applied Math Seminar Winter 22 Typeset by FoilTEX Background O Doherty and Anstey (Geophysical Prospecting,
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 informationPHY102 Electricity Topic 3 (Lectures 4 & 5) Gauss s Law
PHY1 Electricity Topic 3 (Lectures 4 & 5) Gauss s Law In this topic, we will cover: 1) Electric Flux ) Gauss s Law, relating flux to enclosed charge 3) Electric Fields and Conductors revisited Reading
More informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationAssessment of fluctuation-induced and wall-induced anomalous electron transport in HET
Assessment of fluctuation-induced and wall-induced anomalous electron transport in HET IEPC-2015-418 Presented at Joint Conference of 30th International Symposium on Space Technology and Science 34th International
More informationWorkshop on Compressible Navier-Stokes Systems and Related Problems (I) March 5-10, 2018 TITLE & ABSTRACT
Workshop on Compressible Navier-Stokes Systems and Related Problems (I) March 5-10, 2018 TITLE & ABSTRACT (Last updated: 6 March 2018) Classification of asymptotic states for radially symmetric solutions
More informationLocal Density Approximation for the Almost-bosonic Anyon Gas. Michele Correggi
Local Density Approximation for the Almost-bosonic Anyon Gas Michele Correggi Università degli Studi Roma Tre www.cond-math.it QMATH13 Many-body Systems and Statistical Mechanics joint work with D. Lundholm
More informationWebster s horn model on Bernoulli flow
Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular
More informationMULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER
More informationOn the leapfrogging phenomenon in fluid mechanics
On the leapfrogging phenomenon in fluid mechanics Didier Smets Université Pierre et Marie Curie - Paris Based on works with Robert L. Jerrard U. of Toronto) CIRM, Luminy, June 27th 2016 1 / 22 Single vortex
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 informationMethod of infinite system of equations for problems in unbounded domains
Method of infinite system of equations for problems in unbounded domains Dang Quang A Institute of Information Technology 18 Hoang Quoc Viet, Cau giay, Hanoi, Vietnam Innovative Time Integration, May 13-16,
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 informationWorked Examples Set 2
Worked Examples Set 2 Q.1. Application of Maxwell s eqns. [Griffiths Problem 7.42] In a perfect conductor the conductivity σ is infinite, so from Ohm s law J = σe, E = 0. Any net charge must be on the
More informationDepartment of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions
Department of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions Friday May 15th, 2014, 14-17h Examiners: Prof. J. Cline, Prof. H. Guo, Prof. G. Gervais (Chair), and Prof. D. Hanna INSTRUCTIONS
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 informationξ,i = x nx i x 3 + δ ni + x n x = 0. x Dξ = x i ξ,i = x nx i x i x 3 Du = λ x λ 2 xh + x λ h Dξ,
1 PDE, HW 3 solutions Problem 1. No. If a sequence of harmonic polynomials on [ 1,1] n converges uniformly to a limit f then f is harmonic. Problem 2. By definition U r U for every r >. Suppose w is a
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 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 informationElectric Flux. To investigate this, we have to understand electric flux.
Problem 21.72 A charge q 1 = +5. nc is placed at the origin of an xy-coordinate system, and a charge q 2 = -2. nc is placed on the positive x-axis at x = 4. cm. (a) If a third charge q 3 = +6. nc is now
More informationA Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More informationChapter 3. Head-on collision of ion acoustic solitary waves in electron-positron-ion plasma with superthermal electrons and positrons.
Chapter 3 Head-on collision of ion acoustic solitary waves in electron-positron-ion plasma with superthermal electrons and positrons. 73 3.1 Introduction The study of linear and nonlinear wave propagation
More informationNeutralized Drift Compression
Neutralized Drift Compression D. R. Welch, D. V. Rose Mission Research Corporation Albuquerque, NM 87110 S. S. Yu Lawrence Berkeley National Laboratory Berkeley, CA C. L. Olson Sandia National Laboratories
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationOn 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 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 informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More information2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves
More informationEnhanced resolution in structured media
Enhanced resolution in structured media Eric Bonnetier w/ H. Ammari (Ecole Polytechnique), and Yves Capdeboscq (Oxford) Outline : 1. An experiment of super resolution 2. Small volume asymptotics 3. Periodicity
More informationGradient estimates and global existence of smooth solutions to a system of reaction-diffusion equations with cross-diffusion
Gradient estimates and global existence of smooth solutions to a system of reaction-diffusion equations with cross-diffusion Tuoc V. Phan University of Tennessee - Knoxville, TN Workshop in nonlinear PDES
More informationMagnetized Target Fusion: Insights from Mathematical Modelling
: Insights from Mathematical Modelling, UBC, PhD Candidate - Applied Math Brian Wetton, UBC, Supervisor Michael Ward, UBC, Committee Member Sandra Barsky, General Fusion, Committee Member November 28,
More informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
More informationConfinement of toroidal non-neutral plasma in Proto-RT
Workshop on Physics with Ultra Slow Antiproton Beams, RIKEN, March 15, 2005 Confinement of toroidal non-neutral plasma in Proto-RT H. Saitoh, Z. Yoshida, and S. Watanabe Graduate School of Frontier Sciences,
More informationCharacteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:
AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,
More informationV. Electrostatics Lecture 24: Diffuse Charge in Electrolytes
V. Electrostatics Lecture 24: Diffuse Charge in Electrolytes MIT Student 1. Poisson-Nernst-Planck Equations The Nernst-Planck Equation is a conservation of mass equation that describes the influence of
More informationConfinement of toroidal non-neutral plasma in Proto-RT
Workshop on Physics with Ultra Slow Antiproton Beams, RIKEN, March 15, 2005 Confinement of toroidal non-neutral plasma in Proto-RT H. Saitoh, Z. Yoshida, and S. Watanabe Graduate School of Frontier Sciences,
More informationAsymptotic Behavior of Marginally Trapped Tubes
Asymptotic Behavior of Marginally Trapped Tubes Catherine Williams January 29, 2009 Preliminaries general relativity General relativity says that spacetime is described by a Lorentzian 4-manifold (M, g)
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 information13.1 Ion Acoustic Soliton and Shock Wave
13 Nonlinear Waves In linear theory, the wave amplitude is assumed to be sufficiently small to ignore contributions of terms of second order and higher (ie, nonlinear terms) in wave amplitude In such a
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More information