2020-5-27有限元从0开始第八天(插值算子)

继续给出参考资料镇楼
Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007.
2021-5-26 第七天 (有限元的基本概念)

昨日内容简要回忆

在前面昨天的文章中, 讲述了有限元的基本概念, 给出了有限元的基本定义, 尤其是给出引理, 给出验证有限元定义第三条的方法, 即怎么去验证所给出的 $\mathcal{N}$, 定义了给定有限元空间 $\mathcal{P}$ 的对偶空间 ${\mathcal{P}}^{{}^{\prime }}$上的一个基函数.

3.3 插值

参考资料的 P71-77 给出了常用的三角形有限元的例子,可自行参考.
本章要进行的事情, 是希望把前面给出的三角形单元通过适当的方法进行拼接, 利用这样的方法来去构造出 Sobolev 空间的子空间, 首先给出局部插值算子的定义
定义: 局部插值算子
考虑有限元 $\left(K,\mathcal{P},\mathcal{N}\right)$, 记集合 { ϕ i ∣ 1 ≤ i ≤ k } ⊂ P \left\lbrace \phi_{i} | 1 \leq i \leq k \right\rbrace \subset \mathcal{P} {ϕi​∣1≤i≤k}⊂P 表示 $\mathcal{N}$ 的对偶基, 假设 $v$ 为函数, 并且对于 $N_i\in \mathcal{N}$, $N_i(v)$ 能够定义, 则可以定义如下的插值算子:
$${\mathcal{I}}_{K}v:=\sum _{i=1}^k{N}_i(v){\varphi }_i$$

例子:
记 $K$ 表示如图所示的三角形, $\mathcal{P}={\mathcal{P}}_1$ 为次数不超过1的多项式空间, N = { N 1 , N 2 , N 3 } \mathcal{N} = \left\lbrace N_{1}, N_{2}, N_{3} \right\rbrace N={N1​,N2​,N3​}, $N_i(v)=v(i)$, 下面考虑 $f(x,y)=e^{xy}$, 写出 ${\mathcal{I}}_{k}f$
$${\mathcal{I}}_{k}f=N_1(f){\varphi }_1+N_2(f){\varphi }_2+N_3(f){\varphi }_3$$
给定点编号顺序: 从左往右,从下往上依次编号为 1,2,3, 则:
$$\begin{array}{rl}& N_1(f)=N_2(f)=N_3(f)=1\\ & {\varphi }_1(x,y)=1-x-y\\ & {\varphi }_2(x,y)=x\\ & {\varphi }_3(x,y)=y\end{array}$$
从而有:
$$\begin{array}{r}{\mathcal{I}}_{k}f=1-x-y+x+y=1\end{array}$$

插值算子性质:

性质1 ${\mathcal{I}}_K$ 为线性算子.

证明: 由定义: ${\mathcal{I}}_{K}v=\sum _{i=1}^k{\varphi }_i{N}_i(v)$, 利用 $N_i$ 的线性性质, 马上可以知道 ${\mathcal{I}}_k$ 为线性算子.

性质2 $N_i\left({\mathcal{I}}_{K}f\right)=N_{i}f,(\text{插值性质}).$

证明:
$$\begin{array}{rl}{N}_i({\mathcal{I}}_{K}f)& =N_i\left(\sum _{j=1}^k{N}_j(f){\varphi }_j\right)\\ & =\sum _{j=1}^k{N}_j(f)N_i({\varphi }_j)=\sum _{j=1}^k{N}_j(f){\delta }_{ij}\\ & =N_i(f)\end{array}$$
性质3 假设, $f\in \mathcal{P}$, 则有${\mathcal{I}}_{K}f=f$, 从而可以知道, ${\mathcal{I}}_k$ 为幂等算子.
证明, 利用 $N_i\left({\mathcal{I}}_{K}f-f\right)=0,\forall i$ 以及 ${\mathcal{I}}_{K}f-f\in \mathcal{P}$, 由前面引理可知:
$${\mathcal{I}}_{k}f=f$$
利用上式, 可知幂等性显然.

下满考虑单元的拼接, 首先给出以下定义:

定义:(区域的划分) 称有限个单元区域 { K i } \left\lbrace K_{i} \right\rbrace {Ki​}是区域 $\Omega$ 的一个划分, 如以下两条成立

• $\text{int}{K}_i\cap \text{int}{K}_j=\varnothing$,
• $\cup K_i=\overline{\Omega }$.

定义:(全局插值算子) 假设 $\Omega$ 为区域, $\mathcal{T}$ 为区域 $\Omega$ 的一个剖分, 假设在每个单元 $K$ 使用形函数空间 $\mathcal{P}$, 以及定义节点变量 $\mathcal{N}$, 使得 $(K,\mathcal{P},\mathcal{N})$ 为构成有限元, 假设 $m$ 是节点变量所具有的最高阶导数, 假设 $f\in C^m(\overline{\Omega })$, 则可以如下定义全局插值:

$${{\mathcal{I}}_{\mathcal{T}}f|}_{K_i}={\mathcal{I}}_{K_i}f$$

