Partial differential equations and the finite element method 2--SECOND-ORDER ELLIPTIC PROBLEMS

This part is devoted to the discussion of linear second-order elliptic problems.

Consider the model equation
(1)

$$-\nabla (a_1\nabla u)+a_0u=f,$$
with homogeneous Dirichlet boundary conditions
(2) $$u(x)=0$$ on $\partial \Omega.$

Assuming $a_{ij}(x)=a_1(x)\delta_{ij}$, and $b(x)=c(x)=0$.

Classical solution: The solution of (1) is a function $u \in C^2(\Omega)\cap C ( \overline \Omega )$ satisfying the equation (1) everywhere in $\Omega$ and fulfilling the boundary condition (2)
at every $x\in \partial \Omega$.

$C ( \overline \Omega )$=$\{\phi\in C(\Omega)$, $\phi$ is bounded and uniformly continuous function $\}$: because $\overline \Omega$ is closed, the functions in closed sets are uniformly continuous.

Find the solution of (1) using the weak formulation:

1. Multiply (1) with a test function $v\in C_0^{\infty}(\Omega)$
   $$-\nabla (a_1\nabla u)v+a_0uv=fv,$$
2. Integrate over $\Omega$
   $$\int_{\Omega}-\nabla (a_1\nabla u)vdx+\int_{\Omega}a_0uv dx=\int_{\Omega}fv dx,$$

3. reduce the maximum order of the partial derivatives

      By Green's first identity:
   

   <f,−Δx>Ω=<∇f,∇x>Ω−<f,∂nx>∂Ω < f , − Δ x > Ω =< ∇ f , ∇ x > Ω − < f , ∂ n x > ∂ Ω
   <script type="math/tex; mode=display" id="MathJax-Element-30"> _{\Omega}=<\nabla f,\nabla x>_{\Omega}- _{\partial\Omega}</script>
   we have
   $$\int_{\Omega}-\nabla (a_1\nabla u)vdx=\int_{\Omega} \nabla a_1 u \nabla v dx-\int_{{\partial\Omega}}v\partial_n (a_1 u) dS$$
   The fact that $v$ vanishes on the boundary ${\partial\Omega}$ removes the boundary term, then
   (3) $$\int_{\Omega} \nabla a_1 u \nabla v dx+\int_{\Omega}a_0uv dx=\int_{\Omega}fv dx,$$

4. Find the largest possible function spaces for $u, v$ s.t.all integrals are finite
   Originally, identity (3) was derived under very strong regularity assumptions $u \in C^2(\Omega)\cap C ( \overline \Omega )$ and $v\in C_0^{\infty}(\Omega)$.
   All integrals in (3) remain finite when these assumptions are weakened to
   $$u,v\in W_0^{1,2}(\Omega), f\in L^2(\Omega)$$
   Similarly the regularity assumptions for the coefficients $a_1$ and $a_0$ can be reduced to
   $$a_1,a_0\in L^{\infty}(\Omega)$$
   Note: $W^{m,p}=\{u\in L^P(\Omega):D^{|a|}u\in L^P(\Omega): \forall 0<|a|\leq m\}$
   $W_0^{m,p}$ contains the functions in $W^{m,p}$ that have compact support in $\Omega$.

The weak form of the problem (1)-(2) is stated as follows: Given $f\in L^2(\Omega)$, find a
function $u\in W_0^{1,2}(\Omega)$ such that

$$\int_{\Omega} \nabla a_1 u \nabla v dx+\int_{\Omega}a_0uv dx=\int_{\Omega}fv dx,\forall v\in W_0^{1,2}(\Omega)$$

Let $V =W_0^{1,2}(\Omega)$. We define a bilinear form $a(., .): V\times V\rightarrow R$ :

$$a(u,v)=\int_{\Omega} \nabla a_1 u \nabla v dx+\int_{\Omega}a_0uv dx$$
and a linear form $l\in V'$
$$l(v)=<l,v>=\int_{\Omega}fv dx$$
Then the weak formulation of the problem (1), ( 2) reads: Find a function $u\in V$ such that
$$a(u,v)=l(v)$$
This notation is common in the study of partial differential equations and finite element methods.

Bilinear forms, energy norm, and energetic inner product

