Numerical methods for a fractional diffusion/anti-diffusion equation

Transcription

1 Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

2 Outline 1 Introduction 2 Finite difference schemes - Collaboration : Pascal Azerad (Univ. Montp2 ) 3 Splitting methods - Collaboration : Rémi Carles (CNRS & Univ. Montp2) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

3 Outline 1 Introduction 2 Finite difference schemes - Collaboration : Pascal Azerad (Univ. Montp2 ) 3 Splitting methods - Collaboration : Rémi Carles (CNRS & Univ. Montp2) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

4 Model for dune morphodynamics A conservative nonlinear model For all t (0,T) et x R, { ( ) u t (t,x)+ u 2 2 (t,x) u xx(t,x)+i[u(t, )](x) = 0, x u(0,x) = u 0 (x). Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

5 Model for dune morphodynamics A conservative nonlinear and nonlocal model (A.C Fowler, Oxford) For all t (0,T) et x R, { ) u t (t,x)+ ( u 2 2 u(0,x) = u 0 (x). x (t,x) u xx(t,x)+ + 0 ξ 1/3 u xx (t,x ξ)dξ = 0, Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

6 Model for dune morphodynamics A nonlinear and nonlocal conservative model (A.C Fowler, Oxford) For all t (0,T) and x R, { ) u t (t,x)+ ( u 2 2 u(0,x) = u 0 (x). x (t,x) u xx(t,x)+ + 0 ξ 1/3 u xx (t,x ξ)dξ = 0, Remark : This nonlocal term also appears in the work of P.-Y. Lagrée (Paris VI). References : P-Y Lagrée, Asymptotic Methods in Fluid Mechanics : Survey and Recent Advances, lecture notes 523, CISM International Centre for Mechanical Sciences Udine, H. STEINRÜCK Ed., Springer, (2010). A.C. Fowler, Mathematics and environment, lecture note, Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

7 Some expressions for I Integral formula : For all ϕ S(R) and all x R I[ϕ](x) = ϕ(x +z) ϕ(x) ϕ (x)z z 7/3 dz Fractional derivative : For causal functions (i.e. ϕ(x) = 0 for x < 0) 1 + Γ(2/3) 0 ϕ (x ξ) ξ 1/3 dξ = d 2/3 dx 2/3ϕ (x) = d4/3 dx 4/3ϕ(x) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

8 Anti-diffusive effect Pseudo-differential formula For all ϕ S(R) and all ξ R F (I[ϕ])(ξ) = (a I ±b I i) ξ 4/3 Fϕ(ξ), Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

9 Anti-diffusive effect Pseudo-differential formula For all ϕ S(R) and all ξ R F (I[ϕ])(ξ) = (a I ± b I i) ξ 4/3 Fϕ(ξ), = I ( ) λ 2 with λ = 4 3. anti-diffusive effect Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

10 Some results Existence, uniqueness and continuous dependence of the solution u w.r.t. initial datum in L 2 (R), Failure of maximum principle, Existence of travelling waves φ C 1 b (R), Global existence of L 2 perturbations of travelling waves, Instability of constant solutions. References : N. Alibaud, P. Azerad, D. Isèbe, A non-monotone conservation law for dune morphodynamics, Differential Integral Equations, B. Alvarez-Samaniego, P. Azerad, Travelling wave solutions of the Fowler equation, Discrete and Continuous Dynamical Systems, B, A.B. On the instability of a nonlocal scalar conservation law, Disc. Cont. Dyn. Syst., Ser. S Vol. 5, no 3 (2012). A.B. Global existence of solutions to the Fowler equation in a neighbourhood of travelling-waves, Int. J. Diff. Eq. (2011). Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

11 Linearized problem { vt (t,x) v xx (t,x)+i[v(t, )](x) = 0, v(0,x) = v 0 (x). = v(t,x) = K(t, ) v 0, Kernel of I 2 xx K(t, ) = F 1( e tψ I) with ψ I (ξ) = 4π 2 ξ 2 a I ξ 4/3 +b I iξ ξ 1/3 Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

12 Linearized problem { vt (t,x) v xx (t,x)+i[v(t, )](x) = 0, v(0,x) = v 0 (x). = v(t,x) = K(t, ) v 0, Properties of K(t, ),t > 0 C 0 -semi-group : Kernel of I 2 xx K(t, ) = F 1( e tψ I) with ψ I (ξ) = 4π 2 ξ 2 a I ξ 4/3 +b I iξ ξ 1/3 K(t) K(s) = K(t +s) u 0 L 2 (R), lim t 0 K(t) u 0 = u 0 Regularity K(t,x) C ((0,+ ) R) Estimates for the gradient : x K(t) L 2 Ct 3/4 x K(t) L 1 Ct 1/2 Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

13 Linearized problem { vt (t,x) v xx (t,x)+i[v(t, )](x) = 0, v(0,x) = v 0 (x). = v(t,x) = K(t, ) v 0, Kernel of I 2 xx K(t, ) = F 1( e tψ I) with ψ I (ξ) = 4π 2 ξ 2 a I ξ 4/3 +b I iξ ξ 1/3 Evolution of K for t = 0.1 and t = 0.5. = K(t, ) is not positive! Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

14 Outline 1 Introduction 2 Finite difference schemes - Collaboration : Pascal Azerad (Univ. Montp2 ) 3 Splitting methods - Collaboration : Rémi Carles (CNRS & Univ. Montp2) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

15 Explicit finite differences u n+1 j t u n j +F(u n j 1,un j,un j+1 ) un j+1 2un j +u n j 1 x 2 +I x [u n ] j = 0 Two discretizations for the nonlocal term : I x 1 [ϕ] + j = x 4/3 l 1/3 (ϕ j l+1 2ϕ j l +ϕ j l 1 ) l=1 I x 2 [ϕ] j = 4 + ( 9 x 4/3 l 7/3 ϕ j l ϕ j + ϕ ) j+1 ϕ j 1 l 2 l=1 Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

16 Linearization For v > 0, we have u n+1 j t u n j +v un j u n j 1 x ǫ un j+1 2un j +u n j 1 x 2 +ηi x [u n ] j = 0 Numerical scheme reads for I 1 x u n+1 j = ǫ t + ( x 2 un j v t x 2ǫ t x 2 η t ( v t x + ǫ t x 2 +(2 2 1/3 ) η t η t x 4/3 + l=2 x 4/3 x 4/3 ) u n j 1 [ (l +1) 1/3 2l 1/3 +(l 1) 1/3] u n j l. ) u n j Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

17 Linearization For v > 0, we have u n+1 j t u n j +v un j u n j 1 x ǫ un j+1 2un j +u n j 1 x 2 +ηi x [u n ] j = 0 Numerical scheme reads for I 1 x u n+1 j Note : u n+1 j = ǫ t + ( x 2 un j v t x 2ǫ t x 2 η t ( v t x + ǫ t x 2 +(2 2 1/3 ) η t η t x 4/3 + l=2 x 4/3 x 4/3 ) u n j 1 [ (l +1) 1/3 2l 1/3 +(l 1) 1/3] u n j l. is not a convex combination of (u n j ) j N ) u n j Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

18 Amplification factor u(t,x) = e ikx+σt is solution of iff dispersion relationship u t +vu x ǫu xx +ηi[u] = 0, σ +iµk +ǫk 2 ηk 4/3 Γ Exact amplification factor ( ) ( i ) 3 = 0 2 ( G cont = e t( ǫk2 +ηk 4/3 1 2 Γ(2 i t µk+ηk 3)) e 4/3 3 2 Γ(2 3) ) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

19 Amplification factor u(t,x) = e ikx+σt is solution of iff dispersion relationship u t +vu x ǫu xx +ηi[u] = 0, σ +iµk +ǫk 2 ηk 4/3 Γ Exact amplification factor ( ) ( i ) 3 = 0 2 ( G cont = e t( ǫk2 +ηk 4/3 1 2 Γ(2 i t µk+ηk 3)) e 4/3 3 2 Γ(2 3) ) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

20 Strategy High frequencies are responsible of numerical instabilities High frequencies are not amplified in Fowler s model Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

21 Strategy High frequencies are responsible of numerical instabilities High frequencies are not amplified in Fowler s model + Definition (High Frequency stability) We say that a numerical scheme is HF-stable if the high frequencies are strongly stable that is to say : 0 < θ 0 < π such that (θ 0,π], g( x, t,θ) < 1, where g is the discrete amplification factor. = Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

22 High frequency stability Definition (High Frequency stability) We say that a numerical scheme is HF-stable if the high frequencies are strongly stable that is to say : 0 < θ 0 < π such that (θ 0,π], g( x, t,θ) < 1, where g is the discrete amplification factor. Von-Neumann method : Input : one single Fourier mode u n j = û n k eikx j Output : û n+1 k = g( x, t,k)û n k Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

23 Stability criteria Proposition A scheme is HF-stable if x and t satisfy the following conditions : For I x 1. v t x + 2ǫ t x 2 +(2 2 1/3 )η t < 1, x4/3 (1 2 1/3 ) η t 2ǫ t x4/3 x 2. For I 2 x. v t x 4 9 where ζ is the Riemann zeta function. t + 2ǫ x ( ζ( 4 ) η t 9 3 ) 1 < 1, x4/3 ( ζ( 7 ) η t 3 ) 1+ζ(4 3 ) 2ǫ t x4/3 x 2. Azerad P., B. A., Finite difference approximations for a diffusion/anti-diffusion equation, submitted Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

24 Numerical experiments Modified CFL µ t x +2ǫ t x 2+(2 2 1/3 )η t x 4/3 < 1 θ = kδx, CFL mod max θ g(θ,δ t) δ t δt max = 0.042,CFL mod 0.99 CFL mod 1.22 Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

25 Numerical experiments Modified CFL For I x 1 : v t x +2ǫ t +(2 2 1/3 )η t < 1 x 2 x 4/3 VS. For I 2 x : v t x +2ǫ t x ( ζ( 4 9 ) 1) η t x 4/3 < 1 Amplification factors for I x 1 (blue line) and I2 x (dashed line) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

26 Outline 1 Introduction 2 Finite difference schemes - Collaboration : Pascal Azerad (Univ. Montp2 ) 3 Splitting methods - Collaboration : Rémi Carles (CNRS & Univ. Montp2) Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

27 Splitting method { ( ) u u t (t,x)+ 2 2 (t,x) u xx(t,x)+i[u(t, )](x) = 0, x u(0,x) = u 0 (x). Notation : u(t, ) = S t u 0 { ( v t + v 2 2 ւ ) ǫv xx = 0 x v(0, ) = v 0, Notation : v(t, ) = Y t v 0 { ց wt +I[w] ηw xx = 0 w(0, ) = w 0, Notation : w(t, ) = X t w 0 Lie method : Z t L = Xt Y t Finite difference method FFT Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

28 Splitting method { ( ) u u t (t,x)+ 2 2 (t,x) u xx(t,x)+i[u(t, )](x) = 0, x u(0,x) = u 0 (x). Notation : u(t, ) = S t u 0 { ( v t + v 2 2 ւ ) ǫv xx = 0 x v(0, ) = v 0, Notation : v(t, ) = Y t v 0 { ց wt +I[w] ηw xx = 0 w(0, ) = w 0, Notation : w(t, ) = X t w 0 Lie method : Z t L = Xt Y t Note : H. Holden, C. Lubich, N.-H. Risebro; Operator splitting for partial differential equations with Burgers nonlinearity, to appear, Math. Comp (2012). Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

29 Splitting method { ( ) u u t (t,x)+ 2 2 (t,x) u xx(t,x)+i[u(t, )](x) = 0, x u(0,x) = u 0 (x). Notation : u(t, ) = S t u 0 { ( v t + v 2 2 ւ ) ǫv xx = 0 x v(0, ) = v 0, Notation : v(t, ) = Y t v 0 { ց wt +I[w] ηw xx = 0 w(0, ) = w 0, Notation : w(t, ) = X t w 0 Lie method : Z t L = Xt Y t Finite difference method FFT Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

30 Expressions for the flows Linear flow X t v 0 = D(t, ) v 0, where D(t, ) = F 1( ) e t φi with φi (ξ) = 4π 2 ηξ 2 a I ξ 4/3 +b I ξ ξ 1/3 Nonlinear flow (viscous Burgers equation) Y t w 0 = G(t, ) w where G is the heat kernel. Exact flow S t u 0 = K(t, ) u Splitting operator Z t Lu 0 = K(t) u t 0 t 0 t 0 x G(t s, ) (Y s w 0 ) 2 ds, x K(t s, ) (S s u 0 ) 2 ds D(t) G(t s) x (Y s u 0 ) 2 ds Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

31 L 2 local error estimate Proposition Let u 0 H 3 (R). There exists C ( u 0 L (R)) 2 such that for all t [0,1], ZL t u 0 S t u 0 L 2 (R) C ( u 0 L (R)) 2 t 2 u 0 2 H 3 (R). Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

32 L 2 local error estimate Proposition Let u 0 H 3 (R). There exists C ( u 0 L (R)) 2 such that for all t [0,1], ZL t u 0 S t u 0 L 2 (R) C ( u 0 L (R)) 2 t 2 u 0 2 H 3 (R). Ingredients of the proof t ZL tu 0 S t u 0 = xk(t s) ( (S s u 0 ) 2 (ZL su 0) 2) ds +R(t), The remainder R(t) is written as R(t) = 1 t 2 0 R 1(s)ds, with R 1 (s) = x K(t s) (ZL s u 0) 2 D(t) x G(t s, ) (Y s u 0 ) 2, and satisfies : R(t) L 2 (R) C ( u 0 L (R)) 2 t 2 u 0 2 H 3 (R). Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

33 Modified fractional Gronwall Lemma Lemma Let φ : [0,T] R + be a bounded measurable function and P be a polynomial with positive coefficients and no constant term. We assume there exists two positive constants C and θ ]0,1[ such that for all t [0,T], 0 φ(t) φ(0)+p(t)+c d θ dt θφ(t). Then there exists C T (θ) such that for all t [0,T], Riemann-Liouville operator : φ(t) C T (θ)φ(0)+c T (θ)p(t). d θ dt θφ(t) = 1 Γ(θ) t 0 (t s) θ 1 φ(s)ds. Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

34 We will also need the fact that the flow map S t is uniformly Lipschitzean on balls of H 2 (R). Proposition Let T,R > 0. There exists K = K(R,T) < such that if u 0 H 2 (R) R, v 0 H 2 (R) R, then S t u 0 S t v 0 L 2 (R) K u 0 v 0 L 2 (R), t [0,T]. Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

35 Convergence Theorem For all u 0 H 3 (R) and for all T > 0, there exist positive constants c 1,c 2 and t 0 such that for all t ]0, t 0 ] and for all n N such that 0 n t T, (Z t L )n u 0 S n t u 0 L 2 (R) c 1 t and (Z t L )n u 0 H 3 (R) c 2. Here, c 1,c 2 and t 0 depend only on T, ρ = max t [0,T] S t u 0 H 2 (R), and u 0 H 3 (R). B.A., Carles R., Splitting methods for the nonlocal Fowler equation, Math. Comp (2012), to appear. Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

36 Proof The proof follows the same idea as in [1,2] for instance. We prove by induction that there exists γ, t 0 such that if 0 < t t 0, for all n N with n t T, (Z t L )n u 0 L 2 (R) 2ρ, (Z t L )n u 0 S n t u 0 L 2 (R) γ t. The triangle inequality yields (Z t L ) n u 0 S n t u 0 L 2 with u k = (Z t L )k n 1 S (n k 1) t ( ZL t ) u k S (n k 1) t ( S t ) u k L, 2 k=0 References : [1] C. Besse, B. Bidegaray, S. Descombes; Order estimates in time of splitting methods for the nonlinear Schrödinger equation, SIAM J. Numer. Anal., 40 (2002). [2] H. Holden, C. Lubich, N.-H. Risebro; Operator splitting for partial differential equations with Burgers nonlinearity, to appear, Math. Comp (2012). Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

37 Proof Lipschitz property of S t yields S (n k 1) t ( ZL t ) u k S (n k 1) t ( S t ) u k L K 2 Z t L u k S t u L k 2 From L 2 local error estimate, we infer S (n k 1) t ( Z t L u k ) S (n k 1) t ( S t u k ) L 2 CK( t)2 u 0 2 H 3, for some constant C. Therefore, (ZL t ) n u 0 S n t u 0 L 2 nck( t) 2 u 0 2 H3 CTK t, which yields the two estimates of the induction, provided one takes γ = CTK, which is uniform in n and t. Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

38 Numerical experiments : initial data Figure: Initial data used for numerical experiments. Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

39 Numerical convergence for Lie operator Figure: Lie method Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

40 Strang method Strang operator Z t S = Xt/2 Y t X t/2 Figure: Strang method Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30

41 Merci, Thanks, Danke! Afaf Bouharguane Numerical schemes for the Fowler equation Berlin, November / 30