Modeling of striations instabilities and numerical simulations
|
|
- Cameron Nicholson
- 5 years ago
- Views:
Transcription
1 1 Modeling of striations instabilities and numerical simulations C. Besse MIP, Université Paul Sabatier, 118 route de Narbonne, Toulouse cedex, France
2 Joint work with 2 Model derivation, turbulence modeling and numerical simulations F. Deluzet, P. Degond (MIP, Toulouse) J. Claudel, G. Gallice and Ch. Tessieras CEA-CESTA, Le Barp, France Nonlinear instability P. Degond, R. Poncet (MIP, Toulouse) H.-J. Hwang (Duke University)
3 Outline 3 Physical motivations Derivation of a hierarchy of models Striation model in uniform and non uniform magnetic field Numerical simulations Instabilities of the Striation model
4 Physical context and motivations 4 Ionosphere long wave UHF VHF AM radio space wave ground wave sky wave FM radio, TV TV, phones short wave sky wave ground wave Ionosphere : 90 to 1500 km Reflects radio waves : useful for AM radio, long range communications... Many irregularities (aurora, solar eruptions, striations... ) : disturbances of high frequency communications earth-satellites. Focus on striations : long time (many hours) and big scale (100km) irregularities. Striations increase the attenuation and scintillation encoutered by satellite to ground systems
5 Physical context and motivations Northern Hemisphere Evolution of a natural or artificial plasma bubble at altitude around 700 km m k km 15 wind Ionosphere Troposphere Sea level Cosmic Ray Time The plasma bubble is stretched ments Southern Hemisphere along the magnetic field lines. In a plane orthogonal, the plasma bubble is under the influence of the E B instability. The plasma is bent in this plane. Ionospheric Plasma Instabilities 6th MAFPD, Kyoto September 2004
6 Hierarchy of models : euler-maxwell model 6 The modeling is based on the Euler-Maxwell equations with simplifications Only two type of charged particles : O + and electrons Euler equations. No gravity, no energy equations (replace by p e,i = p e,i (n e,i )), no chemical reactions, low electron-ion collisions. Magnetic field B and electric field (dynamo effect) E Maxwell equations. Some parameters are small : m e /m i, u n /c... rescaling Scaling of the Euler-Maxwell equations six dimensionless parameters remain ε α κ τ η β Electron to ion mass ratio squared reciprocal of light speed number of e-n or i-n collisions per rotation period in B-field Mean-time between i-n collisions (dimensionless) Measure of the thermal energy Drift energy relative to magnetic energy measures the strength of the magnetic field perturbation
7 Hierarchy of models : scaled euler-maxwell system 7 t n e + (n e u e ) = 0, τε (Inertia ek ) = η xk p e κ 1 n e (E k + (u e B) k ) ν e n e (u ek u nk ), t n i + (n i u i ) = 0, τ (Inertia ik ) = η xk p i + κ 1 n i (E k + (u i B) k ) ν i (u ik u nk ). α t E B = βj, t B + E = 0, κα β E = ρ, B = 0, ρ = n i n e, κj = n i u i n e u e. Typical values of the parameters for density n i,e = m 3 ε = 10 4, τ = 10 1, η = 10 1, κ = 10 4 α = 10 12, β =
8 Hierarchy of models : scaled euler-maxwell system 8 Euler-Maxwell ε 0, α 0 MHD hierarchy Hall-MHD Dynamo hierarchy κ 0 τ 0 Finite conductivity-mhd Massless Hall-MHD τ 0 β 0 Massless MHD Dynamo β 0 κ 0 Striation MHD hierarchy : valid for large density perturbations Dynamo hierarchy : valid for standard situations
9 MHD hierarchy : massless mhd 9 Plasma is quasi-neutral, ν = ν e + ν i, u := u i = u e, n := n i = n e n t + (nu) = 0, η p(n) = j B νn(u u n ), E + u B = 0, B = βj, B t + E = 0, B = 0. E, u and j can be eliminated from the system gives rise to a system of two equations for n and B.
10 MHD hierarchy : massless mhd 10 Massless MHD : equivalent formulation n t + (nu) = 0, B + (ub Bu) = 0, t with the constraint B = 0 and u given by 1 β (BB) + ( ηp(n) + 1 β ) B 2 2 = νn(u u n ).
11 Dynamo hierarchy : dynamo model 11 Plasma is quasi-neutral n := n i = n e B verifies B = 0 and B = 0 B = B earth B/ t = 0 n t + (nu i) = 0, [ E + u i B = κ ν i (u i u n ) + η ] n p i(n). [ E + u e B = κ ν e (u e u n ) + η ] n p e(n). j = 0, E = 0, κj = n(u i u e ).
12 Multi-Layer Striation model : uniform B field : 12 MHD or Dynamo model always depend on κ or β : are always 3D : limits κ 0 or β 0 Uniform B : B = B x 3 x = (x 1, x 2 ), = ( x1, x2 ), A = (A 1, A 2, A 3 ), A = (A 1, A 2 ). n, u depend on the 3-D coord x n = n(x, t), u = u(x, t) E is orthogonal to B and derives from a potential V : E = (E, 0), E = V V depends on the 2-D coord x V = V (x, t)
13 Multi-Layer Striation model : uniform B field : 13 Suppose p = 0 (for simplicity) n t + (nu) = 0 u = E B ( ) B B + u B 2 n B B ; E = V (x) J = 0, J = 1 B ( σ(x) V + U 2 n B), σ(x) = nν dx 3, U n = nνu n dx 3. σ/ B 2 : field-integrated Pedersen conductivity 3D transport, 2D elliptic
14 Uniform B field : proof 14 Equivalent formulation of the Dynamo model t n + (nu e ) = 0, t n + (nu i ) = 0, u e = M e ( E + κ(ν e u n η p e n )), u i = M i (E + κ(ν i u n η p i n )), E = V, j = 0, κj = n(u i u e ), M e = µ P e,i = µ P e µ H e 0 µ P i µ H i 0 µ H e µ P e 0, M i = µ H i µ P i µ e 0 0 µ i κν e,i (κν e,i ) 2 + B 2, µh e,i = µ e,i = 1 κν e,i. B (κν e,i ) 2 + B 2, M e (resp. M i ) : electron (resp. ion) mobility matrix, We also have j B = n(ν i u i + ν e u e ) nu n (ν i + ν e ).
15 Uniform B field : proof 15 The limit κ 0 : κj = n(u i u e ) u i = u e lim µ P i = lim µ P e = 0, lim µ H i = lim µ H e = 1 κ 0 κ 0 κ 0 κ 0 B, u 1 = x1 V/ B, u 2 = x2 V/ B. third component of the relation leads to kν i u i3 = B x3 V + κν i u n3 x3 V = 0 and u 3 = u n,3. ν = ν i + ν e j B = nν(u u n ) gives the components of j orthogonal to B j = 1 B nν x 1 V + u n2 B 2 x2 V u n1 B
16 Uniform B field : proof 16 j = 0 j 1 x 1 + j 2 x 2 + j 3 x 3 = 0. j 1 and j 2 known a first order differential equation for j 3 integration on the bounded interval [x 3,min, x 3,max ] x 3,min and x 3,max taken inside the neutral atmosphere j 3 vanishes
17 Monolayer striation model in uniform B : 17 Suppose all quantities only depend only on x Suppose u n3 = 0 2D model n t + ( (nu) ) = 0, n V = B 2 ( n u ) n B B 2. with u = E B B 2, E = V u = 0.
18 Non uniform B field : 18 A uniform B field corresponds to cut the earth and make an artificial rectification of a magnetic field tube Extension to a non uniform B differential geometry
19 Non uniform B field : 19 In spherical coordinates, axisymmetric B = (B r (r, ϕ), 0, B ϕ (r, ϕ)) Since B = B = 0, we have (rb ϕ ) r (B r) ϕ = 0 and (r2 sin ϕb r ) r + (r sin ϕb ϕ ) ϕ = 0. β(r, ϕ), γ(r, ϕ), s.t. r β = r sin ϕb ϕ, ϕ β = r 2 sin ϕb r, r γ = B r, ϕ γ = rb ϕ. Transform to (α, β, γ) coordinate system α = θ Orthogonal curvilinear local coordinate system associated to B y z B B ϕ ϕ r B r earth x L earth x r θ y z FIG. 1: The spherical coordinates associated to the earth surrounded by a magnetic field line. (α, β, γ) are known as Euler potentials
20 Non uniform B field : 20 4 β=constant γ=constant 3 2 ts 1 α β γ FIG. 2: Generalized coordinates associated to the magnetic field. FIG. 3: Local coordinates associated to a magnetic field tube.
21 In the (α, β, γ) coordinate system, we have Non uniform B field : 21 ds 2 = 1 B 2 (r2 sin 2 ϕ B 2 dα 2 + dβ2 r 2 sin 2 ϕ + dγ2 ), ( ) 1 f f =, r sin ϕ B f r sin ϕ α β, B f, γ ( ( ) A = B 2 A α + ( ) r sin ϕ A α r sin ϕ B 2 β + ( )) 1 β B γ B A γ ( ( ) fγ r sin ϕ ( )) f β β B γ r sin ϕ B ( A = 1 r sin ϕ B γ (r sin ϕf α) ( )) fγ ( ( α B ) 1 f β ( )) fα B α r sin ϕ B β B,
22 Non uniform B field : 22 γ corresponds to x 3 B-field lines = {(α, β) = Constant } n, u depend on the 3-D coord (α, β, γ) n = n(α, β, γ, t), u = u(α, β, γ, t) E is orthogonal to B and derives from a potential V : E γ = 0, E = V V is constant along magnetic field lines : V = V (α, β, t)
23 Non uniform B field : striation model 23 n t + (nu) = 0, u α = r sin ϕ V β, u β = 1 r sin ϕ B V α, u γ = u nγ, J α α + J β β = 0, J α = A α V α U nβ, J β = A β V β + U nα, A α = U nα = γmax γ min γmax γ min nν dγ r 2 sin 2 ϕ B 4, A β = nνu nα r sin ϕdγ B 2, U nβ = γmax γ min γmax nν r2 sin 2 ϕdγ B 2, γ min nνu nβ dγ r sin ϕ B 3.
24 About the models : 24 Massless MHD : valid for large density perturbations Striation model : Valid for standard situations Striation model : widely used by physicists e.g. Ronchi, Similon, Farley (1989) Zalesak, Ossakov, Chaturvedi (1982) Extension to non uniform B-fields possible
25 Numerical simulations : euler potentials 25 Computations of the Euler potentials Earth dipole field B = µ 0 M (2 cos(ϕ)ˆr + sin(ϕ) ˆϕ), 4π r3 with M = A m 2, µ 0 the vacuum permeability β = µ 0M sin 2 ϕ, γ = µ 0M cos ϕ 4π( r ) 4π r 2 2 ( B 2 µ0 M 4γ 2 = 4π r + β ) 2 r 5.
26 Numerical simulations : scheme 26 Let c = ( α, β ) t and assume u nγ = 0. The multi-layer striation model can be rewritten as : n t + v c n = ns n, c (A c V ) = c J n, with V v = β V α A = A α 0, J n = 0 A β, S n = B 2 β ( 1 U nβ U nα B 2. ) V α,
27 Numerical simulations : scheme 27 n t + v c n = ns n, c (A c V ) = c J n. linear part : a classical convection velocities computed with the electric potential V S n : source term including the magnetic field curvature effects If ν > 0, elliptic equation for the electric potential V with a source term.
28 Numerical simulations : scheme 28 Discretization of the transport equation : Strang splitting method in time and directional splitting method in space On a half time step t n + v c n = 0, On a time step t n = ns n, On a half time step t n + v c n = 0, to ensure the positivity of the density n k+1 = n ( k 1 + t Sn) k, if S k n > 0, n k+1 = n [ k 1 t Sn] k 1, otherwise. linear transport equation : a TVD upwind flux limiter scheme (Ultrabee) Computation of the potential : big linear system to solve
29 Numerical simulations : data 29 computational domain : magnetic field tube immersed in the earth ionosphere the tube is divided into 61 layers central layer located on the equator spreads over 1000 km in the α direction lower and higher altitudes : 700 km and 1900 km each plan is meshed with a Cartesian grid : nodes The neutral wind blows eastward in the α direction with a constant velocity of 45 m s 1.
30 Numerical simulations : data 30 ts Ionosphere upper layer Ionosphere lower layer Magnetic field tube Computational domain FIG. 4: Configuration of a test case. FIG. 5: Mesh of a field tube.
31 Numerical simulations : data 31 eplacements x PSfrag replacements Layer 57 Layer 39 Layer 25 Layer 1 Layers 2, 4, FIG. 6: Slices of a 61 layers discrete field tube. Height (from the earth ground in km) Day plasma density Night plasma density x Plasma density (m 3 ) FIG. 7: Ionospheric ambient plasma.
32 Numerical simulations : data 32 FIG. 8: Initial datum
33 Instabilities of the striation model : 33 Stability analysis of the striation model : phenomenological point of view We consider the monolayer striation model in a uniform B field, B = 1 n t + (nu) = 0 u = V, V = ( x2 V, x1 V ) (n V ) = (nu n ) v = u n = neutral wind (given) after π/2 rotation (n V ) = (nu n ) (nh) = 0 with h = V u n n, x 2 < 0, discontinuous density n(x) = n > n, x 2 > 0, Stable state : V = 0, u n = (0, u n2 ), h = (u n2, 0) x 2 is defined as the direction of the gradient of n.
34 Instabilities of the striation model : 34 Recall that u = V and h = V u n. x 2 u n u n n > x 2 u n u n n > u E u E E u E E u x 1 n < Since u = E B B 2 u E, striations = E B instability. u x 1 n <
35 Instabilities of the striation model : 35
36 Some mathematical results : 36 Local existence in H k, k > 3 Th. Let n 0 H k (R 2 ), n 0 κ > 0. T, a solution (n, u) in C([0, T ], H k ) Lip([0, T ], H k 1 ) Tools : Banach fixed point theorem Energy estimates Linear stability analysis Assumption : u n = (U(x), 0) U of constant sign, smooth, bounded Stationary solution ( n(x), ū = 0) n smooth, bounded, n κ > 0 Instability hypothesis : x 0 such that (U x n) x0 > 0
37 Linearized system Some mathematical results : 37 n = n + η, u = ψ t η + x n y ψ = 0 u = ψ ( n ψ) = U y η Plane wave solutions η(x, y) = η(x) exp(iky) exp(λ k t) ψ(x, y) = ψ(x) exp(iky) exp(λ k t) Growth rate of linear system The growth rate λ k is bounded by a constant Λ. Th. Λ controls the growth rate of solutions of the linearized system : η(t) H s + ψ(t) H s C η(0) H s (1 + t s )e Λt
38 Some mathematical results : 38 Nonlinear instability k 3, ( n(x), 0) steady state. v = (U(x), 0) s.t. x 0, (U x n) x0 > 0 ε 0, δ > 0 small, (n δ, V δ ) solution, T δ > 0 s.t. n δ (0) n H k δ and n δ (T δ ) n L 2,L + V δ (T δ ) L 2,L ε 0 Proof : Follows Grenier, Guo-Strauss, Hwang-Guo,...
Asymptotic-Preserving Schemes for anisotropic diffusion problems related to magnetized plasma simulation
Asymptotic-Preserving Schemes for anisotropic diffusion problems related to magnetized plasma simulation Fabrice Deluzet 1 joint work with C. Besse 2, S. Brull 3, A. Crestetto 1, P. Degond 1, A. Lozinski
More informationQuantum Hydrodynamics models derived from the entropy principle
1 Quantum Hydrodynamics models derived from the entropy principle P. Degond (1), Ch. Ringhofer (2) (1) MIP, CNRS and Université Paul Sabatier, 118 route de Narbonne, 31062 Toulouse cedex, France degond@mip.ups-tlse.fr
More 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 informationPlasma Interactions with Electromagnetic Fields
Plasma Interactions with Electromagnetic Fields Roger H. Varney SRI International June 21, 2015 R. H. Varney (SRI) Plasmas and EM Fields June 21, 2015 1 / 23 1 Introduction 2 Particle Motion in Fields
More informationMODELING THE EARTH S IONOSPHERE: SAMI2 AND SAMI3
MODELING THE EARTH S IONOSPHERE: SAMI2 AND SAMI3 J.D. Huba and G. Joyce Plasma Physics Division Naval Research Laboratory Washington, DC CEDAR Workshop Boulder, CO June, 2010 Icarus Research, Inc. (acknowledge
More informationA Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model
A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program University of Washington IEEE International Conference
More informationOn the occurrence of the equatorial F-region irregularities during solar minimum using radio occultation measurements
On the occurrence of the equatorial F-region irregularities during solar minimum using radio occultation measurements B. A. Carter1, K. Zhang1, R. Norman1, V. V. Kumar2 and S. Kumar3 RMIT University, Australia
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 informationInvestigating COSMIC GPS Radio Occultation Observables as Diagnostics for Ionospheric HF Heating Experiments
Investigating COSMIC GPS Radio Occultation Observables as Diagnostics for Ionospheric HF Heating Experiments ChengYung Huang Institute for Scientific Research, Boston College Chin S. Lin, E. Mishin, and
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationKonvektion und solares Magnetfeld
Vorlesung Physik des Sonnensystems Univ. Göttingen, 2. Juni 2008 Konvektion und solares Magnetfeld Manfred Schüssler Max-Planck Planck-Institut für Sonnensystemforschung Katlenburg-Lindau Convection &
More informationHeliophysics in Atmospheres
Heliophysics in Atmospheres Thermosphere-Ionosphere Response to Geomagnetic Storms Tim Fuller-Rowell NOAA Space Weather Prediction Center and CIRES University of Colorado Atmospheres Gravitationally bound
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 informationA New Equatorial Plasma Bubble Prediction Capability
A New Equatorial Plasma Bubble Prediction Capability Brett A. Carter Institute for Scientific Research, Boston College, USA, http://www.bc.edu/research/isr/, RMIT University, Australia, www.rmit.edu.au/space
More informationGeneral review: - a) Dot Product
General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,
More informationA Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere
A Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere P. Song Center for Atmospheric Research University of Massachusetts Lowell V. M. Vasyliūnas Max-Planck-Institut
More informationKinetic, Fluid & MHD Theories
Lecture 2 Kinetic, Fluid & MHD Theories The Vlasov equations are introduced as a starting point for both kinetic theory and fluid theory in a plasma. The equations of fluid theory are derived by taking
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationOn the sources of day-to-day variability in the occurrence of equatorial plasma bubbles: An analysis using the TIEGCM
On the sources of day-to-day variability in the occurrence of equatorial plasma bubbles: An analysis using the TIEGCM Brett A. Carter, RMIT University, Australia, www.rmit.edu.au/space Institute for Scientific
More informationMHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field
MHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field Y. Matsui, T. Yokoyama, H. Hotta and T. Saito Department of Earth and Planetary Science, University of Tokyo,
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 informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + FV εx, u = 0 is considered in R n. For small ε > 0 it is shown
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 informationSimulating the Ionosphere, one electron at a time.
Simulating the Ionosphere, one electron at a time. Meers Oppenheim CEDAR June 2016 Research supported by NSF, NASA, AFRL, and DOE Grants What? Plasma Physics Particle-in-Cell Simulations Two Examples:
More informationModelling the zonal drift of equatorial plasma irregularities and scintillation. Chaosong Huang Air Force Research Laboratory
Modelling the zonal drift of equatorial plasma irregularities and scintillation Chaosong Huang Air Force Research Laboratory 14 th International Ionospheric Effects Symposium Alexandria, Virginia May 12-14,
More informationChapter 5. Shallow Water Equations. 5.1 Derivation of shallow water equations
Chapter 5 Shallow Water Equations So far we have concentrated on the dynamics of small-scale disturbances in the atmosphere and ocean with relatively simple background flows. In these analyses we have
More informationin Electromagnetics Numerical Method Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD
2141418 Numerical Method in Electromagnetics Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD ISE, Chulalongkorn University, 2 nd /2018 Email: charusluk.v@chula.ac.th Website: Light
More informationSub-Auroral Electric Fields: An Inner Magnetosphere Perspective
Sub-Auroral Electric Fields: An Inner Magnetosphere Perspective Bob Spiro Rice University 2005 GEM/CEDAR Tutorial 1 Introduction/Outline Introduction/Outline Importance of Sub-Auroral E-Fields Early Models
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 informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationInteraction energy between vortices of vector fields on Riemannian surfaces
Interaction energy between vortices of vector fields on Riemannian surfaces Radu Ignat 1 Robert L. Jerrard 2 1 Université Paul Sabatier, Toulouse 2 University of Toronto May 1 2017. Ignat and Jerrard (To(ulouse,ronto)
More informationSolving a magnetic diffusion equation by the Mixed Hybrid Finite Element method.
Solving a magnetic diffusion equation by the Mixed Hybrid Finite Element method. Corinne Aymard, Christophe Fochesato CEA, DAM, DIF, F-91297 Arpajon CEA/DAM/DIF 6 juillet 2010 1 Outline Objectives Model
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 informationCollisions and transport phenomena
Collisions and transport phenomena Collisions in partly and fully ionized plasmas Typical collision parameters Conductivity and transport coefficients Conductivity tensor Formation of the ionosphere and
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 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 informationFluctuation dynamo amplified by intermittent shear bursts
by intermittent Thanks to my collaborators: A. Busse (U. Glasgow), W.-C. Müller (TU Berlin) Dynamics Days Europe 8-12 September 2014 Mini-symposium on Nonlinear Problems in Plasma Astrophysics Introduction
More informationModeling and numerical simulations of swimmers
Plafrim, may 31, 2011 p. 1 Modeling and numerical simulations of swimmers Michel Bergmann INRIA Bordeaux Sud-Ouest, project-team MC2 Institut de Mathématiques Appliquées de Bordeaux 33405 TALENCE cedex,
More informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationNeutral Winds in the Upper Atmosphere. Qian Wu National Center for Atmospheric Research
Neutral Winds in the Upper Atmosphere Qian Wu National Center for Atmospheric Research Outline Overview of the upper atmosphere. Ozone heating. Neutral wind tides (the strongest dynamic feature). Why do
More information1 Energy dissipation in astrophysical plasmas
1 1 Energy dissipation in astrophysical plasmas The following presentation should give a summary of possible mechanisms, that can give rise to temperatures in astrophysical plasmas. It will be classified
More informationAtmospheric delay. X, Y, Z : satellite cartesian coordinates. Z : receiver cartesian coordinates. In the vacuum the signal speed c is constant
Atmospheric delay In the vacuum the signal speed c is constant c τ = ρ = ( X X ) + ( Y Y ) + ( Z Z ) S S S 2 S 2 S 2 X, Y, Z : receiver cartesian coordinates S S S X, Y, Z : satellite cartesian coordinates
More informationSpace Physics. An Introduction to Plasmas and Particles in the Heliosphere and Magnetospheres. May-Britt Kallenrode. Springer
May-Britt Kallenrode Space Physics An Introduction to Plasmas and Particles in the Heliosphere and Magnetospheres With 170 Figures, 9 Tables, Numerous Exercises and Problems Springer Contents 1. Introduction
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationIdeal Magnetohydrodynamics (MHD)
Ideal Magnetohydrodynamics (MHD) Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 1, 2016 These lecture notes are largely based on Lectures in Magnetohydrodynamics
More informationVariations of Ion Drifts in the Ionosphere at Low- and Mid- Latitudes
Variations of Ion Drifts in the Ionosphere at Low- and Mid- Latitudes Edgardo E. Pacheco Jicamarca Radio Observatory Jul, 2014 Outline Motivation Introduction to Ionospheric Electrodynamics Objectives
More informationThe Structure of the Magnetosphere
The Structure of the Magnetosphere The earth s magnetic field would resemble a simple magnetic dipole, much like a big bar magnet, except that the solar wind distorts its shape. As illustrated below, the
More informationElectromagnetic Wave Propagation Lecture 1: Maxwell s equations
Electromagnetic Wave Propagation Lecture 1: Maxwell s equations Daniel Sjöberg Department of Electrical and Information Technology September 2, 2014 Outline 1 Maxwell s equations 2 Vector analysis 3 Boundary
More informationDifferent types of phase transitions for a simple model of alignment of oriented particles
Different types of phase transitions for a simple model of alignment of oriented particles Amic Frouvelle CEREMADE Université Paris Dauphine Joint work with Jian-Guo Liu (Duke University, USA) and Pierre
More informationIdeal MHD Equilibria
CapSel Equil - 01 Ideal MHD Equilibria keppens@rijnh.nl steady state ( t = 0) smoothly varying solutions to MHD equations solutions without discontinuities conservative or non-conservative formulation
More informationElectrostatic Interchange Instabilities of a Rotating, High-Temperature Plasma Confined by a Dipole Magnet: Experiment and Theory
Electrostatic Interchange Instabilities of a Rotating, High-Temperature Plasma Confined by a Dipole Magnet: Experiment and Theory Mike Mauel Columbia University, New York, NY mailto: mauel@columbia.edu
More informationTheoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions
Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates
More informationSome improvements of Xfem for cracked domains
Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,
More informationResearch supported by the NSF Aeronomy & CEDAR Grants
Large-Scale Simulations of Farley- Buneman Turbulence in 2D and 3D and Hybrid Gradient Drift Simulations by Meers Oppenheim, Yakov Dimant, Yann Tambouret Center for Space Physics, Boston University, Boston,
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationGoing with the flow: A study of Lagrangian derivatives
1 Going with the flow: A study of Lagrangian derivatives Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ 12 February
More informationThe Virial Theorem, MHD Equilibria, and Force-Free Fields
The Virial Theorem, MHD Equilibria, and Force-Free Fields Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 10 12, 2014 These lecture notes are largely
More informationAn asymptotic preserving scheme in the drift limit for the Euler-Lorentz system. Stéphane Brull, Pierre Degond, Fabrice Deluzet, Marie-Hélène Vignal
1 An asymptotic preserving scheme in the drift limit for the Euler-Lorentz system. Stéphane Brull, Pierre Degond, Fabrice Deluzet, Marie-Hélène Vignal IMT: Institut de Mathématiques de Toulouse 1. Introduction.
More informationFundamentals of Turbulence
Fundamentals of Turbulence Stanislav Boldyrev (University of Wisconsin - Madison) Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas What is turbulence? No exact definition.
More informationChapter 5 MAGNETIZED PLASMAS. 5.1 Introduction. 5.2 Diamagnetic current
Chapter 5 MAGNETIZED PLASMAS 5.1 Introduction We are now in a position to study the behaviour of plasma in a magnetic field. In the first instance we will re-examine particle diffusion and mobility with
More informationPLASMA ASTROPHYSICS. ElisaBete M. de Gouveia Dal Pino IAG-USP. NOTES: (references therein)
PLASMA ASTROPHYSICS ElisaBete M. de Gouveia Dal Pino IAG-USP NOTES:http://www.astro.iag.usp.br/~dalpino (references therein) ICTP-SAIFR, October 7-18, 2013 Contents What is plasma? Why plasmas in astrophysics?
More informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
More informationJ07M.1 - Ball on a Turntable
Part I - Mechanics J07M.1 - Ball on a Turntable J07M.1 - Ball on a Turntable ẑ Ω A spherically symmetric ball of mass m, moment of inertia I about any axis through its center, and radius a, rolls without
More informationMagnetospheric Currents at Quiet Times
Magnetospheric Currents at Quiet Times Robert L. McPherron Institute of Geophysics and Planetary Physics University of California Los Angeles Los Angeles, CA 90095-1567 e-mail: rmcpherron@igpp.ucla.edu
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 informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationTowards general-relativistic pulsar magnetospheres
Towards general-relativistic pulsar magnetospheres Jérôme Pétri Observatoire Astronomique de Strasbourg, Université de Strasbourg, France. Physics of Neutron Stars, Saint-Petersbourg, 29/7/2014 Jérôme
More informationcos 6 λ m sin 2 λ m Mirror Point latitude Equatorial Pitch Angle Figure 5.1: Mirror point latitude as function of equatorial pitch angle.
Chapter 5 The Inner Magnetosphere 5.1 Trapped Particles The motion of trapped particles in the inner magnetosphere is a combination of gyro motion, bounce motion, and gradient and curvature drifts. In
More informationAsymptotic stability for solitons of the Gross-Pitaevskii and Landau-Lifshitz equations
Asymptotic stability for solitons of the Gross-Pitaevskii and Landau-Lifshitz equations Philippe Gravejat Cergy-Pontoise University Joint work with F. Béthuel (Paris 6), A. de Laire (Lille) and D. Smets
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n
More informationEquatorial and Low Latitude Scintillation Initiated From Low Altitude Forcing via Hurricanes/Typhoons
Equatorial and Low Latitude Scintillation Initiated From Low Altitude Forcing via Hurricanes/Typhoons Dr. Rebecca Bishop The Aerospace Corporation PSL/SSAL 30 March 2012 The Aerospace Corporation 2012
More informationLecture 14: Solar Cycle. Observations of the Solar Cycle. Babcock-Leighton Model. Outline
Lecture 14: Solar Cycle Outline 1 Observations of the Solar Cycle 2 Babcock-Leighton Model Observations of the Solar Cycle Sunspot Number 11-year (average) cycle period as short as 8 years as long as 15
More informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationAncient solutions to geometric flows
Columbia University Banff May 2014 Outline We will discuss ancient or eternal solutions to geometric flows, that is solutions that exist for all time < t < T where T (, + ]. Such solutions appear as blow
More information2. Outline of the MRI-EPS
2. Outline of the MRI-EPS The MRI-EPS includes BGM cycle system running on the MRI supercomputer system, which is developed by using the operational one-month forecasting system by the Climate Prediction
More informationIntroduction to the School
Lucio Crivellari Instituto de Astrofísica de Canarias D.pto de Astrofísica, Universidad de La Laguna & INAF Osservatorio Astronomico di Trieste (Italy) Introduction to the School 10/11/17 1 Setting the
More informationThermospheric Winds. Astrid Maute. High Altitude Observatory (HAO) National Center for Atmospheric Science (NCAR) Boulder CO, USA
Thermospheric Winds Astrid Maute High Altitude Observatory (HAO) National Center for Atmospheric Science (NCAR) Boulder CO, USA High Altitude Observatory (HAO) National Center for Atmospheric Research
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
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 informationE.M.WAVES 1. Taller the antenna longer is the coverage of television broadcast. Justify this statement with the help of a figure. 2.If v g, v x v m represents the speed of gamma rays, X-rays microwaves
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 informationConvergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions
Convergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions Marie Hélène Vignal UMPA, E.N.S. Lyon 46 Allée d Italie 69364 Lyon, Cedex 07, France abstract. We are interested
More informationParity of solar global magnetic field determined by turbulent diffusivity
First Asia-Pacific Solar Physics Meeting ASI Conference Series, 2011, Vol. 1, pp 117 122 Edited by Arnab Rai Choudhuri & Dipankar Banerjee Parity of solar global magnetic field determined by turbulent
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationPhysical Processes in Astrophysics
Physical Processes in Astrophysics Huirong Yan Uni Potsdam & Desy Email: hyan@mail.desy.de 1 Reference Books: Plasma Physics for Astrophysics, Russell M. Kulsrud (2005) The Physics of Astrophysics, Frank
More informationJournal of Computational and Applied Mathematics. Determination of a material constant in the impedance boundary condition for electromagnetic fields
Journal of Computational and Applied Mathematics 234 (2010) 2062 2068 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
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 informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More 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 informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationFinite-time singularity formation for Euler vortex sheet
Finite-time singularity formation for Euler vortex sheet Daniel Coutand Maxwell Institute Heriot-Watt University Oxbridge PDE conference, 20-21 March 2017 Free surface Euler equations Picture n N x Ω Γ=
More informationDifferent types of phase transitions for a simple model of alignment of oriented particles
Different types of phase transitions for a simple model of alignment of oriented particles Amic Frouvelle Université Paris Dauphine Joint work with Jian-Guo Liu (Duke University, USA) and Pierre Degond
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More information