Every bilinear form $a : V \times V\rightarrow R$ in a Banach space $V$ is associated with a unique linear operator $A : V \rightarrow V’$ defined by

(4)

$$(A(u))v=<Au,v>=a(u,v)$$
Note: (4) defines a one-to-one correspondence between continuous bilinear forms $a : V \times V\rightarrow R$ and linear continuous operators $A : V \rightarrow V’$.

Definition 1 (Energetic inner product, energy norm) Let $V$ be a Hilbert space and
$a : V \times V\rightarrow R$ a bounded symmetric V-elliptic bilinear form. The bilinear form defines
an inner product
(5)

$$(u,v)_e=a(u,v)$$
in $V$ , called energetic inner product. The norm induced by the energetic inner product,
(6) $$\parallel u\parallel_e=\sqrt{a(u,u)}$$
is called energy norm.

Theorem 1.5 (Lax-Milgram lemma) Let $V$ be a Hilbert space, $a : V \times V\rightarrow R$ a bounded
V-elliptic bilinear form and $l\in V'$. Then there exists a unique solution to the problem

$$a(u,v)=l(v).$$

Note: the unique solution is the unique represent of the linear form $l\in V'$ with respect to the energetic inner product$(u,v)_e=a(u,v)$. In this sense the Lax-Milgram lemma is a special case of the Riesz representation theorem.

---

Riesz representation theorem.

This theorem establishes an important connection between a Hilbert space and its (continuous) dual space.

Let $H$ be a Hilbert space, and let $H^*$ denote its dual space, consisting of all continuous linear functionals from $H$ into the field ${\mathbb {R} }$ or $\mathbb {C}$ . If $x$ is an element of $H$, then the function $\varphi _{x}$, for all $y$ in $H$ defined by:

$${\displaystyle \varphi _{x}(y)=\left\langle y,x\right\rangle ,}$$

where ${\displaystyle \langle \cdot ,\cdot \rangle }$ denotes the inner product of the Hilbert space, is an element of $H^*$.

Note: The Riesz representation theorem states that every element of $H^*$ can be written uniquely in this form. Given any continuous linear functional $g \in H^*$, the corresponding element ${\displaystyle x_{g}\in H}$ can be constructed uniquely by ${\displaystyle x_{g}=g(e_{1})e_{1}+g(e_{2})e_{2}+...}$ , where ${\displaystyle \{e_{i}\}}$ is an orthonormal basis of $H$, and the value of${\displaystyle x_{g}}$ does not vary by choice of basis. Thus, if ${\displaystyle y\in H,y=a_{1}e_{1}+a_{2}e_{2}+...}$, then ${\displaystyle g(y)=a_{1}g(e_{1})+a_{2}g(e_{2})+...=\langle x_{g},y\rangle }$.

---

Lemma 1. Assume that $a_1(x)\geq C_{min}>0$ and $a_0(x)\geq 0$ a.e. in $\Omega$. Then the weak
problem has a unique solution $u\in V$ .

Nonhomogeneous Dirichlet boundary conditions

More general Dirichlet boundary conditions

$$u(x)=g(x), x\in\partial \Omega$$
and $g\in C(\partial\Omega)$. Consider a new function $G\in C^2(\Omega)\cap C(\overline\Omega)$, so that $g=G$ on $\partial \Omega$.

Note: $G$ is not unique, but we will show later that the solution is invariant under its choice.

Denote u=G+U, then the problem (1) can be rewrited as finding $U\in C_0^2(\Omega)$

$$-\nabla (a_1\nabla (G+U))+a_0(G+U)=f, x\in \Omega$$
s.t.
$$G+U=g, x\in \partial \Omega$$
or, equivalently,
$$-\nabla (a_1\nabla U)+a_0U=f+\nabla (a_1\nabla G)-a_0G, x\in \Omega$$
s.t.
$$U=0, x\in \partial \Omega$$

Independence of the solution $u = U + G$ on the Dirichlet lift $G$:
Assume that
$U_1 + G_I = u_1 \in W_0^{1,2}(\Omega)$ and $U_2 + G_2 = u_2 \in W_0^{1,2}(\Omega)$ are two weak solutions.
So,
a(u1,v)−a(u2,v)=0

