A PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM
|
|
- Ezra Wheeler
- 6 years ago
- Views:
Transcription
1 A PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Abstract. In this paper a singularly perturbed reaction diffusion partial differential equation in two space dimensions is examined. By means of an appropriate decomposition, we describe the asymptotic behaviour of the solution of problems of this kind. A central finite difference scheme is constructed for this problem which involves an appropriate Shishkin mesh. We prove that the numerical approximations are almost second order uniformly convergent (in the maximum norm) with respect to the singular perturbation parameter. Some numerical experiments are given which illustrate in practice the theoretical order of convergence established for the numerical method.. Introduction In this work we consider the two dimensional Dirichlet boundary value reactiondiffusion problem (.) Lu = f, (x, y) Ω, u = g, (x, y) Ω, where the differential operator is defined by Lω = ε(ω xx + ω yy ) + b(x, y)ω. The diffusion parameter satisfies < ε and it can be arbitrarily small. The domain Ω = (, ) 2 is the unit square, the reaction term satisfies b(x, y) 2β > and we assume that f, b C 4,α (Ω) and g C 4,α ( Ω). We also assume that there are sufficient compatibility conditions (see [4]), in order that u C 4,α (Ω). This type of problem is characterised by the presence of a regular exponential layer in a neighbourhood of Ω of width O( ε). It is well-known [3] that uniform meshes are an inappropriate discretization of the domain when these boundary layers are present in the solution. It is known that the differential operator L satisfies a comparison principle: For any y C 2 (Ω), if y on Ω and Ly in Ω then y for all points in the closed domain Ω. To simplify the notation, we consider the following subsets of the boundary Ω Γ = {(x, ) x }, Γ 2 = {(, y) y }, Γ 3 = {(x, ) x }, Γ 4 = {(, y) y }, which are the edges of Ω and the four corners of the domain are denoted by c = Γ Γ 2, c 2 = Γ 2 Γ 3, c 3 = Γ 3 Γ 4, c 4 = Γ 4 Γ. 99 Mathematics Subject Classification. Primary 65N6, 65N2, 65N5; Secondary 35J25. Key words and phrases. reaction-diffusion, uniform convergence, Shihskin mesh, second order. This research was partially supported by the project MCYT/FEDER BFM2 252.
2 2 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN We adopt the following notation for the boundary conditions: g(x, y) = g i (x), (x, y) Γ i, i =, 3; g(x, y) = g i (y), (x, y) Γ i, i = 2, 4. There is an extensive literature on numerical methods for singularly perturbed reaction diffusion problems (see, for example, [,, 2, 7] and the references therein). Our interest lies in examining parameter uniform numerical methods [3, 8] for singularly perturbed problems. That is, we are interested in numerical methods for which the following error bound can be theoretically established u Ū N CN p, p >, where N is the number of mesh elements employed in each coordinate direction, Ū N is a polynomial interpolant generated by the numerical method, is the global pointwise maximum norm and C is a constant independent of ε and N. In general, the gradients of the solution of (.) become unbounded as ε ; however, parameter uniform numerical methods guarantee that the error in the numerical approximation is controlled solely by the size of N. Li and Navon [6] studied (.) within a finite element framework and established a convergence result in the L 2 -norm. Note that the boundary layer function e x/ ε is of measure O(ε.25 ) in both the L 2 -norm and the standard energy norm f 2 = ε f x 2 + ε f y 2 + f 2. Hence, the size of the boundary layer function is negligible in these norms. In this paper, we establish a convergence result in the pointwise norm, where e x/ ε =. In the original work of Shishkin [3], parameter uniform numerical methods were established for an extensive class of linear singularly perturbed differential equations in n-dimensions, including convection diffusion ( ε u + a u + bu, a >, b ) and reaction diffusion ( ε u + bu, b > ). In [3], the author demonstrated the extent of the class of singularly perturbed problems for which the method (simple finite difference operator combined with an appropriate tensor product of piecewise uniform meshes) could be applied, by assuming minimal regularity on the data. As in the case of non-singularly perturbed problems, it is possible to obtain parameter uniform numerical methods with a higher order of parameter uniform convergence by restricting the class of problems to problems with smoother data. In one dimension, standard finite difference operators on an appropriate layeradopted mesh (Shishkin-type, Bakhvalov or variants [7, 8, ]) yield, up to a possible logarithmic factor, parameter-uniform first order convergence for convection diffusion problems ( εu + au + bu, a >, b ). For the one-dimensional singularly perturbed reaction diffusion problem ( εu + bu, b > ) parameter-uniform second order, again up to a possible logarithmic factor, is possible [9, 5]. In the case of convection diffusion problems in two space dimensions, first order schemes have been examined by several authors (for example, see [8, 3]). In this paper we present a parameter-uniform second order scheme for the reaction-diffusion problem in two space dimensions. The analytic properties of the solution to the two-dimensional reaction diffusion problem has been studied by [2] and [4]. In this paper, we take the approach of Shishkin [3] to establish parameter-explicit a priori bounds on the derivatives of the solutions, which are central in deriving our theoretical asymptotic error bound on the numerical approximations. The paper is structured as follows. In section 2 we give the bounds of the exact solution and its derivatives showing their asymptotic behaviour with respect to the
3 A ROBUST METHOD FOR REACTION-DIFFUSION PROBLEM 3 singular perturbation parameter. Moreover, we give an appropriate decomposition of this exact solution. In section 3 we analyze the standard central finite difference scheme constructed on a special mesh of Shishkin type, proving its ε-uniform convergence of second order. Finally in section 4 we show some numerical results, which illustrate the analytical results previously proved. Notation. Throughout this paper, C denotes a generic positive constant which is independent of the diffusion parameter ε and the discretization parameter N. We also use the following notation for the partial derivatives f (k,j) = k+j f x k y j. 2. Decomposition and a priori bounds In this section we examine the asymptotic behaviour of the solution of (.) with respect to the singular perturbation parameter ε. This behaviour will be used later in the analysis of the uniform convergence of the finite difference approximations defined in 3. From Han and Kellogg [4] we have that, if f, b C 2,α (Ω), g s C 4,α ([, ]), s =, 2, 3, 4 and compatibility conditions of second level are satisfied at each of the four corners, that is, (2.a) (2.b) (2.c) g () = g 2 (), g () g 2 () + b(, )g () = f(, ), g (4) () g(4) 2 () + f (2,) (, ) f (,2) (, ) b(, )g () 2b (,) (, )g () b (2,) (, )g () + +b(, )g 2 () + 2b (,) (, )g 2() + b (,2) (, )g 2 () =, and similarly for the other corners, then u C 4,α (Ω). Using a stretching argument and classical results from [5, 6] we can establish crude bounds on the derivatives of the solution of the form (2.2) u (k,j) Cε k/2 j/2, k + j 4, where is the maximum norm. Bounds (2.2) are not sufficient to analyze the uniform convergence of the numerical scheme studied in this paper, because they do not explicitly show the presence of boundary layers in Ω. Below we present a decomposition of u and appropriate bounds of its derivatives with respect to ε, which will be used in the error analysis in 3. This decomposition was first established by Shishkin [3], where the idea of extending the domain was introduced so that a decomposition of the solution into regular and singular components could be effected, without imposing additional artificial compatibility conditions on the data. For the sake of completeness, we present a detailed derivation of these bounds here. Let Ω = ( a, +a) ( a, +a), a > be an extended domain, which contains Ω as a subset. Define smooth extensions b, f and g s of the functions b, f and g s to Ω and [ a, + a], respectively. Note that f Ω = f. Let v = v + εv where v is the solution of the extended reduced problem b v = f,
4 4 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN and v is the solution of the problem (2.3) L v = v, (x, y) Ω, v =, (x, y) Ω. Since v C 4,α (Ω ), then v C 2,α (Ω ). The extensions of all the functions are taken such that the compatibility conditions at the corners of Ω of up to second order (see [4]) are satisfied. Hence v C 4,α (Ω ). The regular solution v is taken to be the solution of the boundary value problem (2.4) Lv = f, (x, y) Ω, v = v, (x, y) Ω. Applying the classical results (2.2) to the extended problem (2.3), we deduce that v Ω = v C 4,α (Ω) and (2.5) v (k,j) C( + ε k/2 j/2 ), k + j 4. The function v is called the regular component of the solution u. Note that by virtue of the extension, it is not necessary to impose compatibility conditions at the corners of Ω so that v C 4,α (Ω). Associated with the bottom edge Γ of the domain we have a boundary layer function w which is defined as follows. The domain Ω is extended in the horizontal direction to a domain Ω = ( a, + a) (, ), a >. The function w is the solution of (2.6a) (2.6b) (2.6c) (2.6d) L w =, (x, y) Ω, w = u v, (x, y) Γ, w(x, ) =, x [ a, + a], w( a, y) = w( + a, y) =, y [, ], and the values of the boundary conditions at the points (x, ) with x ( a, ) (, + a) are constructed so that w C 4,α (Ω ). So, for k 4 the extensions on ( a, ) (, + a) are constructed so that (w ) (k,) (, ) = (u v) (k,) (, ), (w ) (k,) (, ) = (u v) (k,) (, ), (w ) (k,) ( a, ) = (w ) (k,) ( + a, ) =. Moreover, we construct the extensions so that (2.7) (b ) (k,) ( a, y) = (b ) (k,) ( + a, y) =, k, for all y [, ]. Using the comparison principle, we deduce that (2.8) w (x, y) C (a + x)( + a x) e β/εy, (x, y) Ω. a( + a) The crude bounds (2.2) on the derivatives also apply to w, that is, (2.9) (w ) (k,j) Cε k/2 j/2, k + j 4. We now sharpen these bounds on the derivatives in the direction orthogonal to the layer. Using (2.8) and the fact that w ( a, y) = w ( + a, y) =, we get that the derivatives on the sides x = a and x = + a satisfy the bounds (w ) (,) ( a, y) Ce β/εy and (w ) (,) ( + a, y) Ce β/εy.
5 A ROBUST METHOD FOR REACTION-DIFFUSION PROBLEM 5 On the other two sides (w ) (,) (x, ) C and (w ) (,) (x, ) =. By differentiating the differential equation (2.6a) w.r.t. x, we get that L (w ) (,) = (b ) (,) (w ), (x, y) Ω, and using the comparison principle, it follows that (2.) (w ) (,) (x, y) Ce β/εy. Note that (w ) (2,) (x, ) C, (w ) (2,) (x, ) =, and from the regularity of w and the fact that ε((w ) (2,) + (w ) (,2) ) = b w in the closed region Ω, we have that (w ) (2,) ( a, y) = (w ) (,2) ( a, y) =, (w ) (2,) (+a, y) = (w ) (,2) (+a, y) =. Hence, after differentiating (2.6a) twice w.r.t. x the maximum principle establishes (2.) (w ) (2,) (x, y) C(a + x)( + a x)e β/εy. Using (2.) and the fact that (w ) (2,) ( a, y) = (w ) (2,) (+a, y) =, we get that which yields (w ) (3,) ( a, y) Ce β/εy and (w ) (3,) ( + a, y) Ce β/εy, (2.2) (w ) (3,) (x, y) Ce β/εy. Differentiating (2.6a) twice w.r.t. y yields the fact that the mixed derivatives satisfy (w ) (2,2) ( a, y) = (w ) (2,2) ( + a, y) =. Combine this with the result of differentiating (2.6a) twice w.r.t. x yields (2.3) (w ) (4,) ( a, y) = (w ) (4,) ( + a, y) =. Hence (2.4) (w ) (4,) (x, y) C. Associated with the bottom edge Γ we define a boundary layer function w to be the solution of the homogeneous problem (2.5a) (2.5b) (2.5c) Lw =, (x, y) Ω, w = u v, (x, y) Γ, w =, (x, y) Γ 3, w (, y) = w (, y), w (, y) = w (, y). In an analogous fashion, we can define boundary layer functions w k, k = 2, 3, 4 associated with the three other edges and the corresponding bounds on the derivatives of these functions will hold. Associated with the corner c we define a corner layer function z such that (2.6a) (2.6b) (2.6c) Lz =, (x, y) Ω, z = w 2, (x, y) Γ, z = w, (x, y) Γ 2, z =, (x, y) Γ 3, z =, (x, y) Γ 4. Note that Lw = Lw 2 = and w, w 2 C 4,α (Ω). Thus, the compatibility conditions up to second order (see [4]) hold at the four corners of the domain, which
6 6 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN implies that z C 4,α (Ω). The maximum principle and the condition b 2β >, result in the bound (2.7) z (x, y) Ce β/εx e β/εy. In an analogous fashion, we can define corner layer functions z k, k = 2, 3, 4 associated with the three other corners and the corresponding bounds hold. Remark 2.. From (2.2), we have the bounds (w ) (,), (z ) (,), (z ) (,) C ε. These bounds can be sharpened using the following argument. function to get that w (x, y) C y ε φ(x, y) = C(e β/εy e β/ε(2 y) ), e β/εy Use the barrier and (w ) (,) (x, ) C ε e β/ε. Using this and the crude bounds on the derivative of w we have that (2.8) (w ) (,) (x, y) C ε e β/εy. Analogous bounds hold for (w 3) (,) (x, y), (w 2) (,) (x, y) and (w 4) (,) (x, y). Now we sharpen the bounds on the first derivatives of the function z. From the above argument, it follows that w (x, y) C ( e β/εy e β/ε(2 y) ), w 2(x, y) C ( e β/εx e β/ε(2 x) ). Use the barrier function ψ(x, y) = C(e β/εx e β/ε(2 x) )(e β/εy e β/ε(2 y) ), with C sufficiently large and the discrete maximum principle to establish that (2.9) (z ) (,) (x, y), (z ) (,) (x, y) C ε e β/εx e β/εy. Analogous bounds hold for the first derivatives of the other three corner layer functions. We summarize this section with the following result: Theorem 2.2. The solution u of (.) may be written as a sum (2.2) u = v + where 4 w i + i= 4 z i, (2.2) Lv = f, Lw i = Lz i =, i =, 2, 3, 4. i=
7 A ROBUST METHOD FOR REACTION-DIFFUSION PROBLEM 7 Boundary conditions for v, w i, z i i =, 2, 3, 4, can be specified so that the following bounds on the derivatives of the components hold: (2.22a) (2.22b) (2.22c) (2.22d) (2.22e) (2.22f) (2.22g) (2.22h) (2.22i) v (k,j) C( + ε k/2 j/2 ), k + j 4, w (x, y) Ce β/εy, w 2 (x, y) Ce β/εx, w 3 (x, y) Ce β/ε( y), w 4 (x, y) Ce β/ε( x), max{ w (k,j) i, z (k,j) i } Cε k/2 j/2, k + j 4, w (k,) i C, i =, 3, w (,k) i C, i = 2, 4, k 4, z (x, y) Ce β/εy e β/εx, z 2 (x, y) Ce β/ε( y) e β/εx, z 3 (x, y) Ce β/ε( y) e β/ε( x), z 4 (x, y) Ce β/εy e β/ε( x). (3.) 3. The discrete problem To discretize problem (.) we use the standard central difference operator L N U N = ε(δ 2 x + δ 2 y)u + bu = f, (x i, y j ) Ω N, U = u, on the boundary Ω N, where the mesh Ω N is the tensor product of two one dimensional piecewise uniform Shishkin meshes, i.e., Ω N = Ω x Ω y, where Ω x (similarly for Ω y ) splits the interval [, ] into three subintervals [, σ x ], [σ x, σ x ] and [ σ x, ]. The mesh distributes N/4 points uniformly within each of the subintervals [, σ x ] and [ σ x, ] and the remaining N/2 mesh points uniformly in the interior subinterval [σ x, σ x ]. To simplify our discussion we take σ x = σ y ; these transition points are defined as σ = σ x = σ y = min { 4, 2 ε β ln N Below we denote by h = 4σ/N, H = 2( 2σ)/N; h i+ = x i+ x i, k i+ = y i+ y i, i N, and h i = (h i+ + h i )/2, k i = (k i+ +k i )/2, i N. It is well known that for the reaction diffusion problem, central differences is an ε-uniformly stable scheme in the maximum norm; that is, for any mesh function W (3.2) W 2β LN W + max Ω N W. From (3.2) we see that in order to prove uniform convergence, we only need to analyze the local truncation error. Using a standard truncation error argument, we can easily obtain (3.3) L N (U u)(x i, y j ) { Cε( h i u (3,) + k j u (,3) ), if x i = σ x, σ x, or y j = σ y, σ y, Cε(h 2 i u(4,) + kj 2 u(,4) ), otherwise. Nevertheless, to find appropriate bounds of this error we need to decompose the discrete solution similarly to the decomposition of the exact solution. The numerical }.
8 8 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN solution can be written in the form (3.4) U = V + where (3.5) (3.6) (3.7) 4 W k + k= 4 Z k, k= { L N V = f, in Ω N, V = v, in Ω N, { L N W k =, in Ω N, k =, 2, 3, 4, W k = w k, in Ω N, { L N Z k =, in Ω N, k =, 2, 3, 4, Z k = z k, in Ω N. From (2.5) and (3.3) a straightforward computation gives { L N C εn, if x i = σ x, σ x, or y j = σ y, σ y, (V v)(x i, y j ) CN 2, otherwise. Following [9], we define the barrier function Ψ(x i, y j ) = C σ xσ y N 2 (θ(x i ) + θ(y j )) + CN 2, ε where θ(z) is a piecewise linear polynomial defined by z σ, z σ, { 2 N θ(z) =, σ z σ, and δzθ(z) = σ, z = σ, σ, z σ, σ z,, otherwise. From the choice of transition points, it follows that Ψ(x i, y j ) C(N ln N) 2, and also we have { L N CσN + (bψ)(x i, y j ), if x i = σ x, σ x, or y j = σ y, σ y, (Ψ)(x i, y j ) = (bψ)(x i, y j ), otherwise. Then, the discrete maximum principle gives us that (3.8) V v C(N ln N) 2, which is the appropriate bound for the error associated with the regular component. To prove ε-uniform bounds of the errors associated with the edge and corner functions, we use an argument based on appropriate barrier functions. As usual for singularly perturbed problems, we consider the barrier functions j ( ), i ( ), + ks β/ε j, + hs β/ε i, B w ;j = B w2 ;i = s= s=, j =,, i =, N ( ), N ( ), + ks β/ε j N, + hs β/ε i N, B w3 ;j = B w4 ;i = s=j+ s=i+, j = N,, i = N.
9 A ROBUST METHOD FOR REACTION-DIFFUSION PROBLEM 9 These functions are first order Taylor approximations of the exponential functions related to the singular behaviour of the solution of problem (.). Note that for all j it holds exp ( j β/εy j ) = exp ( β/εk s ) B w;j, and for σ <.25 and N/4 j N we have ( (3.9) B w;j B w ;N/4 = + 8 ln N ) N/4 CN 2, N s= (3.) L N B w;j (b(x i, y j ) 2β)B w;j. Analogous bounds hold for the other three edge functions. Proposition 3.. If w k and W k are the solutions of (2.5) and (3.6) respectively, then for k =, 2, 3, 4 (3.) w k (x i, y j ) W k (x i, y j ) C(N ln N) 2, (x i, y j ) Ω N. Proof. We assume throughout that σ <.25. The case of σ =.25 is dealt with in a classical fashion by noting that ε C(ln N) 2 in this case. We only give the details for the edge layer function w. The argument is analogous for the other three boundary layer functions. From (3.6) and Theorem 2.2 it follows that on the boundary Ω N (3.2) W (x i, y j ) = w (x i, y j ) C exp ( β/εy j ) CB w;j, (x i, y j ) Ω N. Also for each internal mesh point (x i, y j ) Ω N, < i, j < N, from (3.6), (3.) and the discrete maximum principle it follows that (3.3) W (x i, y j ) B w;j. Therefore, using Theorem 2.2 and (3.3), we deduce w (x i, y j ) W (x i, y j ) w (x i, y j ) + W (x i, y j ) CB w;j. Hence, for the mesh points that are not close to the edge Γ, from (3.9) we have (3.4) w (x i, y j ) W (x i, y j ) CN 2, i N, N/4 j N. To prove similar bounds of the error in the region Ω N = {(x i, y j ) < i < N, < j < N/4}, we proceed as follows. Using Taylor expansions we obtain { L N Cε( h i w (3,) + kj 2 [W (x i, y j ) w (x i, y j )] w(,4) ), i = N/4, 3N/4, Cε(h 2 i w(4,) + kj 2 w(,4) ), otherwise. From Theorem 2.2, it follows that { L N C((N ln N) 2 + N ε), i = N/4, 3N/4, [W (x i, y j ) w (x i, y j )] C(N ln N) 2, otherwise. Then, similar to the analysis for the regular component, we use the barrier function Ψ(x i, y j ) = C σ x ε N 2 θ(x i ) + C(N ln N) 2, and the discrete maximum principle, now applied only on Ω N, to get (3.5) w (x i, y j ) W (x i, y j ) C(N ln N) 2, (x i, y j ) Ω N.
10 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN From (3.4) and (3.5) the result follows. Proposition 3.2. If z k and Z k are the solutions of (2.6) and (3.7) respectively, then for k =, 2, 3, 4 (3.6) z k (x i, y j ) Z k (x i, y j ) C(N ln N) 2, (x i, y j ) Ω N. Proof. Again we only give the proof of (3.6) for the corner layer function z and in the case of σ <.25. In a similar way to the argument given in Proposition 3., we obtain and Z (x i, y j ) C min{b w2 ;i, B w ;j}, if (x i, y j ) Ω N, z (x i, y j ) Z (x i, y j ) C min{b w2 ;i, B w ;j}, if (x i, y j ) Ω N, < i, j < N. Then, using (3.9) we deduce that z (x i, y j ) Z (x i, y j ) CN 2, (x i, y j ) Ω N \Ω N,2, where Ω N,2 = {(x i, y j ) < i, j < N/4}. Finally, in Ω N,2 the truncation error satisfies L N [Z (x i, y j ) z (x i, y j )] Cεh 2 ( z (4,) + z (,4) ) C(N ln N) 2, where we have used Theorem 2.2. Considering the barrier function Ψ(x i, y j ) = C(N ln N) 2, the discrete maximum principle, used on Ω N,2, proves the required result. Therefore, from the bounds (3.8) and Propositions 3. and 3.2 we deduce the following result of uniform convergence. Theorem 3.3. Let u be the solution of problem (.) and U the numerical solution of (3.) defined on the piecewise uniform Shihskin mesh. Then, the error at the mesh points satisfies (3.7) (u U)(x i, y j ) C(N ln N) 2, (x i, y j ) Ω N. A further property of Shishkin meshes is parameter-uniform interpolation. Let N w I (x, y) = w(x i, y j )φ i (x)φ j (y), i,j= where φ i (x) is the standard piecewise linear basis function associated with the interval [x i, x i+ ]. Define a global bilinear approximation to the solution u of (.) as N Ū N (x, y) = U N (x i, y j )φ i (x)φ j (y). i,j= Note that from the previous theorem u Ū N u u I + C(N ln N) 2. From the argument in [4] and the bounds on the components in the decomposition given in Theorem 2.2 and the sharper bounds on the first derivatives of the layer
11 A ROBUST METHOD FOR REACTION-DIFFUSION PROBLEM components given in (2.8), (2.9), we have the global parameter-uniform error bound (3.8) u Ū N C(N ln N) Numerical experiments In this section we present numerical results obtained by applying the numerical method described in 3 to two particular problems of the form (.). To estimate the maximum errors we use a variant of the double mesh principle. The two-mesh difference is calculated using ( ) Dε N 2N = max Ũ2i,2j Ui,j N, i,j N where {Ũ 2N } is the numerical solution on a mesh which contains the mesh points (x i, y j ) of Ω N and also the midpoints x i+/2 = (x i + x i+ )/2, y i+/2 = (y i + y i+ )/2, i =,,..., N. From these values we define the ε-uniform differences by D N = max Dε N, the numerical orders of convergence are calculated by p N ε = ε log (Dε N /Dε 2N )/ log 2 and the ε-uniform order is given by p N uni = log (DN /D 2N )/ log 2. The first test problem is (4.) ε u + ( + x 2 y 2 )u =, (x, y) (, ) 2, u(x, ) = ( x) 2, x, u(, y) = y, y, u(x, ) = u(, y) =, x, y. Figure shows the numerical solution for ε = 6 and N = 64. From it we see that the solution has a corner layer at (,) and two boundary layers near the edges x = and y = of the unit square. Note that the reduced solution is identically zero and that u = Z X Y Figure. Numerical solution of problem (4.) (ε = 6, N = 64)
12 2 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Table displays the maximum differences Dε N and the numerical orders of convergence p N ε for ε =, 2 2, 2 4,,, These results indicate that the uniform order of convergence is in agreement with Theorem 3.3, even though we do not have sufficiently compatibility conditions (for example in the corner (,) only the condition of order zero holds) for problem (4.). Table. Maximum differences Dε N convergence p N ε for problem (4.) and the numerical orders of ε N = 32 N = 64 N = 28 N = 256 ε = 3.96E E E E ε = E E E E ε = E E 5.78E 5 4.3E ε = E E E 5.989E ε = E 3.38E E E ε = 2.74E E 3.36E E ε = E 2.25E E 3.334E ε = E 2.35E E 3.345E ε = E 2.4E E 3.35E ε = E 2.42E E 3.354E ε = E 2.43E 2 4.E 3.355E ε = E 2.44E 2 4.3E 3.356E ε = E 2.44E 2 4.4E 3.356E ε = E 2.44E 2 4.5E 3.356E D N 2.555E 2.44E 2 4.5E 3.356E 3 p N uni Figure 2 shows the approximate pointwise errors for ε = 6 calculated by comparing the numerical solution obtained with N = 32 with the numerical solution obtained using a finer mesh having N = 52, in such way that the transition points of the mesh Ω 32 and the finer mesh Ω 52 coincide. The finer mesh Ω 52 is obtained by dividing each subinterval of the coarse mesh Ω 32 into sixteen subintervals of equal length. We see that the maximum error occurs in the corner layer region. Finally, following [3], we estimate the ε uniform error constant. From the value p = min N pn uni,
13 A ROBUST METHOD FOR REACTION-DIFFUSION PROBLEM Z X Y Figure 2. Approximate errors U 32 Ũ 52 for problem (4.) (ε = 6 ) we calculate (see table 2) Cp N = DN N p. 2 p Then, the ε uniform error constant is defined by Table 2. Values of Cp N for problem (4.) N = 32 N = 64 N = 28 N = C p = max N CN p, which in this case takes the value As in [3], we propose the following parameter uniform error estimate for the numerical approximations U N of (3.) to the solution u of problem (.) (4.2) u U N 3.79N.29, N 32. The second test problem that we consider is ε u + ( + x) 2 u = (x 2 + 2x) 2, (x, y) (, ) 2, u(x, ) = + 3( x) 2, x, u(, y) = 4( 2y( y)) 2, y, u(x, ) = + 3( x), x, u(, y) = y( y), y. Now the solution has four boundary and corner layers as is shown in Figure 3. Note that the reduced solution is not identically zero and that u = 4. Table 3 displays the differences and the order in this case for the same values of ε as in the previous
14 4 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Z X Y Figure 3. Numerical solution of problem (4.2) (ε = 6, N = 64) example. Again we only have zero order compatibility conditions; nevertheless we observe orders of convergence tending towards two. Figure 4 shows the pointwise errors for ε = 6 using the same technique as before. Again we see that the maximum error occurs in the layer regions. Finally, Z. Y X Figure 4. Approximate errors U 32 Ũ 52 for problem (4.2) (ε = 6 ) Table 4 displays the values of Cp N constant for problem (4.2). and in bold we indicate the ε uniform error
15 A ROBUST METHOD FOR REACTION-DIFFUSION PROBLEM 5 Table 3. Maximum differences Dε N convergence p N ε for problem (4.2) and the numerical orders of ε N = 32 N = 64 N = 28 N = 256 ε =.445E E E E ε = E E E 5.675E ε = E E 4.5E E ε = E 3.78E E 4.2E ε = E E 3.724E E ε = E E E 3.74E ε = E 2 4.7E 2.73E E ε = E E 2.739E E ε = 2 6.E 4.52E 2.743E 2 6.4E ε = 2 8.7E 4.65E 2.745E E ε = 2 2.2E 4.7E 2.746E E ε = E 4.75E 2.746E E ε = E 4.76E 2.747E E ε = E 4.77E 2.747E E D N.22E 4.77E 2.747E E 3 p N uni Table 4. Values of Cp N for problem (4.2) N = 32 N = 64 N = 28 N = References. T. Apel and G. Lube, Anisotropic mesh refinement for singularly perturbed reaction diffusion problems, Appl. Numer. Math. 26 (998), V. F. Butuzov, The asymptotic properties of solutions of the equation µ 2 u k 2 (x, y)u = f(x, y) in a rectangle, Differential nye Uravneniya 9 (973), P.A. Farrell, A.F. Hegarty, J.J.H. Miller, E. O Riordan and G.I. Shishkin, Robust Computational Techniques for Boundary Layers, Chapman and Hall/CRC Press, Boca Raton, U.S.A., H. Han and R.B. Kellogg, Differentiability properties of solutions of the equation ε 2 u + ru = f(x, y) in a square, SIAM J. Math. Anal. 2 (99),
16 6 C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN 5. O.A. Ladyzhenskaya and N. N. Ural tseva, Linear and Quasilinear Elliptic Equations, Academic Press, New York and London, J. Li and I. M. Navon, Uniformly convergent finite element methods for singularly elliptic boundary value problems I: reaction diffusion type, Computers Math. Applic 35 (998), T. Linß, Layer-adapted meshes for convection diffusion problems, Comput. Methods Appl. Mech. Engrg. 92 (23), J. J. H. Miller, E. O Riordan and G. I. Shishkin, Fitted numerical methods for singular perturbation problems, World-Scientific, (Singapore), J.J.H. Miller, E. O Riordan, G.I. Shishkin and L.P. Shishkina, Fitted mesh methods for problems with parabolic boundary layers, Mathematical Proceedings of the Royal Irish Academy 98A, (998), E. O Riordan M. and Stynes, A uniformly accurate finite element method for a singularly perturbed one dimensional reaction diffusion problem, Math. Comput. 47 (986), H.-G. Roos, M. Stynes and L. Tobiska, Numerical Methods for Singularly Perturbed Differential Equations. Convection Diffusion and Flow Problems, Springer Verlag, New York, A. H. Schatz and L. B. Wahlbin, On the finite element method for singularly perturbed reaction-diffusion problems in two and one dimensions, Math. Comp. 4 (983), G.I. Shishkin, Discrete approximation of singularly perturbed elliptic and parabolic equations, Russian Academy of Sciences, Ural section, Ekaterinburg, 992. (In Russian) 4. M. Stynes and E. O Riordan, A uniformly convergent Galerkin method on a Shishkin mesh for a convection-diffusion problem, J. Math. Anal. Appl. 24 (997), G. Sun and M. Stynes, Finite element methods for singularly high order elliptic two point boundary value problems I: reaction diffusion type problems, IMA Journal of Numerical Analysis 5 (995), E. A. Volkov, On the differential properties of solutions of boundary value problems for the Laplace and Poisson equations on a rectangle, Trudy Mat. Inst. Syeklov 77 (965), C. Xenophontos and S. R. Fulton, Uniform approximation of singularly perturbed reactiondiffusion problems by the finite element method on a Shishkin mesh, Numer. Methods Partial Differential Eq. 9 (23), 89. Departamento de Matemática Aplicada, Universidad de Zaragoza, Zaragoza, Spain. address: clavero@unizar.es Departamento de Matemática Aplicada, Universidad de Zaragoza, Teruel, Spain. address: jlgracia@unizar.es School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland. address: eugene.oriordan@dcu.ie
Interior Layers in Singularly Perturbed Problems
Interior Layers in Singularly Perturbed Problems Eugene O Riordan Abstract To construct layer adapted meshes for a class of singularly perturbed problems, whose solutions contain boundary layers, it is
More informationA coupled system of singularly perturbed semilinear reaction-diffusion equations
A coupled system of singularly perturbed semilinear reaction-diffusion equations J.L. Gracia, F.J. Lisbona, M. Madaune-Tort E. O Riordan September 9, 008 Abstract In this paper singularly perturbed semilinear
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationParameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition
Parameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition P. A. Farrell 1,A.F.Hegarty 2, J. J. H. Miller 3, E. O Riordan 4,andG.I. Shishkin 5 1
More informationSuperconvergence analysis of the SDFEM for elliptic problems with characteristic layers
Superconvergence analysis of the SDFEM for elliptic problems with characteristic layers S. Franz, T. Linß, H.-G. Roos Institut für Numerische Mathematik, Technische Universität Dresden, D-01062, Germany
More informationA Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer
A Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer E. O Riordan, J. Quinn School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland. Abstract
More informationIterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations
Iterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations P.A. Farrell 1, P.W. Hemker 2, G.I. Shishkin 3 and L.P. Shishkina 3 1 Department of Computer Science,
More informationA uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction diffusion problems
IMA Journal of Numerical Analysis (2003 23, 627 644 A uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction diffusion problems NIALL MADDEN Department of
More informationA Linearised Singularly Perturbed Convection- Diffusion Problem with an Interior Layer
Dublin Institute of Technology ARROW@DIT Articles School of Mathematics 201 A Linearised Singularly Perturbed Convection- Diffusion Problem with an Interior Layer Eugene O'Riordan Dublin City University,
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 informationA first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh
A first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh James Adler, Department of Mathematics Tufts University Medford, MA 02155 Scott MacLachlan Department
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationParameter robust methods for second-order complex-valued reaction-diffusion equations
Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,
More informationQuintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems
Stud. Univ. Babeş-Bolyai Math. 63(2018, No. 1, 141 151 DOI: 10.24193/subbmath.2018.1.09 Quintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems Ram
More informationFinite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems
Finite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 Received 19 July 2000; accepted
More informationSUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1)
Journal of Computational Mathematics, Vol.5, No., 007, 185 00. SUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1) Ziqing Xie (College of Mathematics and Computer Science,
More informationAn exponentially graded mesh for singularly perturbed problems
An exponentially graded mesh for singularly perturbed problems Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with P. Constantinou (UCY), S. Franz (TU Dresden)
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 informationACCURATE NUMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRANSFER
FDS-2000 Conference, pp. 1 10 M. Sapagovas et al. (Eds) c 2000 MII ACCURATE UMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRASFER B. GAHA 1, J. J. H. MILLER 2, and G. I. SHISHKI
More informationBoundary layers in a two-point boundary value problem with fractional derivatives
Boundary layers in a two-point boundary value problem with fractional derivatives J.L. Gracia and M. Stynes Institute of Mathematics and Applications (IUMA) and Department of Applied Mathematics, University
More informationError estimates for moving least square approximations
Applied Numerical Mathematics 37 (2001) 397 416 Error estimates for moving least square approximations María G. Armentano, Ricardo G. Durán 1 Departamento de Matemática, Facultad de Ciencias Exactas y
More informationA uniformly convergent difference scheme on a modified Shishkin mesh for the singularly perturbed reaction-diffusion boundary value problem
Journal of Modern Methods in Numerical Mathematics 6:1 (2015, 28 43 A uniformly convergent difference scheme on a modified Shishkin mesh for the singularly perturbed reaction-diffusion boundary value problem
More informationConvergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements
Convergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements Sebastian Franz a,1, Gunar Matthies b, a Department of Mathematics and Statistics,
More informationSuperconvergence Using Pointwise Interpolation in Convection-Diffusion Problems
Als Manuskript gedruckt Technische Universität Dresden Herausgeber: Der Rektor Superconvergence Using Pointwise Interpolation in Convection-Diffusion Problems Sebastian Franz MATH-M-04-2012 July 9, 2012
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationA COMPUTATIONAL TECHNIQUE FOR SOLVING BOUNDARY VALUE PROBLEM WITH TWO SMALL PARAMETERS
Electronic Journal of Differential Equations, Vol. 2013 (2013), No. 30, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu A COMPUTATIONAL
More informationCholesky factorisations of linear systems coming from a finite difference method applied to singularly perturbed problems
Cholesky factorisations of linear systems coming from a finite difference method applied to singularly perturbed problems Thái Anh Nhan and Niall Madden The Boundary and Interior Layers - Computational
More informationAnalysis of a Streamline-Diffusion Finite Element Method on Bakhvalov-Shishkin Mesh for Singularly Perturbed Problem
Numer. Math. Theor. Meth. Appl. Vol. 10, No. 1, pp. 44-64 doi: 10.408/nmtma.017.y1306 February 017 Analysis of a Streamline-Diffusion Finite Element Method on Bakhvalov-Shishkin Mesh for Singularly Perturbed
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskunde en Informatica Modelling, Analysis and Simulation Modelling, Analysis and Simulation High-order time-accurate schemes for singularly perturbed parabolic convection-diffusion problems
More informationA New Improvement of Shishkin Fitted Mesh Technique on a First Order Finite Difference Method
A New Improvement of Shishkin Fitted Mesh Technique on a First Order Finite Difference Method With applications on singularly perturbed boundary value problem ODE Rafiq S. Muhammad College of Computers
More informationUNIFORMLY CONVERGENT 3-TGFEM VS LSFEM FOR SINGULARLY PERTURBED CONVECTION-DIFFUSION PROBLEMS ON A SHISHKIN BASED LOGARITHMIC MESH
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 4, Number 4, Pages 35 334 c 23 Institute for Scientific Computing and Information UNIFORMLY CONVERGENT 3-TGFEM VS LSFEM FOR SINGULARLY
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
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 informationFitted operator finite difference scheme for a class of singularly perturbed differential difference turning point problems exhibiting interior layer
Fitted operator finite difference scheme for a class of singularly perturbed differential difference turning point problems exhibiting interior layer Dr. Pratima Rai Department of mathematics, University
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
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 informationMultigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
More informationFinite Difference, Finite element and B-Spline Collocation Methods Applied to Two Parameter Singularly Perturbed Boundary Value Problems 1
European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol. 5, no. 3-4, 2011, pp. 163-180 ISSN 1790 8140
More informationMAXIMUM NORM A POSTERIORI ERROR ESTIMATE FOR A 3D SINGULARLY PERTURBED SEMILINEAR REACTION-DIFFUSION PROBLEM
MAXIMUM NORM A POSTERIORI ERROR ESTIMATE FOR A 3D SINGULARLY PERTURBED SEMILINEAR REACTION-DIFFUSION PROBLEM NARESH M. CHADHA AND NATALIA KOPTEVA Abstract. A singularly perturbed semilinear reaction-diffusion
More informationarxiv:math/ v2 [math.ap] 3 Oct 2006
THE TAYLOR SERIES OF THE GAUSSIAN KERNEL arxiv:math/0606035v2 [math.ap] 3 Oct 2006 L. ESCAURIAZA From some people one can learn more than mathematics Abstract. We describe a formula for the Taylor series
More informationA UNIFORMLY ACCURATE FINITE ELEMENTS METHOD FOR SINGULAR PERTURBATION PROBLEMS. Mirjana Stojanović University of Novi Sad, Serbia and Montenegro
GLASNIK MATEMATIČKI Vol. 38(58(2003, 185 197 A UNIFORMLY ACCURATE FINITE ELEMENTS METHOD FOR SINGULAR PERTURBATION PROBLEMS Mirjana Stojanović University of Novi Sad, Serbia and Montenegro Abstract. We
More informationAvailable online at ScienceDirect. Procedia Engineering 127 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 127 (2015 ) 424 431 International Conference on Computational Heat and Mass Transfer-2015 Exponentially Fitted Initial Value
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 informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationNumerical Studies of Adaptive Finite Element Methods for Two Dimensional Convection-Dominated Problems
DOI 10.1007/s10915-009-9337-6 Numerical Studies of Adaptive Finite Element Methods for Two Dimensional Convection-Dominated Problems Pengtao Sun Long Chen Jinchao Xu Received: 17 April 2009 / Revised:
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationLINK-CUTTING BUBBLES FOR THE STABILIZATION OF CONVECTION-DIFFUSION-REACTION PROBLEMS
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 445 461 c World Scientific Publishing Company LINK-CUTTING BUBBLES FOR THE STABILIZATION OF CONVECTION-DIFFUSION-REACTION PROBLEMS
More informationNUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE
IJRRAS 5 () November NUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE. K. Sharath Babu,* & N. Srinivasacharyulu Research Scholar, Department of Mathematics, National Institute
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 informationNUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE
TWMS Jour. Pure Appl. Math., V.3, N.1, 2012, pp.11-21 NUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE R. MOHAMMADI 1 Abstract. We use adaptive
More informationON MIXED BOUNDARY VALUE PROBLEMS FOR PARABOLIC EQUATIONS IN SINGULAR DOMAINS
Azzam, A. Osaka J. Math. 22 (1985), 691-696 ON MIXED BOUNDARY VALUE PROBLEMS FOR PARABOLIC EQUATIONS IN SINGULAR DOMAINS ALI AZZAM (Received February 28, 1984) 1. Introduction. In this paper we continue
More informationFinite volume element method for analysis of unsteady reaction diffusion problems
Acta Mech Sin (2009) 25:481 489 DOI 10.1007/s10409-009-0237-7 RESEARCH PAPER Finite volume element method for analysis of unsteady reaction diffusion problems Sutthisak Phongthanapanich Pramote Dechaumphai
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 41, pp. 376-395, 2014. Copyright 2014,. ISSN 1068-9613. A ROBUST NUMERICAL SCHEME FOR SINGULARLY PERTURBED DELAY PARABOLIC INITIAL-BOUNDARY-VALUE PROBLEMS
More informationCholesky factorisation of linear systems coming from finite difference approximations of singularly perturbed problems
Cholesky factorisation of linear systems coming from finite difference approximations of singularly perturbed problems Thái Anh Nhan and Niall Madden Abstract We consider the solution of large linear systems
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationSingular Perturbation on a Subdomain*
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS, 937 997 ARTICLE NO. AY97544 Singular Perturbation on a Subdomain* G. Aguilar Departamento de Matematica Aplicada, Centro Politecnico Superior, Uniersidad
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationA uniformly convergent numerical method for singularly perturbed nonlinear eigenvalue problems
A uniformly convergent numerical method for singularly perturbed nonlinear eigenvalue problems Weizhu Bao and Ming-Huang Chai Abstract In this paper we propose a uniformly convergent numerical method for
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationThree remarks on anisotropic finite elements
Three remarks on anisotropic finite elements Thomas Apel Universität der Bundeswehr München Workshop Numerical Analysis for Singularly Perturbed Problems dedicated to the 60th birthday of Martin Stynes
More informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
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 informationFINITE DIFFERENCE APPROXIMATIONS FOR THE FIRST-ORDER HYPERBOLIC PARTIAL DIFFERENTIAL EQUATION WITH POINT-WISE DELAY
International Journal of Pure and Applied Mathematics Volume 67 No. 2, 49-67 FINITE DIFFERENCE APPROXIMATIONS FOR THE FIRST-ORDER HYPERBOLIC PARTIAL DIFFERENTIAL EQUATION WITH POINT-WISE DELAY Parameet
More informationThe continuity method
The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial
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 informationA SINGULARLY PERTURBED CONVECTION DIFFUSION PROBLEM WITH A MOVING INTERIOR LAYER
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 4, Pages 83 843 c 01 Institute for Scientific Computing and Information A SINGULARLY PERTURBED CONVECTION DIFFUSION PROBLEM WITH
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationSYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 12, 1998, 47 59 SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS M. Grossi S. Kesavan F. Pacella M. Ramaswamy
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in -D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf, and Matthias K. Gobbert 1 Department of Mechanical Engineering, University of Maryland, Baltimore County
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationOblique derivative problems for elliptic and parabolic equations, Lecture II
of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationA note on W 1,p estimates for quasilinear parabolic equations
200-Luminy conference on Quasilinear Elliptic and Parabolic Equations and Systems, Electronic Journal of Differential Equations, Conference 08, 2002, pp 2 3. http://ejde.math.swt.edu or http://ejde.math.unt.edu
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationNumerical integration of singularly perturbed delay differential equations using exponential integrating factor
MATHEMATICAL COMMUNICATIONS 51 Math. Commun. 017), 51 64 Numerical integration of singularly perturbed delay differential equations using exponential integrating factor Lakshmi Sirisha Challa and Yanala
More informationRegularity and nonexistence results for anisotropic quasilinear elliptic equations in convex domains
Regularity and nonexistence results for anisotropic quasilinear elliptic equations in convex domains Ilaria FRAGALÀ Filippo GAZZOLA Dipartimento di Matematica del Politecnico - Piazza L. da Vinci - 20133
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationThomas Apel 1, Ariel L. Lombardi 2 and Max Winkler 1
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANISOTROPIC MESH REFINEMENT IN POLYHEDRAL DOMAINS: ERROR ESTIMATES WITH DATA IN
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationHIGH ORDER VARIABLE MESH EXPONENTIAL FINITE DIFFERENCE METHOD FOR THE NUMERICAL SOLUTION OF THE TWO POINTS BOUNDARY VALUE PROBLEMS
JOURNAL OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN (p) 2303-4866, ISSN (o) 2303-4947 www.imvibl.org /JOURNALS / JOURNAL Vol. 8(2018), 19-33 DOI: 10.7251/JMVI1801019P Former BULLETIN OF THE
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationReview Article Robust Numerical Methods for Singularly Perturbed Differential Equations: A Survey Covering
International Scholarly Research Network ISRN Applied Mathematics Volume 2012, Article ID 379547, 30 pages doi:10.5402/2012/379547 Review Article Robust Numerical Methods for Singularly Perturbed Differential
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More information