Quasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths

Size: px
Start display at page:

Download "Quasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths"

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

A note on the plasma sheath and the Bohm Criterion

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

Magnetohydrodynamic waves in a plasma

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

Plasma waves in the fluid picture I

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

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

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy

More information

Figure 1.1: Ionization and Recombination

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

Global existence for the ion dynamics in the Euler-Poisson equations

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

Fundamentals of Plasma Physics

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

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

Fundamentals of Plasma Physics Basic concepts

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

A Variational Analysis of a Gauged Nonlinear Schrödinger Equation

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

Lecture 2. Introduction to plasma physics. Dr. Ashutosh Sharma

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

Chapter 1. Introduction to Nonlinear Space Plasma Physics

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

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

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

Characteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges

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

Waves in plasma. Denis Gialis

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

Bielefeld Course on Nonlinear Waves - June 29, Department of Mathematics University of North Carolina, Chapel Hill. Solitons on Manifolds

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

Minimization formulation of a bi-kinetic sheath

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

Space Plasma Physics Thomas Wiegelmann, 2012

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

The ideal Maxwellian plasma

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

Two-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation

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

xkcd.com It IS about physics. It ALL is.

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

A Kinetic Theory of Planar Plasma Sheaths Surrounding Electron Emitting Surfaces

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

ON THE GLOBAL REGULARITY OF SUB-CRITICAL EULER-POISSON EQUATIONS WITH PRESSURE

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

Numerical Methods for Modern Traffic Flow Models. Alexander Kurganov

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

Dispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson

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

Introduction. Chapter Plasma: definitions

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

Recapitulation: Questions on Chaps. 1 and 2 #A

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

xkcd.com It IS about physics. It ALL is.

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

Boundary Conditions for the Child Langmuir Sheath Model

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

Nonlinear stabilization via a linear observability

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

AST 553. Plasma Waves and Instabilities. Course Outline. (Dated: December 4, 2018)

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

The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.

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

Global existence and asymptotic behavior of the solutions to the 3D bipolar non-isentropic Euler Poisson equation

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

MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT

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

20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:

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

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

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

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

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

Nonlinear electrostatic structures in unmagnetized pair-ion (fullerene) plasmas

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

Analytic Expression For the Electric Potential in the Plasma Sheath

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

Salmon: Lectures on partial differential equations

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

Electromagneic Waves in a non- Maxwellian Dusty Plasma

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

0.2. CONSERVATION LAW FOR FLUID 9

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

MODERN PHYSICS OF PLASMAS (19 lectures)

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

Australian Journal of Physics

Australian 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:

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

Sensors Plasma Diagnostics

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

On Chern-Simons-Schrödinger equations including a vortex point

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

General relativistic computation of shocks in accretion disc and bipolar jets

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

Introduction to Plasma Physics

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

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

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

Electromagnetic Field Theory (EMT)

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

Cluster fusion in a high magnetic field

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

Simple examples of MHD equilibria

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

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

A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin

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

OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES

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

NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE

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

Nonlinear Wave Propagation in 1D Random Media

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

Physique des plasmas radiofréquence Pascal Chabert

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

PHY102 Electricity Topic 3 (Lectures 4 & 5) Gauss s Law

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

TURBULENT TRANSPORT THEORY

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

Euler Equations: local existence

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

Assessment of fluctuation-induced and wall-induced anomalous electron transport in HET

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

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

Local Density Approximation for the Almost-bosonic Anyon Gas. Michele Correggi

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

Webster s horn model on Bernoulli flow

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

MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran

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

On the leapfrogging phenomenon in fluid mechanics

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

Advection, Conservation, Conserved Physical Quantities, Wave Equations

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

Method of infinite system of equations for problems in unbounded domains

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

One dimensional hybrid Maxwell-Boltzmann model of shearth evolution

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

Worked Examples Set 2

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

Department of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions

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

Angular momentum preserving CFD on general grids

Angular 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ξ,

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

Gyrokinetic simulations of magnetic fusion plasmas

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

Beams and magnetized plasmas

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

Electric Flux. To investigate this, we have to understand electric flux.

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

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

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

Neutralized Drift Compression

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

Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations

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

On the Boltzmann equation: global solutions in one spatial dimension

On the Boltzmann equation: global solutions in one spatial dimension On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators

More information

Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas

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

The Euler Equation of Gas-Dynamics

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

Homogenization and error estimates of free boundary velocities in periodic media

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

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

Enhanced resolution in structured media

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

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

Magnetized Target Fusion: Insights from Mathematical Modelling

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

Homogenization of micro-resonances and localization of waves.

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

Confinement of toroidal non-neutral plasma in Proto-RT

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

Characteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:

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

V. Electrostatics Lecture 24: Diffuse Charge in Electrolytes

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

Confinement of toroidal non-neutral plasma in Proto-RT

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

Asymptotic Behavior of Marginally Trapped Tubes

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

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

13.1 Ion Acoustic Soliton and Shock Wave

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

Formation 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 ) 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