Numerical Solutions of PDE by FEM

# **Numerical Solutions of** 
# **Partial Differential Equations by** 
# **Finite Element Method**

</div>
</div>

# **Ch. 1, Sec. 2**
# FEM for the model problem 
# with piecewise linear functions

</div>
</div>

## 回顾

在 1.1 节中，我们已经讨论了如下的模型问题：

$$\begin{aligned}
-u''(x) = f(x)&, \quad x \in (0,1), \\\\
u(0) = 0&, \quad u(1) = 0.
\end{aligned}$$

在接下来的部分中，我们将构造有限的子空间来近似该问题的解（Ch. 1, Sec. 2, 本节），并分析所得近似解的误差（Ch. 1, Sec. 3）。

## 分区与节点

我们将区间 $[0,1]$ 分成 $M+1$ 个子区间：

$$0 = x_0 < x_1 < x_2 < \cdots < x_{M+1} = 1,$$

其中 $x_i$ 称为节点，区间为$I_{j}=(x_{j-1}, x_j)$，$h_j = x_j - x_{j-1}$ 称为第 $j$ 个子区间的长度。

定义最大区间长度为
$$h = \max_{1 \leq i \leq M} h_i.$$

我们使用这个参数 $h$ 来衡量分区的精细程度的。

## 分段线性函数空间

在1.1节中，我们定义了如下的函数空间 $V$：

