Electromagnetic wave propagation in complex dispersive media
|
|
- Joan Anthony
- 5 years ago
- Views:
Transcription
1 Electromagnetic wave propagation in complex dispersive media Associate Professor Department of Mathematics Workshop on Quantification of Uncertainties in Material Science January 15, 2016 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
2 A Polynomial Chaos Method for Microscale Modeling Associate Professor Department of Mathematics Workshop on Quantification of Uncertainties in Material Science January 15, 2016 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
3 Outline 1 Preliminaries (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
4 Outline 1 Preliminaries 2 Maxwell-Debye (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
5 Outline 1 Preliminaries 2 Maxwell-Debye 3 Maxwell-Random Debye (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
6 Outline 1 Preliminaries 2 Maxwell-Debye 3 Maxwell-Random Debye 4 Conclusions (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
7 Acknowledgments Collaborators H. T. Banks (NCSU) V. A. Bokil (OSU) W. P. Winfree (NASA) Students Karen Barrese and Neel Chugh (REU 2008) Anne Marie Milne and Danielle Wedde (REU 2009) Erin Bela and Erik Hortsch (REU 2010) Megan Armentrout (MS 2011) Brian McKenzie (MS 2011) Duncan McGregor (MS 2013, PhD 2016?) (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
8 Preliminaries Preliminaries Havriliak-Negami (H-N) A general (phenomenological) material model that arises from considerations of the multi-scale nature of the spatial microstructure of a broad class of materials (e.g., glassy, soils, biological tissues, and amorphous polymers). a a See A complex plane representation of dielectric and mechanical relaxation processes in some polymers, J. Polym. Sci, Corresponds to a fractional psuedo-differential equation model in the time domain, not suitable for efficient simulation. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
9 Maxwell-Debye Maxwell s Equations Maxwell s Equations B t D t + E = 0, in (0, T ) D (Faraday) + J H = 0, in (0, T ) D (Ampere) D = B = 0, in (0, T ) D E(0, x) = E 0 ; H(0, x) = H 0, in D E n = 0, on (0, T ) D (Poisson/Gauss) (Initial) (Boundary) E = Electric field vector D = Electric flux density H = Magnetic field vector B = Magnetic flux density J = Current density n = Unit outward normal to Ω (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
10 Maxwell-Debye Maxwell s Equations Constitutive Laws Maxwell s equations are completed by constitutive laws that describe the response of the medium to the electromagnetic field. D = ɛe + P B = µh + M J = σe + J s P = M = J s = Polarization Magnetization Source Current ɛ = Electric permittivity µ = Magnetic permeability σ = Electric Conductivity where ɛ = ɛ 0 ɛ and µ = µ 0 µ r. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
11 Maxwell-Debye Dispersive Media Complex permittivity We can usually define P in terms of a convolution P(t, x) = g E(t, x) = t 0 g(t s, x; q)e(s, x)ds, where g is the dielectric response function (DRF). In the frequency domain ˆD = ɛê + ĝê = ɛ 0 ɛ(ω)ê, where ɛ(ω) is called the complex permittivity. ɛ(ω) described by the polarization model (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
12 Maxwell-Debye Dry skin data 10 3 True Data Debye Model Cole Cole Model ε f (Hz) Figure : Real part of ɛ(ω), ɛ, or the permittivity [GLG96]. Note: up to 10% spread in measurements was observed, but only averages were published. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
13 Maxwell-Debye Dry skin data σ True Data Debye Model Cole Cole Model f (Hz) Figure : Imaginary part of ɛ(ω)/ω, σ, or the conductivity. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
14 Maxwell-Debye Dispersive Media Polarization Models P(t, x) = g E(t, x) = t Debye model [1929] q = [ɛ, ɛ d, τ] or 0 g(t s, x; q)e(s, x)ds, g(t, x) = ɛ 0 ɛ d /τ e t/τ τṗ + P = ɛ 0ɛ d E or ɛ(ω) = ɛ + ɛ d 1 + iωτ with ɛ d := ɛ s ɛ and τ a relaxation time. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
15 Maxwell-Debye Dispersive Media Polarization Models P(t, x) = g E(t, x) = t Debye model [1929] q = [ɛ, ɛ d, τ] or 0 g(t s, x; q)e(s, x)ds, g(t, x) = ɛ 0 ɛ d /τ e t/τ τṗ + P = ɛ 0ɛ d E or ɛ(ω) = ɛ + ɛ d 1 + iωτ with ɛ d := ɛ s ɛ and τ a relaxation time. Cole-Cole model [1936] (heuristic generalization) q = [ɛ, ɛ d, τ, α] ɛ d ɛ(ω) = ɛ (iωτ) α (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
16 Maxwell-Debye Dispersive Media Polarization Models P(t, x) = g E(t, x) = t Debye model [1929] q = [ɛ, ɛ d, τ] 0 ɛ(ω) = ɛ + Cole-Cole model [1941] q = [ɛ, ɛ d, τ, α] ɛ(ω) = ɛ + g(t s, x; q)e(s, x)ds, ɛ d 1 + iωτ ɛ d 1 + (iωτ) α Havriliak-Negami model [1967] q = [ɛ, ɛ d, τ, α, β] ɛ(ω) = ɛ + ɛ d (1 + (iωτ) α ) β (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
17 Maxwell-Random Debye Distribution of Parameters Distributions of Relaxation Times The macroscopic Debye polarization model can be derived from microscopic dipole formulations by passing to a limit over the molecular population [see, Elliot1993]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
18 Maxwell-Random Debye Distribution of Parameters Distributions of Relaxation Times The macroscopic Debye polarization model can be derived from microscopic dipole formulations by passing to a limit over the molecular population [see, Elliot1993]. In 1907, von Schweidler observed the need for multiple relaxation times. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
19 Maxwell-Random Debye Distribution of Parameters Distributions of Relaxation Times The macroscopic Debye polarization model can be derived from microscopic dipole formulations by passing to a limit over the molecular population [see, Elliot1993]. In 1907, von Schweidler observed the need for multiple relaxation times. In 1913, Wagner proposed a (continuous) distribution of relaxation times. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
20 Maxwell-Random Debye Distribution of Parameters Distributions of Relaxation Times The macroscopic Debye polarization model can be derived from microscopic dipole formulations by passing to a limit over the molecular population [see, Elliot1993]. In 1907, von Schweidler observed the need for multiple relaxation times. In 1913, Wagner proposed a (continuous) distribution of relaxation times. Empirical measurements suggest a log-normal or Beta distribution [Bottcher-Bordewijk1978]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
21 Maxwell-Random Debye Distribution of Parameters Distributions of Relaxation Times The macroscopic Debye polarization model can be derived from microscopic dipole formulations by passing to a limit over the molecular population [see, Elliot1993]. In 1907, von Schweidler observed the need for multiple relaxation times. In 1913, Wagner proposed a (continuous) distribution of relaxation times. Empirical measurements suggest a log-normal or Beta distribution [Bottcher-Bordewijk1978]. One can show that the S-N (and Cole-Cole) model corresponds to a continuous distribution... it is possible to calculate the necessary distribution function by the method of Fuoss and Kirkwood. [Cole-Cole1941]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
22 Maxwell-Random Debye Distribution of Parameters Distributions of Relaxation Times The macroscopic Debye polarization model can be derived from microscopic dipole formulations by passing to a limit over the molecular population [see, Elliot1993]. In 1907, von Schweidler observed the need for multiple relaxation times. In 1913, Wagner proposed a (continuous) distribution of relaxation times. Empirical measurements suggest a log-normal or Beta distribution [Bottcher-Bordewijk1978]. One can show that the S-N (and Cole-Cole) model corresponds to a continuous distribution... it is possible to calculate the necessary distribution function by the method of Fuoss and Kirkwood. [Cole-Cole1941]. Continuous spectrum relaxation functions are also common in viscoelastic models. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
23 Maxwell-Random Debye Fit to dry skin data with uniform distribution 10 3 True Data Debye (27.79) Cole Cole (10.4) Uniform (13.60) ε f (Hz) Figure : Real part of ɛ(ω), ɛ, or the permittivity [REU2008]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
24 Maxwell-Random Debye Fit to dry skin data with uniform distribution True Data Debye (27.79) Cole Cole (10.4) Uniform (13.60) 10 1 σ f (Hz) Figure : Imaginary part of ɛ(ω)/ω, σ, or the conductivity [REU2008]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
25 Maxwell-Random Debye Distribution of Parameters Distributions of Parameters To account for the effect of possible multiple parameter sets q, consider the following polydispersive DRF h(t, x; F ) = g(t, x; q)df (q), where Q is some admissible set and F P(Q). Then the polarization becomes: P(t, x; F ) = t 0 Q h(t s, x; F )E(s, x)ds. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
26 Maxwell-Random Debye Random Polarization Alternatively we can define the random polarization P(t, x; τ) to be the solution to τṗ + P = ɛ 0 ɛ d E where τ is a random variable with PDF f (τ), for example, f (τ) = 1 τ b τ a for a uniform distribution. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
27 Maxwell-Random Debye Random Polarization Alternatively we can define the random polarization P(t, x; τ) to be the solution to τṗ + P = ɛ 0 ɛ d E where τ is a random variable with PDF f (τ), for example, f (τ) = 1 τ b τ a for a uniform distribution. The electric field depends on the macroscopic polarization, which we take to be the expected value of the random polarization at each point (t, x) P(t, x; F ) = τb τ a P(t, x; τ)f (τ)dτ. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
28 Maxwell-Random Debye Polynomial Chaos Polynomial Chaos Apply Polynomial Chaos (PC) method to approximate each spatial component of the random polarization τṗ + P = ɛ 0 ɛ d E, τ = τ(ξ) = τ r ξ + τ m, ξ F resulting in or (τ r M + τ m I ) α + α = ɛ 0 ɛ d Eê 1 A α + α = f. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
29 Maxwell-Random Debye Polynomial Chaos Polynomial Chaos Apply Polynomial Chaos (PC) method to approximate each spatial component of the random polarization resulting in or τṗ + P = ɛ 0 ɛ d E, τ = τ(ξ) = τ r ξ + τ m, ξ F (τ r M + τ m I ) α + α = ɛ 0 ɛ d Eê 1 A α + α = f. The electric field depends on the macroscopic polarization, the expected value of the random polarization at each point (t, x), which is P(t, x; F ) = E[P] α 0 (t, x). (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
30 Maxwell-Random Debye Polynomial Chaos 10 2 Maximum Difference Calculated for different values of p and r r = 1.00τ r = 0.75τ r = 0.50τ r = 0.25τ Maximum Error p (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
31 Maxwell-Random Debye Inverse Problem Numerical Results Comparison, noise = 0.1, refinement = 1, perturb = 0.8 Data Initial J= Optimal J= Actual J= E t x 10 9 Comparison of simulations to data [Armentrout-G., 2011]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
32 Maxwell-Random Debye Inverse Problem Numerical Results 6 x 1011 Distributions, noise = 0.1, refinement = 1, perturb = 0.8 Initial J= Optimal J= Actual J= Comparison of initial to final distribution [Armentrout-G., 2011]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
33 Maxwell-Random Debye Inverse Problem Numerical Results Distributions, noise = 0.1, refinement = 1,test = 5 3 Initial J= Optimal J= Actual J= (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
34 Maxwell-Random Debye Maxwell-Random Debye system In a polydispersive Debye material, we have with µ 0 H t = E, (1a) ɛ 0 ɛ E t = H P t J τ P t + P = ɛ 0ɛ d E P(t, x; F ) = τb τ a P(t, x; τ)df (τ). (1b) (1c) (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
35 Maxwell-Random Debye Random Debye Dispersion Relation Theorem (G., 2015) The dispersion relation for the system (1) is given by ω 2 ɛ(ω) = k 2 c2 where the expected complex permittivity is given by [ ] 1 ɛ(ω) = ɛ + ɛ d E. 1 + iωτ Where k is the wave vector and c = 1/ µ 0 ɛ 0 is the speed of light. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
36 Maxwell-Random Debye Random Debye Dispersion Relation Theorem (G., 2015) The dispersion relation for the system (1) is given by ω 2 ɛ(ω) = k 2 c2 where the expected complex permittivity is given by [ ] 1 ɛ(ω) = ɛ + ɛ d E. 1 + iωτ Where k is the wave vector and c = 1/ µ 0 ɛ 0 is the speed of light. Note: for a uniform distribution on [τ a, τ b ], this has an analytic form since [ ] 1 1 E = [arctan(ωτ) + i iωτ ω(τ b τ a ) ln ( 1 + (ωτ) 2)] τ=τ a τ=τ b. The exact dispersion relation can be compared with a discrete dispersion relation to determine the amount of dispersion error. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
37 Maxwell-Random Debye FDTD Finite Difference Time Domain (FDTD) We now choose a discretization of the Maxwell-PC Debye model. Note that any scheme can be used independent of the spectral approach in random space employed here. The Yee Scheme (FDTD) This gives an explicit second order accurate scheme in time and space. It is conditionally stable with the CFL condition ν := c t h 1 d where ν is called the Courant number and c = 1/ µ 0 ɛ 0 ɛ is the fastest wave speed and d is the spatial dimension, and h is the (uniform) spatial step. The Yee scheme can exhibit numerical dispersion and dissipation. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
38 Maxwell-Random Debye FDTD Yee Scheme for Maxwell-Debye System (in 1D) become H n+1 H n j+ µ 1 j t E n+ 1 2 j E n 1 2 j ɛ 0 ɛ t 1 τ Pn+ 2 H µ 0 t = E z E ɛ 0 ɛ t = H z P t τ P t = ɛ 0ɛ d E P = E n+ 1 2 j+1 E n+ 1 2 j z = H n H n j+ 1 j z j P n 1 2 j E n+ 1 2 j + E n 1 2 j = ɛ 0 ɛ d t 2 1 Pn+ 2 j P n 1 2 j t Pn+ 1 2 j + P n 1 2 j. 2 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
39 Maxwell-Random Debye FDTD Discrete Debye Dispersion Relation (Petropolous1994) showed that for the Yee scheme applied to the Maxwell-Debye, the discrete dispersion relation can be written ω 2 c 2 ɛ (ω) = K 2 where the discrete complex permittivity is given by ( ) 1 ɛ (ω) = ɛ + ɛ d 1 + iω τ with discrete (mis-)representations of ω and τ given by ω = sin (ω t/2), τ = sec(ω t/2)τ. t/2 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
40 Maxwell-Random Debye FDTD Discrete Debye Dispersion Relation (cont.) The quantity K is given by K = sin (k z/2) z/2 in 1D and is related to the symbol of the discrete first order spatial difference operator by ik = F(D 1, z ). In this way, we see that the left hand side of the discrete dispersion relation ω 2 c 2 ɛ (ω) = K 2 is unchanged when one moves to higher order spatial derivative approximations [Bokil-G,2012] or even higher spatial dimension [Bokil-G,2013]. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
41 Maxwell-Random Debye PC-Debye FDTD Dispersion Relation Theorem (G., 2015) The discrete dispersion relation for the Maxwell-PC Debye FDTD scheme is given by ω 2 c 2 ɛ (ω) = K 2 where the discrete expected complex permittivity is given by ɛ (ω) := ɛ + ɛ d ê T 1 (I + iω A ) 1 ê 1 and the discrete PC matrix is given by A := sec(ω t/2)a. The definitions of the parameters ω and K are the same as before. Recall the exact complex permittivity is given by [ ] 1 ɛ(ω) = ɛ + ɛ d E 1 + iωτ. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
42 Maxwell-Random Debye PC-Debye FDTD Dispersion Analysis Dispersion Error We define the phase error Φ for a scheme applied to a model to be Φ = k EX k k EX, (2) where the numerical wave number k is implicitly determined by the corresponding dispersion relation and k EX is the exact wave number for the given model. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
43 Maxwell-Random Debye PC-Debye FDTD Dispersion Analysis Dispersion Error We define the phase error Φ for a scheme applied to a model to be Φ = k EX k k EX, (2) where the numerical wave number k is implicitly determined by the corresponding dispersion relation and k EX is the exact wave number for the given model. We wish to examine the phase error as a function of ω t in the range [0, π]. t is determined by h τ τ m, while x = y determined by CFL condition. We note that ω t = 2π/N ppp, where N ppp is the number of points per period, and is related to the number of points per wavelength as, N ppw = ɛ νn ppp. We assume a uniform distribution and the following parameters which are appropriate constants for modeling aqueous Debye type materials: ɛ = 1, ɛ s = 78.2, τ m = sec, τ r = 0.5τ m. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
44 Maxwell-Random Debye PC-Debye FDTD Dispersion Analysis 0.4 PC Debye dispersion for FD with h τ =0.01, r=0.5τ, θ=0 0.4 PC Debye dispersion for FD with h τ =0.01, r=0.9τ, θ= Φ M=0 M=1 M=2 M= ω t Φ M=0 M=1 0.1 M=2 M=3 M= M=5 M= ω t Figure : Plots of phase error at θ = 0 for (left column) τ r = 0.5τ m, (right column) τ r = 0.9τ m, using h τ = (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
45 Maxwell-Random Debye PC-Debye FDTD Dispersion Analysis 0.5 PC Debye dispersion for FD with h τ =0.001, r=0.5τ, θ=0 0.5 PC Debye dispersion for FD with h τ =0.001, r=0.9τ, θ=0 Φ M=0 M= M=2 M= ω t Φ M=0 M=1 M=2 M=3 M=4 M=5 M= ω t Figure : Plots of phase error at θ = 0 for (left column) τ r = 0.5τ m, (right column) τ r = 0.9τ m, using h τ = (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
46 Maxwell-Random Debye PC-Debye FDTD Dispersion Analysis PC Debye dispersion for FD with h τ =0.01, r=0.5τ, ωτ µ =1 PC Debye dispersion for FD with h τ =0.01, r=0.9τ, ωτ µ = M=0 M=1 M=2 M= M=0 M=1 M=2 M=3 M=4 M=5 M= Figure : Log plots of phase error versus θ with fixed ω = 1/τ m for (left column) τ r = 0.5τ m, (right column) τ r = 0.9τ m, using h τ = Legend indicates degree M of the PC expansion. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
47 Maxwell-Random Debye PC-Debye FDTD Dispersion Analysis PC Debye dispersion for FD with h τ =0.001, r=0.5τ, ωτ µ =1 PC Debye dispersion for FD with h τ =0.001, r=0.9τ, ωτ µ = M=0 M=1 M=2 M= M=0 M=1 M=2 M=3 M=4 M=5 M= Figure : Log plots of phase error versus θ with fixed ω = 1/τ m for (left column) τ r = 0.5τ m, (right column) τ r = 0.9τ m, using h τ = Legend indicates degree M of the PC expansion. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
48 Conclusions Conclusions We have presented a random ODE model for polydispersive Debye media 1 Gibson, N. L., A Polynomial Chaos Method for Dispersive Electromagnetics, Comm. in Comp. Phys., 2015 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
49 Conclusions Conclusions We have presented a random ODE model for polydispersive Debye media We described an efficient numerical method utilizing polynomial chaos (PC) and finite difference time domain (FDTD) 1 Gibson, N. L., A Polynomial Chaos Method for Dispersive Electromagnetics, Comm. in Comp. Phys., 2015 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
50 Conclusions Conclusions We have presented a random ODE model for polydispersive Debye media We described an efficient numerical method utilizing polynomial chaos (PC) and finite difference time domain (FDTD) Exponential convergence in the number of PC terms was demonstrated 1 Gibson, N. L., A Polynomial Chaos Method for Dispersive Electromagnetics, Comm. in Comp. Phys., 2015 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
51 Conclusions Conclusions We have presented a random ODE model for polydispersive Debye media We described an efficient numerical method utilizing polynomial chaos (PC) and finite difference time domain (FDTD) Exponential convergence in the number of PC terms was demonstrated We have proven (conditional) stability of the scheme via energy decay (not shown) 1 1 Gibson, N. L., A Polynomial Chaos Method for Dispersive Electromagnetics, Comm. in Comp. Phys., 2015 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
52 Conclusions Conclusions We have presented a random ODE model for polydispersive Debye media We described an efficient numerical method utilizing polynomial chaos (PC) and finite difference time domain (FDTD) Exponential convergence in the number of PC terms was demonstrated We have proven (conditional) stability of the scheme via energy decay (not shown) 1 We have derived a discrete dispersion relation and computed phase errors 1 Gibson, N. L., A Polynomial Chaos Method for Dispersive Electromagnetics, Comm. in Comp. Phys., 2015 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
53 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
54 Appendix Polynomial Chaos Polynomial Chaos: Simple example Consider the first order, constant coefficient, linear IVP with ẏ + ky = g, y(0) = y 0 k = k(ξ) = ξ, ξ N (0, 1), g(t) = 0. We can represent the solution y as a Polynomial Chaos (PC) expansion in terms of (normalized) orthogonal Hermite polynomials H j : y(t, ξ) = α j (t)φ j (ξ), j=0 φ j (ξ) = H j (ξ). Substituting into the ODE we get j=0 α j (t)φ j (ξ) + α j (t)ξφ j (ξ) = 0. j=0 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
55 Appendix Polynomial Chaos Triple recursion formula α j (t)φ j (ξ) + α j (t)ξφ j (ξ) = 0. j=0 j=0 We can eliminate the explicit dependence on ξ by using the triple recursion formula for Hermite polynomials ξh j = jh j 1 + H j+1. Thus j=0 α j (t)φ j (ξ) + α j (t)(jφ j 1 (ξ) + φ j+1 (ξ)) = 0. j=0 (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
56 Appendix Polynomial Chaos Galerkin Projection onto span({φ i } p i=0 ) In order to approximate y we wish to find a finite system for at least the first few α i. We take the weighted inner product with the ith basis, i = 0,..., p, j=0 α j (t) φ j, φ i W + α j (t)(j φ j 1, φ i W + φ j+1, φ i W ) = 0, where f (ξ), g(ξ) W := f (ξ)g(ξ)w (ξ)dξ. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
57 Appendix Polynomial Chaos Galerkin Projection onto span({φ i } p i=0 ) In order to approximate y we wish to find a finite system for at least the first few α i. We take the weighted inner product with the ith basis, i = 0,..., p, j=0 α j (t) φ j, φ i W + α j (t)(j φ j 1, φ i W + φ j+1, φ i W ) = 0, where f (ξ), g(ξ) W := f (ξ)g(ξ)w (ξ)dξ. By orthogonality, φ j, φ i W = φ i, φ i W δ ij, we have α i φ i, φ i W + (i + 1)α i+1 φ i, φ i W + α i 1 φ i, φ i W = 0, i = 0,..., p. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
58 Appendix Polynomial Chaos Deterministic ODE system Let α represent the vector containing α 0 (t),..., α p (t). Assuming α 1 (t), α p+1 (t), etc., are identically zero, the system of ODEs can be written α + M α = 0, with M = p 1 0 The degree p PC approximation is y(t, ξ) y p (t, ξ) = p j=0 α j(t)φ j (ξ). The mean value E[y(t, ξ)] E[y p (t, ξ)] = α 0 (t). The variance Var(y(t, ξ)) p j=1 α j(t) 2. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
59 Appendix Polynomial Chaos Figure : Convergence of error with Gaussian random variable by Hermitian-chaos. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
60 Appendix Polynomial Chaos Generalizations Consider the non-homogeneous IVP ẏ + ky = g(t), y(0) = y 0 with then k = k(ξ) = σξ + µ, ξ N (0, 1), α i + σ [(i + 1)α i+1 + α i 1 ] + µα i = g(t)δ 0i, i = 0,..., p, or the deterministic ODE system is α + (σm + µi ) α = g(t) e 1. Note that the initial condition for the PC system is α(0) = y 0 e 1. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
61 Appendix Polynomial Chaos Generalizations For any choice of family of orthogonal polynomials, there exists a triple recursion formula. Given the arbitrary relation ξφ j = a j φ j 1 + b j φ j + c j φ j+1 (with φ 1 = 0) then the matrix above becomes b 0 a 1 c 0 b 1 a 2 M = ap c p 1 b p (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
62 Appendix Polynomial Chaos Generalized Polynomial Chaos Table : Popular distributions and corresponding orthogonal polynomials. Distribution Polynomial Support Gaussian Hermite (, ) gamma Laguerre [0, ) beta Jacobi [a, b] uniform Legendre [a, b] Note: lognormal random variables may be handled as a non-linear function (e.g., Taylor expansion) of a normal random variable. (N.L. Gibson, OSU) Maxwell-PC Debye NIST-UQ4Mat / 35
Polynomial Chaos Approach for Maxwell s Equations in Dispersive Media
Polynomial Chaos Approach for Maxwell s Equations in Dispersive Media Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 15, 2013 Prof. Gibson (OSU) PC-FDTD
More informationElectromagnetic Relaxation Time Distribution Inverse Problems in the Time-domain
Electromagnetic Relaxation Time Distribution Inverse Problems in the Time-domain Prof Nathan L Gibson Department of Mathematics Joint Math Meeting Jan 9, 2011 Prof Gibson (OSU) Inverse Problems for Distributions
More informationElectromagnetic Modeling and Simulation
Electromagnetic Modeling and Simulation Erin Bela and Erik Hortsch Department of Mathematics Research Experiences for Undergraduates April 7, 2011 Bela and Hortsch (OSU) EM REU 2010 1 / 45 Maxwell s Equations
More informationStability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media
Contemporary Mathematics Volume 586 013 http://dx.doi.org/.90/conm/586/11666 Stability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media V. A. Bokil and N.
More informationParameter Estimation Versus Homogenization Techniques in Time-Domain Characterization of Composite Dielectrics
Parameter Estimation Versus Homogenization Techniques in Time-Domain Characterization of Composite Dielectrics H. T. Banks 1 V. A. Bokil and N. L. Gibson, 3 Center For Research in Scientific Computation
More informationSIMULATING POLYDISPERSE MATERIALS WITH DISTRIBUTIONS OF THE DEBYE MODEL
SIMULATING POLYDISPERSE MATERIALS WITH DISTRIBUTIONS OF THE DEBYE MODEL DANIELLE WEDDE AND MARIE MILNE ADVISOR: PROFESSOR GIBSON OREGON STATE UNIVERSITY ABSTRACT. We investigate models which approximate
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationParameter estimation versus homogenization techniques in time-domain characterization of composite dielectrics
c de Gruyter 7 J. Inv. Ill-Posed Problems 15 (7), 117 135 DOI 1.1515 / JIIP.7.6 Parameter estimation versus homogenization techniques in time-domain characterization of composite dielectrics H. T. Banks,
More informationNumerical methods for non-standard fractional operators in the simulation of dielectric materials
Numerical methods for non-standard fractional operators in the simulation of dielectric materials Roberto Garrappa Department of Mathematics - University of Bari - Italy roberto.garrappa@uniba.it Fractional
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume Number 4 ages 657 687 c 04 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume XX Number 0 ages 45 c 03 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND LORNTZ
More informationCONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING FDTD METHODS FOR THE 3D MAXWELL S EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5 Number 6 Pages 747 784 c 08 Institute for Scientific Computing and Information CONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING
More informationFinding damage in Space Shuttle foam
Finding damage in Space Shuttle foam Nathan L. Gibson gibsonn@math.oregonstate.edu In Collaboration with: Prof. H. T. Banks, CRSC Dr. W. P. Winfree, NASA OSU REU p. 1 Motivating Application The particular
More informationA thesis submitted to the University of Manchester for the degree of Doctor of Philosophy in the Faculty of Science and Engineering
INTRODUCTION OF THE DEBYE MEDIA TO THE FILTERED FINITE-DIFFERENCE TIME-DOMAIN METHOD WITH COMPLEX-FREQUENCY-SHIFTED PERFECTLY MATCHED LAYER ABSORBING BOUNDARY CONDITIONS A thesis submitted to the University
More informationarxiv: v1 [math.na] 3 Oct 2018
DISPERSION ANALYSIS OF FINITE DIFFERENCE AND DISCONTINUOUS GALERKIN SCHEMES FOR MAXELL S EQUATIONS IN LINEAR LORENTZ MEDIA YAN JIANG, PUTTHA SAKKAPLANGKUL, VRUSHALI A. BOKIL, YINGDA CHENG, AND FENGYAN
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationPolynomial chaos expansions for sensitivity analysis
c DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for sensitivity analysis B. Sudret Chair of Risk, Safety & Uncertainty
More information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationAntennas and Propagation. Chapter 2: Basic Electromagnetic Analysis
Antennas and Propagation : Basic Electromagnetic Analysis Outline Vector Potentials, Wave Equation Far-field Radiation Duality/Reciprocity Transmission Lines Antennas and Propagation Slide 2 Antenna Theory
More informationStructural Engineering Research Group (SERG) April 10, 2013 Dielectric Modeling of Hydrated Cement Paste Panels
Structural Engineering Research Group (SERG) April 10, 2013 Dielectric Modeling of Hydrated Cement Paste Panels Hao Liu Advisor: Dr. Tzu-Yang Yu Department of Civil and Environmental Engineering University
More informationElectromagnetic (EM) Waves
Electromagnetic (EM) Waves Short review on calculus vector Outline A. Various formulations of the Maxwell equation: 1. In a vacuum 2. In a vacuum without source charge 3. In a medium 4. In a dielectric
More informationHigh-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials
Progress In Electromagnetics Research Letters, Vol. 77, 103 107, 01 High-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials Wei-Jun Chen
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationElectromagnetic Waves Across Interfaces
Lecture 1: Foundations of Optics Outline 1 Electromagnetic Waves 2 Material Properties 3 Electromagnetic Waves Across Interfaces 4 Fresnel Equations 5 Brewster Angle 6 Total Internal Reflection Christoph
More informationPolynomial chaos expansions for structural reliability analysis
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for structural reliability analysis B. Sudret & S. Marelli Incl.
More informationElectrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Electro Dynamic
Electrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Name Electro Dynamic Instructions: Use SI units. Short answers! No derivations here, just state your responses clearly. 1. (2) Write an
More informationTheory and Applications of Dielectric Materials Introduction
SERG Summer Seminar Series #11 Theory and Applications of Dielectric Materials Introduction Tzuyang Yu Associate Professor, Ph.D. Structural Engineering Research Group (SERG) Department of Civil and Environmental
More informationElectromagnetic Wave Propagation Lecture 2: Time harmonic dependence, constitutive relations
Electromagnetic Wave Propagation Lecture 2: Time harmonic dependence, constitutive relations Daniel Sjöberg Department of Electrical and Information Technology September 2, 2010 Outline 1 Harmonic time
More informationEnergy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media
Energy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media Yingda Cheng Michigan State University Computational Aspects of Time Dependent Electromagnetic Wave Problems
More informationDIELECTRIC SPECTROSCOPY. & Comparison With Other Techniques
DIELECTRIC SPECTROSCOPY & Comparison With Other Techniques DIELECTRIC SPECTROSCOPY measures the dielectric and electric properties of a medium as a function of frequency (time) is based on the interaction
More informationElectromagnetically Induced Flows in Water
Electromagnetically Induced Flows in Water Michiel de Reus 8 maart 213 () Electromagnetically Induced Flows 1 / 56 Outline 1 Introduction 2 Maxwell equations Complex Maxwell equations 3 Gaussian sources
More informationPolynomial Chaos Expansions for Random Ordinary Differential Equations
Polynomial Chaos Expansions for Random Ordinary Differential Equations Brian McKenzie Advisor: Dr. Nathan Gibson January 6, 2012 Abstract We consider numerical methods for finding approximate solutions
More informationReport submitted to Prof. P. Shipman for Math 540, Fall 2009
Dynamics at the Horsetooth Volume 1, 009. Three-Wave Interactions of Spin Waves Aaron Hagerstrom Department of Physics Colorado State University aaronhag@rams.colostate.edu Report submitted to Prof. P.
More informationElectromagnetic Wave Propagation Lecture 3: Plane waves in isotropic and bianisotropic media
Electromagnetic Wave Propagation Lecture 3: Plane waves in isotropic and bianisotropic media Daniel Sjöberg Department of Electrical and Information Technology September 2016 Outline 1 Plane waves in lossless
More informationTypical anisotropies introduced by geometry (not everything is spherically symmetric) temperature gradients magnetic fields electrical fields
Lecture 6: Polarimetry 1 Outline 1 Polarized Light in the Universe 2 Fundamentals of Polarized Light 3 Descriptions of Polarized Light Polarized Light in the Universe Polarization indicates anisotropy
More informationGeneral review: - a) Dot Product
General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,
More informationFor the magnetic field B called magnetic induction (unfortunately) M called magnetization is the induced field H called magnetic field H =
To review, in our original presentation of Maxwell s equations, ρ all J all represented all charges, both free bound. Upon separating them, free from bound, we have (dropping quadripole terms): For the
More informationin Electromagnetics Numerical Method Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD
2141418 Numerical Method in Electromagnetics Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD ISE, Chulalongkorn University, 2 nd /2018 Email: charusluk.v@chula.ac.th Website: Light
More informationELECTROMAGNETISM SUMMARY
Review of E and B ELECTROMAGNETISM SUMMARY (Rees Chapters 2 and 3) The electric field E is a vector function. E q o q If we place a second test charged q o in the electric field of the charge q, the two
More informationOptics and Optical Design. Chapter 5: Electromagnetic Optics. Lectures 9 & 10
Optics and Optical Design Chapter 5: Electromagnetic Optics Lectures 9 & 1 Cord Arnold / Anne L Huillier Electromagnetic waves in dielectric media EM optics compared to simpler theories Electromagnetic
More informationBulk permittivity of a composite with coated spheroidal filler particles
JOURNAL OF MATERIALS SCIENCE 5 (2000)5809 586 Bulk permittivity of a composite with coated spheroidal filler particles N. HARFIELD Center for Nondestructive Evaluation, Iowa State University, 95 Scholl
More informationModule I: Electromagnetic waves
Module I: Electromagnetic waves Lecture 3: Time dependent EM fields: relaxation, propagation Amol Dighe TIFR, Mumbai Outline Relaxation to a stationary state Electromagnetic waves Propagating plane wave
More informationMSE 383, Unit 3-3. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.
Dynamic Mechanical Behavior MSE 383, Unit 3-3 Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Scope Why DMA & TTS? DMA Dynamic Mechanical Behavior (DMA) Superposition Principles
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 9
WiSe 202 20.2.202 Prof. Dr. A-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg
More informationMultilayer Reflectivity
Multilayer Reflectivity John E. Davis jed@jedsoft.org January 5, 2014 1 Introduction The purpose of this document is to present an ab initio derivation of the reflectivity for a plane electromagnetic wave
More informationComparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer
Comparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer V. A. Bokil a and M. W. Buksas b Center for Research in Scientific Computation a North Carolina
More informationAn inverse problem for eddy current equations
An inverse problem for eddy current equations Alberto Valli Department of Mathematics, University of Trento, Italy In collaboration with: Ana Alonso Rodríguez, University of Trento, Italy Jessika Camaño,
More informationDielectric Spectroscopy Analysis of Recycled PET with Different Synthetic Polymer Blends
Dielectric Spectroscopy Analysis of Recycled PET with Different Synthetic Polymer Blends Aradoaei S. 1 Constantinescu 1 G. Ciobanu R. 1 1 Technical University Gh. Asachi Iasi D. Mangeron no. 5, tel:+/7868,
More informationElectromagnetic Waves
Physics 8 Electromagnetic Waves Overview. The most remarkable conclusion of Maxwell s work on electromagnetism in the 860 s was that waves could exist in the fields themselves, traveling with the speed
More informationON FRACTIONAL RELAXATION
Fractals, Vol. 11, Supplementary Issue (February 2003) 251 257 c World Scientific Publishing Company ON FRACTIONAL RELAXATION R. HILFER ICA-1, Universität Stuttgart Pfaffenwaldring 27, 70569 Stuttgart,
More informationElectromagnetic Wave Propagation Lecture 2: Time harmonic dependence, constitutive relations
Electromagnetic Wave Propagation Lecture 2: Time harmonic dependence, constitutive relations Daniel Sjöberg Department of Electrical and Information Technology September 2015 Outline 1 Harmonic time dependence
More informationExponential Recursion for Multi-Scale Problems in Electromagnetics
Exponential Recursion for Multi-Scale Problems in Electromagnetics Matthew F. Causley Department of Mathematics Kettering University Computational Aspects of Time Dependent Electromagnetic Wave Problems
More informationCHAPTER 9 ELECTROMAGNETIC WAVES
CHAPTER 9 ELECTROMAGNETIC WAVES Outlines 1. Waves in one dimension 2. Electromagnetic Waves in Vacuum 3. Electromagnetic waves in Matter 4. Absorption and Dispersion 5. Guided Waves 2 Skip 9.1.1 and 9.1.2
More information22 Phasor form of Maxwell s equations and damped waves in conducting media
22 Phasor form of Maxwell s equations and damped waves in conducting media When the fields and the sources in Maxwell s equations are all monochromatic functions of time expressed in terms of their phasors,
More informationJoel A. Shapiro January 20, 2011
Joel A. shapiro@physics.rutgers.edu January 20, 2011 Course Information Instructor: Joel Serin 325 5-5500 X 3886, shapiro@physics Book: Jackson: Classical Electrodynamics (3rd Ed.) Web home page: www.physics.rutgers.edu/grad/504
More informationElectrical and optical properties of materials
Electrical and optical properties of materials John JL Morton Part 4: Mawell s Equations We have already used Mawell s equations for electromagnetism, and in many ways they are simply a reformulation (or
More information(a) Show that the amplitudes of the reflected and transmitted waves, corrrect to first order
Problem 1. A conducting slab A plane polarized electromagnetic wave E = E I e ikz ωt is incident normally on a flat uniform sheet of an excellent conductor (σ ω) having thickness D. Assume that in space
More informationElectrodynamics Qualifier Examination
Electrodynamics Qualifier Examination August 15, 2007 General Instructions: In all cases, be sure to state your system of units. Show all your work, write only on one side of the designated paper, and
More informationMetamaterials. Peter Hertel. University of Osnabrück, Germany. Lecture presented at APS, Nankai University, China
University of Osnabrück, Germany Lecture presented at APS, Nankai University, China http://www.home.uni-osnabrueck.de/phertel Spring 2012 are produced artificially with strange optical properties for instance
More informationUniform Plane Waves. Ranga Rodrigo. University of Moratuwa. November 7, 2008
Uniform Plane Waves Ranga Rodrigo University of Moratuwa November 7, 2008 Ranga Rodrigo (University of Moratuwa) Uniform Plane Waves November 7, 2008 1 / 51 Summary of Last Week s Lecture Basic Relations
More informationAccurate Computation of Vector Potentials in Lossy Media
Accurate Computation of Vector Potentials in Lossy Media Swagato Chakraborty and Vikram Jandhyala {swagato,jandhyala}@ee.washington.edu Dept of EE, University of Washington Seattle WA, 98195-500 UWEE Technical
More informationA Polynomial Chaos Approach to Robust Multiobjective Optimization
A Polynomial Chaos Approach to Robust Multiobjective Optimization Silvia Poles 1, Alberto Lovison 2 1 EnginSoft S.p.A., Optimization Consulting Via Giambellino, 7 35129 Padova, Italy s.poles@enginsoft.it
More informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
More information1 Electromagnetic concepts useful for radar applications
Electromagnetic concepts useful for radar applications The scattering of electromagnetic waves by precipitation particles and their propagation through precipitation media are of fundamental importance
More informationRadio Propagation Channels Exercise 2 with solutions. Polarization / Wave Vector
/8 Polarization / Wave Vector Assume the following three magnetic fields of homogeneous, plane waves H (t) H A cos (ωt kz) e x H A sin (ωt kz) e y () H 2 (t) H A cos (ωt kz) e x + H A sin (ωt kz) e y (2)
More information9 The conservation theorems: Lecture 23
9 The conservation theorems: Lecture 23 9.1 Energy Conservation (a) For energy to be conserved we expect that the total energy density (energy per volume ) u tot to obey a conservation law t u tot + i
More informationANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA
Progress In Electromagnetics Research Letters, Vol. 22, 83 93, 2011 ANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA X. Ai Science and Technology
More informationCausality. but that does not mean it is local in time, for = 1. Let us write ɛ(ω) = ɛ 0 [1 + χ e (ω)] in terms of the electric susceptibility.
We have seen that the issue of how ɛ, µ n depend on ω raises questions about causality: Can signals travel faster than c, or even backwards in time? It is very often useful to assume that polarization
More informationAdvanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III
Advanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III 1 11. Yee s discrete algorithm Maxwell s equations are discretized using central FDs. We set the magnetic loss equal to
More informationELE3310: Basic ElectroMagnetic Theory
A summary for the final examination EE Department The Chinese University of Hong Kong November 2008 Outline Mathematics 1 Mathematics Vectors and products Differential operators Integrals 2 Integral expressions
More informationMicroscopic-Macroscopic connection. Silvana Botti
relating experiment and theory European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre for Computational
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationLecture notes for ELECTRODYNAMICS.
Lecture notes for 640-343 ELECTRODYNAMICS. 1 Summary of Electrostatics 1.1 Coulomb s Law Force between two point charges F 12 = 1 4πɛ 0 Q 1 Q 2ˆr 12 r 1 r 2 2 (1.1.1) 1.2 Electric Field For a charge distribution:
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationPHYS4210 Electromagnetic Theory Quiz 1 Feb 2010
PHYS4210 Electromagnetic Theory Quiz 1 Feb 2010 1. An electric dipole is formed from two charges ±q separated by a distance b. For large distances r b from the dipole, the electric potential falls like
More informationMaxwell s Equations. In the previous chapters we saw the four fundamental equations governging electrostatics and magnetostatics. They are.
Maxwell s Equations Introduction In the previous chapters we saw the four fundamental equations governging electrostatics and magnetostatics. They are D = ρ () E = 0 (2) B = 0 (3) H = J (4) In the integral
More informationH. T. Banks and V. A. Bokil. Center for Research in Scientific Computation North Carolina State University Raleigh, N.C
Parameter Identification for Dispersive Dielectrics Using Pulsed Microwave Interrogating Signals and Acoustic Wave Induced Reflections in Two and Three Dimensions H. T. Banks and V. A. Bokil Center for
More informationIntroduction to Uncertainty Quantification in Computational Science Handout #3
Introduction to Uncertainty Quantification in Computational Science Handout #3 Gianluca Iaccarino Department of Mechanical Engineering Stanford University June 29 - July 1, 2009 Scuola di Dottorato di
More informationPlane Waves and Planar Boundaries in FDTD Simulations
Plane Waves and Planar Boundaries in FDTD Simulations School of Electrical Engineering and Computer Science John B. Schneider Elec. Eng. & Comp. Sci. Washington State Univ. Pullman, WA Robert J. Kruhlak
More informationChapter Three: Propagation of light waves
Chapter Three Propagation of Light Waves CHAPTER OUTLINE 3.1 Maxwell s Equations 3.2 Physical Significance of Maxwell s Equations 3.3 Properties of Electromagnetic Waves 3.4 Constitutive Relations 3.5
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationJust solve this: macroscopic Maxwell s equations
ust solve this: macroscopic Maxwell s equations Faraday: E = B t constitutive equations (here, linear media): Ampere: H = D t + (nonzero frequency) Gauss: D = ρ B = 0 D = ε E B = µ H electric permittivity
More informationElectromagnetic Wave Propagation Lecture 8: Propagation in birefringent media
Electromagnetic Wave Propagation Lecture 8: Propagation in birefringent media Daniel Sjöberg Department of Electrical and Information Technology September 27, 2012 Outline 1 Introduction 2 Maxwell s equations
More informationPart E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann MAXWELL Grid Equations in time domain d 1 h() t MC e( t) dt d 1 e() t M Ch() t j( t) dt Transient Fields system of 1 st order
More informationB. H. Jung Department of Information and Communication Engineering Hoseo University Asan, Chungnam , Korea
Progress In Electromagnetics Research, PIER 77, 111 120, 2007 ANALYSIS OF TRANSIENT ELECTROMAGNETIC SCATTERING WITH PLANE WAVE INCIDENCE USING MOD-FDM B. H. Jung Department of Information and Communication
More informationImaginary. Axis. Real. Axis
Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 25 3:3-6:3 p.m. ffl Print your name and sign the honor code statement ffl You may use your course notes,
More informationGeneralized Analysis of Stability and Numerical Dispersion in the Discrete-Convolution FDTD Method
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 48, NO. 6, JUNE 2000 887 Generalized Analysis of Stability and Numerical Dispersion in the Discrete-Convolution FDTD Method William A. Beck, Member,
More informationPhysics 506 Winter 2004
Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationAn iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems
An iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems K. Hertel 1,2, S. Yan 1,2, C. Pflaum 1,2, R. Mittra 3 kai.hertel@fau.de 1 Lehrstuhl für Systemsimulation
More informationNonstandard Finite Difference Time Domain Algorithm for Berenger s Perfectly Matched Layer
ACES JOURNAL, VOL. 6, NO., FEBRUARY 011 153 Nonstandard Finite Difference Time Domain Algorithm for Berenger s Perfectly Matched Layer Naoki Okada and James B. Cole Graduate School of Systems and Information
More informationFinal Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications
Final Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications PI: George Em Karniadakis Division of Applied Mathematics, Brown University April 25, 2005 1 Objectives
More informationE & M Qualifier. January 11, To insure that the your work is graded correctly you MUST:
E & M Qualifier 1 January 11, 2017 To insure that the your work is graded correctly you MUST: 1. use only the blank answer paper provided, 2. use only the reference material supplied (Schaum s Guides),
More informationFinite Difference Solution of Maxwell s Equations
Chapter 1 Finite Difference Solution of Maxwell s Equations 1.1 Maxwell s Equations The principles of electromagnetism have been deduced from experimental observations. These principles are Faraday s law,
More informationChapter 2. Dielectric Theories
Chapter Dielectric Theories . Dielectric Theories 1.1. Introduction Measurements of dielectric properties of materials is very important because it provide vital information regarding the material characteristics,
More informationThe dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum
Home Search Collections Journals About Contact us My IOPscience The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues This article has been downloaded
More informationWaves. Daniel S. Weile. ELEG 648 Waves. Department of Electrical and Computer Engineering University of Delaware. Plane Waves Reflection of Waves
Waves Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Waves Outline Outline Introduction Let s start by introducing simple solutions to Maxwell s equations
More informationImaginary. Axis. Real. Axis
Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 2005 3:30-6:30 p.m. ffl Print your name and sign the honor code statement ffl You may use your course
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More information