Homogenization of micro-resonances and localization of waves.
|
|
- Maximilian Marsh
- 5 years ago
- Views:
Transcription
1 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 (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
2 Outline: - Classical homogenization for waves ( = low frequencies) - Higher frequencies + higher contrasts ( degeneracies ) resonant homogenization - Effects: (frequency/wavenumber) band gaps, dispersion, negative materials, etc. - Partial degeneracies and resonances (more of effects; general theory) - Photonic Crystal Fibers as an example of partial degeneracies Band gaps in PCFs. (Cooper, Kamotski, V.S., 2012). Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
3 Classical Periodic Homogenisation for waves (= low frequencies) ( anti-plane elastodynamics/ TM electrodynamics, for simplicity) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
4 Classical Periodic Homogenisation for waves (= low frequencies) ( anti-plane elastodynamics/ TM electrodynamics, for simplicity) ρ ε (x)u tt div (a ε (x) u ε )) = f (x, t) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
5 Classical Periodic Homogenisation for waves (= low frequencies) ( anti-plane elastodynamics/ TM electrodynamics, for simplicity) ρ ε (x)u tt div (a ε (x) u ε )) = f (x, t) f (x, t) 0, t 0; u(x, t) 0, t 0. a ε (x) = a(x/ε), ρ ε = ρ(x/ε) a stiffness, ρ density a(y), ρ(y) Q-periodic in y Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
6 Classical Periodic Homogenisation for waves (= low frequencies) ( anti-plane elastodynamics/ TM electrodynamics, for simplicity) ρ ε (x)u tt div (a ε (x) u ε )) = f (x, t) f (x, t) 0, t 0; u(x, t) 0, t 0. a ε (x) = a(x/ε), ρ ε = ρ(x/ε) a stiffness, ρ density a(y), ρ(y) Q-periodic in y Asymptotic expansion u ε (x, t) u 0 (x, x/ε, t) + εu 1 (x, x/ε, t) +... u 0 (x, y, t), u 1 (x, y, t) Q-periodic in y Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
7 Classical Periodic Homogenisation ε 2 : u 0 = u(x, t) ( long wavelengths low frequencies ) alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
8 Classical Periodic Homogenisation ε 2 : u 0 = u(x, t) ( long wavelengths low frequencies ) ε 1 : u 1 = j N j (y) u0 x j (x, t) N solutions of cell problem div y ( a(y)(e j + y N j ) ) = 0 Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
9 Classical Periodic Homogenisation ε 2 : u 0 = u(x, t) ( long wavelengths low frequencies ) ε 1 : u 1 = j N j (y) u0 x j (x, t) ( N solutions of cell problem div y a(y)(e j + y N j ) ) = 0 ( ) ε 0 : ρ hom utt 0 div A hom u 0 ) = f (x), homogenized eqn ρ hom = ρ(y) y A hom = a(y) (I + y N) y homogenised (averaged) density homogenised stiffness tensor Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
10 Classical Periodic Homogenisation ε 2 : u 0 = u(x, t) ( long wavelengths low frequencies ) ε 1 : u 1 = j N j (y) u0 x j (x, t) ( N solutions of cell problem div y a(y)(e j + y N j ) ) = 0 ( ) ε 0 : ρ hom utt 0 div A hom u 0 ) = f (x), homogenized eqn ρ hom = ρ(y) y A hom = a(y) (I + y N) y homogenised (averaged) density homogenised stiffness tensor 0 < νi A hom, ρ hom ν 1 I (uniform positivity inherited) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
11 Classical Periodic Homogenisation ε 2 : u 0 = u(x, t) ( long wavelengths low frequencies ) ε 1 : u 1 = j N j (y) u0 x j (x, t) ( N solutions of cell problem div y a(y)(e j + y N j ) ) = 0 ( ) ε 0 : ρ hom utt 0 div A hom u 0 ) = f (x), homogenized eqn ρ hom = ρ(y) y A hom = a(y) (I + y N) y homogenised (averaged) density homogenised stiffness tensor 0 < νi A hom, ρ hom ν 1 I (uniform positivity inherited) The Homogenisation Theorem C > 0 independent of ε such that u ε (u 0 + εu 1 ) H 1 Cε 1/2 Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
12 High-contrast (= micro-resonant ) homogenization and non-classical two-scale limits (Zhikov 2000, 2004) A ε u = div (a ε (x) u ε ) { a ε ε 2 on Ω (x) = ε 0 ( soft phase) 1 on Ω ε 1 ( stiff phase) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
13 High-contrast (= micro-resonant ) homogenization and non-classical two-scale limits (Zhikov 2000, 2004) A ε u = div (a ε (x) u ε ) { a ε ε 2 on Ω (x) = ε 0 ( soft phase) 1 on Ω ε 1 ( stiff phase) Contrast δ ε 2 is a critical scaling giving rise to non-classical effects (Khruslov 1980s; Arbogast, Douglas, Hornung 1990; Panasenko 1991; Allaire 1992; Sandrakov 1999; Brianne 2002; Bourget, Mikelic, Piatnitski 2003; Bouchitte & Felbaq 2004,...): elliptic, spectral, parabolic, hyperbolic, nonlinear, non-periodic/ random,.... Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
14 High-contrast (= micro-resonant ) homogenization and non-classical two-scale limits (Zhikov 2000, 2004) A ε u = div (a ε (x) u ε ) { a ε ε 2 on Ω (x) = ε 0 ( soft phase) 1 on Ω ε 1 ( stiff phase) Contrast δ ε 2 is a critical scaling giving rise to non-classical effects (Khruslov 1980s; Arbogast, Douglas, Hornung 1990; Panasenko 1991; Allaire 1992; Sandrakov 1999; Brianne 2002; Bourget, Mikelic, Piatnitski 2003; Bouchitte & Felbaq 2004,...): elliptic, spectral, parabolic, hyperbolic, nonlinear, non-periodic/ random,.... WHY? Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
15 Two-scale formal asymptotic expansion: div (a ε (x) u ε ) + ρω 2 u ε = 0 (timeharmonicwaves) A ε u ε = λu ε, λ = ρω 2 (spectral problem). Seek u ε (x) u 0 (x, x/ε) + εu 1 (x, x/ε) +... u j (x, y) Q periodic in y. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
16 Two-scale formal asymptotic expansion: div (a ε (x) u ε ) + ρω 2 u ε = 0 (timeharmonicwaves) A ε u ε = λu ε, λ = ρω 2 (spectral problem). Seek u ε (x) u 0 (x, x/ε) + εu 1 (x, x/ε) +... u j (x, y) Q periodic in y. THEN: Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
17 Two-scale limit problem (Zhikov 2000, 2004) Then u 0 (x, y) = u 0 (x) in Q 1 (still low frequency) w(x,y) in Q 0 ( resonance frequency) alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
18 Two-scale limit problem (Zhikov 2000, 2004) Then u 0 (x, y) = u 0 (x) in Q 1 (still low frequency) w(x,y) in Q 0 ( resonance frequency) (u 0, w), w(x, y) := u 0 (x) + v(x, y), solves the two-scale limit spectral problem: div x (a hom x u(x)) = λu(x) + λ < v > (x) in Ω div y (a (0) y v(x, y)) = λ(u(x) + v(x, y)) in Q 0 v(x, y) = 0 on Q 0 alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
19 Two-scale limit problem (Zhikov 2000, 2004) Then u 0 (x, y) = u 0 (x) in Q 1 (still low frequency) w(x,y) in Q 0 ( resonance frequency) (u 0, w), w(x, y) := u 0 (x) + v(x, y), solves the two-scale limit spectral problem: div x (a hom x u(x)) = λu(x) + λ < v > (x) in Ω div y (a (0) y v(x, y)) = λ(u(x) + v(x, y)) in Q 0 v(x, y) = 0 on Q 0 alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
20 Two-scale limit problem (Zhikov 2000, 2004) Then u 0 (x, y) = u 0 (x) in Q 1 (still low frequency) w(x,y) in Q 0 ( resonance frequency) (u 0, w), w(x, y) := u 0 (x) + v(x, y), solves the two-scale limit spectral problem: div x (a hom x u(x)) = λu(x) + λ < v > (x) in Ω div y (a (0) y v(x, y)) = λ(u(x) + v(x, y)) in Q 0 v(x, y) = 0 on Q 0 Decouple it alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
21 Two-scale limit spectral problem Decouple by choosing v(x, y) = λu(x)b(y) div y (a (0) y b(y)) λu = 1 in Q 0 b(y) = 0 on Q 0 div x (a hom u(x)) = β(λ)u(x), in Ω, where β(λ) = λ + λ 2 φ j 2 y λ j λ, (λ j, φ j ) Dirichlet eigenvalues/functions of inclusion Q 0 ( = micro-resonances : β < 0 negative denstity/magnetism etc) j=1 Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
22 Rigorous analysis: Two-scale Convergence Definition 1. Let u ε (x) be a bounded sequence in L 2 (Ω). We say (u ε ) weakly two-scale converges to u 0 (x, y) L 2 2 (Ω Q), denoted by u ε u 0, if for all φ C0 (Ω), ψ C # (Q) ( x ) u ε (x)φ(x)ψ dx u 0 (x, y)φ(x)ψ(y) dxdy Ω ε Ω Q as ε We say (u ε ) strongly two-scale converges to u 0 L 2 (Ω Q), denoted 2 2 by u ε u 0, if for all v ε v 0 (x, y), u ε (x)v ε (x)dx u 0 (x, y)v 0 (x, y) dxdy Ω as ε 0. (implies convergence of norms upon sufficient regularity) Ω Q Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves.july 13, / 23
23 Rigorous analysis (Zhikov 2000, 2004): 1. Two-scale resolvent convergence: α > 0, A ε u ε + αu ε = f ε L 2 (Ω); u ε H0 1(Ω). If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y). (If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y).) Here u 0 solves two-scale limit problem A 0 u 0 + αu 0 = f 0. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
24 Rigorous analysis (Zhikov 2000, 2004): 1. Two-scale resolvent convergence: α > 0, A ε u ε + αu ε = f ε L 2 (Ω); u ε H0 1(Ω). If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y). (If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y).) Here u 0 solves two-scale limit problem A 0 u 0 + αu 0 = f Spectral band gaps: (Let Ω = R d.) A 0 self-adjoint in H L 2 (R Q 0 ), with a band-gap spectrum σ(a 0 ). σ(a ε ) σ(a 0 ) in the sense of Hausdorff. (Hence a Band-gap effect: For small enough ε waves of certain frequencies do not propagate, in any direction). The proof is based on a two-scale spectral compactness. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
25 Rigorous analysis (Zhikov 2000, 2004): 1. Two-scale resolvent convergence: α > 0, A ε u ε + αu ε = f ε L 2 (Ω); u ε H0 1(Ω). If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y). (If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y).) Here u 0 solves two-scale limit problem A 0 u 0 + αu 0 = f Spectral band gaps: (Let Ω = R d.) A 0 self-adjoint in H L 2 (R Q 0 ), with a band-gap spectrum σ(a 0 ). σ(a ε ) σ(a 0 ) in the sense of Hausdorff. (Hence a Band-gap effect: For small enough ε waves of certain frequencies do not propagate, in any direction). The proof is based on a two-scale spectral compactness. Limit band gaps ={λ : β(λ) < 0}. Interpretation : β(λ) = µ(ω) < 0 negative magnetism/density (Bouchitte & Felbacq, 2004). Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
26 Rigorous analysis (Zhikov 2000, 2004): 1. Two-scale resolvent convergence: α > 0, A ε u ε + αu ε = f ε L 2 (Ω); u ε H0 1(Ω). If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y). (If f ε 2 f 0 (x, y) then u ε 2 u 0 (x, y).) Here u 0 solves two-scale limit problem A 0 u 0 + αu 0 = f Spectral band gaps: (Let Ω = R d.) A 0 self-adjoint in H L 2 (R Q 0 ), with a band-gap spectrum σ(a 0 ). σ(a ε ) σ(a 0 ) in the sense of Hausdorff. (Hence a Band-gap effect: For small enough ε waves of certain frequencies do not propagate, in any direction). The proof is based on a two-scale spectral compactness. Limit band gaps ={λ : β(λ) < 0}. Interpretation : β(λ) = µ(ω) < 0 negative magnetism/density (Bouchitte & Felbacq, 2004). β(λ) > 0 - propagation with (high) dispersion (λ λ j ), due to coupled resonances (V.S. & P. Kuchment, 2007). Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
27 Frequency vs directional gaps and partial degeneracies Cherednichenko, V.S., Zhikov (2006): spatial nonlocality for homogenised limit with highly anisotropic fibers. 1 in Q 1 (matrix) a ε (x) = ε 2 in Q 0 across fibers 1 in Q 0 along fibers alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
28 Frequency vs directional gaps and partial degeneracies Cherednichenko, V.S., Zhikov (2006): spatial nonlocality for homogenised limit with highly anisotropic fibers. 1 in Q 1 (matrix) a ε (x) = ε 2 in Q 0 across fibers 1 in Q 0 along fibers V.S. (2009): a directional localization (via formal asymptotic expansions): for certain frequencies waves can propagate in some directions (e.g. along the fibers above) but cannot in others (e.g. orthogonal to the fibers). Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
29 Frequency vs directional gaps and partial degeneracies Cherednichenko, V.S., Zhikov (2006): spatial nonlocality for homogenised limit with highly anisotropic fibers. 1 in Q 1 (matrix) a ε (x) = ε 2 in Q 0 across fibers 1 in Q 0 along fibers V.S. (2009): a directional localization (via formal asymptotic expansions): for certain frequencies waves can propagate in some directions (e.g. along the fibers above) but cannot in others (e.g. orthogonal to the fibers). Notice, in the above fibres, a ε (x) = a (1) (x/ε) + ε 2 a (0) (x/ε), where a (1) = 0 0 0, i.e. is partially degenerate Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
30 Photonic crystal fibers: A partially degenerate problem (Cooper, I. Kamotski, V.S. 2012); cf scalar prototype problem I.Kamotski V.S. 2006; M. Cherdantsev 2009 ( full contrast) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
31 Photonic crystal fibers: Problem Formulation E = iωµh, H = iωɛe ɛ = ɛ 0 χ 0 (x/ε) + ɛ 1 χ 1 (x/ε) ɛ 0 > ɛ 1, µ constant E=exp(ikx 3 )E(x 1, x 2 ), H=exp(ikx 3 )H(x 1, x 2 ) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
32 Photonic crystal fibers: Problem Formulation E = iωµh, H = iωɛe ɛ = ɛ 0 χ 0 (x/ε) + ɛ 1 χ 1 (x/ε) ɛ 0 > ɛ 1, µ constant E=exp(ikx 3 )E(x 1, x 2 ), H=exp(ikx 3 )H(x 1, x 2 ) In each phase E 3 and H 3 satisfy the following equations E 3 + a ε E 3 = 0, H 3 + a ε H 3 = 0 where a ε = ω 2 µɛ(x/ε) - k 2. E 3 and H 3 coupled across interface Γ ε : [ ɛ ] [ ] [ ] 1 1 [ µ ] ω a ε E 3 n = k a ε H 3 n, k a ε E 3 n = ω a ε H 3 n Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
33 Photonic crystal fibers: A partially degenerate problem Consider the problem in its weak form: ( ωɛ ) ( ωɛ ) ( ) ( ) k k 1 a ε E 3,1 + 2 a ε E 3,2 + 1 a ε H 3,2 2 a ε H 3,1 = ωɛe 3 ( ) ( ) k k ( ωµ ) ( ωµ ) 1 a ε E 3,2 2 a ε E 3,1 1 a ε H 3,1 2 a ε H 3,2 = ωµh 3. alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
34 Photonic crystal fibers: A partially degenerate problem Consider the problem in its weak form: ( ωɛ ) ( ωɛ ) ( ) ( ) k k 1 a ε E 3,1 + 2 a ε E 3,2 + 1 a ε H 3,2 2 a ε H 3,1 = ωɛe 3 ( ) ( ) k k ( ωµ ) ( ωµ ) 1 a ε E 3,2 2 a ε E 3,1 1 a ε H 3,1 2 a ε H 3,2 = ωµh 3. For u = (E 3, H 3 ), find u such that ω ( ) k ({ } { }) ɛ u1 R 2 a ε φ 1 + µ u 2 φ 2 + φ1 a ε, u 2 + u1, φ 2 dx = ωρ(x/ε)u φ φ C0 (R 2 ) R 2 ( ) ɛ(y) 0 {u, v} := u,1 v,2 v,1 u,2 (Poisson bracket); ρ(y) = 0 µ Form positive if k 2 < ω 2 µɛ 1. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
35 Photonic crystal fibers: A partially degenerate problem Consider the problem in its weak form: ( ωɛ ) ( ωɛ ) ( ) ( ) k k 1 a ε E 3,1 + 2 a ε E 3,2 + 1 a ε H 3,2 2 a ε H 3,1 = ωɛe 3 ( ) ( ) k k ( ωµ ) ( ωµ ) 1 a ε E 3,2 2 a ε E 3,1 1 a ε H 3,1 2 a ε H 3,2 = ωµh 3. For u = (E 3, H 3 ), find u such that ω ( ) k ({ } { }) ɛ u1 R 2 a ε φ 1 + µ u 2 φ 2 + φ1 a ε, u 2 + u1, φ 2 dx = ωρ(x/ε)u φ φ C0 (R 2 ) R 2 ( ) ɛ(y) 0 {u, v} := u,1 v,2 v,1 u,2 (Poisson bracket); ρ(y) = 0 µ Form positive if k 2 < ω 2 µɛ 1. Consider a near critical k: k 2 = ω 2 µ(ɛ 1 ε 2 ), ω 2 µ = λ Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
36 Photonic crystal fibers: a partially degenerate problem An emergent high contrast: A ε (x) u φ dx = λ R 2 ρ(x/ε)u φdx R 2 φ C0 (R 2 ) Here A ε (x) = A(x/ε) where A(y) = a (1) (y) + ε 2 a (0) (y) + O(ε 4 ). a (1) 0 BUT a (1) (y) + a (0) (y) > νi, ν > 0 a (1) (y) u u = χ 1 (y) ( u 1,1 + u 2,2 2 + u 1,2 u 2,1 2) (partially) DEGENERATES for u s.t. RHS zero. i.e. satisfies Cauchy-Riemann type equations in matrix phase. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
37 Photonic crystal fibers: a partially degenerate problem An emergent high contrast: A ε (x) u φ dx = λ R 2 ρ(x/ε)u φdx R 2 φ C0 (R 2 ) Here A ε (x) = A(x/ε) where A(y) = a (1) (y) + ε 2 a (0) (y) + O(ε 4 ). a (1) 0 BUT a (1) (y) + a (0) (y) > νi, ν > 0 a (1) (y) u u = χ 1 (y) ( u 1,1 + u 2,2 2 + u 1,2 u 2,1 2) (partially) DEGENERATES for u s.t. RHS zero. i.e. satisfies Cauchy-Riemann type equations in matrix phase. Moral: This appears a particular case of homogenization for partially degenerating PDE systems Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
38 General Partial Degeneracies (I. Kamotski and V.S. 2012) Consider a resolvent problem: Ω R d, α > 0, - div (a ε (x) u ε ) + αρ ε u ε = f ε L 2 (Ω), u ε ( H 1 0 (Ω)) n, n 1. A general degeneracy a ε (x) = a (1) ( ) x ε + ε 2 a (0) ( ) x ε, ( ) n d n d a (l) L # (Q), a (1) 0, a (1) + a (0) > 0, Weak formulation: Ω [ ( a (1) x ) ( u φ(x) + ε 2 a (0) x ) ] u φ(x) + α ρ ε (x) u φ(x) dx ε ε = f ε (x) φ(x) dx, φ ( H0 1 (Ω) ) d. Ω alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
39 General Partial Degeneracies (I. Kamotski and V.S. 2012) Consider a resolvent problem: Ω R d, α > 0, - div (a ε (x) u ε ) + αρ ε u ε = f ε L 2 (Ω), u ε ( H 1 0 (Ω)) n, n 1. A general degeneracy a ε (x) = a (1) ( ) x ε + ε 2 a (0) ( ) x ε, ( ) n d n d a (l) L # (Q), a (1) 0, a (1) + a (0) > 0, Weak formulation: Ω [ ( a (1) x ) ( u φ(x) + ε 2 a (0) x ) ] u φ(x) + α ρ ε (x) u φ(x) dx ε ε = f ε (x) φ(x) dx, φ ( H0 1 (Ω) ) d. A priori estimates: Ω alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
40 Weak two-scale limits. Key assumption on the degeneracy Introduce V := { v ( ) } n H# 1 (Q) a (1) (y) y v = 0. (subspace { of microscopic oscillations ), and ) n d ( W:= ψ (L 2 (a # (Q) divy (1) (y) ) ) 1/2 ψ(y) = 0 in ( microscopic fluxes ) Then, up to a subsequence, u ε 2 u 0 (x, y) L 2 (Ω; V ) ε u ε 2 y u 0 (x, y) ξ ε (x) := ( a (1) (x/ε) ) 1/2 u ε 2 ξ0 (x, y) L 2 (Ω; W ). ( H 1 # (Q) ) n } Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
41 Weak two-scale limits. Key assumption on the degeneracy Introduce V := { v ( ) } n H# 1 (Q) a (1) (y) y v = 0. (subspace { of microscopic oscillations ), and ) n d ( W:= ψ (L 2 (a # (Q) divy (1) (y) ) ) 1/2 ψ(y) = 0 in ( microscopic fluxes ) Then, up to a subsequence, u ε 2 u 0 (x, y) L 2 (Ω; V ) ε u ε 2 y u 0 (x, y) ξ ε (x) := ( a (1) (x/ε) ) 1/2 u ε 2 ξ0 (x, y) L 2 (Ω; W ). ( H 1 # (Q) ) n } Key assumption ) n There exists a constant C > 0 such that for all v (H# 1 (Q) P V v (H 1 # (Q)) n C a (1) (y) y v L 2 Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
42 The two-scale Limit Operator Let Ω be bounded Lipschitz, or Ω = R d. Introduce U L 2 (Ω; V ): { U := u(x, y) L 2 (Ω; V ) ξ(x, y) L2 (Ω; W ) s.t., Ψ(x, y) C ( ( ) ) 1/2 ξ(x, y) Ψ(x, y)dxdy = u(x, y) x a (1) (y) Ψ(x, y) Ω Q Define T : U L 2 by Tu = ξ. Strong two-scale resolvent convergence: Let f ε 2 f 0 (x, y). Then u ε 2 u 0 (x, y) solving: Find u 0 U such that φ U Ω { Tu 0 (x, y) T φ 0 (x, y) + a (0) (y) y u 0 (x, y) y φ 0 (x, y) + Ω Q } + α ρ(y)u 0 (x, y) φ 0 (x, y) dy dx = f 0 (x, y) φ 0 (x, y) dy dx. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23 Q Ω Q
43 Back to Photonic Crystals { V θ := v ( } Hθ 1(Q)) n a(1) (y) y v = 0. v V θ iff v is θ-quasi-periodic, and v 1,1 + v 2,2 = 0, v 1,2 v 2,1 = 0 in Q 1. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
44 Back to Photonic Crystals { V θ := v ( } Hθ 1(Q)) n a(1) (y) y v = 0. v V θ iff v is θ-quasi-periodic, and v 1,1 + v 2,2 = 0, v 1,2 v 2,1 = 0 in Q 1. Theorem 1 (Key assumption holds) There exists a constant c > 0 such that for any u H 1 θ (Q) P V θ u H 1 (Q) c ( u 1,1 + u 2,2 L 2 (Q 0 ) + u 1,2 u 2,1 L 2 (Q 0 )). Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
45 Back to Photonic Crystals { V θ := v ( } Hθ 1(Q)) n a(1) (y) y v = 0. v V θ iff v is θ-quasi-periodic, and v 1,1 + v 2,2 = 0, v 1,2 v 2,1 = 0 in Q 1. Theorem 1 (Key assumption holds) There exists a constant c > 0 such that for any u H 1 θ (Q) P V θ u H 1 (Q) c ( u 1,1 + u 2,2 L 2 (Q 0 ) + u 1,2 u 2,1 L 2 (Q 0 )). Second Result The generalised flux ξ =: Tu is zero. i.e. Tu = 0 u U. (Two-scale limit solution u 0 determined by miscroscopic behaviour only) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
46 Limit Spectral problem: Find u V θ such that a (0) (y) y u(y) y (φ(y)) dy = λ Q Q ρ(y)u(y) φ(y) dy φ V θ. alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
47 Limit Spectral problem: Find u V θ such that a (0) (y) y u(y) y (φ(y)) dy = λ Q Q ρ(y)u(y) φ(y) dy φ V θ. Theorem 2 (Spectral compactness: not trivial) Let λ ε σ(a ε ) (Bloch spectrum). Let u ε be associated normalized Bloch s waves: A ε u ε = λ ε u ε, u ε (x) = e iθε x/ε v ε (x/ε), v ε (y) L 2 (Q) = 1, θ ε ( π, π] d. Let λ λ 0 and θ ε θ 0. Then λ 0 σ 0 (A 0, θ 0 ) (spectrum of the Limit operator), and u ε u 0 (y), eigenfunction of A 0 (θ). Hence the spectra converge. Implication: If the limit problem displays a band-gap, the original problem must also have a gap for small enough ε. alery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
48 Example: Band gaps in ARROW fibres An extreme problem: Q 0 a circle of small radius δ Theorem: There exist constants c 1, c 2 > 0 independent of δ, such that, for all quasimomenta θ 0, λ 2 (θ 0 ) c 1 δ 2 ln δ, λ 3(θ 0 ) c 2 δ 2. This implies that, for small enough δ, there is a wide spectral gap in the limit spectrum, and therefore also for small enough ε for the original problem by the spectral compactness. + Higher gaps: λ λ D j (Q δ 0 ) δ 2 (re micro-resonances ) Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
49 The band gaps in a Photonic Crystal Fiber: Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
50 Summary: - Homogenization for a critical high contrast scaling δ ε 2 gives rise to numerous non-classical effects described by two-scale limit problems due to micro-resonances. - Partial degeneracies often happen in physical problems, and give rise to more of such effects (e.g. Band Gaps in Photonic Crystal Fibers). These however have to be analysed in a new way. - A general two-scale homogenization theory can be constructed for such partial degeneracies, under a generically held decomposition condition. Resulting limit (homogenized) operator is generically two-scale (and non-local ). - Associated two-scale operator and spectral convergence and compactness are held generally or for particular physical examples. - In principle, some of this applies also to nonlinear (cf Cherednichenko & Cherdantsev 2011), discrete-to-continuous, as well as non-periodic problems. Valery Smyshlyaev (University College London, Homogenization UK) of micro-resonances and localization of waves. July 13, / 23
Durham Research Online
Durham Research Online Deposited in DRO: 23 April 2018 Version of attached le: Accepted Version Peer-review status of attached le: Peer-reviewed Citation for published item: Cooper, S. (2018 'Quasi-periodic
More informationResolvent estimates for high-contrast elliptic problems with periodic coefficients
Resolvent estimates for high-contrast elliptic problems with periodic coefficients Joint work with Shane Cooper (University of Bath) 25 August 2015, Centro de Ciencias de Benasque Pedro Pascual Partial
More informationResolvent estimates for high-contrast elliptic problems with periodic coefficients
Resolvent estimates for high-contrast elliptic problems with periodic coefficients K. D. Cherednichenko and S. Cooper July 16, 2015 Abstract We study the asymptotic behaviour of the resolvents A ε + I
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 informationApplications of the periodic unfolding method to multi-scale problems
Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56
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 informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
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 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 informationQuantitative Homogenization of Elliptic Operators with Periodic Coefficients
Quantitative Homogenization of Elliptic Operators with Periodic Coefficients Zhongwei Shen Abstract. These lecture notes introduce the quantitative homogenization theory for elliptic partial differential
More informationPeriodic homogenization and effective mass theorems for the Schrödinger equation
Periodic homogenization and effective mass theorems for the Schrödinger equation Grégoire Allaire September 5, 2006 Abstract The goal of this course is to give an introduction to periodic homogenization
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 informationConvergence of Adaptive Finite Element Methods for Elliptic Eigenvalue Problems with Applications to Photonic Crystals
Convergence of Adaptive Finite Element Methods for Elliptic Eigenvalue Problems with Applications to Photonic Crystals submitted by Stefano Giani for the degree of Doctor of Philosophy of the University
More informationCurriculum Vitae. Personal Details. Current Position. Previous Posts. Education. Research Interests. Grants. S. Cooper p. 1
Personal Details Name: Address: Curriculum Vitae Shane Cooper 8 Bamburgh Road, Newton Hall, Durham, DH1 5NW, UK. salcoops@gmail.com http://salcooper.co.uk Email: Webpage: Date of Birth: 28 th December
More informationEnhanced resolution in structured media
Enhanced resolution in structured media Eric Bonnetier w/ H. Ammari (Ecole Polytechnique), and Yves Capdeboscq (Oxford) Outline : 1. An experiment of super resolution 2. Small volume asymptotics 3. Periodicity
More informationON THE SPECTRUM OF THE DIRICHLET LAPLACIAN IN A NARROW STRIP
ON THE SPECTRUM OF THE DIRICHLET LAPLACIAN IN A NARROW STRIP LEONID FRIEDLANDER AND MICHAEL SOLOMYAK Abstract. We consider the Dirichlet Laplacian in a family of bounded domains { a < x < b, 0 < y < h(x)}.
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
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 informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
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 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 informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationInverse Scattering with Partial data on Asymptotically Hyperbolic Manifolds
Inverse Scattering with Partial data on Asymptotically Hyperbolic Manifolds Raphael Hora UFSC rhora@mtm.ufsc.br 29/04/2014 Raphael Hora (UFSC) Inverse Scattering with Partial data on AH Manifolds 29/04/2014
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationThis is an Open Access document downloaded from ORCA, Cardiff University's institutional repository:
This is an Open Access document downloaded from ORCA, Cardiff University's institutional repository: http://orca.cf.ac.uk/93/ This is the author s version of a work that was submitted to / accepted for
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 informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationTOPICS IN NONLINEAR ANALYSIS AND APPLICATIONS. Dipartimento di Matematica e Applicazioni Università di Milano Bicocca March 15-16, 2017
TOPICS IN NONLINEAR ANALYSIS AND APPLICATIONS Dipartimento di Matematica e Applicazioni Università di Milano Bicocca March 15-16, 2017 Abstracts of the talks Spectral stability under removal of small capacity
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 informationTWO-SCALE CONVERGENCE ON PERIODIC SURFACES AND APPLICATIONS
TWO-SCALE CONVERGENCE ON PERIODIC SURFACES AND APPLICATIONS Grégoire ALLAIRE Commissariat à l Energie Atomique DRN/DMT/SERMA, C.E. Saclay 91191 Gif sur Yvette, France Laboratoire d Analyse Numérique, Université
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationBielefeld Course on Nonlinear Waves - June 29, Department of Mathematics University of North Carolina, Chapel Hill. Solitons on Manifolds
Joint work (on various projects) with Pierre Albin (UIUC), Hans Christianson (UNC), Jason Metcalfe (UNC), Michael Taylor (UNC), Laurent Thomann (Nantes) Department of Mathematics University of North Carolina,
More informationAsymptotic Analysis of the Approximate Control for Parabolic Equations with Periodic Interface
Asymptotic Analysis of the Approximate Control for Parabolic Equations with Periodic Interface Patrizia Donato Université de Rouen International Workshop on Calculus of Variations and its Applications
More informationOn the well-posedness of the Prandtl boundary layer equation
On the well-posedness of the Prandtl boundary layer equation Vlad Vicol Department of Mathematics, The University of Chicago Incompressible Fluids, Turbulence and Mixing In honor of Peter Constantin s
More informationHOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM
Homogenization of a heat transfer problem 1 G. Allaire HOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM Work partially supported by CEA Grégoire ALLAIRE, Ecole Polytechnique
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 Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationPeriodic Homogenization of Elliptic Problems (Draft)
Periodic Homogenization of Elliptic Problems (Draft) Zhongwei Shen 1 Department of Mathematics University of Kentucky 1 Supported in part by NSF grant DMS-0855294 2 Contents 1 Elliptic Operators with Periodic
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 informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationanalysis 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 informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationSOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD. 1. Dirichlet integral
SOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD CRISTIAN E. GUTIÉRREZ FEBRUARY 3, 29. Dirichlet integral Let f C( ) with open and bounded. Let H = {u C ( ) : u = f on } and D(u) = Du(x) 2 dx.
More informationLecture 2: A Strange Term Coming From Nowhere
Lecture 2: A Strange Term Coming From Nowhere Christophe Prange February 9, 2016 In this lecture, we consider the Poisson equation with homogeneous Dirichlet boundary conditions { u ε = f, x ε, u ε = 0,
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationOn some weighted fractional porous media equations
On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationSolutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness
Solutions in the sense of distributions Definition, non uniqueness 1. Notion of distributions In order to build weak solutions to the Hopf equation, we need to define derivatives of non smooth functions,
More informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
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 informationRadiation Integrals and Auxiliary Potential Functions
Radiation Integrals and Auxiliary Potential Functions Ranga Rodrigo June 23, 2010 Lecture notes are fully based on Balanis [?]. Some diagrams and text are directly from the books. Contents 1 The Vector
More informationLong time existence of space periodic water waves
Long time existence of space periodic water waves Massimiliano Berti SISSA, Trieste IperPV2017, Pavia, 7 Settembre 2017 XVII Italian Meeting on Hyperbolic Equations Water Waves equations for a fluid in
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationReconstruction Scheme for Active Thermography
Reconstruction Scheme for Active Thermography Gen Nakamura gnaka@math.sci.hokudai.ac.jp Department of Mathematics, Hokkaido University, Japan Newton Institute, Cambridge, Sept. 20, 2011 Contents.1.. Important
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 informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationAsymptotic behavior of the degenerate p Laplacian equation on bounded domains
Asymptotic behavior of the degenerate p Laplacian equation on bounded domains Diana Stan Instituto de Ciencias Matematicas (CSIC), Madrid, Spain UAM, September 19, 2011 Diana Stan (ICMAT & UAM) Nonlinear
More informationTWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES
TWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES JACK ARBUNICH 1. Setting of Periodic Structures Our aim is to study an application of Bloch-Floquet Theory in the multiscale
More informationConstruction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem
Proceedings of 4th Workshop Group Analysis of Differential Equations & Integrability 2009, 1 23 Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem George W. BLUMAN,
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH College of Informatics and Electronics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS425 SEMESTER: Autumn 25/6 MODULE TITLE: Applied Analysis DURATION OF EXAMINATION:
More informationLecture No 2 Degenerate Diffusion Free boundary problems
Lecture No 2 Degenerate Diffusion Free boundary problems Columbia University IAS summer program June, 2009 Outline We will discuss non-linear parabolic equations of slow diffusion. Our model is the porous
More informationLocal smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping
Local smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping H. Christianson partly joint work with J. Wunsch (Northwestern) Department of Mathematics University of North
More informationGradient estimates for eigenfunctions on compact Riemannian manifolds with boundary
Gradient estimates for eigenfunctions on compact Riemannian manifolds with boundary Xiangjin Xu Department of athematics Johns Hopkins University Baltimore, D 21218 Abstract The purpose of this paper is
More informationA note on W 1,p estimates for quasilinear parabolic equations
200-Luminy conference on Quasilinear Elliptic and Parabolic Equations and Systems, Electronic Journal of Differential Equations, Conference 08, 2002, pp 2 3. http://ejde.math.swt.edu or http://ejde.math.unt.edu
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationElliptic Kirchhoff equations
Elliptic Kirchhoff equations David ARCOYA Universidad de Granada Sevilla, 8-IX-2015 Workshop on Recent Advances in PDEs: Analysis, Numerics and Control In honor of Enrique Fernández-Cara for his 60th birthday
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationarxiv: v3 [math.ap] 13 Nov 2009
Sub-Wavelength Plasmonic Crystals: Dispersion Relations and Effective Properties By Santiago P. Fortes, Robert P. Lipton and Stephen P. Shipman Louisiana State University, Baton Rouge, LA, USA arxiv:0910.0889v3
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationExistence of (Generalized) Breathers in Periodic Media
Existence of (Generalized) Breathers in Periodic Media Guido Schneider Lehrstuhl für Analysis und Modellierung www.iadm.uni-stuttgart.de/lstanamod/schneider/ Collaborators:. Martina Chirilus-Bruckner,
More informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationFree boundaries in fractional filtration equations
Free boundaries in fractional filtration equations Fernando Quirós Universidad Autónoma de Madrid Joint work with Arturo de Pablo, Ana Rodríguez and Juan Luis Vázquez International Conference on Free Boundary
More informationA Recursion Formula for the Construction of Local Conservation Laws of Differential Equations
A Recursion Formula for the Construction of Local Conservation Laws of Differential Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada December
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
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 informationOn the Standard Linear Viscoelastic model
On the Standard Linear Viscoelastic model M. Pellicer (Universitat de Girona) Work in collaboration with: J. Solà-Morales (Universitat Politècnica de Catalunya) (bounded problem) B. Said-Houari (Alhosn
More informationWaves in Honeycomb Structures
Waves in Honeycomb Structures Michael I. Weinstein Columbia University Nonlinear Schrödinger Equations: Theory and Applications Heraklion, Crete / Greece May 20-24, 2013 Joint work with C.L. Fefferman
More informationarxiv:submit/ [math-ph] 15 Jul 2010
Convergent ower Series for Fields in ositive or Negative High-Contrast eriodic Media Santiago. Fortes, Robert. Lipton, and Stephen. Shipman Department of Mathematics, Louisiana State University arxiv:submit/0076540
More informationNON-EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION SYSTEM WITH NONLOCAL SOURCES
Electronic Journal of Differential Equations, Vol. 2016 (2016, No. 45, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu NON-EXTINCTION OF
More informationRemarks on decay of small solutions to systems of Klein-Gordon equations with dissipative nonlinearities
Remarks on decay of small solutions to systems of Klein-Gordon equations with dissipative nonlinearities Donghyun Kim (joint work with H. Sunagawa) Department of Mathematics, Graduate School of Science
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationThe Brézis-Nirenberg Result for the Fractional Elliptic Problem with Singular Potential
arxiv:1705.08387v1 [math.ap] 23 May 2017 The Brézis-Nirenberg Result for the Fractional Elliptic Problem with Singular Potential Lingyu Jin, Lang Li and Shaomei Fang Department of Mathematics, South China
More informationQuasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths
in the presence of plasma sheaths Department of Mathematical Sciences Ulsan National Institute of Science and Technology (UNIST) joint work with Masahiro Suzuki (Nagoya) and Chang-Yeol Jung (Ulsan) The
More informationRegularity of Weak Solution to Parabolic Fractional p-laplacian
Regularity of Weak Solution to Parabolic Fractional p-laplacian Lan Tang at BCAM Seminar July 18th, 2012 Table of contents 1 1. Introduction 1.1. Background 1.2. Some Classical Results for Local Case 2
More informationStrichartz Estimates for the Schrödinger Equation in Exterior Domains
Strichartz Estimates for the Schrödinger Equation in University of New Mexico May 14, 2010 Joint work with: Hart Smith (University of Washington) Christopher Sogge (Johns Hopkins University) The Schrödinger
More informationMultiscale Computations for Flow and Transport in Porous Media
Multiscale Computations for Flow and Transport in Porous Media Thomas Y. Hou Applied Mathematics, 217-5, Caltech, Pasadena, CA 91125, USA. Email: hou@ama.caltech.edu. Abstract. Many problems of fundamental
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
More informationQuasi-conformal minimal Lagrangian diffeomorphisms of the
Quasi-conformal minimal Lagrangian diffeomorphisms of the hyperbolic plane (joint work with J.M. Schlenker) January 21, 2010 Quasi-symmetric homeomorphism of a circle A homeomorphism φ : S 1 S 1 is quasi-symmetric
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More information