Full-wave simulation of tokamak plasma heating
|
|
- Curtis Lloyd
- 5 years ago
- Views:
Transcription
1 Full-wave simulation of tokamak plasma heating Institut Elie Cartan de Nancy, UMR 7502, Université de Lorraine, CNRS & INRIA (project-team CALVI). B.P , Vandœuvre les Nancy Cedex, France. Université Pierre et Marie Curie, 21 juin 2012
2 Joint work with: Jean R. Roche (IECN) Takashi Hattori (IECN) Christoph Kirsch (ex-iecn) Pierre Bertrand (Institut Jean Lamour, U. Lorraine) Jean-Hugues Chatenet (Institut Jean Lamour, U. Lorraine) Yves Peysson (CEA Cadarache)
3 Outline
4 Outline 1 Physical and mathematical modelling 1 Physical assumptions 2 Mixed and augmented formulations 3 Spectral properties of the dielectric tensor 4 Well-posedness of the variational formulations 5 Discretisation (Fourier Finite element) 6 A first numerical test
5 Outline 1 Physical and mathematical modelling 1 Physical assumptions 2 Mixed and augmented formulations 3 Spectral properties of the dielectric tensor 4 Well-posedness of the variational formulations 5 Discretisation (Fourier Finite element) 6 A first numerical test 2 The issue of accessibility 1 Electromagnetic dispersion relation 2 Accessibility criteria 3 More physical test cases
6 Outline 1 Physical and mathematical modelling 1 Physical assumptions 2 Mixed and augmented formulations 3 Spectral properties of the dielectric tensor 4 Well-posedness of the variational formulations 5 Discretisation (Fourier Finite element) 6 A first numerical test 2 The issue of accessibility 1 Electromagnetic dispersion relation 2 Accessibility criteria 3 More physical test cases 3 In progress and to do
7 Part One: Physical and mathematical modelling
8 The name of the game
9 The name of the game (simplified)
10 Physical assumptions Strong external, time-invariant magnetic field B ext. Harmonic time dependence of self-consistent electromagnetic fields (of the form e iωt, with ω > 0).
11 Physical assumptions Strong external, time-invariant magnetic field B ext. Harmonic time dependence of self-consistent electromagnetic fields (of the form e iωt, with ω > 0). Plasma response described by the matrix: S id 0 K = id S 0 + i γ e ε 0 0 P 0 ω }{{}}{{} ε r Landau term (1) Stix frame: third coordinate B ext.
12 Physical assumptions Strong external, time-invariant magnetic field B ext. Harmonic time dependence of self-consistent electromagnetic fields (of the form e iωt, with ω > 0). Plasma response described by the matrix: S id 0 K = id S 0 + i γ e ε 0 0 P 0 ω }{{}}{{} ε r Landau term (1) Stix frame: third coordinate B ext. Cold plasma approximation of ε r is used, valid in the lower hybrid (LH) regime ω ci ω ω ce. (More on this later.)
13 Physical assumptions Strong external, time-invariant magnetic field B ext. Harmonic time dependence of self-consistent electromagnetic fields (of the form e iωt, with ω > 0). Plasma response described by the matrix: S id 0 K = id S 0 + i γ e ε 0 0 P 0 ω }{{}}{{} ε r Landau term (1) Stix frame: third coordinate B ext. Cold plasma approximation of ε r is used, valid in the lower hybrid (LH) regime ω ci ω ω ce. (More on this later.) Both Landau and collisional damping present.
14 Electromagnetic waves: model problem The domain Ω is a torus (tokamak plasma volume). curl curl E (ω/c) 2 KE = 0 in Ω, (2) div(ke) = 0 in Ω. (3)
15 Electromagnetic waves: model problem The domain Ω is a torus (tokamak plasma volume). curl curl E (ω/c) 2 KE = 0 in Ω, (2) div(ke) = 0 in Ω. (3) Boundary condition: Ω = Γ A Γ C = antenna + the rest E n = 0 on Γ C. (Phoney?) (4)
16 Electromagnetic waves: model problem The domain Ω is a torus (tokamak plasma volume). curl curl E (ω/c) 2 KE = 0 in Ω, (2) div(ke) = 0 in Ω. (3) Boundary condition: Ω = Γ A Γ C = antenna + the rest E n = 0 on Γ C. (Phoney?) (4) On Γ A several possible models... Neumann: curl E n = iωµ 0 J s (5) Dirichlet: E n = E A n (6) Silver Müller: (E + i (c/ω) curl E n) n = E n (7) where J s = J s (0, 0, 1) is a current source equivalent to the electric field E A excited at the antenna.
17 Variational Formulations (1) Neumann BC on Γ A (fss) variational space for (2): H C 0 (curl; Ω) = {v H(curl; Ω) : v n = 0 on Γ C} Eq. (3): constraint mixed (unaugmented) formulation: Find (E, p) H C 0 (curl; Ω) H1 0 (Ω) or H1 0,C (Ω) such that a pl (E, F ) + b pl (F, p) = L(F ), F H C 0 (curl; Ω); (8) b pl (E, q) = 0, q H0,(C) 1 (Ω), (9) where H0,C 1 (Ω) = {w H1 (Ω) : w = 0 on Γ C } and a pl (E, F ) := (curl E curl F ) (ω/c) 2 (KE F ), b pl (F, p) := (KF grad p), L(F ) := iωµ 0 J s, F H 1/2 (Γ A ),H 1/2 (Γ A ). Above, ( ) = scalar product in L 2 (Ω).
18 Variational Formulations (2) Other approaches: consider the space X C N (K; Ω) := HC 0 (curl; Ω) H(div K; Ω) Mixed augmented formulation: Find (E, p) X C N (K; Ω) L2 (Ω) such that a s (E, F ) + b(f, p) = L(F ), F X C N (K; Ω); (10) b(e, q) = 0, q L 2 (Ω), (11) where a s (E, F ) := (curl E curl F ) + s (div KE div KF ) (ω/c) 2 (KE F ), b(f, p) := (div KF p).
19 Variational Formulations (2) Other approaches: consider the space X C N (K; Ω) := HC 0 (curl; Ω) H(div K; Ω) Mixed augmented formulation: Find (E, p) X C N (K; Ω) L2 (Ω) such that a s (E, F ) + b(f, p) = L(F ), F X C N (K; Ω); (10) b(e, q) = 0, q L 2 (Ω), (11) where a s (E, F ) := (curl E curl F ) + s (div KE div KF ) (ω/c) 2 (KE F ), b(f, p) := (div KF p). Augmented formulation: Find E X C N (K; Ω) such that a s (E, F ) = L(F ), F X C N (K; Ω). (12)
20 More on plasma response (1) Experimentalists / coupling with other codes maps of external magnetic field B ext (x), density n s (x) and temperature T s (x) of each species s (electrons e, various ions).
21 More on plasma response (1) Experimentalists / coupling with other codes maps of external magnetic field B ext (x), density n s (x) and temperature T s (x) of each species s (electrons e, various ions). Compute maps of plasma and cyclotron frequencies: n s qs 2 ω ps =, ω cs = q s B ext, ε 0 m s m s then plasma parameter and collision frequency Λ = 4π 3 n e ( ) ε0 k B T 3/2 e n e qe 2 ; ν c = ω pe ln Λ. Λ
22 More on plasma response (2) Setting α(x) = ω + iν c (x), the entries of ε r are: S = 1 α ωps 2 ω α s 2 ωcs 2 (13) D = 1 ω cs ωps 2 ω α s 2 ωcs 2 (14) P = 1 1 ωps 2 (15) ωα s
23 More on plasma response (2) Setting α(x) = ω + iν c (x), the entries of ε r are: S = 1 α ωps 2 ω α s 2 ωcs 2 (13) D = 1 ω cs ωps 2 ω α s 2 ωcs 2 (14) P = 1 1 ωps 2 (15) ωα Landau damping rate (phenomenological): ( ) π ω 2 n e qe 2 m 1/2 e γ e = 4(k B T e ) 3/2 k 3 exp ω2 m e 4k B T e k 2. (16) Relies on plane wave analysis: k = k B ext / B ext. s
24 More on plasma response (2) Setting α(x) = ω + iν c (x), the entries of ε r are: S = 1 α ωps 2 ω α s 2 ωcs 2 (13) D = 1 ω cs ωps 2 ω α s 2 ωcs 2 (14) P = 1 1 ωps 2 (15) ωα Landau damping rate (phenomenological): ( ) π ω 2 n e qe 2 m 1/2 e γ e = 4(k B T e ) 3/2 k 3 exp ω2 m e 4k B T e k 2. (16) Relies on plane wave analysis: k = k B ext / B ext. If ν c > 0 or γ e > 0, K is not hermitian. s
25 More on plasma response (2) Setting α(x) = ω + iν c (x), the entries of ε r are: S = 1 α ωps 2 ω α s 2 ωcs 2 (13) D = 1 ω cs ωps 2 ω α s 2 ωcs 2 (14) P = 1 1 ωps 2 (15) ωα Landau damping rate (phenomenological): ( ) π ω 2 n e qe 2 m 1/2 e γ e = 4(k B T e ) 3/2 k 3 exp ω2 m e 4k B T e k 2. (16) Relies on plane wave analysis: k = k B ext / B ext. If ν c > 0 or γ e > 0, K is not hermitian. Actually, a relief! s
26 Spectral properties of K (1) Eigenvalues and uniform (in x), unitary eigenvector basis λ 1 = S + D, λ 2 = S D, λ 3 = P + i γ e ε 0 ω := P L (17) i u 1 = 1 2, u 2 = i, u 3 = 0. (18)
27 Spectral properties of K (1) Eigenvalues and uniform (in x), unitary eigenvector basis λ 1 = S + D, λ 2 = S D, λ 3 = P + i γ e ε 0 ω := P L (17) i u 1 = 1 2, u 2 = i, u 3 = 0. (18) There always holds: Im λ i (x) > 0, hence: Im(v K(x)v) ζ(x) v v, v C 3. Suitable bounds on n s, T s, B ext : v C 3, x Ω, v K(x)v Im(v K(x)v) ζ v 2. (19)
28 Spectral properties of K (2) On the contrary, in the LH frequency range: Re λ 1 (x) > 0, Re λ 2,3 (x) < 0. Without the imaginary parts (damping), the well-posedness proof breaks down!!! (Still OK if ν c = 0 or γ e = 0, but not both.)
29 Well-posedness of the variational formulations (1) Lemma (Sébelin, Bénilan 97) Let V H two Hilbert spaces, with continuous embedding. Let a(, ) be a sesquilinear form on V satisfying the inequalities: Re a(v, v) α v 2 V λ v 2 H, Im a(v, v) γ v 2 H, (20) for some strictly positive constants α, γ, λ and all v V. Then a(, ) is V -elliptic, i.e., a(v, v) α v 2 V.
30 Well-posedness of the variational formulations (1) Lemma (Sébelin, Bénilan 97) Let V H two Hilbert spaces, with continuous embedding. Let a(, ) be a sesquilinear form on V satisfying the inequalities: Re a(v, v) α v 2 V λ v 2 H, Im a(v, v) γ v 2 H, (20) for some strictly positive constants α, γ, λ and all v V. Then a(, ) is V -elliptic, i.e., a(v, v) α v 2 V. With the bound (19): Im(Kv v) ζ v L 2, v L 2 (Ω). a pl (u, v) = (curl u curl v) (ω/c) 2 (Ku v) satisfies (20) on V = H(curl; Ω), with H = L 2 (Ω);
31 Well-posedness of the variational formulations (1) Lemma (Sébelin, Bénilan 97) Let V H two Hilbert spaces, with continuous embedding. Let a(, ) be a sesquilinear form on V satisfying the inequalities: Re a(v, v) α v 2 V λ v 2 H, Im a(v, v) γ v 2 H, (20) for some strictly positive constants α, γ, λ and all v V. Then a(, ) is V -elliptic, i.e., a(v, v) α v 2 V. With the bound (19): Im(Kv v) ζ v L 2, v L 2 (Ω). a pl (u, v) = (curl u curl v) (ω/c) 2 (Ku v) satisfies (20) on V = H(curl; Ω), with H = L 2 (Ω); a s (u, v) = a pl (u, v) + s (div Ku div Kv) satisfies (20) on V = H(curl, div K; Ω) if Re s > 0 and Im s 0.
32 Well-posedness of the variational formulations (1) Lemma (Sébelin, Bénilan 97) Let V H two Hilbert spaces, with continuous embedding. Let a(, ) be a sesquilinear form on V satisfying the inequalities: Re a(v, v) α v 2 V λ v 2 H, Im a(v, v) γ v 2 H, (20) for some strictly positive constants α, γ, λ and all v V. Then a(, ) is V -elliptic, i.e., a(v, v) α v 2 V. With the bound (19): Notice: Im(Kv v) ζ v L 2, v L 2 (Ω). a pl (u, v) = (curl u curl v) (ω/c) 2 (Ku v) satisfies (20) on V = H(curl; Ω), with H = L 2 (Ω); a s (u, v) = a pl (u, v) + s (div Ku div Kv) satisfies (20) on V = H(curl, div K; Ω) if Re s > 0 and Im s 0. coercivity ( Fredholm alternative); no BC required.
33 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram.
34 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram. (K grad u grad v) elliptic on H0 1 (Ω) well-posedness of K u := div(k grad u) = f H 1 (Ω).
35 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram. (K grad u grad v) elliptic on H0 1 (Ω) well-posedness of K u := div(k grad u) = f H 1 (Ω). "Warped" Hodge decomposition for all u L 2 (Ω): u = grad ϕ + u T, with: ϕ H 1 0 (Ω), div Ku T = 0.
36 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram. (K grad u grad v) elliptic on H0 1 (Ω) well-posedness of K u := div(k grad u) = f H 1 (Ω). "Warped" Hodge decomposition for all u L 2 (Ω): u = grad ϕ + u T, with: ϕ H 1 0 (Ω), div Ku T = 0. inf-sup condition for b pl (v, p) = (Kv grad p) on H 0 (curl; Ω) H 1 0 (Ω)
37 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram. (K grad u grad v) elliptic on H0 1 (Ω) well-posedness of K u := div(k grad u) = f H 1 (Ω). "Warped" Hodge decomposition for all u L 2 (Ω): u = grad ϕ + u T, with: ϕ H 1 0 (Ω), div Ku T = 0. inf-sup condition for b pl (v, p) = (Kv grad p) on H 0 (curl; Ω) H 1 0 (Ω) and HC 0 (curl; Ω) H1 0,(C) (Ω)
38 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram. (K grad u grad v) elliptic on H0 1 (Ω) well-posedness of K u := div(k grad u) = f H 1 (Ω). "Warped" Hodge decomposition for all u L 2 (Ω): u = grad ϕ + u T, with: ϕ H 1 0 (Ω), div Ku T = 0. inf-sup condition for b pl (v, p) = (Kv grad p) on H 0 (curl; Ω) H 1 0 (Ω) and HC 0 (curl; Ω) H1 0,(C) (Ω) inf-sup condition for b(v, p) = (div Kv p) on X N (K; Ω) L 2 (Ω) and X C N (K; Ω) L2 (Ω)
39 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram. (K grad u grad v) elliptic on H0 1 (Ω) well-posedness of K u := div(k grad u) = f H 1 (Ω). "Warped" Hodge decomposition for all u L 2 (Ω): u = grad ϕ + u T, with: ϕ H 1 0 (Ω), div Ku T = 0. inf-sup condition for b pl (v, p) = (Kv grad p) on H 0 (curl; Ω) H 1 0 (Ω) and HC 0 (curl; Ω) H1 0,(C) (Ω) inf-sup condition for b(v, p) = (div Kv p) on X N (K; Ω) L 2 (Ω) and X C N (K; Ω) L2 (Ω) Mixed unaugmented (8) (9) and mixed augmented (10) (11) formulations well-posed
40 Well-posedness of the variational formulations (2) Augmented formulation (12) well-posed on the closed subspace X C N (K; Ω) by Lax Milgram. (K grad u grad v) elliptic on H0 1 (Ω) well-posedness of K u := div(k grad u) = f H 1 (Ω). "Warped" Hodge decomposition for all u L 2 (Ω): u = grad ϕ + u T, with: ϕ H 1 0 (Ω), div Ku T = 0. inf-sup condition for b pl (v, p) = (Kv grad p) on H 0 (curl; Ω) H 1 0 (Ω) and HC 0 (curl; Ω) H1 0,(C) (Ω) inf-sup condition for b(v, p) = (div Kv p) on X N (K; Ω) L 2 (Ω) and X C N (K; Ω) L2 (Ω) Mixed unaugmented (8) (9) and mixed augmented (10) (11) formulations well-posed Equivalence of various formulations
41 Dimension reduction (1)
42 Dimension reduction (2) The domain Ω is described by cylindrical coordinates (R, Z, φ), which are related to the cartesian ones (x, y, z) by x = R cos φ, y = R sin φ, z = Z. The 3D problem is reduced to a series of 2D ones by expanding all fields u(r, Z, φ) as Fourier series in φ: u(r, Z, φ) = ν Z u ν (R, Z ) e iνφ, (21) where the coefficients u ν (R, Z ) are defined on a meridian section Ω := {(R, Z ) R 2 (R, Z, π/2) Ω}. Furthermore: curl u = ν Z (curl ν u ν ) e iνφ, div u = ν Z(div ν u ν ) e iνφ, (22) where curl ν, div ν are obtained by setting φ iν in the expressions of curl ν, div ν in cylindrical coordinates.
43 Dimension reduction (3) Simplifying assumptions: External magnetic field B ext (x) = B ext (R, Z ) e φ. Particles constrained to move and diffuse freely along field lines of B ext
44 Dimension reduction (3) Simplifying assumptions: External magnetic field B ext (x) = B ext (R, Z ) e φ. Particles constrained to move and diffuse freely along field lines of B ext plasma properties averaged in φ n s, T s, hence S, D, P, are independent of φ.
45 Dimension reduction (3) Simplifying assumptions: External magnetic field B ext (x) = B ext (R, Z ) e φ. Particles constrained to move and diffuse freely along field lines of B ext plasma properties averaged in φ n s, T s, hence S, D, P, are independent of φ. The r.h.s. J s has a single Fourier mode: J s = J ν 0 s assign a definite value to k (x) = ν 0 /R in (16) the whole matrix K is independent of φ. As a consequence: e iν 0φ div Ku = ν Z(div ν Ku ν ) e iνφ, and, by orthogonality: a s (u, v) = ν Z a s,(ν) (u ν, v ν ), b(v, p) = ν Z b (ν) (v ν, p ν ).
46 Discretisation: Fourier Finite Element Method A 2D problem is obtained for each pair of Fourier coefficients (E ν, p ν ), ν Z, set in the spaces X C N (K; Ω) L 2 ( Ω), traces on Ω of the 3D spaces: a s,(ν) (E ν, F ) + b (ν) (F, p ν ) = L ν (F ), F X C N (K; Ω); (23) b (ν) (E ν, q) = 0, q L 2 ( Ω). (24) The r.h.s. and the solution are zero except when ν = ν 0.
47 Discretisation: Fourier Finite Element Method A 2D problem is obtained for each pair of Fourier coefficients (E ν, p ν ), ν Z, set in the spaces X C N (K; Ω) L 2 ( Ω), traces on Ω of the 3D spaces: a s,(ν) (E ν, F ) + b (ν) (F, p ν ) = L ν (F ), F X C N (K; Ω); (23) b (ν) (E ν, q) = 0, q L 2 ( Ω). (24) The r.h.s. and the solution are zero except when ν = ν 0. We use a Taylor-Hood P 2 -iso-p 1 finite element (continuous fields and multipliers) to discretize the MAVF (23) (24). ( As,(ν) B ) ( ) ( ) (ν) E h Lν ν B (ν) 0 p ν 2h =, (25) 0 where E ν h and p ν 2h nodal values of E ν and p ν. Notice that the matrix of the system (25) is complex-valued and non-hermitian.
48 Essential boundary conditions In a local basis {n, τ, e φ = n τ }, we have: E 1 B 1 0 E n = EB 2 0 = E 3 EB 3 B. 0 EB 2 The boundary condition E n = 0 on Γ C becomes: (EB 2, E B 3 ) = (0, 0). (26) Denoting P = transformation matrix between {e R, e Z, e φ } and {n, τ, n τ } on Γ C nodes, identity elsewhere: v M u = v f v B (P MP) u B = v B (P f ). Thus, one sets (E 2 B, E 3 B ) = (0, 0) in P MP and P f. Generalisation to E n = E A n on Γ A obvious.
49 Numerical test We consider a 2D test problem analogous to (23) (24), with K = I (and s = 1 in the augmented formulation): curl ν curl ν E η 2 E = f in Ω, (27) div ν E = g in Ω, (28) E n = 0 on Γ C, (29) curl ν E n = h on Γ A. (30) The r.h.s. f, g, h are chosen so that the exact solution is E(R, Z ) = (R R 0, Z, 0) H 0 (curl, Ω) H(div, Ω) and ν = 220. This gives rise to six independent tests, for the real and imaginary parts of the three coordinates of E.
50 Figure: The real part of the test solution E.
51 Figure: The imaginary part of the test solution E.
52 Part Two: The issue of accessibility
53 Dispersion relation (1) We look for plane wave solutions of the form Ee i(k x ωt) to Maxwell s equations in the plasma. we set E = Ee ik x in the propagation equation (2): curl curl E (ω/c) 2 KE = 0. No BC is assumed (propagation in an infinite medium).
54 Dispersion relation (1) We look for plane wave solutions of the form Ee i(k x ωt) to Maxwell s equations in the plasma. we set E = Ee ik x in the propagation equation (2): curl curl E (ω/c) 2 KE = 0. No BC is assumed (propagation in an infinite medium). Therefore: ik (ik E) (ω/c) 2 K E = 0. Multiplying by (c/ω) 2, we get the dispersion relation n (n E) + KE = 0. (31) where n = c ω k is the vector index of refraction.
55 Dispersion relation (2) Take a (local) Stix frame (e 1, e 2, e 3 ), with e 3 in the direction of B ext and k e 2 = 0. Thus: n = (n, 0, n ).
56 Dispersion relation (2) Take a (local) Stix frame (e 1, e 2, e 3 ), with e 3 in the direction of B ext and k e 2 = 0. Thus: n = (n, 0, n ). The dispersion relation (31) is rewritten as: ME = 0, (32) where: M = S n 2 id n n id S n 2 n2 0 n n 0 P L n 2.
57 Dispersion relation (2) Take a (local) Stix frame (e 1, e 2, e 3 ), with e 3 in the direction of B ext and k e 2 = 0. Thus: n = (n, 0, n ). The dispersion relation (31) is rewritten as: ME = 0, (32) where: M = S n 2 id n n id S n 2 n2 0 n n 0 P L n 2. To have a non-trivial solution, there must hold det M = 0: Sn 4 [(S + P L)(S n 2 ) D2 ]n 2 + P L[(S n 2 )2 D 2 ] = 0. (33)
58 Accessibility (1) To simplify the discussion, we assume there is no damping (ν c = γ e = 0). Effect of actual damping on this analysis: A propagating wave will still propagate (with a little attenuation). An evanescent wave will not become propagating. Therefore, the entries S, D, P L = P of K are real.
59 Accessibility (1) To simplify the discussion, we assume there is no damping (ν c = γ e = 0). Effect of actual damping on this analysis: A propagating wave will still propagate (with a little attenuation). An evanescent wave will not become propagating. Therefore, the entries S, D, P L = P of K are real. Under this assumption, we say there is strict accessibility if n R 3 solution to (33); vague accessibility if n C 3 solution to (33), with Im n Re n.
60 Accessibility (1) To simplify the discussion, we assume there is no damping (ν c = γ e = 0). Effect of actual damping on this analysis: A propagating wave will still propagate (with a little attenuation). An evanescent wave will not become propagating. Therefore, the entries S, D, P L = P of K are real. Under this assumption, we say there is Recall: strict accessibility if n R 3 solution to (33); vague accessibility if n C 3 solution to (33), with Im n Re n. in the LH regime, there holds: P < 0 < S < D and S, D P. (34)
61 Accessibility (2) We fix n R (Phenomenology...) and write r = S n 2. Then, det M is a second degree polynomial in n 2 : det M = Q(n 2 ) = a r n 4 b r n 2 + c r, where: a r = S, b r = (S + P)r D 2, c r = P(r 2 D 2 ).
62 Accessibility (2) We fix n R (Phenomenology...) and write r = S n 2. Then, det M is a second degree polynomial in n 2 : det M = Q(n 2 ) = a r n 4 b r n 2 + c r, where: a r = S, b r = (S + P)r D 2, c r = P(r 2 D 2 ). Let r := b 2 r 4 a r c r be the discriminant of Q. 1 If r 0 and Q has a non-negative root, there is strict accessibility. 2 If r 0 and both roots are negative, there is strict inaccessibility (all waves are evanescent). 3 If r < 0, there cannot be strict accessibility, but there may be vague accessibility.
63 Strict Accessibility One calculates: r = (S P) 2 r 2 2D 2 (S + P) r + D 2 (D 2 + 4SP). 1 r 0 r / (r, r + ), with: r ± := D2 (S + P) ± δ (S P) 2, δ := 4D 2 SP(D 2 (S P) 2 ). Actually, in the LH regime (34), δ > 0.
64 Strict Accessibility One calculates: r = (S P) 2 r 2 2D 2 (S + P) r + D 2 (D 2 + 4SP). 1 r 0 r / (r, r + ), with: r ± := D2 (S + P) ± δ (S P) 2, δ := 4D 2 SP(D 2 (S P) 2 ). Actually, in the LH regime (34), δ > 0. 2 Both roots negative (strict inaccessibility) iff b r a r < 0 and c r a r > 0 r > r := D2 S + P and r D.
65 The case r < 0 We have vague accessibility iff n C s.t. det M = 0 and Re(n ) Im(n ).
66 The case r < 0 We have vague accessibility iff n C s.t. det M = 0 and Re(n ) Im(n ). We calculate: n 2 = (b r ± i r )/(2a r ) = Re(n ) 2 = b r + br 2 r, Im(n ) 2 = b r + b 2 r r. 4a r 4a r Necessary condition: b r 0.
67 The case r < 0 We have vague accessibility iff n C s.t. det M = 0 and Re(n ) Im(n ). We calculate: n 2 = (b r ± i r )/(2a r ) = Re(n ) 2 = b r + br 2 r, Im(n ) 2 = b r + b 2 r r. 4a r 4a r Necessary condition: b r 0. In this case: Re(n ) Im(n ) = κ κ, where: κ := b2 r r. Thus: vague accessibility b r 0 and κ large enough
68 The case r < 0 We have vague accessibility iff n C s.t. det M = 0 and Re(n ) Im(n ). We calculate: n 2 = (b r ± i r )/(2a r ) = Re(n ) 2 = b r + br 2 r, Im(n ) 2 = b r + b 2 r r. 4a r 4a r Necessary condition: b r 0. In this case: Re(n ) Im(n ) = κ κ, where: κ := b2 r r. Thus: vague accessibility b r 0 and κ large enough r r and b2 r r large enough.
69 The case r < 0 (continued) We investigate the condition ɛ 2 b 2 r + r 0, for ɛ > 0 given. The l.h.s. is a second-degree polynomial in r, it is positive outside its roots: ρ ɛ ± = (ɛ 2 + 1)(S + P)D 2 ± δ ɛ (ɛ 2 + 1)(S 2 + P 2 ) 2(ɛ 2 1)SP, with: δ ɛ = 4D 2 SP[(ɛ 2 + 1)(D 2 S 2 P 2 ) 2(ɛ 2 1)SP]. In the LH regime (34), δ ɛ > 0 and ρ ɛ r ρ ɛ +.
70 Accessibility: conclusions Proposition 1 Assume ω ci ω ω ce and ν c = γ e = 0. Then for ɛ small enough. r < ρ ɛ < r < ρ ɛ + < r + < S, 2 Let λ = D2 P. One has strict accessibility iff n > S r S + 2 Sλ + λ. One has vague accessibility ( Re(n ) / Im(n ) > K ) iff S ρ ɛ < n < S r, with ɛ = 2K /(K 2 1).
71 "Physical" test case We consider an electron-proton plasma with densities n e = n i = m 3 Plasma and cyclotron frequencies (rad/s): (constant). ω pe , ω pi , ω ce , ω ci , for an external magnetic field of intensity B ext = 3.2 T. Lower hybrid frequency ω LH = ω ci ω ce rad/s. Electron temperature: k B T e = 1 kev ν c, γ e.
72 "Physical" test case We consider an electron-proton plasma with densities n e = n i = m 3 Plasma and cyclotron frequencies (rad/s): (constant). ω pe , ω pi , ω ce , ω ci , for an external magnetic field of intensity B ext = 3.2 T. Lower hybrid frequency ω LH = ω ci ω ce rad/s. Electron temperature: k B T e = 1 kev ν c, γ e. We present a series of tests on the physical problem (2) (5) and its 2D VF (23) (24), for the values ω = ω LH, 2ω LH or ω LH /2 of the frequency, and several values of the Fourier index ν.
73 Influence of parameters Parallel wave vector: k (x) = ν R or n (x) = νc Rω.
74 Influence of parameters Parallel wave vector: k (x) = ν R or n (x) = νc Rω. ω = const., ν = { more accessible γ e 0 for ν large enough,
75 Influence of parameters Parallel wave vector: k (x) = ν R or n (x) = νc Rω. ω = const., ν = ν = const., ω = { more accessible γ e 0 for ν large enough, { less accessible γ e 0 for ω large enough. Note: k ν/r 0 in γ e = K = constant.
76 Influence of parameters Parallel wave vector: k (x) = ν R or n (x) = νc Rω. ω = const., ν = ν = const., ω = { more accessible γ e 0 for ν large enough, { less accessible γ e 0 for ω large enough. Note: k ν/r 0 in γ e = K = constant. On the other hand we should have an a priori estimate E ν C ν α rhs (in the suitable spaces). So, for a fixed ω, there is an optimal ν where the wave emitted by the antenna penetrates most the plasma.
77 Figure: The real part of E for ω = ω LH and ν = 200
78 Figure: The real part of E for ω = ω LH and ν = 300
79 Figure: The real part of E for ω = ω LH and ν = 400
80 Figure: The real part of E for ω = ω LH /2 and ν = 100
81 Figure: The real part of E for ω = ω LH /2 and ν = 200
82 Figure: The real part of E for ω = ω LH /2 and ν = 300
83 Figure: The real part of E for ω = 2 ω LH and ν = 400
84 Figure: The real part of E for ω = 2 ω LH and ν = 600
85 Figure: The real part of E for ω = 2 ω LH and ν = 700
86 Figure: The real part of E for ω = 2 ω LH and ν = 800
87 Provisional conclusions
88 In progress and to do Feasibility of simulations incorporating most of the physics
89 In progress and to do Feasibility of simulations incorporating most of the physics (to do: variable plasma properties) Simulations confirm accessibility analysis
90 In progress and to do Feasibility of simulations incorporating most of the physics (to do: variable plasma properties) Simulations confirm accessibility analysis Non-hermitianness crucial for well-posedness, but precludes the use of many numerical methods
91 In progress and to do Feasibility of simulations incorporating most of the physics (to do: variable plasma properties) Simulations confirm accessibility analysis Non-hermitianness crucial for well-posedness, but precludes the use of many numerical methods To do: comparison of antenna boundary conditions
92 In progress and to do Feasibility of simulations incorporating most of the physics (to do: variable plasma properties) Simulations confirm accessibility analysis Non-hermitianness crucial for well-posedness, but precludes the use of many numerical methods To do: comparison of antenna boundary conditions More realistic test-cases: transversal wave length domain size Domain decomposition
93 In progress and to do Feasibility of simulations incorporating most of the physics (to do: variable plasma properties) Simulations confirm accessibility analysis Non-hermitianness crucial for well-posedness, but precludes the use of many numerical methods To do: comparison of antenna boundary conditions More realistic test-cases: transversal wave length domain size Domain decomposition Non-definite and non-hermitian operator make subdomains small enough to use a direct method
WaFu Notes Discussions around the cold plasma model
WaFu Notes Discussions around the cold plasma model Lise-Marie Imbert-Gérard Summer 7 These notes correspond - more or less - to the presentation I gave at the WaFu summer school on July 6th, 7 in Paris.
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 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 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 DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
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 informationPart VIII. Interaction with Solids
I with Part VIII I with Solids 214 / 273 vs. long pulse is I with Traditional i physics (ICF ns lasers): heating and creation of long scale-length plasmas Laser reflected at critical density surface Fast
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 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 informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
More informationSpectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem
Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v
More informationPropagation of Radio Frequency Waves Through Density Filaments
PSFC/JA-15-13 Propagation of Radio Frequency Waves Through Density Filaments A. K. Ram and K. Hizanidis a May 015 a National Technical University of Athens (part of HELLAS) School of Electrical and Computer
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationWaveguide Propagation Modes and Quadratic Eigenvalue Problems
Waveguide Propagation Modes and Quadratic Eigenvalue Problems Prof. André Nicolet, Institut Fresnel, UMR CNRS 6133, Université Paul Cézanne, Marseille, France. Prof. Christophe Geuzaine, Mathematics Department,
More informationRELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE
RELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE Abhay K. Ram Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA 02139. U.S.A. Joan Decker
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 informationTwo ion species studies in LAPD * Ion-ion Hybrid Alfvén Wave Resonator
Two ion species studies in LAPD * Ion-ion Hybrid Alfvén Wave Resonator G. J. Morales, S. T. Vincena, J. E. Maggs and W. A. Farmer UCLA Experiments performed at the Basic Plasma Science Facility (BaPSF)
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 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 informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationVolume and surface integral equations for electromagnetic scattering by a dielectric body
Volume and surface integral equations for electromagnetic scattering by a dielectric body M. Costabel, E. Darrigrand, and E. H. Koné IRMAR, Université de Rennes 1,Campus de Beaulieu, 35042 Rennes, FRANCE
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 informationMicroscopic-Macroscopic connection. Silvana Botti
relating experiment and theory European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre for Computational
More informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationTrefftz-Discontinuous Galerkin Methods for Maxwell s Equations
Trefftz-Discontinuous Galerkin Methods for Maxwell s Equations Ilaria Perugia Faculty of Mathematics, University of Vienna Vienna P D E Ralf Hiptmair (ETH Zürich), Andrea Moiola (University of Reading)
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
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 informationNumerical Analysis of Nonlinear Multiharmonic Eddy Current Problems
Numerical Analysis of Nonlinear Multiharmonic Eddy Current Problems F. Bachinger U. Langer J. Schöberl April 2004 Abstract This work provides a complete analysis of eddy current problems, ranging from
More informationarxiv: v1 [math.ap] 3 Feb 2015
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher arxiv:1502.00954v1 [math.ap] 3 Feb 2015 ELECTROMAGNETIC WAVE PROPAGATION AND ABSORPTION
More informationThe Interior Transmission Eigenvalue Problem for Maxwell s Equations
The Interior Transmission Eigenvalue Problem for Maxwell s Equations Andreas Kirsch MSRI 2010 epartment of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research
More informationIntegration of Fokker Planck calculation in full wave FEM simulation of LH waves
Integration of Fokker Planck calculation in full wave FEM simulation of LH waves O. Meneghini S. Shiraiwa R. Parker 51 st DPP APS, Atlanta November 4, 29 L H E A F * Work supported by USDOE awards DE-FC2-99ER54512
More informationNumerical analysis of problems in electromagnetism
Numerical analysis of problems in electromagnetism ALBERTO VALLI Department of Mathematics, University of Trento Numerical analysis of problems in electromagnetism p.1/195 Finite elements A finite element
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationEstimates for the resolvent and spectral gaps for non self-adjoint operators. University Bordeaux
for the resolvent and spectral gaps for non self-adjoint operators 1 / 29 Estimates for the resolvent and spectral gaps for non self-adjoint operators Vesselin Petkov University Bordeaux Mathematics Days
More informationLecture11: Plasma Physics 1. APPH E6101x Columbia University
Lecture11: Plasma Physics 1 APPH E6101x Columbia University 1 Last Lecture Introduction to plasma waves Basic review of electromagnetic waves in various media (conducting, dielectric, gyrotropic, ) Basic
More informationMULTIGRID METHODS FOR MAXWELL S EQUATIONS
MULTIGRID METHODS FOR MAXWELL S EQUATIONS A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements
More informationTrefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations
Trefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with
More informationThe mathematics of scattering by unbounded, rough, inhomogeneous layers
The mathematics of scattering by unbounded, rough, inhomogeneous layers Simon N. Chandler-Wilde a Peter Monk b Martin Thomas a a Department of Mathematics, University of Reading, Whiteknights, PO Box 220
More informationE E D E=0 2 E 2 E (3.1)
Chapter 3 Constitutive Relations Maxwell s equations define the fields that are generated by currents and charges. However, they do not describe how these currents and charges are generated. Thus, to find
More informationBasics of electromagnetic response of materials
Basics of electromagnetic response of materials Microscopic electric and magnetic field Let s point charge q moving with velocity v in fields e and b Force on q: F e F qeqvb F m Lorenz force Microscopic
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationMixed exterior Laplace s problem
Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationMaxwell s equations in Carnot groups
Maxwell s equations in Carnot groups B. Franchi (U. Bologna) INDAM Meeting on Geometric Control and sub-riemannian Geometry Cortona, May 21-25, 2012 in honor of Andrey Agrachev s 60th birthday Researches
More informationThe Factorization Method for a Class of Inverse Elliptic Problems
1 The Factorization Method for a Class of Inverse Elliptic Problems Andreas Kirsch Mathematisches Institut II Universität Karlsruhe (TH), Germany email: kirsch@math.uni-karlsruhe.de Version of June 20,
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationElectron-Acoustic Wave in a Plasma
Electron-Acoustic Wave in a Plasma 0 (uniform ion distribution) For small fluctuations, n ~ e /n 0
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
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 informationA Generalization for Stable Mixed Finite Elements
A Generalization for Stable Mixed Finite Elements Andrew Gillette joint work with Chandrajit Bajaj Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin,
More informationFull-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 informationA computer-assisted Band-Gap Proof for 3D Photonic Crystals
A computer-assisted Band-Gap Proof for 3D Photonic Crystals Henning Behnke Vu Hoang 2 Michael Plum 2 ) TU Clausthal, Institute for Mathematics, 38678 Clausthal, Germany 2) Univ. of Karlsruhe (TH), Faculty
More informationElectromagnetic (EM) Waves
Electromagnetic (EM) Waves Short review on calculus vector Outline A. Various formulations of the Maxwell equation: 1. In a vacuum 2. In a vacuum without source charge 3. In a medium 4. In a dielectric
More informationElectromagnetic Waves Across Interfaces
Lecture 1: Foundations of Optics Outline 1 Electromagnetic Waves 2 Material Properties 3 Electromagnetic Waves Across Interfaces 4 Fresnel Equations 5 Brewster Angle 6 Total Internal Reflection Christoph
More informationDiscontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints
J A N U A R Y 0 1 8 P R E P R N T 4 7 4 Discontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints gor Voulis * and Arnold Reusken nstitut für Geometrie und Praktische
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
More informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
More information2D full wave analysis of wave structure by TASK/WF2
17th NEXT (Numerical Experiment of Tokamak) Meeting, University of Tokyo, Kashiwa, Japan, March 15-16, 2012 2D full wave analysis of wave structure by TASK/WF2 Y. Maruyama, A. Fukuyama Department of Nuclear
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationDerivation of Appleton s equation Aug 2010 W. Gekelman
Derivation of Appleton s equation Aug 010 W. Gekelman This is a derivation of Appleton s equation, which is the equation for the index of refraction of a cold plasma for whistler waves. The index of refraction
More informationPlasma heating in stellarators at the fundamental ion cyclotron frequency
PHYSICS OF PLASMAS VOLUME 7, NUMBER FEBRUARY 000 Plasma heating in stellarators at the fundamental ion cyclotron frequency V. A. Svidzinski and D. G. Swanson Department of Physics, Auburn University, Auburn,
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationIntroduction to Polarization
Phone: Ext 659, E-mail: hcchui@mail.ncku.edu.tw Fall/007 Introduction to Polarization Text Book: A Yariv and P Yeh, Photonics, Oxford (007) 1.6 Polarization States and Representations (Stokes Parameters
More information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationEDP with strong anisotropy : transport, heat, waves equations
EDP with strong anisotropy : transport, heat, waves equations Mihaï BOSTAN University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Nachos team INRIA Sophia Antipolis, 3/07/2017 Main goals Effective
More informationNumerical Methods and Algorithms for High Frequency Wave Scattering Problems in Homogeneous and Random Media
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Doctoral Dissertations Graduate School 8-014 Numerical Methods and Algorithms for High Frequency Wave Scattering Problems
More informationMethod of Homogenization for the Study of the Propagation of Electromagnetic Waves in a Composite Part 2: Homogenization
, July 5-7, 2017, London, U.K. Method of Homogenization for the Study of the Propagation of Electromagnetic Waves in a Composite Part 2: Homogenization Helene Canot, Emmanuel Frenod Abstract In this paper
More informationWater-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes
Water-bag reduced gyrokinetic model for the study of the spatial structure of ITG instability linear modes PhD work directed by N. Besse with invaluable input from P. Bertrand, E. Gravier, P. Morel,R.
More informationInverse source estimation problems in magnetostatics
Inverse source estimation problems in magnetostatics ERNSI Workshop, Lyon, 2017 Juliette Leblond INRIA Sophia Antipolis (France), Team APICS Joint work with: Laurent Baratchart, Sylvain Chevillard, Doug
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationH ( E) E ( H) = H B t
Chapter 5 Energy and Momentum The equations established so far describe the behavior of electric and magnetic fields. They are a direct consequence of Maxwell s equations and the properties of matter.
More informationSmoothed projections in finite element exterior calculus
Smooted projections in finite element exterior calculus Ragnar Winter CMA, University of Oslo Norway based on joint work wit: Douglas N. Arnold, Minnesota, Ricard S. Falk, Rutgers, and Snorre H. Cristiansen,
More informationThe Interaction of Light and Matter: α and n
The Interaction of Light and Matter: α and n The interaction of light and matter is what makes life interesting. Everything we see is the result of this interaction. Why is light absorbed or transmitted
More informationBreakdown of the KLN Theorem for Charged Particles in Condensed Matter
Breakdown of the KLN Theorem for Charged Particles in Condensed Matter Y. N. Srivastava, and A. Widom Department of Physics & INFN, University of Perugia, Perugia, Italy Department of Physics, Northeastern
More informationWave propagation in an inhomogeneous plasma
DRAFT Wave propagation in an inhomogeneous plasma Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX NP, UK This version is of 7 February 208. Introduction In
More informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
More informationII Theory Of Surface Plasmon Resonance (SPR)
II Theory Of Surface Plasmon Resonance (SPR) II.1 Maxwell equations and dielectric constant of metals Surface Plasmons Polaritons (SPP) exist at the interface of a dielectric and a metal whose electrons
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationFaraday Transport in a Curved Space-Time
Commun. math. Phys. 38, 103 110 (1974) by Springer-Verlag 1974 Faraday Transport in a Curved Space-Time J. Madore Laboratoire de Physique Theorique, Institut Henri Poincare, Paris, France Received October
More informationA Variational Analysis of a Gauged Nonlinear Schrödinger Equation
A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationCoercivity of high-frequency scattering problems
Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 9
WiSe 202 20.2.202 Prof. Dr. A-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg
More informationNonlinear wave-wave interactions involving gravitational waves
Nonlinear wave-wave interactions involving gravitational waves ANDREAS KÄLLBERG Department of Physics, Umeå University, Umeå, Sweden Thessaloniki, 30/8-5/9 2004 p. 1/38 Outline Orthonormal frames. Thessaloniki,
More informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More informationDarwin and higher order approximations to Maxwell s equations in R 3. Sebastian Bauer Universität Duisburg-Essen
Darwin and higher order approximations to Maxwell s equations in R 3 Sebastian Bauer Universität Duisburg-Essen in close collaboration with the Maxwell group around Dirk Pauly Universität Duisburg-Essen
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 informationIntegrated Full Wave Analysis of RF Heating and Current Drive in Toroidal Plasmas
Integrated Full Wave Analysis of RF Heating and Current Drive in Toroidal Plasmas IAEA Fusion Energy Conference Chengdu, China 2006/10/20 A. Fukuyama, S. Murakami, A. Sonoda, M. Honda// Department of Nuclear
More informationRelativistic description of electron Bernstein waves
PSFC/JA-6-7 Relativistic description of electron Bernstein waves J. Decker and A.K. Ram October 6 Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge MA 39 USA This work was
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationElectrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Electro Dynamic
Electrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Name Electro Dynamic Instructions: Use SI units. Short answers! No derivations here, just state your responses clearly. 1. (2) Write an
More informationCoupling of eddy-current and circuit problems
Coupling of eddy-current and circuit problems Ana Alonso Rodríguez*, ALBERTO VALLI*, Rafael Vázquez Hernández** * Department of Mathematics, University of Trento ** Department of Applied Mathematics, University
More informationStrong-Field QED and High-Power Lasers
LC2006 16-05-2006 with: O. Schröder (UoP science + computing, Tübingen) B. Liesfeld, K.-U. Amthor, H. Schwörer and A. Wipf (FSU Jena) R. Sauerbrey, FZ Rossendorf Outline Introduction 1. Introduction QED
More information