> $V=${$v: v$在$[0,1]$上是连续的函数，$v'$在$[0,1]$上是分段连续且有界的，并且$v(0)=v(1)=0$}

</div>
<div class="col">

我们现在构造一个有限维子空间 $V_h$，它由分段线性函数组成，并且是上述空间 $V$ 的子空间。$V_h$ 为如下的函数集合：在每个子区间 $I_j$ 上，$v$ 是线性的；在区间 $[0,1]$ 上，$v$ 是连续的；并且满足边界条件 $v(0) = v(1) = 0$。

这就是满足给定边界条件的分段线性函数空间。如果去除边界条件，则可以得到在 $[0,1]$ 上的给定分段的分段线性函数空间。

</div>
</div>

## 分段线性函数的基函数

在 $V_h$ 中，我们可以定义如下的基函数 $\\{\phi_i\\}_{i=1}^{M}$ ：

$$\phi_i(x) = \begin{cases}
\frac{x - x_{i-1}}{h_i}, & x \in [x_{i-1}, x_i], \\\\
\frac{x_{i+1} - x}{h_{i+1}}, & x \in [x_i, x_{i+1}], \\\\
0, & \text{otherwise}.
\end{cases}$$

这些基函数 $\phi_i$ 在节点 $x_i$ 处取值为 $1$，在其他节点处取值为 $0$。因此，任意的 $v \in V_h$ 都可以表示为
$$v(x) = \sum_{i=1}^{M} v(x_i) \phi_i(x).$$

## 基函数的简单示例

![](CODE/FEM_01_Raw_Basis.png)

</div>

如图所示，基函数在每个节点处呈现出帽子形状。不相邻的节点处，基函数的值为零。

## FEM 近似解的表示

所以，对于 $v \in V_h$，我们有

$$v(x) = \sum_{i=1}^{M} \eta_{i} \phi_i(x).$$

其中，$$\eta_i = v(x_i).$$

每一个 $v \in V_h$ 都可以唯一地表示为基函数 $\phi_i$ 的线性组合。因此，$V_h$ 是一个维数为 $M$ 的线性空间，其基为 $\{\phi_i\}_{i=1}^{M}$。

## FEM 近似解的表示

![加权基函数](CODE/FEM_02_Weighted_Basis.png)

</div>

## FEM 近似解的表示

![最终结果](CODE/FEM_03_Final_Result.png)

</div>

## FEM 近似解的表示

于是，针对给定边界问题的有限元法可以改写成如下的形式：

> 对于函数空间 $V$ 和 $V_h$，我们有如下的问题：
> $$
> \begin{aligned}
> (M\ )&: Find \ u \in V \ such \ that \ F(u)\le F(v)&, \ &\forall \ v \in V. \\\\
> (M_h)&: Find \ u_h \in V_h \ such \ that \ F(u_h)\le F(v)&, \ &\forall \ v \in V_h.
> \end{aligned}
> $$
>

> 或者等价地，可以写成如下的弱形式：
> $$
> \begin{aligned}
> (V\ )&: Find \ u \in V \ such \ that \ (u',v')=(f,v)&, \ &\forall \ v \in V. \\\\
> (V_h)&: Find \ u_h \in V_h \ such \ that \ (u_h',v')=(f,v)&, \ &\forall \ v \in V_h.
> \end{aligned}
> $$
>

## FEM 近似解的表示

我们接下来针对 $(V_h)$ 问题进行讨论。如之前所示，由于 $V_h$ 的维数为 $M$，我们可以将 $u_h$ 表示为基函数的线性组合：
$$u_h(x) = \sum_{i=1}^{M} \eta_i \phi_i(x).$$
将 $u_h$ 代入 $(V_h)$ 问题中，并选择测试函数 $v$ 为每一个基函数 $\phi_j$，我们可以得到如下的线性方程组：
$$\sum_{i=1}^{M} \eta_i (\phi_i', \phi_j') = (f, \phi_j), \quad j = 1, 2, \ldots, M.$$

## FEM 近似解的表示

于是这个线性方程组可以写成矩阵形式为
$$A \eta = b,$$
其中矩阵 $A$ 的元素为 $A_{ji} = (\phi_i',   \phi_j')$，向量 $b$ 的元素为 $b_j = (f, \phi_j)$，向量 $\eta$ 包含了未知系数 $\eta_i$。

具体地说，线性方程组 $A\eta=b$ 形如：
$$
A\eta=\begin{bmatrix}
(\phi_1', \phi_1') & (\phi_2', \phi_1') & \cdots & (\phi_M', \phi_1') \\\\
(\phi_1', \phi_2') & (\phi_2', \phi_2') & \cdots & (\phi_M', \phi_2') \\\\
\vdots & \vdots & \ddots & \vdots \\\\
(\phi_1', \phi_M') & (\phi_2', \phi_M') & \cdots & (\phi_M', \phi_M')
\end{bmatrix}
\begin{bmatrix}
\eta_1 \\\\
\eta_2 \\\\
\vdots \\\\
\eta_M
\end{bmatrix}=
\begin{bmatrix}
(f, \phi_1) \\\\
(f, \phi_2) \\\\
\vdots \\\\
(f, \phi_M)
\end{bmatrix}
=b.
$$

在多维情形下，FEM 近似解的表示也是类似的。
$A_{ij}=(\nabla\phi_i, \nabla\phi_j)$，$b_j=(\phi_j,f)$，并且 $u_h$ 仍然可以表示为基函数的线性组合。

## 矩阵 $A$ 的结构

对于这个给定的问题，矩阵 $A$ 是三对角矩阵的，也就是说，只有主对角线和相邻的两条对角线上的元素是非零的。这是因为基函数 $\phi_i$ 和 $\phi_j$ 只有在它们对应的节点 $x_i$ 和 $x_j$ 相邻时，它们的导数才会有重叠的支持。而且，矩阵 $A$ 是对称的。

对于这个给定的问题，

$$
A_{ij} = \begin{cases}
\frac{1}{h_i} + \frac{1}{h_{i+1}}, & i = j, \\\\
-\frac{1}{h_{i+1}}, & j = i + 1, \\\\
-\frac{1}{h_i}, & j = i - 1, \\\\
0, & \text{otherwise}.
\end{cases}
$$

</div>
<div class="col">

均匀划分的情况下，矩阵 $A$ 的形式为：
$$
A = \frac{1}{h}\begin{bmatrix}
2 & -1 & 0 & \cdots & 0 \\\\
-1 & 2 & -1 & \cdots & 0 \\\\
0 & -1 & 2 & \cdots & 0 \\\\
\vdots & \vdots & \vdots & \ddots & \vdots \\\\
0 & 0 & 0 & \cdots & 2
\end{bmatrix}.
$$

再利用 $b_j = (f, \phi_j)$，我们就可以组装出完整的线性方程组 $A\eta = b$，并最终求出数值解。

</div>
</div>

<span style="color: blue;">在下一节中，我们将讨论有限元法近似解的误差问题。</span>

> - 在给定的 $V_h$ 中，我们可不可以找到更好的近似解？
> - 有限元法的近似解与真实解之间的误差有多大？
> - 随着分区的细化，误差会如何变化？
> - $\cdots$

</div>
</div>

# **Ch. 1, Sec. 3**
# An error estimate for FEM
# for the model problem

</div>
</div>

## 回顾

在上一节中，我们讨论了如下的模型问题：

$$\begin{aligned}
-u''(x) = f(x)&, \quad x \in (0,1), \\\\
u(0) = 0&, \quad u(1) = 0.
\end{aligned}$$

并给出了该问题的有限元法近似解的表示方法，即
$$u_h(x) = \sum_{i=1}^{M} \eta_i \phi_i(x),$$
其中 $\phi_i$ 是分段线性基函数，系数 $\eta_i$ 通过解线性方程组 $A\eta = b$ 得到。

> 在本节中，我们将分析有限元法近似解 $u_h$ 在其空间 $V_h$ 中的最佳性，并给出误差估计。

## 有限元法近似解的最佳性

我们已经知道，有限元法近似解 $u_h$ 满足如下的弱形式：

$$(u_h', v') = (f, v), \quad \forall v \in V_h.$$

同时，真实解 $u$ 满足：

$$(u', v') = (f, v), \quad \forall v \in V.$$

因此，对于任意的 $v \in V_h$，我们有：

$$(u' - u_h', v') = 0.$$

这表明误差 $u - u_h$ 在 $V_h$ 中是正交于所有测试函数 $v$ 的。

## 在给定函数空间中的最佳性

> 对于任意的 $v \in V_h$，我们有：$ \|\|(u-u_{h})'\|\|\leq \|\|(u-v)'\|\|. $
>

**证明：**

$$
\begin{aligned}
\|\|(u - v)'\|\|^2 &= ( (u - v)', (u - v)' ) \\\\
&= ( (u - u_h + u_h - v)', (u - u_h + u_h - v)' ) \\\\
&= \|\|(u - u_h)'\|\|^2 + 2( (u - u_h)', (u_h - v)' ) + \|\|(u_h - v)'\|\|^2 \\\\
&= \|\|(u - u_h)'\|\|^2 + \|\|(u_h - v)'\|\|^2 \\\\
&\geq \|\|(u - u_h)'\|\|^2.
\end{aligned}
$$

于是，我们得到了有限元法近似解在其空间 $V_h$ 中的最佳性。也就是说在 $V_h$ 中，没有其他函数能比 $u_h$ 更接近真实解 $u$。

## 与真解之间的误差分析

通过估计 $\|\|(u - v)'\|\|$，我们可以得到有限元法近似解 $u_h$ 与真实解 $u$ 之间的误差估计。我们选择$u_{h} \in V_h$为 $u$ 的插值函数，即在每个节点 $x_i$ 处，$u_h(x_i) = u(x_i)$。

如果我们采用这样的线性差值函数 $u_h$，则可以证明如下的误差估计：

> 对于给定的 $[0,1]$ 区间划分，我们有这样的误差估计：
> 
> - $\displaystyle\|u'(x)-u_{h}'(x)\| \leq h\max_{0\leq x \leq 1}\|u'(x)\|,$
> - $\displaystyle\|u(x)-u_{h}(x)\| \leq \frac{h^{2}}{8}\max_{0\leq x \leq 1}\|u''(x)\|,$
>
> 证明在数学分析课程中已经给出了。
>

于是可以推出 $\displaystyle\|\|u'(x)-u_{h}'(x)\|\| \leq h\max_{0\leq x \leq 1}\|u''(x)\|$ .

## 与真解之间的误差分析

结合有限元法近似解的最佳性，我们可以得到如下的误差估计：

> 对于有限元法近似解 $u_h$，我们有如下的误差估计：
> - $\displaystyle\|u'(x)-u_{h}'(x)\| \leq h\max_{0\leq x \leq 1}\|u''(x)\|,$
> - $\displaystyle\|u(x)-u_{h}(x)\| \leq \frac{h^{2}}{8}\max_{0\leq x \leq 1}\|u''(x)\|,$
> - $\displaystyle\|\|u'(x)-u_{h}'(x)\|\| \leq h\max_{0\leq x \leq 1}\|u''(x)\|.$
> 
> 这表明随着分区的细化（即 $h$ 的减小），有限元法近似解 $u_h$ 会更接近真实解 $u$。
>

这给我们提供了有限元法在数值求解偏微分方程中的有效性和收敛性的理论基础。我们并不需要解出真实解 $u$，就可以通过误差估计来评估有限元法近似解的质量。

# **Ch. 1, Sec. 4**
# TODO

</div>
</div>

# **Ch. 1, Sec. 5**
# TODO

</div>
</div>

# 敬请老师和同学们批评指正！

<hr/>

参考资料：

[Numerical Solution of Partial Differential Equations<br>by the Finite Element Method, Claes Johnson](https://books.google.com/books?id=0IFCAwAAQBAJ&lpg=PP1&ots=q3tIHf9A8I&dq=info%3AzjpNBiIw5coJ%3Ascholar.google.com&lr&pg=PP1#v=onepage&q&f=false)

我的部分联系方式：

[hsingyu.yen@icloud.com](mailto:hsingyu.yen@icloud.com)

</div>
</div>

## 参考书目

Numerical Solution of Partial Differential Equations by the Finite Element Method, Claes Johnson: