Development of a stable coupling of the Yee scheme with linear current
|
|
- Cora Sullivan
- 5 years ago
- Views:
Transcription
1 Development of a stable coupling of the Yee scheme with linear current Martin Campos Pinto (LJLL), Bruno Després (LJLL) Stéphane Heuraux (IJL), Filipe Da Silva (IPFN+IST) 15 octobre 2013 Munich: Pereverzev legacy 13/10/2013 p. 1 / 30
2 Outline Munich: Pereverzev legacy 13/10/2013 p. 2 / 30
3 The equations No magnetisation (µ = µ 0 ). Maxwell equations with a linear current derive from the linearization µ 0 H B 0 of the Vlasov-Maxwell system (for electrons) around a strong magnetic field B 0 : ε 0 t E + H = q e N e (x)u e, µ 0 t H + E = 0, m e t u e = q e (E + B 0 (x) u e ) νm e u e. Or, writing J = q e N e (x)u e, ε 0 t E = H J, µ 0 t H = E, t J = ε 0 ωpe 2 + ω c b J with ω p (x) = q 2 e N e(x) mε 0, ω c (x) = qe B 0(x) m e and b(x) = B 0(x) B 0 (x). Munich: Pereverzev legacy 13/10/2013 p. 3 / 30
4 Direct simulation of reflectometry configuration The domain is a parallepiped ( 1500 cells in x direction) with an antenna on the side : pulsation ω Cut-off : in O mode (TM), waves propagate if ω ω p (x). Cyclotron resonance : ω = ω c Hybrid resonance : ω 2 = ω p (x) 2 + ω 2 c Munich: Pereverzev legacy 13/10/2013 p. 4 / 30
5 Standard scheme of Xu-Yuan (2006) Based on the Yee scheme for the (E, H) field : general form is ε 0 t (E n+1 E n ) = RH n+ 1 2 J n+ 1 2 µ 0 t (Hn+ 3 2 H n+ 1 2 ) = R t E n+1 1 t (Jn+ 3 2 J n+ 1 2 ) = ε 0 ωpe 2 n+1 + ω c b 1 2 (Jn J n+ 1 2 ). Need to specify the operator h on the Yee grid Munich: Pereverzev legacy 13/10/2013 p. 5 / 30
6 X-mode equations X-mode=Transverse electric (O-mode not discussed in this talk). ε 0 t E x + y H z = J x, J x = en e u x, ε 0 t E y x H z = J y, J y = en e u y, µ 0 t H z + x E y y E x = 0, m e t u x = ee x + eu y Bz 0, m e t u y = ee y eu x Bz 0. Call VLC external Munich: Pereverzev legacy 13/10/2013 p. 6 / 30
7 Main difficulty In fusion plasmas, N e (x) has huge fluctuations along the main axis n e [ m -3 ] n e w/o frm n e w/ frm Negative density (num. or measurement artifact) induces automatically an instability, as well as strong spatial gradient at the plasma edge (phys.) or inside the plasma (phys.). Munich: Pereverzev legacy 13/10/2013 p. 7 / 30
8 Example of unstability (for large times) Magnetic field 20 log 10 H z (where H z = H z L ) vs. time step and level of fluctuations Expertise from F. Da Silva and S. Heuraux. Munich: Pereverzev legacy 13/10/2013 p. 8 / 30
9 Some references Xu-Yuan, FDTD Formulations for Scattering From 3-D Anisotropic Magnetized Plasma Objects, IEEE-2006 Bohner, Simulation of Microwave Propagation in a Fusion Plasma, Diploma Thesis, MaxPlanck, 2011 An unconditionally stable (?) time-domain discretization on cartesian meshes for the simulation of nonuniform magnetized cold plasma, JCP-2012, Tierens-Zutter Yu-Simpson, An E-J Collocated 3-D FDTD Model of Electromagnetic Wave Propagation in Magnetized Cold Plasma, Ieee Transations on antennas and propagation, 58-2, 2010 Smithe, Finite-difference time-domain simulation of fusion plasmas at radiofrequency time scales, Physics of plasmas, 2007 da Silva-Heuraux-Ribeiro-Scott, (2013). Development of a 2D full-wave JE-FDTD Maxwell X-mode code for reflectometry simulation (pp. 16). Presented at the 9th International Reflectometry Workshop. Munich: Pereverzev legacy 13/10/2013 p. 9 / 30
10 Outline Munich: Pereverzev legacy 13/10/2013 p. 10 / 30
11 Energy conservation at the continuous level For simplicity : constant density profile N e (x, t) = N e (x). Inside the computational domain (no boundaries assumed), the total energy is conserved in time, ( d ε0 E 2 + H 2 + m en e (x) u e 2 ) dv = 0. dt Ω 2 2µ 0 2 Using normalized variables Ê := 1 c E, Ĥ := µ 0H and Ĵ := 1 ω pcε 0 J, we have ( ) d Ê 2 dt 2 + Ĥ Ĵ 2 dv = 0. 2 Munich: Pereverzev legacy 13/10/2013 p. 11 / 30
12 Classical stability analysis for the Yee scheme With normalized variables, the Yee scheme (J = 0) reads 1 t (Ê n+1 Ê n ) = crĥ n t (Ĥn+ 1 2 Ĥ n 1 2 ) = cr t Ê n where {Ê := 1 c E Ĥ := µ 0 H. In particular, the energy Ê n := Ê n 2 h + Ĥn h satisfies Ê n+1 Ê n = c t ( RĤ n+ 1 2, Ê n+1 + Ê n R t Ê n, Ĥ n Ĥ n 1 2 ) hence E n := Ê n c t Ê n, RĤn 1 2 is constant. Moreover, Ê n, RĤn R Ên = Ê n (1 c t 2 R ) En = Stability in the energy norm : for c t < 2/ R = h/ 3 Munich: Pereverzev legacy 13/10/2013 p. 12 / 30
13 Stability analysis for an abstract Yee+J scheme With Ê, Ĥ and Ĵ := 1 ω pcε 0 J, the abstract Yee+J scheme is 1 t (Ê n+1 Ê n ) = crĥ n+ 1 2 ω p Ĵ n t (Ĥn+ 1 2 Ĥn 1 2 ) = cr t Ê n 1 t (Ĵn+ 1 2 Ĵ n 1 2 ) = ω p Ê n Ĵ + ω c b n Ĵ n 1 2 h 2 Here the energy Ê n := Ê n 2 + Ĥ n Ĵ n satisfies t ( ω p Ĵ n+ 1 2, Ê n+1 + Ê n ω p Ê n, Ĵ n Ĵ n 1 2 ) provided V, b h V = 0 for all V. Stability in the energy norm : for t 2 (12c2 h 2 + ω p 2 L ) 1 2 < 1. Munich: Pereverzev legacy 13/10/2013 p. 13 / 30
14 Remark on average cross products One can use local averages to define a 2nd order cross product, (b h V ) x := b y {V z } b z {V y } (b h V ) y := (b h V ) z := Then if b(x) = B 0(x) B 0 is uniform, V, b h V = 0 holds for all V previous analysis applies. If b(x) is not uniform this is not so clear... Munich: Pereverzev legacy 13/10/2013 p. 14 / 30
15 Improved stability for a new Yee+J scheme Discretizing the current on t n, t n+1,... yields a new scheme 1 t (Ê n+1 Ê n ) = crĥ n+ 1 Ĵ n ω Ĵn p 2 1 t (Ĥn+ 1 2 Ĥn 1 2 ) = cr t Ê n 1 t (Ĵn+1 Ĵn 1 ) = ω p {Ê} n+ 1 Ĵ n ω c b Ĵn h The energy Ê n satisfies Ê n+1 Ê n =c t ( RĤ n+ 1 2, Ê n+1 + Ê n R t Ê n, Ĥ n Ĥ n ) 1 2 t ( ω p {Ĵ} n+ 1 2, 2{Ê} n+ 1 2 ω p {Ê} n+ 1 2, 2{Ĵ} n+ ) 1 2 once again provided V, b h V = 0 for all V. Stability in the energy norm : for c t < 2/ R = h/ Munich: Pereverzev legacy 13/10/2013 p. 15 / 30
16 Outline Munich: Pereverzev legacy 13/10/2013 p. 16 / 30
17 Example of the Xu-Yuan scheme Based on the Yee scheme for the (E, H) field : general form is ε 0 t (E n+1 E n ) = RH n+ 1 2 J n+ 1 2 µ 0 t (Hn+ 3 2 H n+ 1 2 ) = R t E n+1 1 t (Jn+ 3 2 J n+ 1 2 ) = ε 0 ωpe 2 n+1 + ω c b 1 2 (Jn J n+ 1 2 ). Munich: Pereverzev legacy 13/10/2013 p. 17 / 30
18 Example of the Xu-Yuan scheme Based on the Yee scheme for the (E, H) field : general form is ε 0 t (E n+1 E n ) = RH n+ 1 2 J n+ 1 2 µ 0 t (Hn+ 3 2 H n+ 1 2 ) = R t E n+1 1 t (Jn+ 3 2 J n+ 1 2 ) = ε 0 ωpe 2 n ω c b h 2 (Jn J n+ 1 2 ). Need an explicit solver with the b h operator. Munich: Pereverzev legacy 13/10/2013 p. 17 / 30
19 Problem with the X-Y approach Consider once again the cross product by local averages (b h V ) x := b y {V z } b z {V y } (b h V ) y := (b h V ) z := The result ε 0 t (E n+1 E n ) = RH n+ 1 2 J n+ 1 2 µ 0 t (Hn+ 3 2 H n+ 1 2 ) = R t E n+1 1 t (Jn+ 3 2 J n+ 1 2 ) = ε 0 ωpe 2 n ω c b h 2 (Jn J n+ 1 2 ). is a global scheme which needs a linear solver to invert the matrice. Munich: Pereverzev legacy 13/10/2013 p. 18 / 30
20 Solution : use clustered cross-products Instead, choose a pattern (α, β, γ) { 1, +1} 3 and define the first order cross product with local clusters : (b h V ) x i+ α 2,j,k := b y V z i,j,k+ γ b z{v y } 2 i,j+ β 2,k (b h V ) y i,j+ β 2,k := b zv x i+ α 2,j,k b x {V z } i,j,k+ γ 2 (b h V ) z i,j,k+ γ := b xv y 2 i,j+ β 2,k b y {V x } i+ α 2,j,k The resulting scheme ε 0 t (E n+1 E n ) = RH n+ 1 2 J n+ 1 2 µ 0 t (Hn+ 3 2 H n+ 1 2 ) = R t E n+1 1 t (Jn+ 3 2 J n+ 1 2 ) = ε 0 ωpe 2 n ω c b h 2 (Jn J n+ 1 2 ). can be solved with a local procedure (i.e. solution is explicit and local). Munich: Pereverzev legacy 13/10/2013 p. 19 / 30
21 Abstract criterion The criterion for explicit scheme writes : (b h ) 4 = (b h ) 2. Indeed one has the implications J αb h J = Z, J α 2 (b h ) 2 J = (I + αb h )Z, (1 + α 2 )(b h ) 2 J = (b h ) 2 (I + αb h )Z, J = (I + αb h )Z + α2 1 + α 2 (b h) 2 (I + αb h )Z. This algebra is enough to compute the solution by means of explicit and local formulas (for MXYK and new Kernel). Munich: Pereverzev legacy 13/10/2013 p. 20 / 30
22 Short summary The coupling of the Yee scheme and a linear current is for : (V, b h V ) = 0 for : (b h ) 4 = (b h ) 2 Solution (so far) is clustered first order product Additional and natural condition is that ω p and ω c are the same within a cluster. Munich: Pereverzev legacy 13/10/2013 p. 21 / 30
23 Outline Munich: Pereverzev legacy 13/10/2013 p. 22 / 30
24 Data N e Zoom around the foot of the ramp Cut of the electronic density in the horizontal direction. An additional kink (in red) is sometimes added at x = 500 to evaluate the effect of an extremely strong gradient. Munich: Pereverzev legacy 13/10/2013 p. 23 / 30
25 H z plot With the kick and 30% noise Without the kick but 40% noise An instability shows up near x = 500 cells on the left, near x = 1000 Munich: Pereverzev legacy 13/10/2013 p. 24 / 30
26 20 log 10 H z L With respect to the time and to the level of noise. With the kick on the left, without the kick on the right. Munich: Pereverzev legacy 13/10/2013 p. 25 / 30
27 With the first order vectorial product 20 log 10 H z L, with respect to the time and to the level of noise. The computation is done We observe unconditional stability, with however more amplitude for a higher level of noise. The number of time steps is much greater than in previous figure to illustrate the long time stability of the method. Munich: Pereverzev legacy 13/10/2013 p. 26 / 30
28 Energy dissipation Initial data is a Dirac mass, at the exact foot of the electronic density ramp. The external magnetic field used in this set of runs was B 0 = 0.95T. The plasma density N e (x) is linear, with its edge at x = 500 grid point. The number of iterations considered is N = 700 (far from PML layer). (slide courtesy of F. Da Silva) Munich: Pereverzev legacy 13/10/2013 p. 27 / 30
29 Additional remarks and Need to use clustered multiplications by scalar fields, consistent with clustered cross products. Counter-intuitive : the stable and explicit scheme is globally first order (and not second order like the standard Yee scheme). Possibility to average in time by alternating the cluster patterns (α, β, γ) in { 1, +1} 3 Work in progress for direct simulation of time-dependent densities N e = N e (x, t) (Doppler reflectometry) A paper is being written Munich: Pereverzev legacy 13/10/2013 p. 28 / 30
30 An open question (is it really?) Look at ε 0 t E + H = q e N e (x)u e, µ 0 t H + E = 0, m e t u e = q e (E + B 0 (x) u e ) νm e u e plus harmonic forcing on the boundary, plus initial condition, plus friction ν > 0. Assume resonance configuration (cyclotron, hybrid,...) : do we have lim lim = lim lim T T ν 0 + ν 0 +? In other words, do we have Limit absorption=limit amplitude? If not, which one is the correct physical solution? Munich: Pereverzev legacy 13/10/2013 p. 29 / 30
31 lim ν 0 + lim T (L.M. Imbert-Grard) Munich: Pereverzev legacy 13/10/2013 p. 30 / 30
Full-Wave Maxwell Simulations for ECRH
Full-Wave Maxwell Simulations for ECRH H. Hojo Plasma Research Center, University of Tsukuba in collaboration with A. Fukuchi, N. Uchida, A. Shimamura, T. Saito and Y. Tatematsu JIFT Workshop in Kyoto,
More informationDevelopment of a 2D full-wave JE-FDTD Maxwell X-mode code for reflectometry simulation
Development of a 2D full-wave JE-FDTD Maxwell X-mode code for reflectometry simulation F. da Silva, S. Heuraux, T. Ribeiro and B. Scott Associação EURATOM/IST Instituto de Plasmas e Fusão Instituto Superior
More informationplasmas Lise-Marie Imbert-Gérard, Bruno Després. July 27th 2011 Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris.
Lise-Marie Imbert-Gérard, Bruno Després. Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris. July 27th 2011 1 Physical and mathematical motivations 2 approximation of the solutions 3 4 Plan
More informationFinite Difference Solution of Maxwell s Equations
Chapter 1 Finite Difference Solution of Maxwell s Equations 1.1 Maxwell s Equations The principles of electromagnetism have been deduced from experimental observations. These principles are Faraday s law,
More informationElectromagnetic waves in magnetized plasma The dispersion relation
Electromagnetic waves in magnetized plasma The dispersion relation Bruno Després (LJLL-UPMC) Electromagnetic waves in magnetized plasma The dispersion relation p. 1 / 32 Vlasov-Maxwell Vectors are in bold.
More informationCONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING FDTD METHODS FOR THE 3D MAXWELL S EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5 Number 6 Pages 747 784 c 08 Institute for Scientific Computing and Information CONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING
More informationFDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
FDTD for 1D wave equation Equation: 2 H = t 2 c2 2 H x 2 Notations: o t = nδδ, x = iδx o n H nδδ, iδx = H i o n E nδδ, iδx = E i discretization H 2H + H H 2H + H n+ 1 n n 1 n n n i i i 2 i+ 1 i i 1 = c
More informationChapter 6: Vector Analysis
Chapter 6: Vector Analysis We use derivatives and various products of vectors in all areas of physics. For example, Newton s 2nd law is F = m d2 r. In electricity dt 2 and magnetism, we need surface and
More informationWaves in plasmas. S.M.Lea
Waves in plasmas S.M.Lea 17 1 Plasma as an example of a dispersive medium We shall now discuss the propagation of electromagnetic waves through a hydrogen plasm an electrically neutral fluid of protons
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 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 informationEnergy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media
Energy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media Yingda Cheng Michigan State University Computational Aspects of Time Dependent Electromagnetic Wave Problems
More informationRadio Propagation Channels Exercise 2 with solutions. Polarization / Wave Vector
/8 Polarization / Wave Vector Assume the following three magnetic fields of homogeneous, plane waves H (t) H A cos (ωt kz) e x H A sin (ωt kz) e y () H 2 (t) H A cos (ωt kz) e x + H A sin (ωt kz) e y (2)
More informationEfficient domain decomposition methods for the time-harmonic Maxwell equations
Efficient domain decomposition methods for the time-harmonic Maxwell equations Marcella Bonazzoli 1, Victorita Dolean 2, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier 4 1 Inria Saclay (Defi
More informationTu Efficient 3D MT Inversion Using Finite-difference Time-domain Modelling
Tu 11 16 Efficient 3D MT Inversion Using Finite-difference Time-domain Modelling S. de la Kethulle de Ryhove* (Electromagnetic Geoservices ASA), J.P. Morten (Electromagnetic Geoservices ASA) & K. Kumar
More informationLecture 4: The particle equations (1)
Lecture 4: The particle equations (1) Presenter: Mark Eric Dieckmann Department of Science and Technology (ITN), Linköping University, Sweden July 17, 2014 Overview We have previously discussed the leapfrog
More informationPHYSICS OF HOT DENSE PLASMAS
Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced
More 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 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 informationComputational Methods in Plasma Physics
Computational Methods in Plasma Physics Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Purpose of Talk Describe use of numerical methods to solve simple problem in plasma
More informationFDFD. The Finite-Difference Frequency-Domain Method. Hans-Dieter Lang
FDFD The Finite-Difference Frequency-Domain Method Hans-Dieter Lang Friday, December 14, 212 ECE 1252 Computational Electrodynamics Course Project Presentation University of Toronto H.-D. Lang FDFD 1/18
More informationExplanation of prompt growth of ECE signal in tokamak runaway electron experiments
Chang Liu et al. 2nd IAEA TM on the Fusion Data Processing, Validation and Analysis 1 Explanation of prompt growth of ECE signal in tokamak runaway electron experiments Chang Liu 1, Lei Shi 2, Eero Hirvijoki
More informationOverview of FRC-related modeling (July 2014-present)
Overview of FRC-related modeling (July 2014-present) Artan Qerushi AFRL-UCLA Basic Research Collaboration Workshop January 20th, 2015 AFTC PA Release# 15009, 16 Jan 2015 Artan Qerushi (AFRL) FRC modeling
More informationAn iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems
An iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems K. Hertel 1,2, S. Yan 1,2, C. Pflaum 1,2, R. Mittra 3 kai.hertel@fau.de 1 Lehrstuhl für Systemsimulation
More informationPlasma turbulence measured by fast sweep reflectometry on TORE SUPRA
Plasma turbulence measured by fast sweep reflectometry on TORE SUPRA F. Clairet &, L. Vermare &, S. Heuraux, G. Leclert # & Association Euratom-CEA sur la fusion, DSM/DRFC/SCCP C.E. Cadarache, 8 Saint-Paul-lès-Durance,
More informationPHY481: Electromagnetism
PHY481: Electromagnetism Vector tools Sorry, no office hours today I ve got to catch a plane for a meeting in Italy Lecture 3 Carl Bromberg - Prof. of Physics Cartesian coordinates Definitions Vector x
More informationA Survey of Computational High Frequency Wave Propagation I
A Survey of Computational High Frequency Wave Propagation I Bjorn Engquist University of Texas at Austin CSCAMM Workshop on High Frequency Wave Propagation, University of Maryland, September 19-22, 2005
More informationMEMORANDUM-4. n id (n 2 + n2 S) 0 n n 0. Det[ɛ(ω, k)]=0 gives the Dispersion relation for waves in a cold magnetized plasma: ω 2 pα ω 2 cα ω2, ω 2
Fundamental dispersion relation MEMORANDUM-4 n S) id n n id n + n S) 0 } n n 0 {{ n P } ɛω,k) E x E y E z = 0 Det[ɛω, k)]=0 gives the Dispersion relation for waves in a cold magnetized plasma: n P ) [
More informationAnisotropic electron distribution functions and the transition between the Weibel and the whistler instabilities
Anisotropic electron distribution functions and the transition between the Weibel and the whistler instabilities F. Pegoraro, L. Palodhi, F. Califano 5 th INTERNATIONAL CONFERENCE ON THE FRONTIERS OF PLASMA
More informationECE 497 JS Lecture - 13 Projects
ECE 497 JS Lecture - 13 Projects Spring 2004 Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jose@emlab.uiuc.edu 1 ECE 497 JS - Projects All projects should be accompanied
More informationHigh-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials
Progress In Electromagnetics Research Letters, Vol. 77, 103 107, 01 High-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials Wei-Jun Chen
More informationOutline. 1 Why PIC Simulations Make Sense. 2 The Algorithm. 3 Examples Surface High Harmonics Generation. 4 Extensions Of The PIC Algorithm
PIC Simulations an Introduction GRK 1203 Meeting February 12-15, 2008, Oelde Outline 1 Simulations Make Sense 2 3 Surface High Harmonics Generation 4 Of PIC Plasma Physics Is Complex Experiment real thing
More informationPHY481: Electromagnetism
PHY481: Electromagnetism Vector tools Lecture 4 Carl Bromberg - Prof. of Physics Cartesian coordinates Definitions Vector x is defined relative to the origin of 1 coordinate system (x,y,z) In Cartsian
More informationarxiv:physics/ v1 19 Feb 2004
Finite Difference Time Domain (FDTD) Simulations of Electromagnetic Wave Propagation Using a Spreadsheet David W. Ward and Keith A. Nelson Department of Chemistry Massachusetts Institute of Technology,
More informationGuangye Chen, Luis Chacón,
JIFT workshop! Oct. 31, 2014 New Orleans, LA.! Guangye Chen, Luis Chacón, CoCoMANs team Los Alamos National Laboratory, Los Alamos, NM 87545, USA gchen@lanl.gov 1 Los Alamos National Laboratory Motivation
More informationLecture10: Plasma Physics 1. APPH E6101x Columbia University
Lecture10: Plasma Physics 1 APPH E6101x Columbia University Last Lecture - Conservation principles in magnetized plasma frozen-in and conservation of particles/flux tubes) - Alfvén waves without plasma
More informationEdge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas
Edge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas by T.M. Wilks 1 with A. Garofalo 2, K.H. Burrell 2, Xi. Chen 2, P.H. Diamond 3, Z.B. Guo 3, X. Xu
More informationAdaptive semi-lagrangian schemes for transport
for transport (how to predict accurate grids?) Martin Campos Pinto CNRS & University of Strasbourg, France joint work Albert Cohen (Paris 6), Michel Mehrenberger and Eric Sonnendrücker (Strasbourg) MAMCDP
More informationFundamentals of wave kinetic theory
Fundamentals of wave kinetic theory Introduction to the subject Perturbation theory of electrostatic fluctuations Landau damping - mathematics Physics of Landau damping Unmagnetized plasma waves The plasma
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 informationANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA
Progress In Electromagnetics Research Letters, Vol. 22, 83 93, 2011 ANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA X. Ai Science and Technology
More informationCOLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE
Progress In Electromagnetics Research M, Vol. 3, 239 252, 213 COLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE Lijuan Shi 1, 3, Lixia Yang 2, *, Hui Ma 2, and Jianning Ding 3 1 School
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationJasmin Smajic1, Christian Hafner2, Jürg Leuthold2, March 23, 2015
Jasmin Smajic, Christian Hafner 2, Jürg Leuthold 2, March 23, 205 Time Domain Finite Element Method (TD FEM): Continuous and Discontinuous Galerkin (DG-FEM) HSR - University of Applied Sciences of Eastern
More informationChapter 9 WAVES IN COLD MAGNETIZED PLASMA. 9.1 Introduction. 9.2 The Wave Equation
Chapter 9 WAVES IN COLD MAGNETIZED PLASMA 9.1 Introduction For this treatment, we will regard the plasma as a cold magnetofluid with an associated dielectric constant. We then derive a wave equation using
More informationA new method for sensitivity analysis of photonic crystal devices
A new method for sensitivity analysis of photonic crystal devices Georgios Veronis, Robert W. Dutton, and Shanhui Fan Department of Electrical Engineering, Stanford University, Stanford, California 9435
More informationFull-wave Electromagnetic Field Simulations in the Lower Hybrid Range of Frequencies
Full-wave Electromagnetic Field Simulations in the Lower Hybrid Range of Frequencies P.T. Bonoli, J.C. Wright, M. Porkolab, PSFC, MIT M. Brambilla, IPP, Garching, Germany E. D Azevedo, ORNL, Oak Ridge,
More informationMODELING THE ELECTROMAGNETIC FIELD
MODELING THE ELECTROMAGNETIC FIELD IN ANISOTROPIC INHOMOGENEOUS MAGNETIZED PLASMA OF ECR ION SOURCES G. Torrisi (1), D. Mascali (1), A. Galatà (2), G. Castro (1), L. Celona (1), L. Neri (1), G. Sorbello
More informationAdvanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III
Advanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III 1 11. Yee s discrete algorithm Maxwell s equations are discretized using central FDs. We set the magnetic loss equal to
More informationFast ion physics in the C-2U advanced, beam-driven FRC
Fast ion physics in the C-2U advanced, beam-driven FRC Richard Magee for the TAE Team 216 US-Japan Workshop on the Compact Torus August 23, 216! High β FRC embedded in magnetic mirror is a unique fast
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationAntennas and Propagation. Chapter 2: Basic Electromagnetic Analysis
Antennas and Propagation : Basic Electromagnetic Analysis Outline Vector Potentials, Wave Equation Far-field Radiation Duality/Reciprocity Transmission Lines Antennas and Propagation Slide 2 Antenna Theory
More informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
More informationComparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses
Journal of the Korean Physical Society, Vol. 58, No. 4, April 2011, pp. 1014 1020 Comparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses Jian-Bo Chen, Yan Shen,
More informationPublication II Wiley Periodicals. Reprinted by permission of John Wiley & Sons.
Publication II Ilkka Laakso and Tero Uusitupa. 2008. Alternative approach for modeling material interfaces in FDTD. Microwave and Optical Technology Letters, volume 50, number 5, pages 1211-1214. 2008
More informationThe importance of including XMHD physics in HED codes
The importance of including XMHD physics in HED codes Charles E. Seyler, Laboratory of Plasma Studies, School of Electrical and Computer Engineering, Cornell University Collaborators: Nat Hamlin (Cornell)
More informationDoppler Reflectometry Simulations for ASDEX Upgrade
Doppler Reflectometry Simulations for ASDEX Upgrade C. Lechte IGVP University of Stuttgart Pfaffenwaldring 31, 70569 Stuttgart Germany Phone +49 711 685 62306 Fax +49 711 685 63102 G. D. Conway, T. Görler,
More informationSIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(0,n) MODES*
25th IEEE International Conference on Plasma Science Raleigh, North Carolina June 1-4, 1998 SIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(,n) MODES* Ron L. Kinder and Mark J.
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationHybrid Simulation Method ISSS-10 Banff 2011
Hybrid Simulation Method ISSS-10 Banff 2011 David Burgess Astronomy Unit Queen Mary University of London With thanks to Dietmar Krauss-Varban Space Plasmas: From Sun to Earth Space Plasma Plasma is (mostly)
More informationEMAFF: MAgnetic Equations with FreeFem++
EMAFF: MAgnetic Equations with FreeFem++ The Grad-Shafranov equation & The Current Hole Erwan DERIAZ Bruno DESPRÉS Gloria FACCANONI š Kirill Pichon GOSTAF Lise-Marie IMBERT-GÉRARD Georges SADAKA Remy
More informationParallel Heating Associated with Interaction of Forward and Backward Electromagnetic Cyclotron Waves
J. Geomag. Geoelectr., 40, 949-961, 1988 Parallel Heating Associated with Interaction of Forward and Backward Electromagnetic Cyclotron Waves Yoshiharu OMURA1, Hideyuki USUI2, and Hiroshi MATSUMOTO1 2Department
More informationGravitational-like interactions in a cloud of cold atoms?
Gravitational-like interactions in a cloud of cold atoms? J. Barré (1), B. Marcos (1), A. Olivetti (1), D. Wilkowski (2), M. Chalony (2) (+ R. Kaiser (2) ) 1 Laboratoire JA Dieudonné, U. de Nice-Sophia
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationNumerical Analysis of Electromagnetic Fields in Multiscale Model
Commun. Theor. Phys. 63 (205) 505 509 Vol. 63, No. 4, April, 205 Numerical Analysis of Electromagnetic Fields in Multiscale Model MA Ji ( ), FANG Guang-You (ྠ), and JI Yi-Cai (Π) Key Laboratory of Electromagnetic
More informationA Brief Revision of Vector Calculus and Maxwell s Equations
A Brief Revision of Vector Calculus and Maxwell s Equations Debapratim Ghosh Electronic Systems Group Department of Electrical Engineering Indian Institute of Technology Bombay e-mail: dghosh@ee.iitb.ac.in
More informationLinear and non-linear evolution of the gyroresonance instability in Cosmic Rays
Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays DESY Summer Student Programme, 2016 Olga Lebiga Taras Shevchenko National University of Kyiv, Ukraine Supervisors Reinaldo
More informationAn Upwind Leapfrog Scheme for Computational Electromagnetics: CL-FDTD
156 Progress In Electromagnetics Research Symposium 2005, Hangzhou, China, August 22-26 An Upwind Leapfrog Scheme for Computational Electromagnetics: CL-FDTD Yong Zhang, K. R. Shao, and X. W. Hu Huazhong
More informationB. H. Jung Department of Information and Communication Engineering Hoseo University Asan, Chungnam , Korea
Progress In Electromagnetics Research, PIER 77, 111 120, 2007 ANALYSIS OF TRANSIENT ELECTROMAGNETIC SCATTERING WITH PLANE WAVE INCIDENCE USING MOD-FDM B. H. Jung Department of Information and Communication
More informationMulti-Domain Approaches for the Solution of High-Frequency Time-Harmonic Propagation Problems
Académie universitaire Wallonie Europe Université de Liège Faculté des Sciences Appliquées Collège de doctorat en Électricité, électronique et informatique Multi-Domain Approaches for the Solution of High-Frequency
More informationSimulation results for magnetized plasmas
Chapter 4 Simulation results for magnetized plasmas In this chapter, we consider the dust charge fluctuation mode and lower hybrid wave damping in a magnetized plasma. Also, we consider plasma instabilities
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 informationProgress In Electromagnetics Research, PIER 52, , 2005 FDTD ANALYSIS OF MICROSTRIP PATCH ANTENNA COVERED BY PLASMA SHEATH
Progress In Electromagnetics Research, PIER 52, 173 183, 25 FDTD ANALYSIS OF MICROSTRIP PATCH ANTENNA COVERED BY PLASMA SHEATH Z. H. Qian and R. S. Chen Department of Communication Engineering Nanjing
More informationElectron-MHD Turbulence in Neutron Stars
Electron-MHD Turbulence in Neutron Stars Toby Wood Newcastle University (previously University of Leeds) Wood, Hollerbach & Lyutikov 2014, Physics of Plasmas 21, 052110 Outline What is a neutron star?
More informationA MATLAB GUI FOR SIMULATING THE PROPAGATION OF THE ELECTROMAGNETIC FIELD IN A 2-D INFINITE SPACE
A MATLAB GUI FOR SIMULATING THE PROPAGATION OF THE ELECTROMAGNETIC FIELD IN A 2-D INFINITE SPACE Ioana SĂRĂCUŢ Victor POPESCU Marina Dana ŢOPA Technical University of Cluj-Napoca, G. Bariţiu Street 26-28,
More informationEXAMINATION QUESTION PAPER
Faculty of Science and Technology EXAMINATION QUESTION PAPER Exam in: Fys-2009 Introduction to Plasma Physics Date: 20161202 Time: 09.00-13.00 Place: Åsgårdvegen 9 Approved aids: Karl Rottmann: Matematisk
More informationarxiv: v1 [physics.plasm-ph] 12 Nov 2017
Astronomy & Astrophysics manuscript no. PIC_DPR9 c ESO 2018 April 24, 2018 Double plasma resonance instability as a source of solar zebra emission J. Benáček 1 and M. Karlický 2 arxiv:1714281v1 [physics.plasm-ph]
More informationEnergy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas
Energy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas Yingda Cheng Andrew J. Christlieb Xinghui Zhong March 18, 2014 Abstract In this paper, we apply our recently developed energy-conserving
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationModel Order Reduction and Stability Enforcement of Finite-Difference Time-Domain Equations Beyond the CFL Limit. Xihao Li
Model Order Reduction and Stability Enforcement of Finite-Difference Time-Domain Equations Beyond the CFL Limit by Xihao Li A thesis submitted in conformity with the requirements for the degree of Master
More informationBerk-Breizman and diocotron instability testcases
Berk-Breizman and diocotron instability testcases M. Mehrenberger, N. Crouseilles, V. Grandgirard, S. Hirstoaga, E. Madaule, J. Petri, E. Sonnendrücker Université de Strasbourg, IRMA (France); INRIA Grand
More informationLower Hybrid Current Drive Experiments on Alcator C-Mod: Comparison with Theory and Simulation
Lower Hybrid Current Drive Experiments on Alcator C-Mod: Comparison with Theory and Simulation P.T. Bonoli, A. E. Hubbard, J. Ko, R. Parker, A.E. Schmidt, G. Wallace, J. C. Wright, and the Alcator C-Mod
More informationOptimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients
Optimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients Victorita Dolean, Martin Gander, Erwin Veneros To cite this version: Victorita Dolean, Martin Gander, Erwin Veneros. Optimized
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 informationSummer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations
1856-30 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The particle-in-cell M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The particle-in-cell
More informationEITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity
EITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity Daniel Sjöberg Department of Electrical and Information Technology Spring 2018 Outline 1 Basic reflection physics 2 Radar cross section definition
More informationHigh-order FDTD methods via derivative matching for Maxwell s. equations with material interfaces. Abstract
High-order FDTD methods via derivative matching for Maxwell s equations with material interfaces Shan Zhao 1 and G. W. Wei 1,2 1 Department of Mathematics, Michigan State University, East Lansing, MI 48824
More informationFINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC
FINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC WAVEGUIDES Chin-ping Yu (1) and Hung-chun Chang (2) (1) Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei,
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 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 informationOptical Properties of Left-Handed Materials by Nathaniel Ferraro 01
Optical Properties of Left-Handed Materials by Nathaniel Ferraro 1 Abstract Recently materials with the unusual property of having a simultaneously negative permeability and permittivity have been tested
More informationA PARTICLE-IN-CELL METHOD FOR THE SIMULATION OF PLASMAS BASED ON AN UNCONDITIONALLY STABLE WAVE EQUATION SOLVER. Eric Matthew Wolf
A PARTICLE-IN-CELL METHOD FOR THE SIMULATION OF PLASMAS BASED ON AN UNCONDITIONALLY STABLE WAVE EQUATION SOLVER By Eric Matthew Wolf A DISSERTATION Submitted to Michigan State University in partial fulfillment
More informationFinite-Difference Time-Domain and Beam Propagation Methods for Maxwell s Equations
Finite-Difference Time-Domain and Beam Propagation Methods for Maxwell s Equations Wolfgang Freude and Jan Brosi Institute of High-Frequency and Quantum Electronics (IHQ), University of Karlsruhe, Germany
More informationScattering of ECRF waves by edge density fluctuations and blobs
PSFC/JA-14-7 Scattering of ECRF waves by edge density fluctuations and blobs A. K. Ram and K. Hizanidis a June 2014 Plasma Science and Fusion Center, Massachusetts Institute of Technology Cambridge, MA
More informationSpherical Waves. Daniel S. Weile. Department of Electrical and Computer Engineering University of Delaware. ELEG 648 Spherical Coordinates
Spherical Waves Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Spherical Coordinates Outline Wave Functions 1 Wave Functions Outline Wave Functions 1
More informationIMPROVED SHIFT-OPERATOR FDTD METHOD FOR ANISOTROPIC MAGNETIZED COLD PLASMAS WITH ARBITRARY MAGNETIC FIELD DECLINATION
Progress In Electromagnetics Research B, Vol. 38, 39 56, 01 IMPROVED SHIFT-OPERATOR FDTD METHOD FOR ANISOTROPIC MAGNETIZED COLD PLASMAS WITH ARBITRARY MAGNETIC FIELD DECLINATION X. Yin *, H. Zhang, H.-Y.
More informationPre-Corrected FFT/AIM Algorithm for. Department of Electrical & Computer Engineering
A Multigrid Enhanced Pre-Corrected FF/AIM Algorithm for Multiscale Integral Equation Analysis K. Yang, F. Wei, and A. E. Yilmaz Department of Electrical & Computer Engineering University of exas at Austin
More informationfield using second order edge elements in 3D
The current issue and full text archive of this journal is available at http://www.emerald-library.com using second order edge elements in 3D Z. Ren Laboratoire de GeÂnie Electrique de Paris, UniversiteÂs
More informationSimulation Study of High-Frequency Magnetosonic Waves Excited by Energetic Ions in Association with Ion Cyclotron Emission )
Simulation Study of High-Frequency Magnetosonic Waves Excited by Energetic Ions in Association with Ion Cyclotron Emission ) Mieko TOIDA 1),KenjiSAITO 1), Hiroe IGAMI 1), Tsuyoshi AKIYAMA 1,2), Shuji KAMIO
More informationProgress In Electromagnetics Research B, Vol. 56, , 2013
Progress In Electromagnetics Research B, Vol. 56, 409 427, 2013 A COMBINED FDTD/TLM TIME DOMAIN METHOD TO SOLVE EFFICIENTLY ELECTROMAGNETIC PROB- LEMS Nathanael Muot 1, Christophe Girard 1, Xavier Ferrieres
More information