a ( u 1 , v ) − a ( u 2 , v ) = 0

$$a(u_1,v)-a(u_2,v)=0$$
set $v=u_1-u_2$, then
$$0=a(u_1-u_2,u_1-u_2)\geq C\parallel u_1-u_2\parallel$$
that is $u_1=u_2$.

Neumann boundary conditions

Neumann boundary conditions of the form

$$\frac{\partial u}{\partial \nu}=g, x \in\partial \Omega$$

1.Assume $u\in$ Multiply (1) with a test function $v\in C_0^{\infty}(\Omega)$

$$-\nabla (a_1\nabla u)v+a_0uv=fv,$$
2. Integrate over $\Omega$
$$\int_{\Omega}-\nabla (a_1\nabla u)vdx+\int_{\Omega}a_0uv dx=\int_{\Omega}fv dx,$$
3. reduce the maximum order of the partial derivatives

       By Green's first identity:

<f,−Δx>Ω=<∇f,∇x>Ω−<f,∂nx>∂Ω

< f , − Δ x > Ω =< ∇ f , ∇ x > Ω − < f , ∂ n x > ∂ Ω

<script type="math/tex; mode=display" id="MathJax-Element-124"> _{\Omega}=<\nabla f,\nabla x>_{\Omega}- _{\partial\Omega}</script>
we have
$$\int_{\Omega}-\nabla (a_1\nabla u)vdx=\int_{\Omega} \nabla a_1 u \nabla v dx-\int_{{\partial\Omega}}v\partial_n (a_1 u) dS$$
As $v$ does not vanish on the boundary ${\partial\Omega}$, we have
$$\int_{\Omega} \nabla a_1 u \nabla v dx+\int_{\Omega}a_0uv dx-\int_{{\partial\Omega}}a_1 v\frac{\partial u }{\partial \nu}dS=\int_{\Omega}fv dx,$$

Using the boundary condition, rewrite the problem with Neumann condition: given $f\in L^2(\Omega)$, $g\in L^2(\partial\Omega)$, find $u$ s.t.

$$\int_{\Omega} \nabla a_1 u \nabla v dx+\int_{\Omega}a_0uv dx=\int_{\Omega}fv dx+\int_{{\partial\Omega}}a_1 vgdS,$$
4. Find the largest possible function spaces for $u, v$ s.t.all integrals are finite

All integrals remain finite when these assumptions are weakened to, as $u,v$ need not to vanish on the boundary

$$u,v\in W^{1,2}(\Omega), f\in L^2(\Omega)$$

Remark (Essential and natural boundary conditions)
1. Dirichler boundary conditions are sometimes called essential since they essentially influence the weak formulation: They determine the function space in which the solution is sought.
2. Neumann boundary conditions do not influence the function space and can be naturally incorporated into the boundary integrals. Therefore they are called natural.

Newton (Robin) boundary conditions

boundary conditions involves a combination of function values and normal derivatives.

---

Energy of elliptic problems

In this part we introduce the explicit form of the abstract energy, at least for symmetric problems. The most important numerical scheme based on the minimization of the abstract energy.

Theorem. Let $V$ be a linear space, $a : V\times V \rightarrow R$ a symmetric V-elliptic bilinear
form and $l\in V’$. Then the functional of abstract energy,
(T.1) $$E(v)=\frac{1}{2}a(v,v)-l(v)$$
attains its minimum in $V$ at an element $u\in V$ if and only if
(T.2) $$a(u,v)= l ( v ),\forall v\in V$$
Moreover; the minimizer $u\in V$ is unique.

Proof: Let (T.2) holds. Then

$$E(u+tv)=a(u+tv,u+tv)/2-l(u+tv)=(u+tv,u+tv)_e/2-l(u+tv)=(u,u)_e/2+t^2(v,v)_e/2+t(u,v)_e-l(u)-tl(v)$$

$$E(u+tv)=E(u)+t^2(v,v)_e/2>E(u)$$ .

if $E$ has a minimum at $u\in V$, then for every $v\in V$ the derivative of the $\phi(t)=E(u+tv)$ must vanish at 0

$$0=\phi'(0)=\frac{E(u+tv)}{\partial t}|_{t=0}=(u,v)_e-l(v)$$

…………