DG

3.1 Legendre polynomials and nodal elements
$\displaystyle x(r)=x^k_l+\frac{1+r}{2}h^k,h^k=x^k_r-x_l^k,r\in [-1,1]$
给出$u(x,t)$的modal 和nodal表达：
$\displaystyle u^k_h(x(r),t)=\sum^{N_p}_{n=1}\hat{u}^k_n(t)\phi_n(r)=\sum^{N_p}_{n=1}u^k_h(x^k_i,t)l^k_i(r)$
注意，第二式称为节点表达式，$x_i^k$称之为节点，就是第$k$个区间里的第$i$个节点。这的$\displaystyle \phi_n(r)=\frac{P_{n-1}(r)}{\sqrt{\gamma_{n-1}}},\gamma_{n}=\frac{2}{2n+1}$,$P_{n}$是Legendre多项式，这样的话就有$\displaystyle \int_{-1}^{1}\phi_i(r)\phi_j(r)dr=\delta_{ij},M=((\phi_i,\phi_j))_{N_p\times N_p}=I$.
$\displaystyle (u(r),\phi(r)_m)=\sum^{N_p}_{n=1}\hat{u}_n(\phi_n,\phi_m) \Rightarrow M\hat{u}=u\Rightarrow \hat{u}=u\Rightarrow \hat{u}_i=(u,\phi_i)=\int_{-1}^1u\phi_idr$,
这里的$\hat{u}=(\hat{u}_1,...\hat{u}_{N_p})$是待求的系数，$u=((u,\phi_1),(u,\phi_2)...(u,\phi_{N_p}))$.
对于这个积分，用高斯数值积分公式$\displaystyle \int_{-1}^1u\phi_idr \simeq \sum_{k=1}^{N_p}u(r_k)\phi_i(r_k)\omega_k$
注意在modal表达式中取节点$\{ \xi_i\}$ :

u(ξi)=∑n=1Npu^nϕn(ξi)

$$u(\xi_i)=\sum_{n=1}^{N_p}\hat{u}_n\phi_n(\xi_i)$$
记$V_{in}=\phi_n(\xi_i),u=(u(\xi_1),u(\xi_2)...u(\xi_{N_p})),\hat{u}=(\hat{u}_i)$,上式可以写成矩阵表达式
$$V\hat{u}=u$$
试着推导关系 $$V^Tl(r)=\bar{P}(r)$$
番外篇：取数值通量$(au)^*=au^-$,$u_i$代表$u$在区间$[x_i,x_{i+1}]$之间的逼近
$$\int_{x_i}^{x_{i+1}}(\frac{\partial u_i}{\partial t}+a \frac{\partial u_i}{\partial x})\phi_n dx=\phi_n(x_i)(u_{i-1}(x_i)-u_i(x_i))$$
练习1：推导格式,令$\displaystyle \rho_i=\sum_{m=0}^p \rho_{i,m}L_m(\tilde{x}),where\ \tilde{x}=\frac{2x-x_i-x_{i+1}}{x_{i+1}-x_i}$
$$\begin{aligned} &\ \sum_{m=0}^p M_{nm}\frac{d \rho_{i,m}}{dt}+a\sum_{m=0}^p D_{nm}\rho_{i,m}\\ &=aL_n(-1)\sum_{m=0}^p \rho_{i-1,m}L_m(1)-aL_n(-1)\sum_{m=0}^p \rho_{i,m}L_m(-1)\\ &=a(-1)^n\sum_{m=0}^p \rho_{i-1,m}-(-1)^{n+m}a\sum_{m=0}^p \rho_{i,m} \end{aligned}$$
写成矩阵的形式就是：
$$M\frac{d\rho_i}{dt}+a D\rho_i=a[(-1)^n]\rho_{i-1}-[(-1)^{n+m}]\rho_i$$ 其中需要注意的是：$[(-1)^n]$是矩阵，它的每一行都是$(-1)^n,n$是行指标。$\rho_i$也是一个列向量，代表在$i-$cell里面的系数。
以下将讨论如何将其化成Runge-Kutta格式来计算：
$$\frac{d\rho_i}{dt}=(\frac{2n+1}{x_{i+1}-x_i})(-a(D+[(-1)^{n+m}]\rho_i)+a\rho_{i-1}[(-1)^n])$$
三阶Runge—Kutta格式：取时间步长为$\delta$,编程时就直接取$\delta =dt$,开始下面的循环：
$$\begin{aligned} \sigma&=\rho_i(t)\\ \sigma&=\rho_i(t)+\frac{\delta}{3}\frac{d}{dt}\sigma\\ \sigma&=\rho_i(t)+\frac{\delta}{2}\frac{d}{dt}\sigma\\ \sigma&=\rho_i(t)+\frac{\delta}{1}\frac{d}{dt}\sigma\\ \rho_i(t+\delta)&=\sigma \end{aligned}$$
练习：从Taylor级数$\displaystyle \rho_i(t+dt)=\rho_i(t)+\rho_i'dt+\rho_i''\frac{dt^2}{2}+\frac{dt^2}{3!}\rho_i'''+o{dt^4}$推导出上面的循环。
回答以下问题：
（1）简要概述RK-DG的算法；
（2）怎样计算$M$?写出程序来；
（3）怎样计算$D$?写出程序来；
（4）怎样实现RK循环？写出程序。
（5）为了回答以上四个问题，你还需要知道怎样生成Legendre基。
Legendre polynomial $$L_{m+1}(x)=\frac{2m+1}{m+1}xL_{m}(x)-\frac{m}{m+1}L_{m-1}(x);L_0=1,L_1=x$$
Orthogonality of the Legendre polynomial $$\int_{-1}^1 L_n(x)L_m(x)dx =0$$
计算矩阵$D$: $\displaystyle \frac{d \rho_i}{d\tilde{x}}=\sum_{m=0}^p \rho_{i,m}^{(1)} L_{m}(\tilde{x}), \rho_{i,m}^{(1)}=(2m+1)\sum_{j=m+1}^p\rho_{i,j}$

$L_2$范数误差估计

error=∫10(uh−u)2dx−−−−−−−−−−−−√

$$error=\sqrt{\int_0^1 (u_h-u)^2dx}$$
$$\begin{aligned} e_h&=Ch^r\\ e_{2h}&=C(2h)^2\\ \frac{e_{2h}}{e_h}&=2^r\\ \frac{ln\frac{e_{2h}}{e_h}}{ln2}&=r \end{aligned}$$

fun=@(x)(sin(x)-x).^2;%必须注意：数值积分的时候^ 使用.^
error2=sqrt(integral(fun,0,1))%被积分函数必须写成@(x)函数

error2 =

    0.5774