Numerical methods for non-standard fractional operators in the simulation of dielectric materials
|
|
- Jennifer Green
- 5 years ago
- Views:
Transcription
1 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 PDEs: Theory, Algorithms and Applications ICERM - Providence, June 18 22, 2018 R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
2 Outline Numerical methods for non-standard fractional operators in the simulation of dielectric materials 1 The problem 2 The operator 3 The numerical method Financial support: EU Cost Action Fractional Systems Istituto Nazionale di Alta Matematica (INdAM-GNCS) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
3 Maxwell s equations: Standard Maxwell s equations: H = ɛ 0 t E H E = µ 0 t E: electric field H: magnetic field Ampere s law Faraday s law R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
4 Maxwell s equations: Standard Maxwell s equations: H = ɛ 0 t E H E = µ 0 t E: electric field H: magnetic field Ampere s law Faraday s law Real world applications design of data and energy storage devices design of antennas medical diagnostic (MRI), cancer therapy,... R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
5 Maxwell s equations: Standard Maxwell s equations with polarization: H = ɛ 0 t E + t P H E = µ 0 t Ampere s law Faraday s law E: electric field H: magnetic field P: polarization R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
6 Maxwell s equations: Standard Maxwell s equations with polarization: H = ɛ 0 t E + t P H E = µ 0 t Ampere s law Faraday s law E: electric field H: magnetic field P: polarization The complex susceptibility The polarization P depends on the electric field E (constitutive law) ˆP = ε 0 ˆχ(ω) Ê ˆχ(ω) is a specific feature of the matter (or system) Simplified notation (just for easy of presentation) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
7 Determining the complex susceptibility ˆχ(ω) How to derive ˆχ(ω) = ˆχ (ω) iˆχ (ω)? R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
8 Determining the complex susceptibility ˆχ(ω) How to derive ˆχ(ω) = ˆχ (ω) iˆχ (ω)? Experimental data (in the frequency domain): ˆP = ˆχ(ω)Ê 10 0 Exper. data 10 1 ˆχ (ω) frequency ω Match experimental data into a mathematical model R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
9 Determining the complex susceptibility ˆχ(ω) 10 0 Exper. data Debye model 10 1 ˆχ (ω) frequency ω The Debye model: ˆχ(ω) = iω (standard materials) Ordinary differential equation: d P(t) + P(t) = E(t) dt R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
10 Determining the complex susceptibility ˆχ(ω) 10 0 Exper. data 10 1 ˆχ (ω) frequency ω Materials with anomalous dielectric properties: amorphous polymers complex systems (biological tissues) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
11 Determining the complex susceptibility ˆχ(ω) 10 0 Exper. data Debye model 10 1 ˆχ (ω) frequency ω Debye model not satisfactory R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
12 Determining the complex susceptibility ˆχ(ω) Exper. data Debye model C C model ˆχ (ω) frequency ω The Cole-Cole model: ˆχ(ω) = Fractional differential equation: Only partially satisfactory (iω) α 0 < α < 1 α d α P(t) + P(t) = E(t) dtα R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
13 Determining the complex susceptibility ˆχ(ω) Exper. data Debye model C C model H N model ˆχ (ω) frequency ω The Havriliak-Negami model: ˆχ(ω) = Better matching thanks to three parameters α, γ and 1 ( 1 + (iω) α ) γ 0 < α, αγ 1 S. Havriliak and S. Negami A complex plane representation of dielectric and mechanical relaxation processes in some polymers. In: Polymer (1967) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
14 Determining the complex susceptibility ˆχ(ω) Exper. data Debye model C C model H N model ˆχ (ω) frequency ω The Havriliak-Negami model: ˆχ(ω) = 1 ( 1 + (iω) α ) γ 0 < α, αγ < 1 Fractional pseudo-differential equation: (1 + α d α ) γ P(t) = E(t)? dt α R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
15 Other models for complex susceptibility ˆχ(ω) Modified Havriliak-Negami or JWS (Jurlewicz, Weron and Stanislavsky) ˆχ JWS (iω) = 1 1 (1 + ( i ω ) α ) γ = 1 (i ω) αγ χ HN (iω) EW: Excess wing (Hilfer, Nigmatullin and others) ˆχ EW (iω) = Multichannel excess wing (Hilfer) 1 + ( 2 iω) α 1 + ( 2 iω) α + 1 iω. 1 ˆξ(s) = [ n 1 + (iω k ) α k k=1 This talk focuses on the Havriliak-Negami model ] 1 Garrappa R., Mainardi F. and Maione G., Models of dielectric relaxation based on completely monotone functions. In: Frac. Calc. Appl. Anal. 19(5) (2016) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
16 Dealing with the Havriliak-Negami model ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) Few contributions on simulation of this constitutive law: C.S.Antonopoulos, N.V.Kantartzis, I.T.Rekanos FDTD Method for Wave Propagation in Havriliak-Negami Media Based on Fractional Derivative Approximation. In: IEEE Trans Magn. 53(6) (2017) P.Bia et al. A novel FDTD formulation based on fractional derivatives for dispersive Havriliak Negami media. In: Signal Processing, 107 (2015) M.F.Causley, P.G.Petropoulos, On the Time-Domain Response of Havriliak-Negami Dielectrics. In: IEEE Trans. Antennas Propag., 61(6) (2013) M.F.Causley, P.G.Petropoulos, and S. Jiang Incorporating the Havriliak-Negami dielectric model in the FD-TD method. In: J. Comput. Phys., 230 (2011), Main problems: Define time-domain operator for HN Discretize the operator for simulations R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
17 Operators in the time domain ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) ( α 0 D α t + 1 ) γ P(t) = E(t)??? How to define the fractional pseudo-differential operator ( α 0D α t + 1 ) γ? R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
18 Operators in the time domain ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) ( α 0 D α t + 1 ) γ P(t) = E(t)??? How to define the fractional pseudo-differential operator ( α 0D α t + 1 ) γ? Combination of fractional operators ( α 0 Dt α + 1 ) γ = exp ( t ) α α 0 Dt 1 α αγ 0D αγ t ( t ) exp α α 0 Dt 1 α Useful for theoretical investigations R.R.Nigmatullin and Y.E.Ryabov Cole Davidson dielectric relaxation as a self similar relaxation process. In: Physics of the Solid State 39.1 (1997) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
19 Operators in the time domain ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) ( α 0 D α t + 1 ) γ P(t) = E(t)??? How to define the fractional pseudo-differential operator ( α 0D α t + 1 ) γ? Expansion in infinite series ( α 0 D α t + 1 ) γ = k=0 ( ) γ α(γ k) 0D α(γ k) t k No satisfactory for error control V.Novikov et al. Anomalous relaxation in dielectrics. Equations with fractional derivatives. In: Mater. Sci. Poland 23.4 (2005) P.Bia et al. A novel FDTD formulation based on fractional derivatives for dispersive Havriliak Negami media. In: Signal Processing, 107 (2015) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
20 Operators in the time domain ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) ( α 0 D α t + 1 ) γ P(t) = E(t)??? How to define the fractional pseudo-differential operator ( α 0D α t + 1 ) γ? Prabhakar derivative ( α 0 D α t + 1 ) γ P(t) = d dt t 0 ( u Prabhakar function: E γ α,β (z) = 1 Γ(γ) ) αγ ( ( E γ u ) α ) α,1 αγ P(t u) du k=0 Γ(γ + k)z k k!γ(αk + β) R.Garra, R.Gorenflo, F.Polito and Z.Tomovski Hilfer-Prabhakar derivatives and some applications. In: Appl. Math. Comput. 242 (2014) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
21 Operators in the time domain ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) ( α 0 D α t + 1 ) γ P(t) = E(t)??? How to define the fractional pseudo-differential operator ( α 0D α t + 1 ) γ? Prabhakar derivative of Caputo type ( α 0 Dt α + 1 ) γ d t ( u ) αγ ( ( P(t) = E γ u ) α ) α,1 αγ P(t u) du dt 0 RL C( α 0 Dt α + 1 ) t ( γ u P(t) = 0 ) αγ ( ( E γ u ) α ) α,1 αγ P (t u) du Caputo C( α 0 D α t + 1 ) γ P(t) = ( α 0 D α t + 1 ) γ ( P(t) P(0 + ) ) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
22 Operators in the time domain ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) ( α 0 D α t + 1 ) γ P(t) = E(t)??? How to define the fractional pseudo-differential operator ( α 0D α t + 1 ) γ? Fractional differences of Grünwald Letnikov ( α 0 Dt α + 1 ) γ 1 P(t) = lim h 0 h αγ Ω k P(t kh) k=0 Ω 0 = 1, Ω k = α α + h α k ( (1 + γ)j j=1 k ) 1 ω (α) j Ω k j, Ω k = ( α + h α)γ Ω k R.Garrappa On Grünwald Letnikov operators for fractional relaxation in Havriliak-Negami models. In: Commun. Nonlinear. Sci. Numer. Simul. 38 (2016) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
23 Operators in the time domain ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) ( α 0 D α t + 1 ) γ P(t) = E(t)??? How to define the fractional pseudo-differential operator ( α 0D α t + 1 ) γ? Fractional differences of Grünwald Letnikov ( α 0 Dt α + 1 ) γ P(tn ) 1 n h αγ Ω k P(t n k ) + O ( h ) C( α 0 D α t + 1 ) γ P(tn ) 1 h αγ k=0 n ( Ω k P(tn k ) P(t 0 ) ) + O ( h ) k=0 Approximation of first order only! R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
24 Integral operator ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) P(t) = ( α 0J α t + 1 ) γ E(t) Prabhakar integral ( α 0J α t P(t) = 1 αγ t ) γ (t u) αγ 1 E γ α,αγ ( ( t u ) α ) E(u) du, R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
25 Integral operator ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) P(t) = ( α 0J α t + 1 ) γ E(t) Prabhakar integral ( α 0J α t P(t) = 1 αγ t ) γ (t u) αγ 1 E γ α,αγ ( ( t u ) α ) E(u) du, Non-local operator Weakly singular integral Series of RL integrals 1 ( α 0J α t + 1 ) γ ( ) γ 1 = ( 1) k k k=0 Jα(k+γ) α(k+γ) 0 1 Giusti A., A comment on some new definitions of fractional derivative. Nonlinear Dyn. (2018) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
26 Integral operator ˆP(ω) = 1 ((iω) α + 1) γ Ê(ω) P(t) = ( α 0J α t + 1 ) γ E(t) Prabhakar integral ( α 0J α t P(t) = 1 αγ t ) γ (t u) αγ 1 E γ α,αγ ( ( t u ) α ) E(u) du, Preferable in numerical simulation Numerical integration is easier than numerical differentiation Approximation by quadrature rules What kind of quadrature rule? R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
27 Discretization of Maxwell s equations H = ɛ 0 t E + t P H E = µ 0 t Ampere s law Faraday s law F = ( Fz y F ) ( y Fx i + z z F ) ( z Fy j + x x F ) x k y The Finite Differences Time Domain (FDTD) method: 1 Introduced by K. Yee in Approximation based on centred finite difference operators 3 Use of staggered grids in space and time R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
28 Discretization of Maxwell s equations - FDTD t n+2 t n+1 t n Ampere s law Faraday s law t n 1 t n 2 x m 2 x m 1 x m x m+1 x m+1 H n+ 1 2,m = H m+ 1 2,n+ 1 2 H m 1 2,n+ 1 2 x + O ( 2 ) x E m+ 1 2,n = E m+1,n H m 1,n + O ( 2 ) x x R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
29 Discretization of Maxwell s equations - FDTD t n+2 t n+1 t n E H t n 1 t n 2 x m 2 x m 1 x m x m+1 x m+1 H n+ 1 2,m = H m+ 1 2,n+ 1 2 H m 1 2,n+ 1 2 x + O ( 2 ) x E m+ 1 2,n = E m+1,n H m 1,n + O ( 2 ) x x R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
30 Discretization of the constitutive law H = ɛ 0 t E + t P E = µ 0 H t Ampere s law Faraday s law P = ( α 0J α t + 1 ) γ E Constitutive law It is necessary a second-order approximation Use of the values of E just on grid points (t n, x m ) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
31 Discretization of the constitutive law H = ɛ 0 t E + t P E = µ 0 H t Ampere s law Faraday s law P = ( α 0J α t + 1 ) γ E Constitutive law It is necessary a second-order approximation Use of the values of E just on grid points (t n, x m ) Convolution quadratures by Lubich [Lubich, 1988] Based on the LT ˆχ(s) = 1 ( α s α + 1 ) γ (i.e. the complex susceptibility) Extension of classical linear multistep methods (LMM) for ODEs Preservation of convergence and stability of the LMM R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
32 Convolution quadrature rules I (t) = t Generalizes LMM for ODEs k ρ j y n j = h j=0 0 f (t u)g(u)du k σ j f (t n j, y n j ), n k. j=0 Based on the generating function of the LMM ρ(ξ) = k ρ j ξ k j σ(ξ) = j=0 I (t n ) = h µ k σ j ξ k j j=0 ν n w n,j g(t j ) + h µ ω n j g(t j ) j=0 } {{ } starting term j=0 } {{ } convolution term δ(ξ) = ρ(1/ξ) σ(1/ξ) Lubich C. Convolution quadrature and discretized operational calculus (I and II), In: Numer. Math. 52, and (1988). R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
33 Convolution quadrature rules I (t) = t 0 ν n f (t u)g(u)du = I (t n ) = h µ w n,j g(t j ) + h µ ω n j g(t j ) j=0 } {{ } starting term j=0 } {{ } convolution term Use of LT F (s) of f (t) F (s) analytic in a sector Σ ϕ,c = { s C : arg(s c) < π ϕ } F (s) M s µ µ > 0 M < ( ) δ(ξ) Weights F = ω n ξ n h n=0 ϕ c Lubich C. Convolution quadrature and discretized operational calculus (I and II), In: Numer. Math. 52, and (1988). R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
34 Convolution quadrature for Havriliak-Negami I (t) = t 0 ν n f (t u)g(u)du = I (t n ) = h µ w n,j g(t j ) + h µ ω n j g(t j ) j=0 } {{ } starting term j=0 } {{ } convolution term ( f (t) = t αγ 1 Eα,αγ γ t α / α) 1 F (s) = ˆχ(s) = ( α s α + 1) γ F (s) analytic outside (, 0] F (s) M s αγ µ = αγ ( ) δ(ξ) Weights ˆχ = ω n ξ n h n=0 R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
35 Convolution quadrature for Havriliak-Negami P = ( α 0J α t + 1 ) γ E P(x, t n ) = h αγ ν n w n,j E(x, t j ) + h αγ ωn HN E(x, t j ) j=0 } {{ } starting term Trapezoidal rule y n+1 y n = h 2 Generating function δ(ξ) = Convolution weights h αγ n=0 2(1 ξ) 1 + ξ j=0 } {{ } convolution term ( f (tn, y n ) + f (t n+1, y n+1 ) ) ( ) 2(1 ξ) ωn HN ξ n = ˆχ h(1 + ξ) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
36 Evaluation of convolution weights h αγ n=0 ) γ ωn HN ξ n = ( h hαγ α (1 ξ)α h 2 α α + (1 + ξ) α, h = 2 1) Recursion, FFT and power of FPS 2) Numerical integration of Cauchy integral R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
37 Evaluation of convolution weights h αγ n=0 ) γ ωn HN ξ n = ( h hαγ α (1 ξ)α h 2 α α + (1 + ξ) α, h = 2 1) Recursion, FFT and power of FPS ( h α + (1 ξ) α = (1 + ξ) α = ˆω n ξ n ˆω 0 = 1 ˆω n = n=0 ( 1 α + 1 n ( ) 1 α ω n ξ n ω 0 = 1 ω n = 1 n n=0 ) ˆω n 1 ω n 1 (1 ξ) α (1 + ξ) α = ω n ξ n FFT of weights of (1 ξ) α and (1 + ξ) α n=0 ) γ ( ) γ ω n ξ n 1 ω n = (1 + h 1 + α ) γ 1 + h ω nξ n unitary FPS α n=0 n=1 R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
38 Evaluation of convolution weights h αγ n=0 ) γ ωn HN ξ n = ( h hαγ α (1 ξ)α h 2 α α + (1 + ξ) α, h = 2 1) Recursion, FFT and power of FPS Miller s Formula: power β C of unitary Formal Power Series ( 1 + a1 ξ + a 2 ξ 2 + a 3 ξ ) β (β) = v 0 + v (β) 1 ξ + v (β) 2 ξ 2 + v (β) 3 ξ n ( ) v (β) (β + 1)j 0 = 1, v n (β) = 1 a j v (β) n j n. j=1 Direct evaluation ω HN 0 = 1 ω HN n = n ( (1 γ)j j=1 n ) 1 ω j ωn j HN (1 ξ) α (1+ξ) α R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
39 Evaluation of convolution weights h αγ n=0 ) γ ωn HN ξ n = ( h hαγ α (1 ξ)α h 2 α α + (1 + ξ) α, h = 2 1) Recursion, FFT and power of FPS Weights evaluated in an exact way Expensive computation O(N 2 ) Instability due to long recurrences R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
40 Evaluation of convolution weights h αγ n=0 ) γ ωn HN ξ n = ( h hαγ α (1 ξ)α h 2 α α + (1 + ξ) α, h = 2 2) Numerical integration of Cauchy integral ω HN n = 1 n! d n dξ n G(ξ) G(ξ) = ( h α + ξ=0 Representation in terms of Cauchy integrals ωn HN = 1 ξ n 1 G(ξ)dξ, 2πi Approximation by quadrature rule C ) γ (1 ξ)α (1 + ξ) α R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
41 Evaluation of convolution weights h αγ n=0 ) γ ωn HN ξ n = ( h hαγ α (1 ξ)α h 2 α α + (1 + ξ) α, h = 2 2) Numerical integration of Cauchy integral Contour ξ = ρ 0 < ρ < 1 0 < ρ < 1 Equispaced quadrature nodes ξ l = ρe 2πil/L, l = 0, 1,..., L Trapezoidal rule ωn,l HN = ρ n L 1 G(ξ l )e 2πinl/L L l=0 Fast computation by FFT O(L log L) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
42 Evaluation of convolution weights h αγ n=0 ) γ ωn HN ξ n = ( h hαγ α (1 ξ)α h 2 α α + (1 + ξ) α, h = 2 2) Numerical integration of Cauchy integral Choice of ρ and L: error analysis 0 < ρ < Discretization error 0 < ρ < r < 1 ω n HN ωn,l HN ρ n (r/ρ) n a C r (r/ρ) L 1 Round-off error ω HN n,l ˆω n,l HN ρ n C ρ ɛ Balancing of the two errors: target accuracy ε > ɛ a Trefethen & Weideman in [SIAM Review 2014] R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
43 Numerical integration of Cauchy integral 10 8 Error α=0.8 γ=0.9 = 10 2 h= n 10 9 Error α=0.7 γ=0.8 = 10 3 h= n R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
44 Evaluation of the starting term P(x, t n ) = h αγ ν n w n,j E(x, t j ) + h αγ ωn HN E(x, t j ) j=0 } {{ } starting term j=0 } {{ } convolution term Handling the singularity at the origin w n,j obtained by imposing the rule exact for {1, t αγ, t} when αγ > 1 2 Solution of small systems of algebraic equations n n αγ E γ α,αγ αγ w ( (t n/) α ) j=0 n,0 w n,1 = γ α n 2αγ E γ α,2αγ+1 ( (t n j n/) α ) w n,2 j=0 n n αγ+1 E γ α,αγ+2 ( (t n/) α ) γ α = Γ(αγ + 1) Simple inversion of the matrix j=0 ω HN n j ω HN n j αγ ω HN n jj R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
45 Evaluation of Prabhakar (three-parameter ML) function Matlab code available Possible evaluation of E γ α,β (z) when arg(z) > απ Also available code for ML function with matrix arguments R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
46 Testing the quadrature rule N N P(t) = ( α 0Jt α + 1) γ E(t) E(t) = a k cos ω k t + b k sin ω k t k=0 k=0 4 α γ Set Set Set α=0.9 γ=0.9 = α=0.8 γ=0.8 =10 0 α=0.8 γ=0.7 = t R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
47 Numerical experiments: errors and EOC N N P(t) = ( α 0Jt α + 1) γ E(t) E(t) = a k cos ω k t + b k sin ω k t k=0 k=0 α = 0.9 γ = 0.9 α = 0.8 γ = 0.8 α = 0.8 γ = 0.7 = 10 2 = 10 0 = 10 1 h Error EOC Error EOC Error EOC ( 3) 6.21( 3) 2.04( 3) ( 3) ( 3) ( 4) ( 4) ( 4) ( 4) ( 5) ( 5) ( 5) ( 5) ( 5) ( 6) ( 6) ( 6) ( 6) EOC = log 2 ( E(h)/E( h 2 )) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
48 Numerical experiments: errors and EOC N N N P(t) = ( α 0Jt α + 1) γ E(t) E(t) = a k t αk + b k t γk + c k t k k=0 k=0 k=0 α = 0.9 γ = 0.9 α = 0.8 γ = 0.8 α = 0.8 γ = 0.7 = 10 2 = 10 0 = 10 1 h Error EOC Error EOC Error EOC ( 4) 1.26( 3) 1.43( 4) ( 5) ( 4) ( 5) ( 6) ( 5) ( 6) ( 6) ( 5) ( 6) ( 7) ( 6) ( 7) ( 7) ( 6) ( 7) EOC = log 2 ( E(h)/E( h 2 )) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
49 The Excess Wing model P(t) = 1 + ( 2 iω) α 1 + ( 2 iω) α + 1 iω E(t) 1) Linear multiterm fractional differential equation D t P(t) + 2 α 0Dt α P(t) + 1 ( 1 P(t) = E(t) + α ) 2 0Dt α E(t) ) Convolution quadrature rule h 1 α n=0 h α + 2 α α (1 ξ) α ωn EW ξ n = h 1 α 2 (1 + ξ) α h + 2 α 2 α (1 ξ)α h1 α (1 + ξ) α + (1 ξ) 1 (1 + ξ) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
50 Concluding remarks The Havriliak-Negami dielectric model can be described in the time-domain by Prabhakar derivatives and integrals Prabhakar derivatives and integrals are non standard fractional-order operators depending on several parameters Discretization of Prbhakar operators by means of Lubich s convolution quadratures: direct use of the complex susceptibility General approach: extension to other models: Excess Wing, multichannel Excess Wing, and others... Further developments Improving computational efficiency (memory term) Validation with experimental data (Technical University of Bari) Releasing robust software for solving Maxwell s systems with anomalous dielectric relaxation R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
51 Some references Garra R., Gorenflo R., Polito F. and Tomovski Z., Hilfer-Prabhakar derivatives and some applications. In: Appl. Math. Comput. 242 (2014) Garra R., Garrappa R., The Prabhakar or three parameter Mittag Leffler function: theory and application. In: CNSNS 56 (2018) Garrappa R., Numerical evaluation of two and three parameter Mittag-Leffler functions. In: SIAM J. Numer. Anal. 53(3) (2015) Garrappa R., Grünwald Letnikov operators for Havriliak Negami fractional relaxation. In: CNSNS 38 (2016) Garrappa R., Mainardi F. and Maione G., Models of dielectric relaxation based on completely monotone functions. In: FCAA 19(5) (2016) Giusti A., Colombaro R., Prabhakar-like fractional viscoelasticity. In: CNSNS 56 (2018) Lubich C. Convolution quadrature and discretized operational calculus (I and II), In: Numer. Math., 52 (1998) Mainardi F. and Garrappa R., On complete monotonicity of the Prabhakar function and non-debye relaxation in dielectrics. In: JCP 293 (2015) R.Garrappa (Univ. of Bari - Italy) Non-standard fractional operators ICERM / 25
WORKSHOP ON FRACTIONAL CALCULUS AND ITS APPLICATIONS Roma Tre, March 11, 2015 PROGRAM
On March 11, 2015, a workshop on fractional calculus was held, at the Department of Mathematics and Physics, Roma Tre University. The purpose of this one-day event was to present few contributions by several
More informationElectromagnetic wave propagation in complex dispersive media
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,
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 informationPolynomial 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 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 informationTrapezoidal methods for fractional differential equations: theoretical and computational aspects
Trapezoidal methods for fractional differential equations: theoretical and computational aspects Roberto Garrappa Università degli Studi di Bari, Dipartimento di Matematica Via E. Orabona n., 71 Bari,
More informationA Survey of Numerical Methods in Fractional Calculus
A Survey of Numerical Methods in Fractional Calculus Fractional Derivatives in Mechanics: State of the Art CNAM Paris, 17 November 2006 Kai Diethelm diethelm@gns-mbh.com, k.diethelm@tu-bs.de Gesellschaft
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 informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationNumerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag Leffler Functions
mathematics Article Numerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag Leffler Functions Marina Popolizio ID Dipartimento di Matematica e Fisica Ennio De Giorgi,
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 informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
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 informationMaximum principle for the fractional diusion equations and its applications
Maximum principle for the fractional diusion equations and its applications Yuri Luchko Department of Mathematics, Physics, and Chemistry Beuth Technical University of Applied Sciences Berlin Berlin, Germany
More informationThe generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations.
The generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations. Johannes Tausch Abstract An extension of the Euler-Maclaurin formula to singular integrals was introduced
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 informationAnomalous relaxation in dielectrics. Equations with fractional derivatives
Materials Science-Poland, Vol. 23, No. 4, 25 Anomalous relaxation in dielectrics. Equations with fractional derivatives V.V. NOVIKOV, K.W. WOJCIECHOWSKI 2*, O.A. KOMKOVA, T. THIEL 3 Odessa National Polytechnical
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 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 informationSIGNALING PROBLEM FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATION IN A HALF-PLANE. Yuriy Povstenko. Abstract
SIGNALING PROBLEM FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATION IN A HALF-PLANE Yuriy Povstenko Abstract The time-fractional diffusion-wave equation is considered in a half-plane. The Caputo fractional derivative
More informationOn the fractional-order logistic equation
Applied Mathematics Letters 20 (2007) 817 823 www.elsevier.com/locate/aml On the fractional-order logistic equation A.M.A. El-Sayed a, A.E.M. El-Mesiry b, H.A.A. El-Saka b, a Faculty of Science, Alexandria
More informationComputing the matrix Mittag Leffler function with applications to fractional calculus
Computing the matrix Mittag Leffler function with applications to fractional calculus arxiv:1804.04883v1 [math.na] 13 Apr 2018 Roberto Garrappa Università degli Studi di Bari, Dipartimento di Matematica,
More informationHigh Order Numerical Methods for the Riesz Derivatives and the Space Riesz Fractional Differential Equation
International Symposium on Fractional PDEs: Theory, Numerics and Applications June 3-5, 013, Salve Regina University High Order Numerical Methods for the Riesz Derivatives and the Space Riesz 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 informationBoundary layers in a two-point boundary value problem with fractional derivatives
Boundary layers in a two-point boundary value problem with fractional derivatives J.L. Gracia and M. Stynes Institute of Mathematics and Applications (IUMA) and Department of Applied Mathematics, University
More informationProperties of the Mittag-Leffler relaxation function
Journal of Mathematical Chemistry Vol. 38, No. 4, November 25 25) DOI: 1.17/s191-5-699-z Properties of the Mittag-Leffler relaxation function Mário N. Berberan-Santos Centro de Química-Física Molecular,
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationarxiv: v1 [math.ap] 26 Mar 2013
Analytic solutions of fractional differential equations by operational methods arxiv:134.156v1 [math.ap] 26 Mar 213 Roberto Garra 1 & Federico Polito 2 (1) Dipartimento di Scienze di Base e Applicate per
More informationHigher order numerical methods for solving fractional differential equations
Noname manuscript No. will be inserted by the editor Higher order numerical methods for solving fractional differential equations Yubin Yan Kamal Pal Neville J Ford Received: date / Accepted: date Abstract
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationIndia
italian journal of pure and applied mathematics n. 36 216 (819 826) 819 ANALYTIC SOLUTION FOR RLC CIRCUIT OF NON-INTGR ORDR Jignesh P. Chauhan Department of Applied Mathematics & Humanities S.V. National
More informationA Spectral Time-Domain Method for Computational Electrodynamics
A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation
More informationWe start with a discussion of the Sampling Theorem.
Math 625 Lecture Notes - March 10, 2014 We start with a discussion of the Sampling Theorem. Discussion 0.1 Assume that the Fourier transform f of a function f satisfies f(w) = 0 for w > Ω, that is, f is
More informationQuadrature Based Approximations of Non-integer Order Integrator on Finite Integration Interval
Quadrature Based Approximations of Non-integer Order Integrator on Finite Integration Interval Jerzy Baranowski Abstract Implementation of non-integer order systems is the subject of an ongoing research.
More informationFractional Derivatives and Fractional Calculus in Rheology, Meeting #3
Fractional Derivatives and Fractional Calculus in Rheology, Meeting #3 Metzler and Klafter. (22). From stretched exponential to inverse power-law: fractional dynamics, Cole Cole relaxation processes, and
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 informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
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 informationdu dθ We wish to approximate the initial value problem for the fractional ordinary differential equation of order β :
Mixed Collocation for Fractional Differential Equations François Dubois and Stéphanie Mengué December 5, 3 This contribution is dedicaced to Claude Brézinski at the occasion of his sixtieth birthday. Abstract
More informationLecture 3 Fiber Optical Communication Lecture 3, Slide 1
Lecture 3 Optical fibers as waveguides Maxwell s equations The wave equation Fiber modes Phase velocity, group velocity Dispersion Fiber Optical Communication Lecture 3, Slide 1 Maxwell s equations in
More informationFRACTIONAL CALCULUS AND WAVES IN LINEAR VISCOELASTICITY
FRACTIONAL CALCULUS AND WAVES IN LINEAR VISCOELASTICITY by Francesco Mainardi (University of Bologna, Italy) E-mail: Francesco.Mainardi@bo.infn.it Imperial College Press, London 2, pp. xx+ 347. ISBN: 978--8486-329-4-8486-329-
More informationAnomalous diffusion equations
Anomalous diffusion equations Ercília Sousa Universidade de Coimbra. p.1/34 Outline 1. What is diffusion? 2. Standard diffusion model 3. Finite difference schemes 4. Convergence of the finite difference
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 informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationUSAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION
October 4-6, 2007 - Chiinu, Rep.Moldova USAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION Ionu- P. NICA, Valeriu Gh. DAVID, /tefan URSACHE Gh. Asachi Technical University Iai, Faculty
More informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
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 informationOptical Solitons. Lisa Larrimore Physics 116
Lisa Larrimore Physics 116 Optical Solitons An optical soliton is a pulse that travels without distortion due to dispersion or other effects. They are a nonlinear phenomenon caused by self-phase modulation
More informationHOMOTOPY PERTURBATION METHOD TO FRACTIONAL BIOLOGICAL POPULATION EQUATION. 1. Introduction
Fractional Differential Calculus Volume 1, Number 1 (211), 117 124 HOMOTOPY PERTURBATION METHOD TO FRACTIONAL BIOLOGICAL POPULATION EQUATION YANQIN LIU, ZHAOLI LI AND YUEYUN ZHANG Abstract In this paper,
More informationRapid Design of Subcritical Airfoils. Prabir Daripa Department of Mathematics Texas A&M University
Rapid Design of Subcritical Airfoils Prabir Daripa Department of Mathematics Texas A&M University email: daripa@math.tamu.edu In this paper, we present a fast, efficient and accurate algorithm for rapid
More informationDevelopment of a stable coupling of the Yee scheme with linear current
Development of a stable coupling of the Yee scheme with linear current Martin Campos Pinto (LJLL), Bruno Després (LJLL) Stéphane Heuraux (IJL), Filipe Da Silva (IPFN+IST) 15 octobre 2013 Munich: Pereverzev
More informationOn The Leibniz Rule And Fractional Derivative For Differentiable And Non-Differentiable Functions
On The Leibniz Rule And Fractional Derivative For Differentiable And Non-Differentiable Functions Xiong Wang Center of Chaos and Complex Network, Department of Electronic Engineering, City University of
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 50, No. 4, pp. 849 860 c 0 Society for Industrial and Applied Mathematics TWO RESULTS CONCERNING THE STABILITY OF STAGGERED MULTISTEP METHODS MICHELLE GHRIST AND BENGT FORNBERG
More informationMODELLING OF RECIPROCAL TRANSDUCER SYSTEM ACCOUNTING FOR NONLINEAR CONSTITUTIVE RELATIONS
MODELLING OF RECIPROCAL TRANSDUCER SYSTEM ACCOUNTING FOR NONLINEAR CONSTITUTIVE RELATIONS L. X. Wang 1 M. Willatzen 1 R. V. N. Melnik 1,2 Abstract The dynamics of reciprocal transducer systems is modelled
More informationScalar electromagnetic integral equations
Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how
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 informationStochastic processes for symmetric space-time fractional diffusion
Stochastic processes for symmetric space-time fractional diffusion Gianni PAGNINI IKERBASQUE Research Fellow BCAM - Basque Center for Applied Mathematics Bilbao, Basque Country - Spain gpagnini@bcamath.org
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 informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
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 informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
More informationarxiv: v1 [math.na] 6 Nov 2013
arxiv:1311.1388v1 math.na 6 Nov 213 Exponential quadrature rules for linear fractional differential equations Roberto Garrappa Università degli Studi di Bari - Dipartimento di Matematica Via E. Orabona
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 informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationarxiv: v1 [math.ca] 28 Feb 2014
Communications in Nonlinear Science and Numerical Simulation. Vol.18. No.11. (213) 2945-2948. arxiv:142.7161v1 [math.ca] 28 Feb 214 No Violation of the Leibniz Rule. No Fractional Derivative. Vasily E.
More informationTHE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY
THE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY RUDOLF GORENFLO Free University of Berlin Germany Email:gorenflo@mi.fu.berlin.de
More information2007 Summer College on Plasma Physics
SMR/1856-10 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Dielectric relaxation and ac universality in materials with disordered structure. J. Juul Rasmussen Risø National Laboratory
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationThe solutions of time and space conformable fractional heat equations with conformable Fourier transform
Acta Univ. Sapientiae, Mathematica, 7, 2 (25) 3 4 DOI:.55/ausm-25-9 The solutions of time and space conformable fractional heat equations with conformable Fourier transform Yücel Çenesiz Department of
More informationarxiv: v1 [math.pr] 19 Nov 2017
Article A note on Hadamard fractional differential equations with varying coefficients and their applications in probability Roberto Garra 1,, Enzo Orsingher 1 and Federico Polito 2 arxiv:1711.07016v1
More informationWave propagation in discrete heterogeneous media
www.bcamath.org Wave propagation in discrete heterogeneous media Aurora MARICA marica@bcamath.org BCAM - Basque Center for Applied Mathematics Derio, Basque Country, Spain Summer school & workshop: PDEs,
More informationFRACTIONAL DIFFERENTIAL EQUATIONS
FRACTIONAL DIFFERENTIAL EQUATIONS An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of their Solution and some of their Applications by Igor Podlubny Technical University
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 informationA novel difference schemes for analyzing the fractional Navier- Stokes equations
DOI: 0.55/auom-207-005 An. Şt. Univ. Ovidius Constanţa Vol. 25(),207, 95 206 A novel difference schemes for analyzing the fractional Navier- Stokes equations Khosro Sayevand, Dumitru Baleanu, Fatemeh Sahsavand
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
More informationz k, z C, (1.1) , for k 2. j! j=0
ON THE CONVERGENCE OF KRYLOV SUBSPACE METHODS FOR MATRIX MITTAG-LEFFLER FUNCTIONS IGOR MORET AND PAOLO NOVATI Abstract. In this paper we analyze the convergence of some commonly used Krylov subspace methods
More informationNonlocal problems for the generalized Bagley-Torvik fractional differential equation
Nonlocal problems for the generalized Bagley-Torvik fractional differential equation Svatoslav Staněk Workshop on differential equations Malá Morávka, 28. 5. 212 () s 1 / 32 Overview 1) Introduction 2)
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
More informationOn the Finite Caputo and Finite Riesz Derivatives
EJTP 3, No. 1 (006) 81 95 Electronic Journal of Theoretical Physics On the Finite Caputo and Finite Riesz Derivatives A. M. A. El-Sayed 1 and M. Gaber 1 Faculty of Science University of Alexandria, Egypt
More informationSuperconvergence analysis of multistep collocation method for delay Volterra integral equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 216, pp. 25-216 Superconvergence analysis of multistep collocation method for delay Volterra integral equations
More informationRepresentations of solutions of general nonlinear ODEs in singular regions
Representations of solutions of general nonlinear ODEs in singular regions Rodica D. Costin The Ohio State University O. Costin, RDC: Invent. 2001, Kyoto 2000, [...] 1 Most analytic differential equations
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More informationProgress In Electromagnetics Research M, Vol. 13, 29 40, 2010
Progress In Electromagnetics Research M, Vol. 13, 9 40, 010 SHIFT-OPERATOR FINITE DIFFERENCE TIME DO- MAIN ANALYSIS OF CHIRAL MEDIUM A. M. Attiya Electrical Engineering Department King Saud University,
More informationJournal of Computational Physics 161, (2000) doi: /jcph , available online at
Journal of Computational Physics 161, 723 727 (2000) doi:10.1006/jcph.2000.6524, available online at http://www.idealibrary.com on NOTE Note on Nonsymmetric Finite Differences for Maxwell s Equations Tobin
More informationQuadratic nonlinear interaction
Nonlinear second order χ () interactions in III-V semiconductors 1. Generalities : III-V semiconductors & nd ordre nonlinear optics. The strategies for phase-matching 3. Photonic crystals for nd ordre
More informationChapter 10 Exercises
Chapter 10 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rl/fdmbook Exercise 10.1 (One-sided and
More informationKinetic models of Maxwell type. A brief history Part I
. A brief history Part I Department of Mathematics University of Pavia, Italy Porto Ercole, June 8-10 2008 Summer School METHODS AND MODELS OF KINETIC THEORY Outline 1 Introduction Wild result The central
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 informationA Note on Hadamard Fractional Differential Equations with Varying Coefficients and Their Applications in Probability
mathematics Article A Note on Hadamard Fractional Differential Equations with Varying Coefficients and Their Applications in Probability Roberto Garra 1, *, Enzo Orsingher 1 and Federico Polito 2 1 Dipartimento
More informationACM 106a: Lecture 1 Agenda
1 ACM 106a: Lecture 1 Agenda Introduction to numerical linear algebra Common problems First examples Inexact computation What is this course about? 2 Typical numerical linear algebra problems Systems of
More informationThe 3-wave PDEs for resonantly interac7ng triads
The 3-wave PDEs for resonantly interac7ng triads Ruth Mar7n & Harvey Segur Waves and singulari7es in incompressible fluids ICERM, April 28, 2017 What are the 3-wave equa.ons? What is a resonant triad?
More informationarxiv:cs/ v1 [cs.na] 23 Aug 2004
arxiv:cs/0408053v1 cs.na 23 Aug 2004 WEIGHTED AVERAGE FINITE DIFFERENCE METHODS FOR FRACTIONAL DIFFUSION EQUATIONS Santos B. Yuste,1 Departamento de Física, Universidad de Extremadura, E-0071 Badaoz, Spain
More informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Using Quadrature-Based Projection Methods Julian Köllermeier, Manuel Torrilhon 29th International Symposium on Rarefied Gas Dynamics, Xi an July 14th, 2014 Julian Köllermeier, Manuel Torrilhon 1 / 12 Outline
More informationA tutorial on meta-materials and THz technology
p.1/49 A tutorial on meta-materials and THz technology Thomas Feurer thomas.feurer@iap.unibe.ch Institute of Applied Physics Sidlerstr. 5, 3012 Bern Switzerland p.2/49 Outline Meta-materials Super-lenses
More informationNonlinear equations. Norms for R n. Convergence orders for iterative methods
Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector
More informationCOLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE
Progress In Electromagnetics Research M, Vol. 3, 239 252, 213 COLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE Lijuan Shi 1, 3, Lixia Yang 2, *, Hui Ma 2, and Jianning Ding 3 1 School
More information