EDP with strong anisotropy : transport, heat, waves equations
|
|
- Georgiana Dickerson
- 6 years ago
- Views:
Transcription
1 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
2 Main goals Effective models when disparate scales occur ; Derivation of the limit models, well-posedness, balances, symmetries, etc ; Convergence, error estimates ; Applications : transport of charged particles (magnetic confinement), heat equations, Maxwell equations.
3 1. Multi-scale analysis for linear first order PDE { t u ε + a y u ε + 1 ε b y u ε = 0, (t, y) R + R m Hypotheses u ε (0, y) = u in (y), y R m. a L 1 loc (R +; W 1, loc (Rm )), b W 1, loc (Rm ) = smooth flows div y a = 0, div y b = 0 = measure preserving flows a(t, y) + b(y) C(1 + y ) = global flows
4 Vlasov equation with strong magnetic field t f ε + 1 ε (v 1 x1 +v 2 x2 )f ε +v 3 x3 f ε + q m E v f ε + qb mε (v 2 v1 v 1 v2 )f ε = 0 m = 6, y = (x, v), a x,v = v 3 x3 + q m E(x) v 1 ε b x,v = 1 ε (v 1 x1 + v 2 x2 ) + ω c ε (v 2 v1 v 1 v2 )
5 Question: behavior when ε 0? Main idea: filtering out the fast oscillations dy ds = b(y (s; y)), Y (0; y) = y, (s, y) R Rm New coordinates Search for a profile z = Y ( t/ε, y) or equivalently y = Y (t/ε, z) u ε (t, y) = v ε (t, Y ( t/ε; y) ) }{{} z
6 Why this change of coordinates? t v ε (t, z) + y Y ( t/ε; Y (t/ε; z))a(t, Y (t/ε; z)) z v ε (t, z) = 0, ( }{{} ϕ(t/ε)a(t) v ε (0, z) = u in (z), z Stability for (v ε ) ε>0 Behavior when ε 0 ϕ(s)a = y Y ( s; Y (s; ))a(y (s; )) y Y ( t/ε; Y (t/ε; z))a(t, Y (t/ε; z)) = ϕ(t/ε)a(t)
7 If involution between a(t) and b [b, a(t)] = 0 = ϕ(s)a(t) = a(t), s R { t v ε (t, z) + a(t, z) z v ε (t, z) = 0, (t, z) R + R m v ε (0, z) = u in (z), v ε (t, z) = u in (Z( t; z)) = v(t, z), dz dt z R m = a(t, Z(t; z)) u ε (t, y) = v(t, Y ( t/ε; y)) = u in (Z( t; Y ( t/ε; y))) Splitting : advection along a and advection along 1 ε b.
8 Two scale approach : t and s = t/ε Use ergodicity lim T + 1 T Key point T 0 y Y ( s; Y (s; ))a(t, Y (s; )) ds = lim T + Emphasize a C 0 -group of unitary transformations and use : von Neumann s Ergodic Mean Theorem 1 T T 0 ϕ(s)a(t) ds Let (G(s)) s R be a C 0 -group of unitary operators on a Hilbert space (H, (, )) and A be the infinitesimal generator of G. Then, for any x H, we have lim T + 1 T r+t uniformly with respect to r R. r G(s)x ds = Proj ker A x, strongly in H
9 Average vector field X Q = {c(y) : R m R m measurable : Q(y) : c(y) c(y) dy < + } R m Q = P 1, P = t P > 0, [b, P] := (b y )P y b P P t y b = 0 (c, d) Q = Q(y) : c(y) d(y) dy, c, d X Q. R m Proposition (ϕ(s)) s R is a C 0 -group of unitary operators on X Q. Theorem We denote by L the infinitesimal generator of the group (ϕ(s)) s R. Then for any vector field a X Q, we have the strong convergence in X Q a := lim T + 1 r+t y Y ( s; Y (s; ))a(y (s; )) ds = Proj T ker L a r uniformly with respect to r R.
10 { t c L 2 c = 0, t R + Theorem (Long time behavior) For any vector field a X Q, we consider the problem c(0, ) = a( ) with Lc = [b, c]. Then the solution c(t) converges weakly in X Q, as t +, toward the orthogonal projection on ker L lim c(t) = Proj ker La, weakly in X Q. t + Moreover, if the range of L is closed, then the previous convergence holds strongly in X Q and has exponential rate.
11 Theorem (Convergence) The family (v ε ) ε>0 converges strongly in L loc (R +; L 2 (R m )) to a weak solution v L (R + ; L 2 (R m )) of the transport problem { t v + a(t, ) z v = 0, (t, z) R + R m v(0, z) = u in (z), z R m. Moreover, if v is smooth enough, we have v ε = v + O(ε) in L loc (R +; L 2 (R m )), as ε 0. Formal proof t v ε + ϕ(t/ε)a(t) z v ε = 0 v ε (t, z) = v(t, s = t/ε, z) + εv 1 (t, s = t/ε, z) +... s v = 0, t v + ϕ(s)a(t) z v + s v 1 = 0 ( 1 T ) t v + lim ϕ(s)a(t) ds z v = 0. T + T 0
12 Rigorous proof Lemma Let c L (R + ; X Q ), d L 1 (R + ; X P ) such that c := lim T + 1 r+t c(s) ds strongly in X Q, uniformly w.r.t. r R. T r Then we have lim ε 0 d(t), c(t/ε) P,Q dt = R + d(t), c P,Q dt R + Application R + t v ε + ϕ(t/ε)a(t) z v ε = 0 R m v ε (t, z)ϕ(t/ε)a(t) z ξ(t, z) dzdt, ξ C 1 c (R + R m ) c(s) = ϕ(s)a, d(t) = v ε (t, ) z ξ(t, ), lim T + 1 T r+t r c(s) ds = a.
13 Non linear multi-scale problems Vlasov-Poisson equations t f ε + v x f ε + q m ( x φ ε + v ε 0 x φ ε = ρ ε := q ) B(t, x) v f = 0 ε f ε (t, x, v) dv f ε : particle presence density in the phase space f ε (t, x, v)dxdv : particle number inside the volume dxdv E ε = x φ ε : self consistent electric field
14 Main purposes Homogenization, two scale convergence Averaging with respect to the fast cyclotronic motion Effective non linear Vlasov equation Strong convergence results for any initial conditions (not necessarily well prepared) Conservation laws, Hamiltonian structure
15 The finite Larmor radius regime The reference time T is much larger than the cyclotronic period (strong magnetic field) i.e., T q Bε m 1, with 0 < ε << 1. ε The kinetic energy is much larger than the potential energy m V 2 qφ 1 ε The Larmor radius is of the same order as the Debye length i.e., λ 2 D = ε 0φ nq ρ2 L.
16 Finite Larmor radius regime strong uniform magnetic field B ε = (0, 0, B/ε), perpendicular to x 1 Ox 2. x = (x 1, x 2 ), v = (v 1, v 2 ), ω c = qb m is the rescaled cyclotronic frequency T c = 2π ω c is the rescaled cyclotronic period the real cyclotronic frequency is ωc ε = ω c /ε the real cyclotronic period is Tc ε = 2π ωc ε References = ε 2π ω c Frénod, Sonnendrücker, Golse, Han-Kwan,... v = (v 2, v 1 ) = εt c
17 Two dimensional setting t f ε + 1 ε (v xf ε +ω c v v f ε ) x φ ε v f ε = 0, (t, x, v) R + R 2 R 2 x φ ε = ρ ε (t, x) := f ε (t, x, v) dv, (t, x) R + R 2 R 2 f ε (0, x, v) = f in (x, v), (x, v) R 2 R 2 Main goal What is the behavior of (f ε, φ ε ) ε>0 when ε 0? What is the limit Vlasov-Poisson system?
18 The characteristic equations The trajectories in (x, v) oscillate at the cyclotronic frequency dx ε dt = V ε (t), ε dv ε dt = ω c ε V ε (t) x φ ε (t, X ε (t)) But X ε (t) = X ε (t) + V ε (t)/ω c, Ṽ ε (t) = R (ω c t/ε) V ε (t) are left invariant with respect to the cyclotronic dynamics d X ε dt x φ ε (t, X ε (t)) =, ω c dṽ ε dt = R (ω c t/ε) x φ ε (t, X ε (t)) Main idea stability for ( X ε, Ṽ ε ) ε>0 d X dt = V(t, X (t), Ṽ (t)), dṽ dt = A(t, X (t), Ṽ (t))
19 Change of phase space coordinates Presence densities in the phase space ( x, ṽ) f ε (t, x, ṽ) = f ε (t, x, v), x = x + v ω c, ṽ = R (ω c t/ε) v t f ε x φ ε x f ε R (ω c t/ε) x φ ε ṽ f ε = 0, (t, x, ṽ) R + R 2 R 2 ω c f ε (0, x, ṽ) = f in ( x We expect that lim ε 0 f ε = f and ṽ ), ṽ, ( x, ṽ) R 2 R 2. ω c t f + V x f + A ṽ f = 0, (t, x, ṽ) R + R 2 R 2
20 The effective trajectories e(z) = 1 2π ln z, z R2 \ {0} The solution of the Poisson equation φ ε (t, x) = e(x y)ρ ε (t, y) dy = R 2 R 2 R 2 e(x y)f ε (t, y, w) dwdy d X ε x φ ε (t, X ε (t)) = dt ω c = 1 ω c = 1 ω c R 2 e (X ε (t) y) f ε (t, y, w) dwdy R 2 R 2 ( e X ε (t) ỹ 1 ) R ( ω c t/ε) (Ṽ ε (t) w) f ε (t, ỹ, w) ω c R 2
21 Average over a cyclotronic period X ε (t + T ε c ) X ε (t) = 1 ω c T ε c Tc ε t+t ε c t R 2 R 2 e ( f ε (τ, ỹ, w) d wdỹdτ = 1 ( 1 2π e ω c R 2 R 2 2π 0 ξ = ω c R 2 ( X ε (τ) ỹ R ω cτ ε X (Ṽ (t) ỹ R(θ) ) ) (Ṽ ε (τ) w) ω c ) (t) w) dθ f d wdỹ ω c R 2 E( X (t) ỹ, Ṽ (t) w) f (t, ỹ, w) d wdỹ + o(1), ε 0 E(ξ, η) = 1 2π 2π 0 ( ) e ξ ωc 1 R(θ) η dθ, (ξ, η) (R 2 R 2 )\{(0, 0)}
22 Effective velocity field d X dt = V[ f (t)]( X (t), Ṽ (t)) ξ V[ f (t)]( x, ṽ) = E( x ỹ, ṽ w) f (t, ỹ, w) d wdỹ, ( x, ṽ) R 4 ω c R 2 R 2 Effective acceleration field dṽ dt = A[ f (t)]( X (t), Ṽ (t)) A[ f (t)]( x, ṽ) =ω c η R 2 R 2 E( x ỹ, ṽ w) f (t, ỹ, w) d wdỹ, ( x, ṽ) R 4
23 Hamiltonian structure d X dt = V[ f (t)]( X (t), Ṽ (t)), dṽ dt = A[ f (t)]( X (t), Ṽ (t)) V[ f (t)]( x, ṽ) = ωc 1 x φ[ f (t)], A[ f (t)]( x, ṽ) = ωc ṽ φ[ f (t)] φ[ f (t)]( x, ṽ) = E( x ỹ, ṽ w) f (t, ỹ, w) d wdỹ R 2 R 2 This characteristic system is Hamiltonian, with respect to the conjugate variables ( x 2, ωc 1 ṽ 1 ) and (ω c x 1, ṽ 2 ) and the Hamiltonian function φ[ f ]. d dt d dt t ( X 2, ωc 1 Ṽ 1 ) = ωc x 1,ṽ 2 φ[ f (t)]( X (t), Ṽ (t)) t (ω c X1, Ṽ2) = x2,ωc 1 ṽ 1 φ[ f (t)]( X (t), Ṽ (t)).
24 The function E E(ξ, η) is the average of the fundamental solution e( ) over the circle of center ξ and radius η / ω c If ξ > η / ω c, the function z e(z) is harmonic in the open set R 2 \ {0}, which contains the disc {z R 2 : z ξ η / ω c } and thus, the mean property applied to the function e( ) and the circle of center ξ and radius η / ω c yields E(ξ, η) = 1 2π 2π 0 e ( ξ R(θ) ) η ω c dθ = e(ξ) = 1 η ln ξ, ξ > 2π ω c.
25 The function E The function E has also the symmetry property E(ω 1 c η, ω c ξ) = E(ξ, η) for any (ξ, η) (R 2 R 2 ) \ {(0, 0)}. If ξ < η / ω c, then Finally E(ξ, η) = E ( η E(ξ, η) = e ω c ( ) ( ) η η, ω c ξ = e ω c ω c ) 1 { ξ η / ωc } + e(ξ)1 { ξ > η / ωc }.
26 First order derivatives of E ξ E(ξ, η) = e(ξ)1 { ξ > η / ωc } ( ) η η E(ξ, η) = ωc 1 e 1 { ξ η / ωc } Second order derivatives of E ω c ξ E = 1 { ξ = η / ω c }dσ(ξ, η) 2π ξ 1 + ω 2 c, η E = 1 { ξ = η / ω c }dσ(ξ, η) 2π η. 1 + ωc 2
27 The velocity and acceleration fields associated to any density f = f ( x, ṽ) write V[ f ]( x, ṽ) = 1 ( x ỹ) 2πω c R 2 R 2 x ỹ f 2 (ỹ, w)1 ṽ w { x ỹ > ωc } d wdỹ A[ f ]( x, ṽ) = ω c 2π R 2 Non linear limit model R 2 (ṽ w) ṽ w 2 f (ỹ, w)1 { x ỹ ṽ w ωc } d wdỹ. t f + V[ f ]( x, ṽ) x f + A[ f ]( x, ṽ) ṽ f = 0.
28 Convergence results f in = f in (x, v) non negative presence density satisfying H1 H2 H3 R 2 R 2 f in (x, v) dvdx < + R 2 R 2 v 2 2 f in (x, v) dvdx < + there is a bounded, non increasing function F in = F in (r) L L 1 (R + ; rdr), such that f in (x, v) F in ( v ), (x, v) R 2 R 2. Let (f ε, φ ε ) ε>0 be the weak solutions for the Vlasov-Poisson system t f ε + 1 ε (v xf ε +ω c v v f ε ) x φ ε v f ε = 0, (t, x, v) R + R 4 x φ ε = ρ ε (t, x) := f ε (t, x, v) dv, (t, x) R + R 2 R 2 f ε (0, x, v) = f in (x, v), (x, v) R 2 R 2
29 Theorem (2D) Let ( f ε ) ε>0 be the densities ( f ε (t, x, ṽ) = f ε t, x R ( ω ct/ε) ṽ, ) R ( ω c t/ε) ṽ ω c Then (ε k ) k 0 s. t. ( f ε k ) k converges strongly in L 2 ([0, T ]; L 2 (R 2 R 2 )), for any T R +, to f t f +V[ f (t)]( x, ṽ) x f +A[ f (t)]( x, ṽ) ṽ f = 0, (t, x, ṽ) R + R 2 R 2 with the initial condition f (0, x, ṽ) = f in ( x ṽ ), ṽ, ( x, ṽ) R 2 R 2 ω c where the velocity and acceleration vector fields V, A are given by V[ f (t)]( x, ṽ) = ωc 1 x φ[ f (t)], φ= 1 { ṽ w ln 2π ω c R 2 R 2 A[ f (t)]( x, ṽ) = ωc ṽ φ[ f (t)] 1 { x ỹ ṽ w ωc } + ln x ỹ 1 { x ỹ > ṽ w ωc } } f
30 Lemma Let U = U(z, t, s) : O R + R s R be a function in L 1 (O R + ; C # (R s )), where O is an open set of R N and C # (R s ) stands for the set of continuous periodic functions of (fixed) period L > 0. Then we have the convergence lim U(z, t, t/ε) dtdz = 1 ε 0 O R + L O R + L 0 U(z, t, s) dsdtdz.
31 Conservations 1. Let f = f (t, x, ṽ) be a solution of the limit problem such that 1, x, ṽ, x 2, ṽ 2 are integrable functions with respect to f (0, x, ṽ)dṽd x = f in ( x ωc 1 ṽ, ṽ)dṽd x. For any t R + we have d {1, x, ṽ, x 2, ṽ 2 } f (t, x, ṽ) dṽd x = 0 dt R 2 R 2 2. Let f = f (t, x, ṽ) be a solution of the limit problem such that R 2 R 2 R 2 R 2 E( x ỹ, ṽ w) f (0, ỹ, w) f (0, x, ṽ) d wdỹ dṽd x < +. The electric energy is preserved in time d 1 dt 2 2 φ[ f (t)]( x, ṽ) f (t, x, ṽ) dṽd x = R 2 R R 2 R 2 φ[ f (t)] t f dṽd x = 0.
32 Lemma Let f = f (t, x, ṽ) be a solution of the limit model and ψ = ψ( x, ṽ) be a smooth integrable function with respect to f (0, x, ṽ)dṽd x = f in ( x ωc 1 ṽ, ṽ)dṽd x. For any t R + we have 2 d ψ( x, ṽ) f (t, x, ṽ) dṽd x = dt 2 f (t, ỹ, w) f (t, x, ṽ) R 2 R 2 R 2 R 2 R 2 R [ 1 ( ỹ ψ(ỹ, w) x ψ( x, ṽ) ) e( x ỹ) 1 ω { x ỹ > ṽ w / ωc } c + ( ṽ ψ( x, ṽ) w ψ(ỹ, w)) e (ṽ ) ] w 1 { x ỹ < ṽ w / ωc } ω c
33 Averaged Poisson equation ρ s (t, x, ṽ) = f (t, x ωc 1 R( ω c s) (ṽ w), w) d w R 2 We introduce the average charge density ρ(t, x, ṽ) = 1 Tc ρ s (t, x, ṽ) ds, (t, x, ṽ) R + R 2 R 2 T c 0 The definition of the function E implies φ[ f (t)]( x, ṽ) = E( x ỹ, ṽ w) f (t, ỹ, w) d wdỹ R 2 R 2 = 1 Tc e( x ỹ ωc 1 R( ω c s) (ṽ w)) f (t, ỹ, w) d wdỹds T c 0 R 2 R 2 = 1 Tc e( x z) f (t, z ωc 1 R( ω c s) (ṽ w), w)d wd zds T c 0 R 2 R 2 = 1 Tc e( x z) ρ s (t, z, ṽ) d zds = e( x z) ρ(t, z, ṽ) d z =(e ρ(t,, ṽ))( T c 0 R 2 R 2
34 2. Strongly anisotropic parabolic problems t u ε div y (D(y) y u ε ) 1 ε div y (b(y) b(y) y u ε ) = 0, (t, y) R + R m Limit model (ε 0) D = (1) u ε (0, y) = u in (y), y R m (2) t u div y ( D y u) = 0, (t, y) R + R m u(0, ) = u in 1 S lim Y ( s; Y (s; ))D(Y (s; )) t Y ( s; Y (s; )) ds. S + S 0
35 Parabolic equations with stiff transport t u ε div y (D(y) y u ε ) + 1 ε b(y) y u ε = 0, (t, y) R + R m u ε (0, y) = u in (y), y R m. Change of variable y = Y (t/ε; z), v ε (t, z) = u ε (t, Y (t/ε; z)), (t, z) R + R m, ε > 0 div y {D y (w Y ( s, ))} = (div z (G(s)D z w)) Y ( s, ) (G(s)D)(z) = Y ( s; Y (s; z)) D(Y (s; z)) t Y ( s; Y (s; z)) { t v ε div z ((G(t/ε)D) z v ε ) = 0, (t, z) R + R m v ε (0, z) = u ε (0, z) = u in (z), z R m, ε > 0
36 Average matrix field { } H Q = A(y) : Q(y)A(y) : A(y)Q(y) dy < + R m (, ) Q : H Q H Q R, (A, B) Q = Q(y)A(y) : B(y)Q(y) dy. R m G(s) : H Q H Q, G(s)A = y Y 1 (s; )A(Y (s; )) t y Y 1 (s; ), s R Proposition (G(s)) s R is a C 0 -group of unitary operators on H Q. Theorem Let L be the infinitesimal generator of the group (G(s)) s R. For any matrix field A H Q we have the strong convergence in H Q A := lim T + 1 T r+t uniformly with respect to r R. r y Y 1 (s; )A(Y (s; )) t y Y 1 (s; ) ds = Proj ker L A
37 { t C L 2 C = 0, t R + Theorem (Long time behavior) For any matrix field A H Q we consider the problem C(0, ) = A( ) where L(C) = [b, C] = (b y )C y b C C t y b. Then C(t) converges weakly in H Q, as t +, toward the orthogonal projection on ker L lim C(t) = Proj ker LA, weakly in H Q. t + Moreover, if the range of L is closed, then the previous convergence holds strongly in H Q, and has exponential rate.
38 3. Strongly anisotropic wave equation 2 t u ε div y (D(y) y u ε ) 1 ε 2 div y (b(y) b(y) y u ε ) = 0, (t, y) R + R m u ε (0, ) = u in ( ), t u ε (0, ) = u in ( ) The limit model is a wave equation associated to the matrix field D Strong convergences for well prepared initial conditions lim ε 0 R m ( b y u ε ) (0, y) 2 dy = 0 ε
39 Finite speed propagation Spectral analysis of D min λ min λ λ SpQ 1/2 DQ 1/2 λ SpQ 1/2 D Q 1/2 max λ max λ λ SpQ 1/2 D Q 1/2 λ SpQ 1/2 DQ 1/2 The limit solution propagates with speed c such that c 2 max max λ λ SpQ 1/2 D Q 1/2 The limit solution propagates with speed c max λ λ SpQ 1/2 DQ 1/2 1/2
40 The Maxwell equations t D roth = 0, t B + rote = 0, divd = 0, divb = 0 D = ɛ 0 ɛ r E, B = µ 0 H Energy balance 1 d {D E + B H} + div(e H) = 0. 2 dt
41 Strongly anisotropic electric permittivity n i = indice propre du milieux ɛ r = diag(n 2 1, n 2 2, n 2 3) ɛ 1/2 r = M + b b ε Maxwell equations [ t D ε rot Bε = 0, t B ε + rot µ 0 1 ɛ 0 divb ε = 0, divd ε = 0. ( M + b b ) ] 2 D ε = 0 ε
42 The case TM B = (0, 0, B), D = (D 1, D 2, 0) M ε = M + b b, M = λ b b ε The limit wave equation for B ( λ t 2 B c0 2 2 ξ 4 ) div ξ 4 (ξ ξ) B = 0 ξ = b b, = the average operator along the flow of ξ.
43 Perspectives 1. limit models with non uniform magnetic fields (magnetic confinement) 2. two scale convergence with non periodic fast variable 3. replace periodicity by ergodicity 4. gyrokinetic collisional models, averaged collision operators, invariants, equilibria (work in progress) 5. wave equation, Maxwell equations (general case) 6. spectral analysis of the average matrix field.
analysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationTransport equations with disparate advection fields. Application to the gyrokinetic models in plasma physics
Transport equations with disparate advection fields. Application to the gyrokinetic models in plasma physics Mihai Bostan (July 22, 2 Abstract The subject matter of this paper concerns the asymptotic regimes
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 informationTHE THREE-DIMENSIONAL FINITE LARMOR RADIUS APPROXIMATION
THE THREE-DIMENSIONAL FINITE LARMOR RADIUS APPROXIMATION DANIEL HAN-KWAN Abstract. Following Frénod and Sonnendrücker ([12], we consider the finite Larmor radius regime for a plasma submitted to a large
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationWhat place for mathematicians in plasma physics
What place for mathematicians in plasma physics Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 15-19 September 2008 Eric Sonnendrücker (U. Strasbourg) Math
More informationON THE VLASOV-MAXWELL SYSTEM WITH A STRONG EXTERNAL MAGNETIC FIELD
ON THE VLASOV-MAXWELL SYSTEM WITH A STONG EXTENAL MAGNETIC FIELD Francis Filbet, Tao Xiong, Eric Sonnendr ucker To cite this version: Francis Filbet, Tao Xiong, Eric Sonnendr ucker. ON THE VLASOV-MAXWELL
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationCourse no. 4. The Theory of Electromagnetic Field
Cose no. 4 The Theory of Electromagnetic Field Technical University of Cluj-Napoca http://www.et.utcluj.ro/cs_electromagnetics2006_ac.htm http://www.et.utcluj.ro/~lcret March 19-2009 Chapter 3 Magnetostatics
More informationGeometry and Motion Selected answers to Sections A and C Dwight Barkley 2016
MA34 Geometry and Motion Selected answers to Sections A and C Dwight Barkley 26 Example Sheet d n+ = d n cot θ n r θ n r = Θθ n i. 2. 3. 4. Possible answers include: and with opposite orientation: 5..
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 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 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 informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationAutomorphic Equivalence Within Gapped Phases
1 Harvard University May 18, 2011 Automorphic Equivalence Within Gapped Phases Robert Sims University of Arizona based on joint work with Sven Bachmann, Spyridon Michalakis, and Bruno Nachtergaele 2 Outline:
More informationLectures on basic plasma physics: Introduction
Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters
More informationFluid Dynamics from Kinetic Equations
Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu
More informationarxiv: v3 [math.ap] 28 Nov 2010
Quasineutral limit of the Vlasov-Poisson system with massless electrons Daniel Han-Kwan arxiv:4.93v3 [math.ap] 8 Nov Abstract In this paper, we study the quasineutral limit in other words the limit when
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationOn perturbations in the leading coefficient matrix of time-varying index-1 DAEs
On perturbations in the leading coefficient matrix of time-varying index-1 DAEs Institute of Mathematics, Ilmenau University of Technology Darmstadt, March 27, 2012 Institute of Mathematics, Ilmenau University
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 informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationSemiclassical limit of the Schrödinger-Poisson-Landau-Lifshitz-Gilbert system
Semiclassical limit of the Schrödinger-Poisson-Landau-Lifshitz-Gilbert system Lihui Chai Department of Mathematics University of California, Santa Barbara Joint work with Carlos J. García-Cervera, and
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationGeometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas
Geometric Gyrokinetic Theory and its Applications to Large-Scale Simulations of Magnetized Plasmas Hong Qin Princeton Plasma Physics Laboratory, Princeton University CEA-EDF-INRIA School -- Numerical models
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationn 2 xi = x i. x i 2. r r ; i r 2 + V ( r) V ( r) = 0 r > 0. ( 1 1 ) a r n 1 ( 1 2) V( r) = b ln r + c n = 2 b r n 2 + c n 3 ( 1 3)
Sep. 7 The L aplace/ P oisson Equations: Explicit Formulas In this lecture we study the properties of the Laplace equation and the Poisson equation with Dirichlet boundary conditions through explicit representations
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More informationGeneralized transition waves and their properties
Generalized transition waves and their properties Henri Berestycki a and François Hamel b a EHESS, CAMS, 54 Boulevard Raspail, F-75006 Paris, France & University of Chicago, Department of Mathematics 5734
More informationKAM for NLS with harmonic potential
Université de Nantes 3rd Meeting of the GDR Quantum Dynamics MAPMO, Orléans, 2-4 February 2011. (Joint work with Benoît Grébert) Introduction The equation : We consider the nonlinear Schrödinger equation
More informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
More informationRégularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen
Régularité des équations de Hamilton-Jacobi du premier ordre et applications aux jeux à champ moyen Daniela Tonon en collaboration avec P. Cardaliaguet et A. Porretta CEREMADE, Université Paris-Dauphine,
More informationCHARGED PARTICLE MOTION IN CONSTANT AND UNIFORM ELECTROMAGNETIC FIELDS
CHARGED PARTICLE MOTION IN CONSTANT AND UNIFORM ELECTROMAGNETIC FIELDS In this and in the following two chapters we investigate the motion of charged particles in the presence of electric and magnetic
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationRenormalization of microscopic Hamiltonians. Renormalization Group without Field Theory
Renormalization of microscopic Hamiltonians Renormalization Group without Field Theory Alberto Parola Università dell Insubria (Como - Italy) Renormalization Group Universality Only dimensionality and
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
More informationSingular Perturbations of Stochastic Control Problems with Unbounded Fast Variables
Singular Perturbations of Stochastic Control Problems with Unbounded Fast Variables Joao Meireles joint work with Martino Bardi and Guy Barles University of Padua, Italy Workshop "New Perspectives in Optimal
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
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 information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationMAE294B/SIOC203B: Methods in Applied Mechanics Winter Quarter sgls/mae294b Solution IV
MAE9B/SIOC3B: Methods in Applied Mechanics Winter Quarter 8 http://webengucsdedu/ sgls/mae9b 8 Solution IV (i The equation becomes in T Applying standard WKB gives ɛ y TT ɛte T y T + y = φ T Te T φ T +
More informationBlow-up on manifolds with symmetry for the nonlinear Schröding
Blow-up on manifolds with symmetry for the nonlinear Schrödinger equation March, 27 2013 Université de Nice Euclidean L 2 -critical theory Consider the one dimensional equation i t u + u = u 4 u, t > 0,
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 informationNosé-Hoover Thermostats
Nosé- Nosé- Texas A & M, UT Austin, St. Olaf College October 29, 2013 Physical Model Nosé- q Figure: Simple Oscillator p Physical Model Continued Heat Bath T p Nosé- q Figure: Simple Oscillator in a heat
More informationWeak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System
Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly
More informationNew ideas in the non-equilibrium statistical physics and the micro approach to transportation flows
New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,
More informationHyperbolic Geometry on Geometric Surfaces
Mathematics Seminar, 15 September 2010 Outline Introduction Hyperbolic geometry Abstract surfaces The hemisphere model as a geometric surface The Poincaré disk model as a geometric surface Conclusion Introduction
More informationFundamental Concepts of Particle Accelerators III : High-Energy Beam Dynamics (2) Koji TAKATA KEK. Accelerator Course, Sokendai. Second Term, JFY2012
.... Fundamental Concepts of Particle Accelerators III : High-Energy Beam Dynamics (2) Koji TAKATA KEK koji.takata@kek.jp http://research.kek.jp/people/takata/home.html Accelerator Course, Sokendai Second
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationHigh-electric-field limit for the Vlasov-Maxwell-Fokker-Planck system
High-electric-field limit for the Vlasov-Maxwell-Fokker-Planck system Mihai Bostan, Thierry Goudon (November 1, 5 Abstract In this paper we derive the high-electric-field limit of the three dimensional
More informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
More informationStatistical mechanics of random billiard systems
Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore
More informationDr. Allen Back. Dec. 3, 2014
Dr. Allen Back Dec. 3, 2014 forms are sums of wedge products of the basis 1-forms dx, dy, and dz. They are kinds of tensors generalizing ordinary scalar functions and vector fields. They have a skew-symmetry
More informationThe KPP minimal speed within large drift in two dimensions
The KPP minimal speed within large drift in two dimensions Mohammad El Smaily Joint work with Stéphane Kirsch University of British Columbia & Pacific Institute for the Mathematical Sciences Banff, March-2010
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationZ. Zhou On the classical limit of a time-dependent self-consistent field system: analysis. computation
On the classical limit of a time-dependent self-consistent field system: analysis and computation Zhennan Zhou 1 joint work with Prof. Shi Jin and Prof. Christof Sparber. 1 Department of Mathematics Duke
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
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 informationOn the leapfrogging phenomenon in fluid mechanics
On the leapfrogging phenomenon in fluid mechanics Didier Smets Université Pierre et Marie Curie - Paris Based on works with Robert L. Jerrard U. of Toronto) CIRM, Luminy, June 27th 2016 1 / 22 Single vortex
More informationOn a general definition of transition waves and their properties
On a general definition of transition waves and their properties Henri Berestycki a and François Hamel b a EHESS, CAMS, 54 Boulevard Raspail, F-75006 Paris, France b Université Aix-Marseille III, LATP,
More informationPhysics GRE: Electromagnetism. G. J. Loges 1. University of Rochester Dept. of Physics & Astronomy. xkcd.com/567/
Physics GRE: Electromagnetism G. J. Loges University of Rochester Dept. of Physics & stronomy xkcd.com/567/ c Gregory Loges, 206 Contents Electrostatics 2 Magnetostatics 2 3 Method of Images 3 4 Lorentz
More informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More informationOn the Boltzmann equation: global solutions in one spatial dimension
On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators
More informationAn inverse scattering problem for short-range systems in a time-periodic electric field. François Nicoleau
An inverse scattering problem for short-range systems in a time-periodic electric field. François Nicoleau Laboratoire Jean Leray UMR CNRS-UN 6629 Département de Mathématiques 2, rue de la Houssinière
More informationLecture 3: Asymptotic Methods for the Reduced Wave Equation
Lecture 3: Asymptotic Methods for the Reduced Wave Equation Joseph B. Keller 1 The Reduced Wave Equation Let v (t,) satisfy the wave equation v 1 2 v c 2 = 0, (1) (X) t2 where c(x) is the propagation speed
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationLectures notes on Boltzmann s equation
Lectures notes on Boltzmann s equation Simone Calogero 1 Introduction Kinetic theory describes the statistical evolution in phase-space 1 of systems composed by a large number of particles (of order 1
More informationEntropic structure of the Landau equation. Coulomb interaction
with Coulomb interaction Laurent Desvillettes IMJ-PRG, Université Paris Diderot May 15, 2017 Use of the entropy principle for specific equations Spatially Homogeneous Kinetic equations: 1 Fokker-Planck:
More informationHARMONIC ANALYSIS. Date:
HARMONIC ANALYSIS Contents. Introduction 2. Hardy-Littlewood maximal function 3. Approximation by convolution 4. Muckenhaupt weights 4.. Calderón-Zygmund decomposition 5. Fourier transform 6. BMO (bounded
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationThe branch of mechanics which deals with the motion of object is called dynamics. It is divided into two branches:
M th KINEMATICS 4 CHAPTER INTRODUCTION The branch of mechanics which deals with the motion of object is called dynamics. It is divided into two branches: Kinematics: (i) Kinematics (ii) Kinetics The branch
More informationLarge deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationIsometric elastic deformations
Isometric elastic deformations Fares Al-Azemi and Ovidiu Calin Abstract. This paper deals with the problem of finding a class of isometric deformations of simple and closed curves, which decrease the total
More informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 La Tremblade, Congrès SMAI 2017 5
More informationPARABOLIC LIMIT AND STABILITY OF THE VLASOV FOKKER PLANCK SYSTEM
Mathematical Models and Methods in Applied Sciences Vol. 10, No. 7 (2000 1027 1045 c World Scientific Publishing Company PARABOLIC LIMIT AND STABILITY OF THE VLASOV FOKKER PLANCK SYSTEM F. POUPAUD Laboratoire
More informationFYS 3120: Classical Mechanics and Electrodynamics
FYS 3120: Classical Mechanics and Electrodynamics Formula Collection Spring semester 2014 1 Analytical Mechanics The Lagrangian L = L(q, q, t), (1) is a function of the generalized coordinates q = {q i
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationWELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2)
WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS We will use the familiar Hilbert spaces H = L 2 (Ω) and V = H 1 (Ω). We consider the Cauchy problem (.1) c u = ( 2 t c )u = f L 2 ((, T ) Ω) on [, T ] Ω u() = u H
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More informationKPP Pulsating Traveling Fronts within Large Drift
KPP Pulsating Traveling Fronts within Large Drift Mohammad El Smaily Joint work with Stéphane Kirsch University of British olumbia & Pacific Institute for the Mathematical Sciences September 17, 2009 PIMS
More informationNew phenomena for the null controllability of parabolic systems: Minim
New phenomena for the null controllability of parabolic systems F.Ammar Khodja, M. González-Burgos & L. de Teresa Aix-Marseille Université, CNRS, Centrale Marseille, l2m, UMR 7373, Marseille, France assia.benabdallah@univ-amu.fr
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationOn the cosmic no-hair conjecture in the Einstein-Vlasov setting
On the cosmic no-hair conjecture in the Einstein-Vlasov setting KTH, Royal Institute of Technology, Stockholm Recent Advances in Mathematical General Relativity Institut Henri Poincaré, Paris September
More informationExercises - Chapter 1 - Chapter 2 (Correction)
Université de Nice Sophia-Antipolis Master MathMods - Finite Elements - 28/29 Exercises - Chapter 1 - Chapter 2 Correction) Exercise 1. a) Let I =], l[, l R. Show that Cl) >, u C Ī) Cl) u H 1 I), u DĪ).
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More informationThe Gaussian free field, Gibbs measures and NLS on planar domains
The Gaussian free field, Gibbs measures and on planar domains N. Burq, joint with L. Thomann (Nantes) and N. Tzvetkov (Cergy) Université Paris Sud, Laboratoire de Mathématiques d Orsay, CNRS UMR 8628 LAGA,
More information