NUMERICAL ANALYSIS OF NONLINEAR SUBDIFFUSION EQUATIONS
|
|
- Clifford Garrett
- 6 years ago
- Views:
Transcription
1 NUMERICAL ANALYSIS OF NONLINEAR SUBDIFFUSION EQUATIONS BANGTI JIN, BUYANG LI, AND ZHI ZHOU Abstract. We present a general framework for the rigorous numerical analysis of time-fractional nonlinear parabolic partial differential equations, with a fractional derivative of order α (, 1) in time. It relies on three technical tools: a fractional version of the discrete Grönwall-type inequality, discrete maximal regularity, and regularity theory of nonlinear equations. We establish a general criterion for showing the fractional discrete Grönwall inequality, and verify it for the L1 scheme and convolution quadrature generated by BDFs. Further, we provide a complete solution theory, e.g., existence, uniqueness and regularity, for a time-fractional diffusion equation with a Lipschitz nonlinear source term. Together with the known results of discrete maximal regularity, we derive pointwise L 2 (Ω) norm error estimates for semidiscrete Galerkin finite element solutions and fully discrete solutions, which are of order O(h 2 ) (up to a logarithmic factor) and O(τ α ), respectively, without any extra regularity assumption on the solution or compatibility condition on the problem data. The sharpness of the convergence rates is supported by the numerical experiments. Keywords: nonlinear fractional diffusion equation, discrete fractional Grönwall inequality, L1 scheme, convolution quadrature, error estimate 1. Introduction. Time-fractional parabolic partial differential equations (PDEs) have been very popular for modeling anomalously slow transport processes in the past two decades. These models are commonly referred to as fractional diffusion or subdiffusion. At a microscopic level, the underlying stochastic process is continuous time random walk [32]. So far they have been successfully applied in a broad range of diversified research areas, e.g., thermal diffusion in fractal domains [35], flow in highly heterogeneous aquifer [6] and single-molecular protein dynamics [2], just to name a few. Hence, the rigorous numerical analysis of such problems is of great practical importance. For the linear problem, various efficient time stepping schemes have been proposed, which include mainly two classes: L1 type schemes and convolution quadrature (CQ). L1 type schemes approximate the fractional derivative by replacing the integrand with its piecewise polynomial interpolation [24, 26, 37, 3] and thus generalize the classical finite difference method. The piecewise linear case has a local truncation error O(τ 2 α ) for sufficiently smooth solution, where τ denotes the time step size. See also [31, 33] for the discontinuous Galerkin method. CQ is a flexible framework introduced by Lubich [27, 28] for constructing high-order time discretization methods for approximating fractional derivatives. It approximates the fractional derivative in the Laplace domain and automatically inherits the stability property of general linear multistep methods. See [1, 39, 4, 16] for CQ type schemes. Optimal error estimates have been derived for both spatially semidiscrete and fully discrete schemes, including problems with nonsmooth data [1, 14, 31, 16]. However, up to now, there has been very few work on the rigorous numerical analysis of nonlinear time fractional diffusion equations. In this paper, we present a general framework for analyzing discretization errors of nonlinear problems. The error of the numerical solution can be split into a linear part and a nonlinear part. While the linear part has been carefully The work of B. Jin is partially supported by UK EPSRC grant EP/M2516/1. The work of B. Li is partially supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region (Project No ). The work of Z. Zhou is partially supported by the AFOSR MURI center for Material Failure Prediction through peridynamics and the ARO MURI Grant W911NF Department of Computer Science, University College London, Gower Street, London, WC1E 6BT, UK (b.jin@ucl.ac.uk, bangti.jin@gmail.com) Department of Applied Mathematics, The Hong Kong Polytechnic University, Kowloon, Hong Kong. (buyang.li@polyu.edu.hk, libuyang@gmail.com) Department of Applied Mathematics, The Hong Kong Polytechnic University, Kowloon, Hong Kong. (zhizhou@polyu.edu.hk, zhizhou125@gmail.com) 1
2 studied, the analysis of the nonlinear part requires different mathematical machineries, in order to derive sharp error estimates. Besides regularity estimates for the nonlinear problem, it requires discrete maximal l p regularity, and a fractional version of the discrete Grönwall s inequality for time stepping schemes. The former gives a bound on the discrete fractional derivative due to the nonlinear part, whereas the latter allows combining the nonlinear part with the linear part to obtain a global error estimate. To the best of our knowledge, a fractional version of discrete Grönwall s inequality for time stepping schemes is still unavailable in the literature. We shall establish such discrete Grönwall s inequality for both L1 scheme and CQs generated by backward difference formulas (BDFs) up to order 6 in Theorem 2.6. Further, in Theorem 2.5, we present a general criterion under which the fractional discrete Grönwall s inequality holds. To illustrate the main idea of this framework, we consider the following nonlinear problem in a bounded convex polygonal domain Ω R d, d 1: t α u u = f(u) in Ω (, T ), (1.1) u = on Ω (, T ), u = u in Ω {}, where u H 1 (Ω) H 2 (Ω) is a given function and f : R R is a Lipschitz continuous function, i.e., f(s) f(t) L s t for all s, t R, and α t u denotes the Caputo fractional derivative of order α (, 1) in time [19, pp. 91] (1.2) α t u(t) := 1 Γ(1 α) (t s) α d ds u(s) ds, with Γ(z) := s z 1 e s ds. Let S h H 1 (Ω) be the continuous piecewise linear finite element space subject to a quasi-uniform shape regular triangulation of Ω, with a mesh size h, and let h : S h S h denote the Galerkin finite element approximation of the Dirichlet Laplacian, defined by ( h w h, v h ) := ( w h, v h ), w h, v h S h. Let = t < t 1 <... < t N = T be a uniform partition of the time interval [, T ], with grid points t n = nτ and step size τ = T/N. Upon rewriting the Caputo derivative t α u as a Riemann-Liouville one [19, pp. 91], we consider a linearized time-stepping scheme: for the given initial value u h = R hu (Ritz projection of u ), find u n h, n = 1, 2,..., N, such that (1.3) α τ (u n h u h) h u n h = P h f(u n 1 h ), where P h denotes the L 2 projection onto the finite element space S h, and α τ u n h denotes either the CQ generated by the backward Euler method or L1 scheme; see (2.8) and (2.9) below. These methods are popular for discretizing the fractional derivative in time. After proving the fractional discrete Grönwall s inequality in Section 2 and the regularity estimate in Section 3, we present an error analysis for the fully discrete scheme (1.3) in Section 4. By introducing an intermediate spatially semidiscrete Galerkin problem (1.4) α t u h (t) h u h (t) = P h f(u h (t)) t (, T ], we split the error into two parts: u(t n ) u n h = (u(t n) u h (t n )) + (u h (t n ) u n h ), and derive the following error estimates for each component in Theorems 4.3 and 4.4: max u(t) u h(t) L t T 2 (Ω) cl 2 hh 2 and max u h(t n ) u n h L 1 n N 2 (Ω) cτ α, 2
3 where l h = log(2 + 1/h). These estimates are sharp with respect to the regularity of the solution in Theorem 3.1 (up to a logarithmic factor l h ), and are confirmed by the numerical experiments in Section 6. Besides, we show how to simplify the analysis of nonlinear problems by applying the fractional-type discrete maximal l p -regularity established in [17], an extension of the discrete maximal l p -regularity of standard parabolic equations [18, 21, 25], which has been applied to numerical analysis of nonlinear parabolic equations in the literature [1, 2, 22]. Last we mention the interesting works [1, 34] on integro-differential equations, where a Riemann-Liouville fractional integral operator appears in front of the Laplacian. These models are closely related to (1.1), but have different smoothing properties. Cuesta et al [1] proposed the CQ generated by the second-order BDF for a semilinear problem, and proved an O(τ 2 ) error bound of the temporal error. In [34], a Crank-Nicolson type method for a semilinear problem with variable time step size was studied. In these works, a variant of the discrete Grönwall s inequality due to Chen et al [8] plays a crucial role, which differs substantially from the discrete Grönwall s inequality we shall establish below. Throughout this paper, the notation c denotes a generic constant, which may vary at different occurrences, but it is always independent of the mesh size h and time step size τ. 2. Discrete Grönwall s inequality for time-fractional diffusion. In this section, we establish a fractional version of Grönwall s inequality and its discrete analogue for time stepping schemes. These inequalities are crucial in analyzing numerical schemes for nonlinear subdiffusion equations, and are of independent interest Continuous Grönwall s inequality. We begin with the continuous Grönwall s inequality for fractional differential equations in a general Banach space setting. Theorem 2.1 (Fractional Grönwall s inequality). Let X be any given Banach space. For α (, 1) and p (1/α, ), if a function u C([, T ]; X) satisfies α t u L p (, T ; X), u() = and (2.1) α t u L p (,s;x) κ u L p (,s;x) + σ, s (, T ], for some positive constants κ and σ, then (2.2) u C([,T ];X) + α t u L p (,T ;X) cσ, where the constant c is independent of σ, u and X, but may depend on α, p, κ and T. Proof. Due to the zero initial condition u() =, the Riemann Liouville and Caputo fractional derivatives coincide. Hence, the function u(t) can be expressed in terms of t α u (cf. [19, pp. 96, Lemma 2.22] and [19, pp. 74, Lemma 2.5]): u(t) = 1 t Γ(α) (t ξ)α 1 ξ α u(ξ) dξ. Since p > 1/α, Hölder s inequality implies (2.3) ( ) p 1 u(t) X c (t ξ) (α 1)p p p 1 dξ α ξ u L p (,t;x) c ξ α u L p (,t;x). Upon taking the supremum with respect to t (, s) for any s (, T ] in (2.3), we obtain u L (,s;x) c α ξ u L p (,s;x) cκ u L p (,s;x) + cσ ɛκ u L (,s;x) + c ɛ κ u L 1 (,s;x) + cσ, s [, T ], where ɛ > can be arbitrary. By choosing ɛ = 1 2κ, the L -norm on the right-hand side can be eliminated by the left-hand side, and the last inequality reduces to u L (,s;x) c κ u L1 (,s;x) + cσ, s [, T ]. 3
4 That is, we have u(s) X c κ s u(ξ) Xdξ+cσ for s (, T ]. Now the standard Grönwall s inequality yields max s [,T ] u(s) X e cκt cσ. Substituting it into (2.1) yields (2.2). The proof of Theorem 2.1 is complete Discrete Grönwall s inequality. In this part, we establish the discrete analogue of the Grönwall s inequality in Theorem 2.1 for time stepping schemes that approximate the fractional derivative t α v(t n ) by a discrete convolution: (2.4) α τ v n := 1 τ α n K n j v j, n =, 1, 2,... j= where v n is an approximation of v(t n ), and K j, j =, 1, 2,..., are the weights independent of the time step size τ. Throughout, we denote by K(ζ) the generating function of the discrete fractional derivative τ α, defined by (2.5) K(ζ) := 1 τ α K j ζ j, which is an analytic function in the (open) unit disk D := {z C : z < 1}, continuously differentiable up to the boundary D\{±1}, except for the two points ±1. Then we have (2.6) K(ζ) v n ζ n = ( τ α v n )ζ n. n= Example 2.1. The CQ generated by the k th -order BDF [27, 1] is given by (2.4), where the coefficients K j, j =, 1,..., are determined by the power series expansion ( k ) α 1 (2.7) (1 ζ)j = K j ζ j. j j=1 The special case k = 1, i.e., the backward Euler CQ, is very popular and commonly known as Grünwald Letnikov approximation, and the coefficients K j, j =, 1, 2,..., are given by (2.8) (1 ζ) α = K j ζ j. Example 2.2. The popular L1 scheme [26] is also of the form (2.4) with [17, pp. 8] (1 ζ) 2 (2.9) ζγ(2 α) Li α 1(ζ) = K j ζ j, where Li p (z) = j=1 zj /j p is the polylogarithmic function, which is well defined for z < 1 and can be analytically continued to the split complex plane C \ [1, ) [11]. Now we turn to the discrete Grönwall s inequality. For 1 p, we denote by l p (X) the space of sequences v n X, n =, 1,..., such that (v n ) n= lp (X) <, where j= n= j= ( τ v n p (v n ) X n= l p (X) := n= j= j= ) 1 p if 1 p <, sup v n X if p =. n 4
5 For a finite sequence v n X, n =, 1,..., m, we denote (v n ) m n= l p (X) := (v n ) n= l p (X), by setting v n = for n > m. The following theorem is a discrete analogue of Theorem 2.1 for the backward Euler CQ. It is foundational to the proof of the discrete Grönwall s inequalities for other time-stepping schemes. Theorem 2.2 (Discrete fractional Grönwall s inequality: backward Euler). Let X be any given Banach space, and let α τ denote the backward Euler CQ given by (2.4) and (2.8). If α (, 1) and p (1/α, ), and a sequence v n X, n =, 1, 2,..., with v =, satisfies (2.1) ( α τ v n ) m n=1 l p (X) κ (v n ) m n=1 l p (X) + σ, m N, for some positive constants κ and σ, then there exists a τ > such that for any τ < τ there holds (2.11) (v n ) N n=1 l (X) + ( α τ v n ) N n=1 l p (X) cσ, where the constants c and τ are independent of σ, τ, N, X and v n, but may depend on α, p, κ and T. To prove Theorem 2.2, we need a technical lemma, which gives a discrete analogue of the Hardy type inequality (2.3). Lemma 2.3 (Discrete Hardy type inequality). Let α (, 1), and X be any given Banach space. If v n X and w n X, n =, 1, 2,..., satisfy ( ) α 1 ζ (2.12) v n ζ n = w n ζ n, τ in the sense that both sides are analytic in D, then for p (1/α, ), there holds (2.13) n= n= (v n ) m n= l (X) c (w n ) m n= lp (X), m N, where the constant c is independent of τ, m, N and X, but may depend on α, p and T. Proof. We define φ n, n =, 1,..., to be the coefficients of the power series expansion (1 ζ) α = φ n ζ n. n= Then direct calculations yield φ = 1 and φ n = n j=1 inequality ln(1 + x) x for x > 1, we have n ( ln φ n = ln 1 + α 1 ) n (α 1) j j=1 j=1 ( 1 + α 1 j ) for n 1. By the trivial j 1 (α 1) ln(n + 1). That is, φ n (n + 1) α 1 for n. It follows from (2.12) that ( ) α τ v n ζ n = w n ζ n = τ α( φ n ζ n)( w n ζ n). 1 ζ n= With p = p p 1, the last identity yields n ( X v n X = τ n ) 1 ( α φ n j w j τ α φ n j p p n w j p X (2.14) j= n= j= τ α 1/p ( n 5 j= n= j= n= ) 1 p ) 1 1 p (w j ) n (j + 1) p (1 α) j= lp (X).
6 If p > 1/α, then < p (1 α) < 1 and so n j= Hence, (2.14) reduces to 1 n+1 (j + 1) p (1 α) ds s p (1 α) (n + 1)1 p (1 α) = 1 p. (1 α) v n X τ α 1/p (n + 1) α 1/p (w j ) n (2T ) α 1/p (1 p (1 α)) j= l 1/p p (X) (w j ) n (1 p (1 α)) j= l 1/p p (X), where we have used the fact τ(n + 1) 2T in the last inequality. Since the last inequality holds for all n =,..., m, it follows that (2.13) holds. Now we are ready to prove Theorem 2.2. Proof of Theorem 2.2. For the backward Euler CQ we have K(ζ) = ( ) 1 ζ α. τ Since, v =, τ α v =, and the identity (2.6) can be written as ( ) 1 ζ α τ n= vn ζ n = n= ( τ α v n )ζ n. Then Lemma 2.3 and (2.1) imply (v n ) m n= l (X) c ( α τ v n ) m n= l p (X) cκ (v n ) m n= l p (X) + cσ ɛκ (v n ) m n= l (X) + c ɛ κ (v n ) m n= l 1 (X) + cσ, 1 m N. By choosing ɛκ = 1/2 and collecting terms, and using the fact v =, we obtain (v n ) m n=1 l (X) c κ (v n ) m n=1 l1 (X) + cσ, 1 m N. That is, v m X c κ τ m n=1 vn X + cσ for 1 m N. Then the standard discrete Grönwall s inequality gives, for sufficiently small step size τ, max 1 n N vn X e cκt cσ. Substituting this into (2.1) yields (2.11). The proof of Theorem 2.2 is complete. To analyze other time-stepping schemes, we shall need the following lemma of discrete Mikhlin multipliers, which is a simple consequence of Blunck s multiplier theorem [7, Theorem 1.3] through the transform ζ = e iθ. Here, a UMD space X denotes a Banach space such f(s) t s ds is bounded on Lp (R; X) for all 1 < p < that the Hilbert transform Hf(t) := R [23]. Examples of UMD spaces include R d, d 1, and L q (Ω), 1 < q <, and their closed subspaces (e.g. the finite element space S h equipped with the L q (Ω) norm). Lemma 2.4 (Discrete Mikhlin multipliers). Let X be a UMD space and let M : D C be an analytic function, continuously differentiable up to D\{±1}, such that the set { } { M(ζ) : ζ D\{±1} (1 ζ)(1 + ζ)m (ζ) : ζ D\{±1} } is bounded, and denote its bound by c R. Then for any 1 < p < and any sequence (f n ) n= l p (X), the coefficients u n X, n =, 1,..., in the power series expansion M(ζ) f n ζ n = u n ζ n, ζ D, n= n= satisfy (u n ) n= l p (X) c p,x c R (f n ) n= l p (X), 6
7 where the constant c p,x is independent of the operators M(ζ), ζ D. Now other time-stepping schemes can be connected to the backward Euler CQ. The next result gives a general criterion for the discrete fractional Grönwall s inequality. Theorem 2.5 (General criterion for discrete fractional Grönwall s inequality). Let X be a UMD space. If the generating function K(ζ) = 1 τ α n= K nζ n satisfies (2.15) K(ζ) 1 c 1 ζ τ α and (1 ζ)(1 + ζ)k (ζ) c K(ζ), ζ D\{±1}, then the discrete fractional Grönwall s inequality holds: if α (, 1) and p (1/α, ), and a sequence v n X, n =, 1, 2,..., with v =, satisfies (2.16) ( α τ v n ) m n=1 l p (X) κ (v n ) m n=1 l p (X) + σ, 1 m N, for some positive constants κ and σ, then there exists a τ > such that for any τ < τ there holds (2.17) (v n ) N n=1 l (X) + ( α τ v n ) N n=1 l p (X) cσ. where the constants c and τ are independent of σ, τ, N and v n, but may depend on α, p, κ, X and T. Proof. First, we note that τ α v n = τ α n j= K jv n j, n =, 1, 2,..., are the coefficients in the power series expansion (2.18) K(ζ) v n ζ n = ( τ α v n )ζ n, n= n= it follows that (2.19) v n ζ n = n= ( τ ) α [ 1 1 ζ K(ζ) ( ) α ] 1 ζ ( τ τ α v n )ζ n = n= where F n, n =, 1,..., are the coefficients in the expansion [ ( ) α ] 1 1 ζ ( K(ζ) τ τ α v n )ζ n = F n ζ n. n= By applying Lemma 2.3 to (2.19), we obtain (2.2) (v n ) m n= l (X) c (F n ) m n= l p (X), 1 m N. n= ( ) α τ F n ζ n, 1 ζ n= Let m be fixed and define Ẽn = α τ v n if n m and Ẽn = if n > m. Let F n be the coefficients of the power series (2.21) F n ζ n = n= [ 1 K(ζ) ( ) α ] 1 ζ Ẽ n ζ n, τ n= then F n = F n for n m. Now the conditions in (2.15) imply ( ) α 1 1 ζ c and K(ζ) τ (1 ζ)(1 + ζ) d [ ( ) α ] 1 1 ζ c, ζ D \ {±1}. dζ K(ζ) τ 7
8 ( By choosing M(ζ) = 1 1 ζ ) α K(ζ) τ and applying Lemma 2.4 to equation (2.21), we obtain which further implies ( F n ) n= lp (X) c (Ẽn ) n= lp (X), (F n ) m n= lp (X) = ( F n ) m n= lp (X) c (Ẽn ) n= lp (X) = c ( α τ v n ) m n= lp (X), where the constant c is independent of m. The last inequality and (2.2) yield (v n ) m n= l (X) c ( α τ v n ) m n= lp (X). Substituting (2.16) into the last inequality gives (2.22) (v n ) m n=1 l (X) cκ (v n ) m n=1 l p (X) + cσ ɛκ (v n ) m n=1 l (X) + c ɛ κ (v n ) m n=1 l 1 (X) + cσ, 1 m N. where ɛ > is arbitrary. By choosing ɛκ = 1/2, we obtain (v n ) m n=1 l (X) c κ (v n ) m n=1 l 1 (X) + cσ, 1 m N. That is, v m X c κ τ m n=1 vn X + cσ for 1 m N. Then the standard discrete Grönwall s inequality gives, for sufficiently small step size τ, max 1 n N vn X e cκt cσ. This together with (2.16) and (2.22) yields (2.17). The proof of Theorem 2.5 is complete. By Theorem 2.5, the discrete fractional Grönwall s inequality can be proved for the L1 scheme and general BDF CQs. Theorem 2.6 (Discrete Grönwall s inequality for L1 scheme and BDF CQ). Let X be a UMD space. For both L1 scheme and CQ generated by the k th -order BDF, with 1 k 6, the discrete fractional Grönwall s inequality holds: if α (, 1) and p (1/α, ), and a sequence v n X, n =, 1, 2,..., with v =, satisfies ( α τ v n ) m n=1 l p (X) κ (v n ) m n=1 l p (X) + σ, 1 m N, for some positive constants κ and σ, then there exists a τ > such that for any τ < τ there holds (v n ) N n=1 l (X) + ( α τ v n ) N n=1 l p (X) cσ, where the constants c and τ are independent of σ, τ, N and v n, but may depend on α, p, κ, X and T. Proof. By Theorem 2.5, it suffices to show that the generating functions K(ζ) of the L1 scheme and CQ satisfy (2.15). We discuss them separately. First, for the L1 scheme, 1 (1 ζ) K(ζ) = 2 Γ(2 α)τ α ζ Li α 1 (ζ) converges for ζ D\{1} and has the following asymptotic expansion (cf. [11, Theorem 1], or [17, equation (4.6)]) If ζ D\{1} is sufficiently close to 1, then τ α K(ζ) = (1 ζ) α + o((1 ζ) α ), as ζ 1. τ α K(ζ) ζ α. 8
9 Meanwhile, we recall the following series expansion (cf. [17, equation (4.5)]) Li α 1 (e iθ ( ) = (2π) α 2 cos ( (2 α)π ) (Aθ + B θ ) i sin ( ) (2 α)π ) (Aθ B θ ), Γ(2 α) 2 2 ( ) k + θ α 2 and Bθ = ( k + 1 θ 2π ) α 2. Thus, if ζ = e iθ is away where A θ = k= 2π k= from 1, then θ is away from and 2π, and thus A θ + B θ c. This shows Li α 1 (e iθ ) > c. Since 1 ζ 2 c 1 ζ α when ζ = e iθ is away from 1, it follows that τ α K(ζ) = Li α 1(ζ) Γ(2 α) 1 ζ 2 c 1 ζ 2 c 1 ζ α. Overall, the first inequality of (2.15) holds for the generating function K(ζ) of the L1 scheme. The second inequality of (2.15) has been proved in [17, Lemma 4.3]. This shows the assertion for the L1 scheme. Next we turn to the CQ. For the CQ generated by the k th -order BDF, the generating function K(ζ) satisfies ( ) α ( k ) α τ 1 K(ζ) = (1 ζ)j 1. 1 ζ j j=1 Since the function k j=1 1 j (1 ζ)j 1 has no root on the unit circle D for 1 k 6 (see [9, Proof of Lemma 2] or [12, pp ]), it follows that ( ) α τ k K(ζ) 1 ζ = 1 (1 ζ)j 1 j j=1 This proves the first inequality of (2.15). Note that (1 + ζ)(1 ζ)k (ζ) = (1 + ζ)(1 ζ) α ( k τ α j=1 ( k = α (1 + ζ)(1 ζ)α τ α j=1 α c. ) α 1 k 1 (1 ζ)j (1 ζ) j 1 j j=1 ) α 1 k 1 1 (1 ζ)j 1 j j= (1 ζ) j, and so for any ζ D\{±1}, there holds ( (1 + ζ)(1 ζ)k k (ζ) K(ζ) = α(1 + ζ) j=1 1 j (1 ζ)j 1) α 1 k 1 j= (1 ζ)j ( k j=1 1 j (1 ζ)j 1) α c. where the last inequality holds, since the denominator ( k j=1 1 j (1 ζ)j 1) α has no root on D. This shows the second part of (2.15), completing the proof of the theorem. Remark 2.1. In Theorems 2.5 and 2.6, if we assume ( α τ v n ) m n=1 l p (X) κ (v n ) m 1 n=1 l p (X) + σ, 1 m N, i.e., the index on the right-hand side is slightly changed, then we have (v n ) N n=1 l (X) + ( α τ v n ) N n=1 l p (X) cσ, without any restriction on the step size τ. 9
10 3. Regularity of the solution. Now we discuss the existence, uniqueness and regularity for the solutions to (1.1) and (1.4). These results are needed in the numerical analysis in Section 4. The main result of this section is the following theorem. Theorem 3.1. Let u H 1 (Ω) H 2 (Ω), and let f : R R be Lipschitz continuous. Then problem (1.1) has a unique solution u such that (3.1) (3.2) u C α ([, T ]; L 2 (Ω)) C([, T ]; H 1 (Ω) H 2 (Ω)), t u(t) L 2 (Ω) and t u(t) L 2 (Ω) ct α 1 for t (, T ]. α t u C([, T ]; L 2 (Ω)), Similarly, problem (1.4) has a unique solution u h such that (3.3) (3.4) u h C α ([,T ];L 2 (Ω)) + h u h C([,T ];L 2 (Ω)) + α t u h C([,T ];L 2 (Ω)) c, t u h (t) L2 (Ω) ct α 1 for t (, T ]. The constant c above is independent of the mesh size h, but may depend on T. Remark 3.1. For smooth initial data and right-hand side, in the absence of extra compatibility conditions, the regularity results (3.1)-(3.2) and the h-independent estimates (3.3)-(3.4) are sharp with respect to the Hölder continuity in time. The regularity (3.1) was shown in [36] for linear subdiffusion equations and in [29] for a semilinear problem with Neumann boundary conditions under certain compatibility conditions. However, we are not aware of any existing results such as (3.2) and (3.3)-(3.4) for semilinear problems without compatibility conditions, which are important for the numerical analysis in Section 4. Remark 3.2. If f is smooth but not Lipschitz continuous, and problems (1.1) and (1.4) have unique bounded solutions, respectively, then f(u), f (u), f(u h ) and f (u h ) are still bounded. In this case, the estimates (3.1)-(3.2) and (3.3)-(3.4) are still valid, which can be seen from the proof of Theorem 3.1. We begin with some preliminary results. Let L 2 h (Ω) be the vector space S h equipped with the norm of L 2 (Ω) and let Hh 2(Ω) be the vector space S h equipped with the norm v h H 2 h (Ω) := v h L2 (Ω) + h v h L2 (Ω), v h S h. To analyze u(t) and u h (t) in a unified way, we consider the following abstract problem: { α t u(t) Au(t) = P f(u(t)) for t (, T ], (3.5) u() = u, where the notation (X, D, A, u, P, u ) denotes either (L 2 (Ω), H 1 (Ω) H 2 (Ω),, u, I, u ) or (L 2 h (Ω), H2 h (Ω), h, u h, P h, R h u ), with I denoting the identity operator. In a bounded convex polygonal domain Ω, the norm of D is equivalent to the graph norm, i.e., (3.6) v D v X + Av X, v D. Let X X be the operator norm on the space X. Then the operator A satisfies the following resolvent estimate [4, Example and Theorem ]: (z A) 1 X X c φ z 1, z Σ φ, φ (, π), where for φ (, π), Σ φ := {z C\{} : arg(z) < φ}. This further implies (3.7) (z α A) 1 X X c φ,α z α, A(z α A) 1 X X c φ,α, 1 z Σ φ, φ (, π), z Σ φ, φ (, π).
11 Let g(t) = P f(u(t)), and w := u u. Then w satisfies the following equation (3.8) α t w(t) Aw(t) = Au + g(t), with w() =. By means of Laplace transform, denoted by, we obtain z α ŵ(z) Aŵ(z) = z 1 Au + ĝ(z), which together with (3.7) implies ŵ(z) = (z α A) 1 (z 1 Au + ĝ(z)). By inverse Laplace transform and convolution rule, the solution w(t) to (3.8) is given by (3.9) w(t) = F (t)au + E(t s)g(s)ds, where the operators F (t) : X X and E(t) : X X are defined by (3.1) F (t) := 1 e zt z 1 (z α A) 1 dz and E(t) := 1 e zt (z α A) 1 dz, 2πi Γ θ,δ 2πi Γ θ,δ respectively. Clearly, we have E(t) = F (t). The contour Γ θ,δ is defined by (3.11) Γ θ,δ = {z C : z = δ, arg z θ} {z C : z = ρe ±iθ, ρ δ}, oriented with an increasing imaginary part, where θ (π/2, π) is fixed. In view of (3.9), u is the solution of problem (3.5) if and only if it is the solution of (3.12) u(t) u = F (t)au + E(t s)p f(u(s))ds. The next lemma summarizes the mapping properties of the operators F and E. These are partially known [36, Section 2] and [3]. We only sketch the proof for completeness. Lemma 3.2. For the operators F and E, the following properties hold. (i) t α F (t) X X + t 1 α F (t) X X + AF (t) X X c, t (, T ], (ii) F (t) : X D is continuous with respect to t [, T ], and AF () =. (iii) t 1 α E(t) X X + t 2 α E (t) X X + t AE(t) X X c, t (, T ]. (iv) E(t) : X D is continuous with respect to t (, T ]. Proof. First, consider (ii) in the case X = L 2 (Ω), D = H 1 (Ω) H 2 (Ω) and A =. By setting f(u(t)) and A = in (3.12), [36, Theorem 2.1] implies that F (t) = F (t) : L 2 (Ω) L 2 (Ω) is continuous with respect to t [, T ]. Thus, F (t) : L 2 (Ω) H 1 (Ω) H 2 (Ω) is continuous with respect to t [, T ]. Then taking t in (3.12) yields F () =. This proves (ii) in the case X = L 2 (Ω), D = H 1 (Ω) H 2 (Ω) and A =. The proof for the case X = L 2 h (Ω), D = H2 h (Ω) and A = h is similar. For any integers k and m =, 1, by choosing δ = t 1 in the contour Γ θ,δ and using the identity A(z α A) 1 = I + z α (z α A) 1, the resolvent estimate (3.7), and change of variables z = s cos ϕ + is sin ϕ, we have (with dz being the arc length element of Γ θ,δ ) dk Am dt k F (t) = 1 e zt z k 1 A m (z α A) 1 dz X X 2πi Γ θ,δ X X c e Re(z)t z k 1+(m 1)α dz Γ θ,δ c cos θ δ ct (m 1)α k. e st cos θ s k 1+(m 1)α ds + c 11 θ θ e cos ϕ δ k+(m 1)α dϕ
12 Since E(t) = F (t), the last inequality yields (i) and (iii). The continuity of F (t) : X D and E(t) : X D for t (, T ] follows from the equivalent norm in (3.6), showing (iv). Now we are ready to present the proof of Theorem 3.1. Proof of Theorem 3.1. The proof is divided into four steps. Step 1: Existence and uniqueness. We denote by C([, T ]; X) λ the function space C([, T ]; X) equipped with the following weighted norm: v λ := max t T e λt v(t) X, v C([, T ]; X), which is equivalent to the standard norm of C([, T ]; X) for any fixed parameter λ >. Then we define a nonlinear map M : C([, T ]; X) λ C([, T ]; X) λ by Mv(t) = u + F (t)au + E(t s)p f(v(s))ds. For any λ >, u C([, T ]; X) is a solution of (3.12) if and only if u is a fixed point of the map M : C([, T ]; X) λ C([, T ]; X) λ. It remains to prove that for some λ >, the map M : C([, T ]; X) λ C([, T ]; X) λ has a unique fixed point. In fact, the definition of M and Lemma 3.2(iii) immediately yield (3.13) e λt (Mv 1 (t) Mv 2 (t)) X = e λt E(t s)(p f(v 1 (s)) P f(v 2 (s)))ds ce λt (t s) α 1 v 1 (s) v 2 (s) X ds c (t s) α 1 e λ(t s) max s [,T ] e λs (v 1 (s) v 2 (s)) X ds ( 1 ) = cλ α (1 θ) α 1 (λt) α e λt(1 θ) dθ v 1 v 2 λ (change of variable s = tθ) c sup ([λt(1 θ)] α/2 e λt(1 θ)) ( 1 ) (t/λ) α/2 (1 θ) α/2 1 dθ v 1 v 2 λ λ>,t t> θ [,1] X c(t/λ) α/2 v 1 v 2 λ, v 1, v 2 C([, T ]; X) λ. By choosing a sufficiently large λ, the last inequality implies e λt (Mv 1 (t) Mv 2 (t)) X 1 2 v 1 v 2 λ, v 1, v 2 C([, T ]; X) λ. Hence, the map M is contractive on the space C([, T ]; X) λ. The Banach fixed point theorem implies that M has a unique fixed point, which is also the unique solution of (3.12). Step 2: C α ([, T ]; X) regularity. Consider the difference quotient for h > (3.14) u(t + h) u(t) h α F (t + h) F (t) = h α Au h α +h t E(s)P f(u(t s))ds P f(u(t + h s)) P f(u(t s)) E(s) h α ds =: 12 3 I i (t, h). i=1
13 A simple consequence of Lemma 3.2(i) is that h α F (t + h) F (t) X X c, which implies I 1 (t, h) X c. By appealing to Lemma 3.2(iii), we have I 2 (t, h) X = 1 t+h E(s)P f(u(t s))ds c 1 h α h α t +h t X s α 1 ds = c (t + h) α t α α h α c. By the Lipschitz continuity of f, we have e λt I 3 (t, h) X = P f(u(s + h)) P f(u(s)) e λt E(t s) h α ds X c 1 e λ(t s) (t s) α 1 e λs u(s + h) u(s) h α ds. X By substituting the estimates of I i (t, h), i = 1, 2, 3, into (3.14) and denoting W h (t) = e λt h α u(t + h) u(t) X, we obtain W h (t) c + c 1 e λ(t s) (t s) α 1 W h (s)ds c + c 1 (T/λ) α 2 max W h(s), s [,T ] where the last inequality can be derived in the same way as (3.13). By choosing a sufficiently large λ and taking maximum of the left-hand side with respect to t [, T ], it implies max t [,T ] W h(t) c, which further yields h α u(t + h) u(t) X ce λt c, where the constant c is independent of h. Thus, we have proved u C α ([,T ];X) c. Step 3: C([, T ]; D) regularity. By applying the operator A to both sides of (3.12) and using the identity AF (t) = AE(t s)ds, cf. Lemma 3.2, we obtain (3.15) Au(t) Au = AF (t)au + = AF (t) (Au + P f(u(t))) + = I 4 (t) + I 5 (t). AE(t s)p f(u(s))ds AE(t s)(p f(u(s) P f(u(t)))ds By Lemma 3.2(iii) and the C α ([, T ]; X) regularity from Step 2, we have I 5 (t) X = AE(t s)(p f(u(s)) P f(u(t)))ds c u(s) u(t) X ds t s X c t s α t s ds ctα, t (, T ]. Lemma 3.2(iv) implies that I 5 (t) is continuous for t (, T ], and the last inequality implies that I 5 (t) is also continuous at t =. Hence I 5 C([, T ]; X). Moreover, Lemma 3.2(ii) gives I 4 C([, T ]; X) and I 4 (t) X c Au + P f(u(t)) X c. 13
14 Substituting the estimates of I 4 (t) and I 5 (t) into (3.15) yields Au C([,T ];X) c, which further implies u C([,T ];D) c. The regularity result u C([, T ]; D) together with (3.5) yields α t u = Au + P f(u) C([, T ]; X). Step 4: Estimate of u (t) X. By differentiating (3.12) with respect to t, we obtain u (t) = F (t)au + E(t)P f(u ) + = E(t)(Au + P f(u )) + By multiplying this equation by t 1 α, we get t 1 α u (t) = t 1 α E(t)(Au + P f(u )) + E(s)P f (u(t s))u (t s)ds E(t s)p f (u(s))u (s)ds. t 1 α s α 1 E(t s)p f (u(s))s 1 α u (s)ds, which together with the L stability of P h [38, Lemma 6.1] directly implies that e λt t 1 α u (t) X e λt t 1 α E(t) X X Au + P f(u ) X + e λ(t s) t 1 α s α 1 (t s) α 1 P f (u(s)) L (Ω)e λs s 1 α u (s) X ds ce λt Au + P f(u ) X + c(t/λ) α 2 max s [,T ] e λs s 1 α u (s) X. where the last line follows similarly as (3.13). By choosing a sufficiently large λ and taking maximum of the left-hand side with respect to t [, T ], it implies max t [,T ] e λt t 1 α u (t) X c, which further yields (3.2). The proof of Theorem 3.1 is complete. 4. Error estimates. Now, we derive error estimates for the numerical solutions of problem (1.1) using the discrete Grönwall s inequality from Section 2 and discrete maximal l p -regularity from [17]. To illustrate the general framework for the numerical analysis of nonlinear time fractional diffusion equations, we focus on the L1 scheme and backward Euler CQ. Other time stepping schemes can be analyzed similarly. The convergence rates we show below are sharp (up to a logarithmic factor) with respect to the solution regularity in Theorem 3.1, and also confirmed by the numerical experiments in Section Preliminaries on the linear problem. First we recall some error estimates for the following linear subdiffusion equation: (4.1) α t v(t) v(t) = g(t), t (, T ], where g is a given function. The semidiscrete FEM for (4.1) seeks v h (t) S h such that (4.2) α t v h (t) h v h (t) = P h g(t), t (, T ], with v h () = R h v(), and the fully discrete scheme seeks v n h S h, n = 1,..., N, such that (4.3) α τ (v n h v h) h v n h = P h g(t n ), with vh = v h(), where τ α vh n denotes either the backward Euler CQ or the L1 scheme. The semidiscrete solution v h satisfies the following error estimate [14, 13, 16]. Lemma 4.1 (Semidiscrete solution of linear problems). For the semidiscrete solution v h to problem (4.2), there holds with l h = log(2 + 1/h) max t [,T ] v h(t) v(t) L2 (Ω) ch 2 v() H2 (Ω) + ch 2 l 2 h g L (,T ;L 2 (Ω)). 14
15 The solution vh n of the fully discrete scheme (4.3) satisfies the following error estimate. For the backward Euler CQ, it was proved in [16, Theorems 3.5 and 3.6], while the proof for the L1 scheme will be given in Section 5. Lemma 4.2 (Fully discrete solutions of linear problems). For the fully discrete solutions to problem (4.3) with the L1 scheme or backward Euler CQ, there holds v n h v h (t n ) vh n L2 (Ω) cτt α 1 n ( v() L2 (Ω) + g() L2 (Ω)) + cτ n (t n s) α 1 g (s) L2 (Ω)ds. Remark 4.1. If 1 d 3 and v() H 1 (Ω) H 2 (Ω), then the error estimates in Lemmas 4.1 and 4.2 are still valid if v h () is the Lagrange interpolation of v(), due to the smoothing property of the solution operator [14, Lemma 3.1]. Consequently, all the results in Section 4.2 remain valid in this case. Lemmas 4.1 and 4.2 will be used below in the analysis of the nonlinear problem Error estimates for the nonlinear problem. Now we can present error estimates for problem (1.1). Like in the linear case, we discuss the spatial error and temporal error separately. First, we derive the spatial discretization error. Theorem 4.3. Let u H 1 (Ω) H 2 (Ω), and f : R R be Lipschitz continuous. Then the semidiscrete problem (1.4) has a unique solution u h C([, T ]; L 2 h (Ω)), which satisfies (4.4) max u(t) u h(t) L2 (Ω) cl 2 hh 2. t T Proof. By Theorem 3.1, the existence and uniqueness of the solution u h hold. It remains to establish the estimate (4.4). To this end, we define v h (t) as the solution of α t v h (t) h v h (t) = P h f(u(t)), with v h () = u h () = R h u. This together with Lemma 4.1 yields the following estimate for t (4.5) (u v h )(t) L2 (Ω) ch 2 u() H2 (Ω) + ch 2 l 2 h f(u) L (,T ;L 2 (Ω)) ch 2 l 2 h. Meanwhile, we note that ρ h := v h u h satisfies the following equation α t ρ h (t) ρ h (t) = P h f(u(t)) P h f(u h (t)), with ρ h () =. Then, by the Lipschitz continuity of f and the maximal L p -regularity of fractional evolution equations [5, Corollary 1], we obtain the following estimate for any p (1, ) α t ρ h Lp (,T ;L 2 (Ω)) c P h f(u) P h f(u h ) Lp (,T ;L 2 (Ω)) c u u h Lp (,T ;L 2 (Ω)) c u v h Lp (,T ;L 2 (Ω)) + c ρ h Lp (,T ;L 2 (Ω)) ch 2 l 2 h + c ρ h L p (,T ;L 2 (Ω)). Then by the fractional Grönwall s inequality in Theorem 2.1, we have max ρ h(t) L 2 (Ω) ch 2 l 2 h. t [,T ] This and (4.5) directly imply the desired result. 15
16 Next we give the temporal discretization error. Theorem 4.4. Let u H 1 (Ω) H 2 (Ω), and f : R R be Lipschitz continuous. Then the fully discrete scheme (1.3), with either the L1 scheme or backward Euler CQ for time discretization, has a unique solution u n h S h, n = 1,..., N, and the solutions satisfy (4.6) max 1 n N u h(t n ) u n h L2 (Ω) cτ α. Proof. For given u h,, un 1 h, (1.3) is essentially a linear system with a symmetric positive definite matrix, and thus it has a unique solution u n h S h. It suffices to establish the estimate (4.6). Like before, we decompose the fully discrete solution u n h into two parts, u n h = vn h + ρn h, where vn h and ρn h respectively satisfy (4.7) (4.8) α τ (v n h v h) h v n h = P h f(u h (t n )), α τ ρ n h h ρ n h = P h f(u n 1 h ) P h f(u h (t n )), with vh = u h() = R h u and ρ h =. Equation (4.7) can be viewed as the time discretization of (1.4), with the right-hand side being a given function. Hence, by Lemma 4.2 and using s u h (s) L2 (Ω) cs α 1 (cf. Theorem 3.1) and Rademacher s theorem, we have ( ) u h (t n ) vh n L2 (Ω) ct α 1 n τ h u h () L2 (Ω) + f(u h ()) L2 (Ω) (4.9) + cτ n ct α 1 n τ + cτ ct α 1 n (t n s) α 1 f (u h (s)) s u h (s) L2 (Ω)ds n τ + ct 2α 1 n τ cτ α. (t n s) α 1 s α 1 ds It remains to estimate ρ n h. By applying the discrete maximal lp -regularity to (4.8) (choosing X = L 2 h (Ω) in [17, Theorems 3.1 and 4.1]), we obtain that for all 1 < p < : ( τ α ρ n h) m n=1 l p (L 2 (Ω)) c (f(u n 1 h ) f(u h (t n ))) m n=1 l p (L 2 (Ω)) c (f(u n 1 h ) f(u h (t n 1 ))) m n=1 lp (L 2 (Ω)) + c (f(u h (t n 1 )) f(u h (t n ))) m n=1 lp (L 2 (Ω)). By the Lipschitz continuity of f and the triangle inequality, we arrive at (f(u n 1 h ) f(u h (t n 1 ))) m n=1 l p (L 2 (Ω)) c (u h (t n 1 ) u n 1 h ) m n=1 l p (L 2 (Ω)) c (u h (t n 1 ) v n 1 h ) m n=1 lp (L 2 (Ω)) + c (ρ n 1 h cτ α + c (ρ n h) m 1 n=1 l p (L 2 (Ω)), ) m n=1 lp (L 2 (Ω)) where the last inequality follows from (4.9). Similarly, by the Lipschitz continuity of f and the a priori estimate u h C α ([,T ];L 2 (Ω)) c (cf. Theorem 3.1), we deduce ( f(u h (t n 1 )) f(u h (t n )) L 2 (Ω)) m n=1 l p c ( u h (t n 1 ) u h (t n ) L 2 (Ω)) m n=1 l p Combining the preceding three estimates yields c (cτ α ) m n=1 l p. ( α τ ρ n h) m n=1 l p (L 2 (Ω)) c (ρ n h) m 1 n=1 l p (L 2 (Ω)) + cτ α. 16
17 By choosing p > 1/α and applying the discrete Grönwall s inequality (with X = L 2 (Ω) in Theorem 2.6), we obtain (4.1) max 1 n N ρn h L2 (Ω) cτ α. In view of the decomposition u h (t n ) u n h = (u h(t n ) vh n) ρn h, the two estimates (4.9) and (4.1) imply (4.6), completing the proof of the theorem. Remark 4.2. If the nonlinear source f is not Lipschitz continuous but problem (1.1) has a unique bounded solution u, then Theorems 4.3 and 4.4 are still valid by proving the boundedness of the semidiscrete solution u h and the fully discrete solution u n h. For simplicity, we have assumed f to be Lipschitz continuous in order to avoid these technicalities. 5. Proof of Lemma 4.2 for the L1 scheme. The L1 scheme was analyzed in [15] only for the homogeneous problem. Below we give a proof for the general case. First, we assume that g is time-independent, i.e., g(t) g(). Then using Laplace transform, one can derive the following error representation (cf. [15, eq. (2.7) and (2.9)]): v h (t n ) vh n = 1 e ztn z 1 (z α h ) 1 ( h v h () + P h g())dz 2πi Γ θ,δ 1 e ztn µ(e zτ ) 1 (β τ (e zτ ) h ) 1 ( h v h () + P h g()) dz, 2πi Γ τ θ,δ where the contour Γ θ,δ is defined in (3.11), Γ τ θ,δ = {z Γ θ,δ : Im(z) 1/τ}, and µ(z) = 1 e zτ τe zτ and β τ (e zτ ) = (1 e zτ ) 2 e zτ τ α Γ(2 α) Li α 1(e zτ ), which satisfy the following estimates (cf. [15, Section 3]): (5.1) (5.2) c z µ(e zτ ) c 1 z and µ(e zτ ) z cτ z 2, z Γ τ θ,δ, β τ (e zτ ) c z τ 1 α and β τ (e zτ ) z α c z 2 τ 2 α, z Γ τ θ,δ. By using (5.1) (5.2), direct calculations yield (5.3) z 1 (z α h ) 1 µ(e zτ ) 1 (β τ (e zτ ) h ) 1 L2 (Ω) L 2 (Ω) c z α τ. Now we split the error v h (t n ) vh n into two components, i.e., v h(t n ) vh n = I 1 + I 2, where I 1 = 1 e ( ztn z 1 (z α h ) 1 µ(e zτ ) 1 (β τ (e zτ ) h ) 1) ( h v h () + P h g()) dz, 2πi Γ τ θ,δ I 2 = 1 e ztn z 1 (z α h ) 1 ( h v h () + P h g())dz. 2πi Γ θ,δ \Γ τ θ,δ By using (5.3) and (3.7), and choosing δ 1/t n, the argument from [15] yields (5.4) I 1 L2 (Ω) + I 2 L2 (Ω) ctn α 1 τ h v h () + P h g() L2 (Ω). Second, we consider the case v() = g() =. Then Taylor s expansion gives (5.5) P h g(t) = P h g() + 1 P h g (t) = 1 P h g (t). 17
18 In view of (3.9), the semidiscrete solution v h (t n ) can be represented by (5.6) Similarly, we have v h (t n ) = (E P h g)(t n ) = (E (1 P h g ))(t n ) = ((E 1) P h g )(t n ). (β τ (ξ) h ) 1 = n= Eτ n ξ n with Eτ n = τ e znτ (β τ (e zτ ) h ) 1 dz. 2πi Γ τ θ,δ Hence the fully discrete solution vh n can be represented by vn h = n second inequality of (5.2) implies (5.7) E n τ L2 (Ω) L 2 (Ω) ct α 1 n τ. j= En j τ P h g(t j ), and the Let E τ,ɛ (t) = n= En τ δ tn ɛ(t), where δ tn ɛ is the Dirac Delta function concentrated at t n ɛ, with ɛ (, τ). Then vh n can be rewritten as (5.8) v n h = lim ɛ (E τ,ɛ P h g)(t n ) = lim ɛ (E τ,ɛ (1 P h g ))(t n ) = (lim ɛ (E τ,ɛ 1) P h g )(t n ). The representations (5.6) and (5.8) yield (5.9) v h (t n ) v n h L 2 (Ω) [lim ɛ ((E E τ,ɛ ) 1) P h g ](t n ) L 2 (Ω). Using Laplace transform and Cauchy s integral formula, we deduce (lim(e E τ,ɛ ) 1)(t n ) = 1 e ztn z 1 (z α h ) 1 dz ɛ 2πi Γ θ,δ 1 e ztn µ(e zτ ) 1 (β τ (e zτ ) h ) 1 dz. 2πi Then using the estimate (5.3) we obtain (5.1) (lim(e E τ,ɛ ) 1)(t n ) L2 (Ω) L 2 (Ω) cτtn α 1. ɛ It remains to prove the following extension of the estimate (5.1): (5.11) (lim ɛ (E E τ,ɛ ) 1)(t) L 2 (Ω) L 2 (Ω) cτt α 1, t (, T ). Γ τ θ,δ Then this and (5.9) yield the second part on the right-hand side of (4.2), and completes the proof of Lemma 4.2. To prove (5.11), we consider the Taylor expansion of (E(t) E τ,ɛ (t)) 1 at t = t n, i.e., (5.12) ((E E τ,ɛ ) 1)(t) = ((E E τ,ɛ ) 1)(t n ) In view of Lemma 3.2 (iii), there holds n E(s) ds c L2 (Ω) L 2 (Ω) Similarly, appealing to (5.7), we have n lim E τ,ɛ (s) ds ɛ t t L 2 (Ω) L 2 (Ω) n t n t (E E τ,ɛ )(s) ds. s α 1 ds cτt α 1. = E n τ L 2 (Ω) L 2 (Ω) ct α 1 n τ. Substituting (5.1) and the last two inequalities into (5.12) yields (5.11). 18
19 6. Numerical experiments. In this section, we present numerical examples to verify the theoretical results in Theorems 4.3 and 4.4. We consider problem (1.1) with a diffusion coefficient.1 in the unit square Ω = (, 1) 2, with the following two sets of problem data: (a) u (x, y) = xy(1 x)(1 y) H 1 (Ω) H 2 (Ω) and f = 1 + u 2 ; (b) u (x, y) = x(1 x) sin(2πy) H 1 (Ω) H 2 (Ω) and f = 1 u 3. In the computation, we divided the domain Ω into regular right triangles with M equal subintervals of length h = 1/M on each side of the domain. The numerical solutions are computed by using the Galerkin FEM in space, and the backward Euler (BE) CQ or the L1 scheme in time. To evaluate the convergence, we compute the spatial error e t and temporal error e s, respectively, defined by e s = max u h(t n ) u(t n ) L2 (Ω) and e t = max 1 n N 1 n N un h u h (t n ) L2 (Ω). Since the exact solution to problem (1.1) is unavailable, we compute reference solutions on a finer mesh, i.e., the continuous solution u(t n ) with a fixed time step τ = 1/1 and mesh size h = 1/128, and the semidiscrete solution u h (t n ) with h = 1/1 and τ = 1/( ). In case (a), since the nonlinearity f is Lipschitz continuous, the theory in Section 4 applies. The numerical results for case (a) are shown in Tables 1 and 2, where the numbers in the bracket in the last column refer to the theoretical predictions from Section 4. We observe an O(h 2 ) rate for the spatial error e s, and an O(τ α ) rate for the temporal error e t for both backward Euler CQ and L1 scheme. These observations fully confirm Theorems 4.3 and 4.4. Table 1 Numerical results for case (a): the spatial error e s with T = 1, with N = 1, h = 1/M. α\m rate e-2 2.e e e e (2.).6 7.6e-2 2.5e e e e (2.) e e-2 5.8e e e (2.) Table 2 Numerical results for case (a): the temporal error e t with T = 1, τ = T/N, N = k 1 4, and h =.1. α k rate.4 BE 1.16e e e e e-4.39 (.4) L1 2.6e e e e e-4.38 (.4).6 BE 1.79e e e-5 5.1e e-5.6 (.6) L1 3.5e-4 2.2e e-4 8.8e e-5.6 (.6).8 BE 1.73e e e e e-6.8 (.8) L1 3.91e e e e e-6.8 (.8) In case (b), the nonlinear source f is not Lipschitz continuous. Nonetheless, one observes an O(h 2 ) and O(τ α ) convergence rate for the spatial and temporal errors, respectively, cf. Tables 3 and 4. This concurs with the discussions in Remarks 3.1 and 4.2. Further, the absolute accuracy of the L1 scheme and backward Euler CQ is comparable with each other for both cases (a) and (b). Interestingly, the spatial error e s increases slightly with the fractional order α, but the temporal error e t decreases with α. Acknowledgements. The authors are grateful to the anonymous referees for their constructive comments, which are very helpful to improve the presentation of the paper. 19
20 Table 3 Numerical results for case (b): the spatial error e s with T = 1, with N = 1, h = 1/M. α\m rate e e e e e-4 2. (2.).6 5.9e e e e e-4 2. (2.) e e e e-3 3.8e (2.) Table 4 Numerical results for case (b): the temporal error e t with T = 1, τ = T/N, N = k 1 4, h =.1. α k rate.4 BE 1.53e e-3 9.7e e e-4.38 (.4) L1 2.73e e e e e-4.38 (.4).6 BE 2.43e-4 1.6e-4 1.5e e e-5.6 (.6) L1 4.14e e e-4 1.2e e-5.6 (.6).8 BE 2,35e e e-6 4.4e e-6.8 (.8) L1 5.3e-5 3.4e e-5 1.e e-6.8 (.8) REFERENCES [1] G. Akrivis and B. Li, Maximum norm analysis of implicit-explicit backward difference formulae for nonlinear parabolic equations, IMA J. Numer. Anal., DOI: 1.193/imanum/drx8. [2] G. Akrivis, B. Li, and C. Lubich, Combining maximal regularity and energy estimates for time discretizations of quasilinear parabolic equations, Math. Comp., 86 (217), pp [3] A. A. Alikhanov, A new difference scheme for the time fractional diffusion equation, J. Comput. Phys., 28 (215), pp [4] W. Arendt, C. J. Batty, M. Hieber, and F. Neubrander, Vector-valued Laplace Transforms and Cauchy Problems, Birkhäuser, Basel, 2nd ed., 211. [5] E. Bazhlekova, Strict L p solutions for fractional evolution equations, Fract. Calc. Appl. Anal., 5 (22), pp [6] B. Berkowitz, J. Klafter, R. Metzler, and H. Scher, Physical pictures of transport in heterogeneous media: Advection-dispersion, random-walk, and fractional derivative formulations, Water Res. Research, 38 (22), pp [7] S. Blunck, Maximal regularity of discrete and continuous time evolution equations, Studia Math., 146 (21), pp [8] C. Chen, V. Thomée, and L. B. Wahlbin, Finite element approximation of a parabolic integrodifferential equation with a weakly singular kernel, Math. Comp., 58 (1992), pp [9] D. M. Creedon and J. J. H. Miller, The stability properties of q-step backward difference schemes, BIT, 15 (1975), pp [1] E. Cuesta, C. Lubich, and C. Palencia, Convolution quadrature time discretization of fractional diffusion-wave equations, Math. Comp., 75 (26), pp [11] P. Flajolet, Singularity analysis and asymptotics of Bernoulli sums, Theoret. Comput. Sci., 215 (1999), pp [12] E. Hairer, S. P. Nørsett, and G. Wanner, Solving Ordinary Differential Equations. I, Springer- Verlag, Berlin, second ed., 21. Nonstiff problems. [13] B. Jin, R. Lazarov, J. Pasciak, and Z. Zhou, Error analysis of semidiscrete finite element methods for inhomogeneous time-fractional diffusion, IMA J. Numer. Anal., 35 (215), pp [14] B. Jin, R. Lazarov, and Z. Zhou, Error estimates for a semidiscrete finite element method for fractional order parabolic equations, SIAM J. Numer. Anal., 51 (213), pp [15], An analysis of the L1 scheme for the subdiffusion equation with nonsmooth data, IMA J. Numer. Anal., 36 (216), pp [16], Two fully discrete schemes for fractional diffusion and diffusion-wave equations with nonsmooth data, SIAM J. Sci. Comput., 38 (216), pp. A146 A17. [17] B. Jin, B. Li, and Z. Zhou, Discrete maximal regularity of time-stepping schemes for fractional evolution equations. Preprint, arxiv: (to appear in Numer. Math.). [18] T. Kemmochi, Discrete maximal regularity for abstract Cauchy problems, Studia Math., 234 (216), pp [19] A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differ- 2
21 ential Equations, Elsevier Science B.V., Amsterdam, 26. [2] S. Kou, Stochastic modeling in nanoscale biophysics: Subdiffusion within proteins, Ann. Appl. Stat., 2 (28), pp [21] B. Kovács, B. Li, and C. Lubich, A-stable time discretizations preserve maximal parabolic regularity, SIAM J. Numer. Anal., 54 (216), pp [22] P. C. Kunstmann, B. Li, and C. Lubich, Runge-Kutta time discretization of nonlinear parabolic equations studied via discrete maximal parabolic regularity, Preprint, arxiv: (to appear in Found. Comput. Math.). [23] P. C. Kunstmann and L. Weis, Maximal L p-regularity for parabolic equations, Fourier multiplier theorems and H -functional calculus, in Functional Analytic Methods for Evolution Equations, vol of Lecture Notes in Math., Springer, Berlin, 24, pp [24] T. A. M. Langlands and B. I. Henry, The accuracy and stability of an implicit solution method for the fractional diffusion equation, J. Comput. Phys., 25 (25), pp [25] D. Leykekhman and B. Vexler, Discrete maximal parabolic regularity for Galerkin finite element methods, Numer. Math., 135 (217), pp [26] Y. Lin and C. Xu, Finite difference/spectral approximations for the time-fractional diffusion equation, J. Comput. Phys., 225 (27), pp [27] C. Lubich, Discretized fractional calculus, SIAM J. Math. Anal., 17 (1986), pp [28], Convolution quadrature and discretized operational calculus. I, Numer. Math., 52 (1988), p- p [29] Y. Luchko, W. Rundell, M. Yamamoto, and L. Zuo, Uniqueness and reconstruction of an unknown semilinear term in a time-fractional reactiondiffusion equation, Inverse Problems, 29 (213), p [3] W. McLean, Regularity of solutions to a time-fractional diffusion equation, ANZIAM J., 52 (21), pp [31] W. McLean and K. Mustapha, Time-stepping error bounds for fractional diffusion problems with non-smooth initial data, J. Comput. Phys., 293 (215), pp [32] R. Metzler and J. Klafter, The random walk s guide to anomalous diffusion: a fractional dynamics approach, Phys. Rep., 339 (2), pp [33] K. Mustapha, B. Abdallah, and K. M. Furati, A discontinuous Petrov-Galerkin method for timefractional diffusion equations, SIAM J. Numer. Anal., 52 (214), pp [34] K. Mustapha and H. Mustapha, A second-order accurate numerical method for a semilinear integrodifferential equation with a weakly singular kernel, IMA J. Numer. Anal., 3 (21), pp [35] R. R. Nigmatulin, The realization of the generalized transfer equation in a medium with fractal geometry, Phys. Stat. Sol. B, 133 (1986), pp [36] K. Sakamoto and M. Yamamoto, Initial value/boundary value problems for fractional diffusionwave equations and applications to some inverse problems, J. Math. Anal. Appl., 382 (211), pp [37] Z.-Z. Sun and X. Wu, A fully discrete scheme for a diffusion wave system, Appl. Numer. Math., 56 (26), pp [38] V. Thomée, Galerkin Finite Element Methods for Parabolic Problems, Springer-Verlag, Berlin, second ed., 26. [39] S. B. Yuste, Weighted average finite difference methods for fractional diffusion equations, J. Comput. Phys., 216 (26), pp [4] F. Zeng, C. Li, F. Liu, and I. Turner, Numerical algorithms for time-fractional subdiffusion equation with second-order accuracy, SIAM J. Sci. Comput., 37 (215), pp. A55 A78. 21
An Analysis of the L1 Scheme for the Subdiffusion Equation with Nonsmooth Data
IMA Journal of Numerical Analysis (214) Page 1 of 23 doi:1.193/imanum/drnxxx An Analysis of the L1 Scheme for the Subdiffusion Equation with Nonsmooth Data BANGTI JIN Department of Computer Science, University
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 information1. Introduction. We consider the model initial boundary value problem for the fractional order parabolic differential equation (FPDE) for u(x, t):
ERROR ESTIMATES FOR A SEMIDISCRETE FINITE ELEMENT METHOD FOR FRACTIONAL ORDER PARABOLIC EQUATIONS BANGTI JIN, RAYTCHO LAZAROV, AND ZHI ZHOU Abstract. We consider the initial boundary value problem for
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationCOMBINING MAXIMAL REGULARITY AND ENERGY ESTIMATES FOR TIME DISCRETIZATIONS OF QUASILINEAR PARABOLIC EQUATIONS
COMBINING MAXIMAL REGULARITY AND ENERGY ESTIMATES FOR TIME DISCRETIZATIONS OF QUASILINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS, BUYANG LI, AND CHRISTIAN LUBICH Abstract. We analyze fully implicit and linearly
More informationA POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS
A POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS D. LEYKEKHMAN AND L. B. WAHLBIN Abstract. This paper considers a posteriori error estimates by averaged gradients
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationCOMBINING MAXIMAL REGULARITY AND ENERGY ESTIMATES FOR TIME DISCRETIZATIONS OF QUASILINEAR PARABOLIC EQUATIONS
COMBINING MAXIMAL REGULARITY AND ENERGY ESTIMATES FOR TIME DISCRETIZATIONS OF QUASILINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS, BUYANG LI, AND CHRISTIAN LUBICH Abstract. We analyze fully implicit and linearly
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationNUMERICAL ANALYSIS OF FRACTIONAL-ORDER DIFFERENTIAL EQUATIONS WITH NONSMOOTH DATA. A Dissertation ZHI ZHOU
NUMERICAL ANALYSIS OF FRACTIONAL-ORDER DIFFERENTIAL EQUATIONS WITH NONSMOOTH DATA A Dissertation by ZHI ZHOU Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial
More informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationMaximum-norm stability of the finite element Ritz projection with mixed boundary conditions
Noname manuscript No. (will be inserted by the editor) Maximum-norm stability of the finite element Ritz projection with mixed boundary conditions Dmitriy Leykekhman Buyang Li Received: date / Accepted:
More informationA POSTERIORI ERROR ESTIMATES FOR THE BDF2 METHOD FOR PARABOLIC EQUATIONS
A POSTERIORI ERROR ESTIMATES FOR THE BDF METHOD FOR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND PANAGIOTIS CHATZIPANTELIDIS Abstract. We derive optimal order, residual-based a posteriori error estimates for
More informationConvergence analysis of a discontinuous Galerkin method for a sub-diffusion equation
Numer Algor (29 52:69 88 DOI 1.17/s1175-8-9258-8 ORIGINAL PAPER Convergence analysis of a discontinuous Galerkin method for a sub-diffusion equation William McLean Kassem Mustapha Received: 16 September
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationTHE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands
THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands 1. Introduction This paper deals with initial value problems for delay
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 informationOn boundary value problems for fractional integro-differential equations in Banach spaces
Malaya J. Mat. 3425 54 553 On boundary value problems for fractional integro-differential equations in Banach spaces Sabri T. M. Thabet a, and Machindra B. Dhakne b a,b Department of Mathematics, Dr. Babasaheb
More informationA NUMERICAL APPROXIMATION OF NONFICKIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXX, 1(21, pp. 75 84 Proceedings of Algoritmy 2 75 A NUMERICAL APPROXIMATION OF NONFICIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA R. E. EWING, Y. LIN and J. WANG
More informationL p MAXIMAL REGULARITY FOR SECOND ORDER CAUCHY PROBLEMS IS INDEPENDENT OF p
L p MAXIMAL REGULARITY FOR SECOND ORDER CAUCHY PROBLEMS IS INDEPENDENT OF p RALPH CHILL AND SACHI SRIVASTAVA ABSTRACT. If the second order problem ü + B u + Au = f has L p maximal regularity for some p
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 informationCONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS
J. Korean Math. Soc. 34 (1997), No. 3, pp. 515 531 CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS S. K. CHUNG, A.K.PANI AND M. G. PARK ABSTRACT. In this paper,
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More 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 informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationBull. Math. Soc. Sci. Math. Roumanie Tome 60 (108) No. 1, 2017, 3 18
Bull. Math. Soc. Sci. Math. Roumanie Tome 6 8 No., 27, 3 8 On a coupled system of sequential fractional differential equations with variable coefficients and coupled integral boundary conditions by Bashir
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 informationON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP
Dedicated to Professor Gheorghe Bucur on the occasion of his 7th birthday ON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP EMIL POPESCU Starting from the usual Cauchy problem, we give
More informationError estimates for a finite volume element method for parabolic equations in convex polygonal domains
Error estimates for a finite volume element method for parabolic equations in convex polygonal domains P Chatzipantelidis, R D Lazarov Department of Mathematics, Texas A&M University, College Station,
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
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 information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationGALERKIN AND RUNGE KUTTA METHODS: UNIFIED FORMULATION, A POSTERIORI ERROR ESTIMATES AND NODAL SUPERCONVERGENCE
GALERKIN AND RUNGE KUTTA METHODS: UNIFIED FORMULATION, A POSTERIORI ERROR ESTIMATES AND NODAL SUPERCONVERGENCE GEORGIOS AKRIVIS, CHARALAMBOS MAKRIDAKIS, AND RICARDO H. NOCHETTO Abstract. We unify the formulation
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationMEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EQUATIONS AND THEIR APPLICATIONS. Yalchin Efendiev.
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX MEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EUATIONS AND THEIR APPLICATIONS
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationIdentification of a memory kernel in a nonlinear parabolic integro-differential problem
FACULTY OF ENGINEERING AND ARCHITECTURE Identification of a memory kernel in a nonlinear parabolic integro-differential problem K. Van Bockstal, M. Slodička and F. Gistelinck Ghent University Department
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (211) 219 223 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Laplace transform and fractional differential
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
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 informationDIfferential equations of fractional order have been the
Multistage Telescoping Decomposition Method for Solving Fractional Differential Equations Abdelkader Bouhassoun Abstract The application of telescoping decomposition method, developed for ordinary differential
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
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 informationProperties of a New Fractional Derivative without Singular Kernel
Progr. Fract. Differ. Appl. 1, No. 2, 87-92 215 87 Progress in Fractional Differentiation and Applications An International Journal http://dx.doi.org/1.12785/pfda/122 Properties of a New Fractional Derivative
More informationExistence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions
Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions Mouffak Benchohra a,b 1 and Jamal E. Lazreg a, a Laboratory of Mathematics, University
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationOn Solutions of Evolution Equations with Proportional Time Delay
On Solutions of Evolution Equations with Proportional Time Delay Weijiu Liu and John C. Clements Department of Mathematics and Statistics Dalhousie University Halifax, Nova Scotia, B3H 3J5, Canada Fax:
More informationInterface Dynamics for Quasi- Stationary Stefan Problem
Interface Dynamics for Quasi- Stationary Stefan Problem R. Andrushkiw Department of Mathematical Sciences New Jersey Institute of Technology Newark, New Jersey 7 USA V. Gafiychuk Institute of Computer
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationMemoirs on Differential Equations and Mathematical Physics
Memoirs on Differential Equations and Mathematical Physics Volume 51, 010, 93 108 Said Kouachi and Belgacem Rebiai INVARIANT REGIONS AND THE GLOBAL EXISTENCE FOR REACTION-DIFFUSION SYSTEMS WITH A TRIDIAGONAL
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationResearch Article A New Fractional Integral Inequality with Singularity and Its Application
Abstract and Applied Analysis Volume 212, Article ID 93798, 12 pages doi:1.1155/212/93798 Research Article A New Fractional Integral Inequality with Singularity and Its Application Qiong-Xiang Kong 1 and
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 informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationComplex Analysis, Stein and Shakarchi Meromorphic Functions and the Logarithm
Complex Analysis, Stein and Shakarchi Chapter 3 Meromorphic Functions and the Logarithm Yung-Hsiang Huang 217.11.5 Exercises 1. From the identity sin πz = eiπz e iπz 2i, it s easy to show its zeros are
More informationCUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION
CUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION HANS CHRISTIANSON Abstract. This paper shows how abstract resolvent estimates imply local smoothing for solutions to the Schrödinger equation.
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationResearch Article Existence and Uniqueness of the Solution for a Time-Fractional Diffusion Equation with Robin Boundary Condition
Abstract and Applied Analysis Volume 211, Article ID 32193, 11 pages doi:1.1155/211/32193 Research Article Existence and Uniqueness of the Solution for a Time-Fractional Diffusion Equation with Robin Boundary
More informationA finite element method for time fractional partial differential equations
CHESTER RESEARCH ONLINE Research in Mathematics and its Applications ISSN 25-661 Series Editors: CTH Baker, NJ Ford A finite element method for time fractional partial differential equations Neville J.
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationFractional Evolution Integro-Differential Systems with Nonlocal Conditions
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 5, Number 1, pp. 49 6 (21) http://campus.mst.edu/adsa Fractional Evolution Integro-Differential Systems with Nonlocal Conditions Amar
More informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
More informationRUNGE-KUTTA TIME DISCRETIZATIONS OF NONLINEAR DISSIPATIVE EVOLUTION EQUATIONS
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 631 640 S 0025-571805)01866-1 Article electronically published on December 19, 2005 RUNGE-KUTTA TIME DISCRETIZATIONS OF NONLINEAR DISSIPATIVE EVOLUTION
More informationPositive solutions for a class of fractional boundary value problems
Nonlinear Analysis: Modelling and Control, Vol. 21, No. 1, 1 17 ISSN 1392-5113 http://dx.doi.org/1.15388/na.216.1.1 Positive solutions for a class of fractional boundary value problems Jiafa Xu a, Zhongli
More informationFractional differential equations with integral boundary conditions
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 (215), 39 314 Research Article Fractional differential equations with integral boundary conditions Xuhuan Wang a,, Liping Wang a, Qinghong Zeng
More informationANALYTIC SEMIGROUPS AND APPLICATIONS. 1. Introduction
ANALYTIC SEMIGROUPS AND APPLICATIONS KELLER VANDEBOGERT. Introduction Consider a Banach space X and let f : D X and u : G X, where D and G are real intervals. A is a bounded or unbounded linear operator
More informationOptimal order a posteriori error estimates for a class of Runge Kutta and Galerkin methods
Numerische Mathematik manuscript No. (will be inserted by the editor) Optimal order a posteriori error estimates for a class of Runge Kutta and Galerkin methods Georgios Akrivis 1, Charalambos Makridakis
More informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
More informationAn Introduction to Complex Analysis and Geometry John P. D Angelo, Pure and Applied Undergraduate Texts Volume 12, American Mathematical Society, 2010
An Introduction to Complex Analysis and Geometry John P. D Angelo, Pure and Applied Undergraduate Texts Volume 12, American Mathematical Society, 2010 John P. D Angelo, Univ. of Illinois, Urbana IL 61801.
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationALMOST PERIODIC SOLUTIONS OF HIGHER ORDER DIFFERENTIAL EQUATIONS ON HILBERT SPACES
Electronic Journal of Differential Equations, Vol. 21(21, No. 72, pp. 1 12. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu ALMOST PERIODIC SOLUTIONS
More informationFinite Difference Method for the Time-Fractional Thermistor Problem
International Journal of Difference Equations ISSN 0973-6069, Volume 8, Number, pp. 77 97 203) http://campus.mst.edu/ijde Finite Difference Method for the Time-Fractional Thermistor Problem M. R. Sidi
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 informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationL p -boundedness of the Hilbert transform
L p -boundedness of the Hilbert transform Kunal Narayan Chaudhury Abstract The Hilbert transform is essentially the only singular operator in one dimension. This undoubtedly makes it one of the the most
More informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More information