1 Introduction
1.1 Terminology and classification of PDEs
A differential eqation is an equation satisfied by a function u, which involves besides u its derivates as well. If u depends only on one variable, e. g., the time t, we call the equation ordinary differential equation (ODE). n If it depends on several variables, e. g., on z = (z1, z2, ..., zn) ∈ R , then there occur partial derivatives
∂u(z) ∂iu(z) := , ∂xi and we call the equation partial differential equation (PDE). n A general second-order, linear PDE for a function u(z), z ∈ R is
n n X X − ∂iaij(z)∂ju + bi(z)∂iu + c(z)u = f(z), (1.1) i,j=1 i=1
In case aij(z), bi(z) and c(z) are independent of z, we have a PDE with constant coeffi- cients. d Note, that for time-dependent problems in R we have z = (x, t) and n = d + 1 where d x ∈ R . For time-independent problems z = x. n If we are searching for u defined in the open set O ⊂ R then we need the regularity 2 1 0 assumption u ∈ C (O), aij ∈ C (O), bi, c, f ∈ C (O) such that all derivatives exist in classical sense. These regularity assumptions will be reduced later when we will formulate the PDE in weak sense. 2 For u ∈ C (O) the partial derivates can be switched, i. e., ∂i∂j = ∂j∂i. If aij are not symmetric we can symmetrize the coefficients by
new orig orig aij := (aij + aji )/2, and by adjusting the remaining coefficients bi such that the form of the PDE remains (1.1). n So, we can assume that the coefficient matrix A(z) = {aij(z)}i,j=1 is symmetric. PDEs are classified into
• In elliptic equations an incident at a point z influences all point in their neigh- bourhood.
5 • In parabolic equations there is in one direction from a point z for which the influence is only for larger values. This is often the time direction, where only later times are influenced.
• In hyberbolic equations there are areas of directions around a point z which are influenced.
The principal part (Hauptteil) of the PDE
n X − ∂iaij(z)∂ju i,j=1 is mainly responsible for their classification.
Definition 1.1 (Classification of second-order PDEs). Consider a second-order, linear PDE of the form (1.1) with the symmetric coefficient matrix A(z).
1. The equation is said to be elliptic at z ∈ O if the eigenvalues of A(z) are all non-zero and of same sign, so e. g., for positive-definite A(z).
2. The equation is said to be parabolic at z ∈ O if one eigenvalue of A(z) vanishes whereas the others are all non-zero and of same sign, e. g., A(z) is positive semi- definite, but not positive definite, and the rank of (A(z), b(z)) is equal to n.
3. The equation is said to be hyberbolic at z ∈ O if A(z) has only non-zero eigenval- ues, whereas n − 1 of one sign and one of the other sign.
A partial differential equation (1.1) is said to be elliptic, parabolic or hyberbolic in a set O of points z, if it has the property for all z ∈ O.
The classification covers most (linear) physical models , but is not complete.
Remark 1.2. Do not mix second-order elliptic PDEs with elliptic bilinear forms, which we will discuss later.
1.2 Partial differential equations as mathematical models
A mathematical model is the description of a system using mathematical language. Mathematical models are obtained by a combination of first principles and constitutive equations. First principles are law of nature like conservative equations, actio = reactio, etc. They are based on basic assumptions. Constitutive equations represent material properties and are based often on measure- ments. They include empirical parameters.
6 1.2.1 Heat transfer equation We have as first principle the convervation of energy which correspond in absence of work to the conservation of temperature ∂u (x, t) + div j(x, t) = f(x, t) (1.2) ∂t where u → temperature [u] = 1K W j → heat flux [j] = 1 m2 W f → heat source/sink [f] = 1 m3 and as constitutive equation Fourier’s law
j(x, t) = −A(x) grad u(x, t), (1.3) saying, that the heat flux is proportional to the temperature gradient. The matrix A represents the heat conductivity tensor (Leitf¨ahigkeitsmatrix). We call the material
homogen : A(x) = A, isotrop : A(x) = α(x) · 1, and inhomogen or anisotrop otherwise. α(x) is the conductivity. The heat equation is parabolic if A is positive definite. For a unique solution we have to complete the system by initial conditions
u(x, 0) = u0(x), and boundary conditions. We may describe the temperature on the boundary of the domain of interest Ω
u(x, t) = g(x, t) for x ∈ ∂Ω or the heat flux
j(x, t) · n(x) = h(x, t) for x ∈ ∂Ω.
Here, n(x) is the normalised outer normal vector.
1.2.2 Diffusion and Poisson equation ∂u The static limit of the heat equation, i. e., with ∂t = 0 in (1.3), we get the elliptic system
j = −A(x) grad u , (FL) div j = f, (EL) which has to be equipped with the above boundary conditions. The same system arises in
7 • Electrostatics: u → electric potential [u] = 1V As j → displacement current (D)[j] = 1 m2 As f → charge density (ρ)[f] = 1 m3 Here, A stands for the dielectric tensor, which is usually designated by . The relationship (EL) is Gauss’ law, and (FL) arises from Faraday’s law curl E = 0 and the linear constitutive law D = E. • Stationary electric currents: u → electric potential [u] = 1V A j → electric current [j] = 1 m2 Inside conductors, where tensor A represents the conductivity. The source term f usually vanishes and excitation is solely provided by non-homogeneous boundary conditions. In this context (FL) arises from Ohm’s and Amp`erescircuit law and (EL) is a consequence of the conservation of charge. A similar elliptic system is j = −A(x) grad u , (FL) div j = f − c(x)u, (EL) which appear in mol • Molecular diffusion: u → concentration [u] = 1 m3 mol j → flux [j] = 1 m2s mol f → production/consumption rate [f] = 1 m3s Here A stands for the diffusion constant and, if non-zero, c denotes a so-called reaction coefficient. The equation (EL) guarantees the conservation of total mass of the relevant species.
1.2.3 Wave equation For the vertical displacement of a thin membran can be modelled by the first principle, mass times acceleration is the force, ∂2u m(x) (x, t) = div σ(x, t) + f(x, t), ∂t2 and Hook’s law as constitutive equation σ(x, t) = A(x) grad u(x, t). Here, we have u → vertical displacement [u] = 1m σ → stress vector [σ] = 1J f → external force [f] = 1N grad u → (local) deformation of the membran ∂u ∂t → (local) speed of the membran ∂2u ∂t2 → (local) acceleration kg and the density of mass m(x) ([m] = 1 m2 ).
8 Equation Modelled system Classification Advection equation Transport of particles in a flow. Hyberbolic Navier-Stokes equations Flow of fluids. Hyberbolic Euler Flow of incompressible fluids. Hyberbolic Shallow water equations Flow in shallow domains (earth atmosphere). Hyberbolic Stokes equations Static limit of flow of viscous flows. Elliptic Elasticity equations Stresses in elastic bodies. Elliptic Maxwell’s equations Electromagnetic fields. Hyberbolic Time-harmonic Maxwell’s equations Electromagnetic fields. Elliptic Eddy current modell Electromagnetic fields in low frequency. Elliptic Schr¨odinger’sequation Electron structure of matters. Elliptic Black-Scholes equation Prices of stocks. Parabolic
Table 1.1: Some examples for systems modelled by PDEs.
The second-order PDE is ∂2u m(x) (x, t) = div(A(x) grad u) + f, (1.4) ∂t2 and for homogeneous isotrope material with A = 1 we have ∂2u m(x) (x, t) = ∆u + f, ∂t2 which is hyperbolic. We have to complete the PDE with initial conditions
u(x, 0) = u0(x), and boundary conditions. We may describe the displacement on the boundary of the membran u(x, t) = g(x, t) for x ∈ ∂Ω or the boundary stress σ(x, t) · n(x) = h(x, t) for x ∈ ∂Ω.
1.2.4 Time-harmonic wave-equation −i ωt −i ωt −i ωt Let f = Re(fbe ) as well as g = Re(gbe ), h = Re(bh e ) for some (angular) + frequency ω ∈ R . With the independance of the coefficient function of t we can write −i ωt u = Re(ube ) and insertion into (1.4) leads to the time-harmonic wave-equation 2 − div(A(x) grad ub) − ω m(x)ub = f,b (1.5) which is a linear elliptic PDE of second order. The boundary condition transfer to ub where g and h are replaced by gb and bh. Remark 1.3. Linear parabolic or hyberbolic PDEs of second-order get elliptic in their static limit or in their time-harmonic form.
9 1.2.5 Further models Further systems modelled by PDEs are given in Table 1.1.
1.3 Well-posedness
Definition 1.4 (Hadamard’s well-posedness). A problem is said to be well-posed (kor- rekt gestellt) if
1. it has a unique solution u,
2. the solution depends continuously on the given data f, i. e.,
kuk ≤ Ckfk with a constant independant of u and f. Otherwise the problem is ill-posed.
Corollary 1.5. Small pertubations in the data of linear well-posed problems leads to small pertubations of the solution.
Proof. Let fe = f + δf with kδfk 1 (very small). Then ue is the solution of the system with data fe and δu = ue − u the solution of the system with data δf. So kδuk ≤ Ckδfk 1.
Example 1.6 (Ill-posed problem). Let Ω = [0, 1]×[0, ∞) and u the solution of the PDE
∆u = 0 1 u(x , 0) = sin(nx ) 1 n 1 ∂2u(x1, 0) = 0, i. e., instead of boundary conditions on all boundaries we prescibe “initial condition” on the boundary x2 = 0. For n → ∞ the boundary data gets smaller and smaller, whereas the solution 1 u(x , x ) = sin(nx ) cosh(nx ) 1 2 n 1 2 explodes for a fixed x2 > 0.
10 1.4 Some norms and spaces
Repeating the definition of a vector space and a norm.
Definition 1.7 (Vector space). A vector space over R or C is a set X, whose elements are called vectors, for which the operations addition and scalar multiplication are defined, and for any vectors x, y ∈ X and λ ∈ R or C it holds
(i) x + y ∈ X, (ii) λx ∈ X, (iii) 0 ∈ X.
Definition 1.8 (Norm). Let X be a real (or complex) vector space. We call k·k : X → R a norm on X if
kxk = 0 ⇐⇒ x = 0, (definiteness), (N1) kxk ≥ 0 ∀x ∈ X, (N2) kλxk = |λ|kxk ∀λ ∈ R(or λ ∈ C), (homogeneity), (N3) kx + yk ≤ kxk + kyk ∀x, y ∈ X, (triangle inequality). (N4)
11 2 Finite differences scheme for second-order elliptic PDEs
2.1 The strong formulation
d We consider the problem the connected bounded open set Ω ⊂ R with boundary ∂Ω and the PDE with f ∈ C(Ω), g ∈ C(∂Ω), 0 ≤ c ∈ C(Ω)
Lu = −∆u + cu = f in Ω, (2.1a) u = g on ∂Ω. (2.1b)
Lemma 2.1 (Basic Maximum (minimum) principle for c = 0). Let u ∈ C2(Ω) ∩ C(Ω) be solution of (2.1) with c = 0. If f ≤ 0 (f ≥ 0) in Ω, then the maximum (minimum) of u in Ω is attained on the boundary ∂Ω. Furthermore, if the maximum (minimum) is attained at an interior point of Ω, then the function u is constant.
Proof. First, we carry out the proof for the stronger assumption f < 0 in Ω. Suppose that there exists some xe ∈ Ω (a maximum point in the interior) such that
u(xe) = sup u(x) > sup u(x). (2.2) x∈Ω x∈∂Ω
With
d 2 (2.1a) X ∂ u 0 > f(x) = (Lu)(x) = −∆u(x) = − (x) e e ∂x2 e i=1 i
12 which is a contradiction to (2.2), as for a maximum point for all 1 ≤ i ≤ d
∂2u 2 (xe) ≤ 0. ∂xi
Now, let f ≤ 0, and xe is again the maximum of u in the interior and (2.2) holds. As the u on the boundary is smaller than the maximal value we can find a sufficiently small β > 0 such that the function
d X 2 w = u + β (xi − xei) i=1 attains its maximum at an interior point x0. Since
Lw = fw = f − dβ < f the maximum of w cannot be attained at an interior point, and we have a contradiction as well, i. e.,(2.2) does not hold and u attains its maximum on the boundary. Let the maximum is denoted by M = maxx∈∂Ω. There might be still points xe in the interior with u(xe) = M. Assume that u is not constant in Ω. So there exists a ball B ⊂ Ω of positive radius R and mid-point x0 with supx∈B2 u(x) < M for any proper subset B2 ( B, but with a boundary point xe ∈ ∂B for which u(xe) = M. Let G := B\B(x0,R/2) the open bounded ring domain with outer boundary Γe and inner boundary Γ. Let furthermore
2 v(x) := e−λ|x−x0| > 0 with λ > 0 large enough such that
2 −λ|x−x0| Lv = λ(d − λ|x − x0|)e < 0 for all x ∈ G. Note, that
v(x) = 0 on Γe, (2.3) and so
wε := u + εv coincides with u on Γ.e Since supx∈Γ u(x) < M by assumption we can choose a ε > 0 so small such that
wε(x) − wε(xe) = wε(x) − M < 0 in G. (2.4)
Since Lwε = f + εLv < 0 the maximum of wε in G is on ∂G, which is by (2.3) and (2.4) only attained in xe. Hence,
grad wε(xe) = 0
13 and as grad v(xe) =6 0 it holds
grad u(xe) =6 0, i. e., xe is not a maximum point which contradicts the assumptions. So there exists no interior maximum point if u is not constant. For f ≥ 0 we consider ue := −u and apply the previous steps. Example 2.2. Consider an open bounded set Ω = (−1, 1)2 and the Poisson equation
−∆u = −4 < 0 in Ω,
2 2 which is solved by the function u = u(x1, x2) = x1 + x2 with appropiate boundary condi- tion. Since f ≤ 0 in Ω the solution u attains its maximum on the boundary ∂Ω. This is indeed true.
Remark 2.3 (Non-local influence of sources). Let f ≤ 0 and g = 0. As the maximum of the solution u of (2.1) is attained on the boundary, there exists no subdomain G ⊂ Ω with |G| > 0 where u ≡ 0. Any point in the domain Ω “feels” the source. Lemma 2.4 (Basic Maximum (minimum) principle for c ≥ 0). Let u ∈ C2(Ω) ∩ C(Ω) be solution of (2.1) with c ≥ 0. Let furthermore f ≤ 0 (f ≥ 0) in Ω. Then if u has a non-negative maximum (non-positive minimum) in Ω is attained on the boundary ∂Ω. Furthermore, if the positive maximum (negative minimum) is attained at an interior point of Ω, then the function u is constant. Corollary 2.5 (Comparison principle). Let u, v ∈ C2(Ω) ∩ C(Ω) solve the equations Lu = fu and Lv = fv, respectively, and
fu ≤ fv in Ω u ≤ v on ∂Ω.
Then, u ≤ v in Ω.
14 Proof. Let w := v − u and so Lw = fv − fu ≤ 0. If the minimum of w in Ω is positive, then w is positive in Ω and the corollary is true. If the minimum of w in Ω is not positive, then by the minimum principle (Lemma 2.4) it is attained on boundary. But w ≥ 0 on the boundary, i. e., the minimal value is not below 0, and so w ≥ 0 in Ω.
Corollary 2.6 (Uniqueness). Let u a solution of the boundary value problem (2.1). Then u is unique. Proof. Assume that v =6 u solves (2.1). Then, w := v − u solves (2.1) with f ≡ 0 and g ≡ 0. Applying the comparison principle we conclude that 0 ≥ w ≥ 0 which contradicts the assumption.
Definition 2.7 (Classical solution). If u ∈ C2(Ω) ∩ C(Ω) satisfy (2.1) (or even (1.1)) in a pointwise sense, and the prescribed boundary conditions, then these functions are called a classical solution of the boundary value problem. In general there does not exist a classical solution. A special case with a classical solution is that of the second-order elliptic PDE with smooth boundary and constant coefficient functions, e. g.,(2.1). Here, we can make even a statement about the regularity of the solution. k ∞ Lemma 2.8 (Elliptic shift theorem). Let f ∈ C (Ω), k ∈ N0, g ∈ C (∂Ω) and ∂Ω ∈ C∞. Let u the unique solution of (2.1). Then u ∈ Ck+2(Ω). Lemma 2.9 (Continuous extension theorem). Let u be the solution of (2.1) and u ≡ 0 in the open bounded subset G ⊂ Ω. Then, u ≡ 0 in Ω. Proof. See [2].
2.2 Approximation by the Finite Difference Method (FDM)
If the solution is smooth enough we can replace the derivatives by difference quotients.
2.2.1 Difference quotients n+1 Let h > 0 and x ∈ R. Let I = [x − h, x + h], u ∈ C for some n ≥ 0. Taylor’s theorem gives 2 3 n 0 h 00 h 000 (±h) (n) u(x ± h) = u(x) ± hu (x) + u (x) ± u (x) + ... + u (x) + Rn(u; x, h) 2 3! n! with the remainder Z x±h n+1 1 n (n+1) (±h) (n+1) Rn(u; x, h) = (x − t) u (t) dt = u (ξ) for some ξ ∈ I. n! x (n + 1)! The forward difference for u ∈ C3(I) is u(x + h) − u(x) h (D+u)(x) := = u0(x) + u00(x) + O(h2), h 2
15 whereas the backward difference is given by
u(x) − u(x − h) h (D−u)(x) := = u0(x) − u00(x) + O(h2). h 2
The central difference for u ∈ C5(I) is
u(x + h) − u(x − h) h2 (D0)(x) := = u0(x) + u000(x) + O(h4). 2h 3! For the second derivative u00(x) we get for u ∈ C4(I):
u(x + h) − 2u(x) + u(x − h) (D+D−u)(x) := = u00(x) + O(h2). h2 In two dimensions we have for u ∈ C4([−h, h] × [−h, h] + x) the five-point-stencil for the Laplacian
u(x1 + h, x2) + u(x1 − h, x2) + u(x1, x2 + h) + u(x1, x2 − h) − 4u(x1, x2) (∆u)(x) = +O(h2). h2 | {z } =: Λu
Let us also specify the non-equidistant finite difference for u ∈ C4(I)
2 u(x + h ) − u(x) u(x) − u(x − h ) h − h u00(x) = R − L − R L u000(x) + O(h2). hL + hR hR hL 3
The Laplacian can be approximated for u ∈ C3([−h, h] × [−h, h] + x) 2 u(x1 + h , x2) − u(x1, x2) u(x1, x2) − u(x1 − h , x2) (∆u)(x) = R − L hL + hR hR hL 2 u(x1, x2 + h ) − u(x1, x2) u(x1, x2) − u(x1, x2 − h ) + T − B + O(h). hB + hT hT hB
∗ We call the difference quotient Λ u, which generalises Λu (special case h = hL = hR = hT = hB).
2.3 The Finite Difference method
The PDE (2.1) is approximated on a uniform grid
Th := {x = (jh, ih): x ∈ Ω}, (2.5) whereas the boundary has the grid
Gh := {x = (jh, x2) or x = (x1, ih) for some x1, x2 ∈ R : x ∈ ∂Ω}.
16 Figure 2.1: Grid Th for a computational domain Ω and boundary grid Gh.
The union is T h := Th ∪ Gh (see Fig. 2.1), and |T h| its cardinality. −2 −1 The cardinalities of Th and Gh behave asymptotically like |Th| ∼ h and |Gh| ∼ h . We call a grid function a function vh ∈ T h. We ask the finite difference approximant uh to u to fulfill
(Lhuh)(x) = f(x) for all x ∈ Th, (2.6a)
uh(x) = g(x) for all x ∈ Gh, (2.6b) with
∗ (Lhuh)(x) := −Λ uh(x) + cuh(x),
∗ and Λ meaning the general difference quotient for ∆ with h ≥ hL, hR, hT , hB > 0 the > > > > minimal values such that x − (hL, 0) , x + (hR, 0) , x + (0, hT ) , x − (0, hB) ∈ T h.
2.3.1 Existence and uniqueness
Lemma 2.10 (Discrete maximum principle). Let f ≤ 0 (f ≥ 0) on Th. If the maximum (minimum) of the solution uh of (2.6) is non-negative (non-positive) it is attained on Gh.
Proof. Let Lhuh ≤ 0 for all x ∈ Th, and xe ∈ Th such that
uh(xe) = max uh(x) x∈T h and uh(xe) ≥ 0.
17 As consequence of the fact that xe is a maximum and an interior point it is 2 u (x1 + h , x2) − u (x1, x2) u (x1, x2) − u (x1 − h , x2) h e R e h e e − h e e h e L e ≤ 0, hL + hR hR hL 2 u (x1, x2 + h ) − u (x1, x2) u (x1, x2) − u (x1, x2 − h ) h e e T h e e − h e e h e e B ≤ 0, hB + hT hT hB and so
∗ 0 ≥ (Lhuh)(xe) = −(Λ uh)(xe) + c(xe)uh(xe) ≥ c(xe)uh(xe) ≥ 0
So, it remains only (Lhuh)(xe) = 0, which is possible for
c(xe) = 0 and uh(x) = uh(xe) or c(xe) > 0 and uh(x) = uh(xe) = 0 for all points x in the stencil of xe. > If xe + (hR, 0) ∈ GT the statement of the lemma is proved, otherwise we repeat the > previous steps for xe + (hR, 0) instead of xe until the (right) boundary of the domain is reached, and the lemma holds. For f ≥ 0 we consider ueh := −uh and apply the previous steps.
Corollary 2.11. The only solution of (2.6) with f ≡ 0 on Th and g ≡ 0 on Gh is uh ≡ 0.
Proof. As f ≤ 0 there can be only a non-negative maximum of uh on the boundary. In view of the boundary conditions uh ≤ 0. Vice-versa with f ≥ 0 we have uh ≥ 0 and the proof is complete.
Corollary 2.12 (Existence and uniqueness). The system of equations (2.6) provides a unique solution uh.
Proof. We have shown in Corollary 2.11 that the null space of the FDM matrix is empty, and as the matrix is finite-dimensional the range is the full and a solutions exists for any f and g.
Lemma 2.13 (Discrete comparison principle). Let uh,1 and uh,2 the solutions of (2.6) with f = f1, g = g1 or f = f2, g = g2, respectively. If for the data it hold the inequalities
|f1(x)| ≤ f2(x) ∀x ∈ Th,
|g1(x)| ≤ g2(x) ∀x ∈ Gh.
Then,
|uh,1(x)| ≤ uh,2(x) ∀x ∈ T h.
18 Proof. Let vh := uh,1 − uh,2, then
Lhvh = f1 − f2 ≤ 0 on Th,
vh = g1 − g2 ≤ 0 on Gh.
If there is a non-negative maximum, it is attained on Gh, and so vh(x) ≤ 0 for all x ∈ T h. Otherwise, the maximum is negative and so even vh(x) < 0 for all x ∈ T h. Repeating the same steps with wh := −uh,1 − uh,2 finishes the proof.
2.3.2 Convergence The finite difference method provides approximation on grid points only, therefore we can also only use norms on the grid to measure the error. Grid functions are measured by the `∞ norm
kvhk∞ := max |vh(x)| x∈T h or the weighted `2 norm defined by 1/2 1 X 2 kvhk2,h := p |vh(x)| ∼ hkvhk2. |Th| x∈T h The weight is to make the norm independent of h, so that we can compare different mesh widths. For example, the constant function 1 has same norm on meshes T h1 and
T h2 even if h1 =6 h2. Note, that this holds also true for meshes of different domains in 2 difference to the L -norm. We use the same norms also on Th. By the Cauchy-Schwarz inequality we have for any grid function vh
kvhk2,h ≤ kvhk∞. (2.7)
0 Let v ∈ C (Ω). Then, we call the grid function (Rhv) ∈ T h its restriction to T h. With the consistency error we measure the accuracy if the exact solution would be inserted in the numerical scheme. Definition 2.14 (Consistency error). Let u be the solution of (2.1). The consistency error is defined as
kLuh − Lhuk = kf − Lhuk with k · k some norm on the mesh Th. The consistency error is of order k if for h → 0 it can be bounded by C hk for some constant C > 0. Lemma 2.15 (Estimate of the consistency error). Let u ∈ C3+`(Ω) with ` = 0, 1 the solution of (2.1). There exist two constants C∞,C2 > 0 independent of h such that on grids Th the consistency error can be bounded as
kf − Lhuk∞ ≤ C∞h, 1+`/2 kf − Lhuk2,h ≤ C2 h .
19 Proof. Since for any x ∈ Th
∗ f − Lhu = −∆u + Λ u, we have to bound the error of the finite difference approximation of the Laplacian. The non-equidistant finite difference quotient Λ∗, which is needed close to the bound- ary, approximates the Laplacian point-wise with O(h), and dominates so the consistency error in the `∞ norm. For u ∈ C4(Ω) we have on O(h−2) points with the difference quotient Λ an point-wise error of O(h2) and on O(h−1) points with the difference quotient Λ∗ a point-wise error of O(h) and so
∗ 2 2 −2 4 −1 2 3 k∆u − Λ uk2,h ≤ O(h ) · O(h )O(h ) + O(h )O(h ) = O(h ).
The error on the boundary dominates again. If u ∈ C3(Ω) the pointwise error is anyway O(h) and we loose half an order.
As the numerical scheme works with finite precision we have to consider the stability of the scheme, i. e., the influence of small pertubation of the data to the solution, as well.
Lemma 2.16 (Stability of the FDM in the `∞-norm). There exists a constant C > 0 only depending on the size of Ω such that for the solution uh of (2.6) it holds
diam(Ω)2 ku k∞ ≤ kfk∞ + kgk∞. h 4
Proof. Without loss of generality let 0 ∈ Ω and let R = diam(Ω). Then, for any constants α, β > 0 the grid function
2 2 2 vh = α(R − x1 − x2) + β (2.8)
∗ is positive for all x ∈ T h. Inserting (2.8) into Λ we have
∗ ∗ 2 2 (Λ vh)(x) = −Λ (x1 + x2) 2 2 2 2 2 (x1 + hR) − x1 x1 − (x1 − hL) = −α 2 − 2 hL + hR hR hL 2 2 2 2 2 (x2 + hT ) − x2 x2 − (x2 − hB) + 2 − 2 hT + hB hT hB = −4 α < 0.
Thus, vh is by Corollary 2.11 the unique solution of