注:如果不对划分作进一步的假设, 则全局插值不一定具有连续性.后面就是去介绍在什么条件下能够获得这种全局的连续性. 下面只考虑二维的三角形单元, 类似的结果可以做相应的推广.

定义:(多边形区域的三角剖分) 多边形区域的三角划分,处理满足以上条件,还需要满足:

• 任意一个三角形, 定点不会落在其它三角形边的内部或者其它三角形的内部.

例子:(全局插值算子) 考虑如下矩形区域的三角剖分
考虑拉格朗日有限元, 可以知道, $T_1$ 的对偶基为: { 1 − x − y , x , y } \left\lbrace 1 - x - y, x, y \right\rbrace {1−x−y,x,y}, 以及 $T_2$ 的对偶基为: { 1 − x , 1 − y , x + y − 1 } \left\lbrace 1-x, 1-y, x+y-1\right\rbrace {1−x,1−y,x+y−1}, 考虑 $f=\sin \left(\pi (x+y)/2\right)$, 计算 $f$ 到区域上的全局插值.

由:

$$\begin{array}{r}f(0,0)=0,f(1,0)=1\\ f(0,1)=1,f(1,1)=0\end{array}$$
从而有:
$$\begin{array}{rl}& {{\mathcal{I}}_{\mathcal{T}}f|}_{T_1}=0\cdot (1-x-y)+1\cdot x+1\cdot y=x+y\\ & {{\mathcal{I}}_{\mathcal{T}}f|}_{T_2}=1\cdot (1-x)+1\cdot (1-y)+0\cdot (x+y-1)=2-x-y\end{array}$$
从而可以知道全局插值算子:
$${\mathcal{I}}_{\mathcal{T}}f=\{\begin{array}{l}x+y\text{on}{T}_1\\ 2-x-y\text{on}{T}_2.\end{array}$$

注: 为了对齐次 Dirichlet 边值问题进行逼近,将使用由符合边界条件的分片多项式来作为有限维空间, 即:
$$V_T=\left\{{\mathcal{I}}_{T}f:f\in C^m\left(\overline{\Omega }\right),{f|}_{\partial \Omega }=0\right\}$$

定义(’’$C^r$’‘有限元空间) 称插值具有 $r$ 阶连续性,简记为 ‘’$C^r$’’ 如果对于一切的 $f\in C^m\left(\overline{\Omega }\right)$, 成立 ${\mathcal{I}}_{T}f\in C^r$, 称空间 V T = { I T f : f ∈ C m } V_{T} = \left\lbrace \mathcal{I}_{T}f : f \in C^{m} \right\rbrace VT​={IT​f:f∈Cm} 为 “$C^r$” 有限元空间.

注释 用于构造 $C^r$ 空间的有限元通常称为是 $C^r$ 连续的, 并不是任何一种选取节点的方式, 构造出来的有限元都能是 $C^r$ 连续的, 节点的选择必须符合一定的规律. 现有的有限元当中, 选取节点的一个关键点是, 节点应该按照不依赖于坐标的, 关于边中点具有对称性的顺序进行配置.

性质: Lagrange 有限元以及 Hernite 有限元都是 $C^0$ 有限元, $Argyris$ 有限元为 $C^1$ 有限, 更准确地来说, 对于一个给定区域 $\Omega$ 的三角形网格 剖分 $\mathcal{T}$, 能够找到合适的节点, 构造出有限元 $(K,\mathcal{P},\mathcal{N})$, $K\in \mathcal{T}$, 使得当 $f\in C^m$ (Lagrange 元 $m=0$, Hermite元, $m=1$, Argyris 元 $m=2$ ) 时, 成立 ${\mathcal{I}}_{\mathcal{T}}f\in C^r$ (Lagrange 以及 Hermite 元, $r=0$, Argyris 元, $r=1$. ) 特别的, 如果取每条边$\overline{\mathbf{x}{\mathbf{x}}^{\prime }}$的节点为: ${\xi }_i({\mathbf{x}}^{\prime }-\mathbf{x})+\mathbf{x}$, 其中 { ξ i = 1 , ⋯   , k − 1 − 2 m } \left\lbrace \xi_{i} = 1, \cdots, k - 1 -2m\right\rbrace {ξi​=1,⋯,k−1−2m} 固定, 并且关于 $\xi =1/2$ 对称, 除此以外, 在这些假设之下, 成立 ${\mathcal{I}}_{\tau }f\in W_{\infty }^{r+1}.$

只需要证明结论在跨越边界的地方成立即可, 考虑两个具有公共边 $e$ 的单元, 由于边上的节点以某种对称的方式不依赖于坐标地进行选取, 因此两个单元在边 $e$ 上公用节点, 记 $w={\mathcal{I}}_{T_1}f-{\mathcal{I}}_{T_2}f$, 则 $w$ 仍然是一个相应次数的多项式, 并且将其限制到边上, 具有一维的 Lagrange, Hermite, Argyris 节点, 并且都为 0, 从而知道 ${w|}_e=0$. 从而跨越每条边时, 插值都是连续的.