Hopf-Cole 變換與高階差分格式

非線性發(fā)展方程的有限差分方法

$1.4$ Hopf-Cole 變換與高階差分格式

$1.4.1$ Hopf-Cole 變換

令
$\begin{aligned} &u(x, t)=-2 \nu \frac{w_{x}(x, t)}{w(x, t)} \end{aligned}\tag{1.66}$

則有

$\begin{aligned} &u_{t}=-2 \nu\left(\frac{w_{x}}{w}\right)_{t}=-2 \nu \frac{w_{x t} w-w_{x} w_{t}}{w^{2}}=-2 \nu\left(\frac{w_{t}}{w}\right)_{x} \\ &u_{x}=-2 \nu\left(\frac{w_{x}}{w}\right)_{x} \\ &u_{x x}=-2 \nu\left(\frac{w_{x}}{w}\right)_{x x} \end{aligned}$

將上式代入 (1.1), 得到

$-2 \nu\left(\frac{w_{t}}{w}\right)_{x}+\left(-2 \nu \frac{w_{x}}{w}\right)\left(-2 \nu\left(\frac{w_{x}}{w}\right)_{x}\right)=\nu\left[-2 \nu\left(\frac{w_{x}}{w}\right)_{x x}\right]$

即

$\left(\frac{w_{t}}{w}\right)_{x}-\nu\left[\left(\frac{w_{x}}{w}\right)^{2}\right]_{x}=\nu\left(\frac{w_{x}}{w}\right)_{x x} \text {, }$

或

$\left[\frac{w_{t}}{w}-\nu\left(\frac{w_{x}}{w}\right)^{2}-\nu\left(\frac{w_{x}}{w}\right)_{x}\right]_{x}=0$

上式又可以寫成

$\left(\frac{w_{t}-\nu w_{x x}}{w}\right)_{x}=0$

因而

$\frac{w_{t}-\nu w_{x x}}{w}=q(t)$

可以將上式寫成

$w_{t}-q(t) w=\nu w_{x x}$

上式兩邊同乘以 $e^{-\int_{0}^{t} q(s) \mathrmbdrx1xh s}$ , 則可得到

$\left[w e^{-\int_{0}^{t} q(s) \mathrmtrjjn3r s}\right]_{t}=\nu\left[w e^{-\int_{0}^{t} q(s) \mathrmnl1d33f s}\right]_{x x}$

$\tilde{w}(x, t)=w(x, t) e^{-\int_{0}^{t} q(s) \mathrmhzxdvdz s}$

則

$-2 \nu \frac{\widetilde{w}_{x}}{\widetilde{w}}=-2 \nu \frac{w_{x}}{w}=u(x, t)$

即對于任意 $q(t)$ , 不影響 $u(x, t)$ . 因而取 $q(t)=0$ . 于是得到如下等價問題

$\begin{aligned} &w_{t}-\nu w_{x x}=0, \quad 0<x<L, 0<t \leqslant T \\ &w(x, 0)=\tilde{\varphi}(x), \quad 0 \leqslant x \leqslant L, \\ &w_{x}(0, t)=0, \quad w_{x}(L, t)=0, \quad 0 \leqslant t \leqslant T \end{aligned}\tag{1.67-1.69}$

其中

$\widetilde{\varphi}(x)=e^{-\frac{1}{2 \nu} \int_{0}^{x} \varphi(s) \mathrm9lhbph3 s} .\tag{1.70}$

稱 (1.66) 為 Hopf-Cole 變換.

$1.4.2$ 差分格式的建立

下面給出幾個帶積分余項的數值微分公式.

引理 $1.4([20]) \quad$ 記 $\alpha(s)=(1-s)^{2}\left[5-(1-s)^{2}\right]$ .

(a) 設 $g(x) \in C^{6}\left[x_{0}, x_{1}\right]$ , 有

$\begin{aligned} & \frac{5}{6} g^{\prime \prime}\left(x_{0}\right)+\frac{1}{6} g^{\prime \prime}\left(x_{1}\right)-\frac{2}{h}\left[\frac{g\left(x_{1}\right)-g\left(x_{0}\right)}{h}-g^{\prime}\left(x_{0}\right)\right] \\ =&-\frac{h}{6} g^{\prime \prime \prime}\left(x_{0}\right)+\frac{h^{3}}{90} g^{(5)}\left(x_{0}\right)+\frac{h^{4}}{180} \int_{0}^{1} g^{(6)}\left(x_{0}+s h\right) \alpha(s) d s \\ =&-\frac{h}{6} g^{\prime \prime \prime}\left(x_{0}\right)+\frac{h^{3}}{90} g^{(5)}\left(x_{0}\right)+\frac{h^{4}}{240} g^{(6)}\left(x_{0}+\theta_{0} h\right), \quad \theta_{0} \in(0,1) \end{aligned}$

(b) 設 $g(x) \in C^{6}\left[x_{m-1}, x_{m}\right]$ , 有

$\begin{aligned} & \frac{1}{6} g^{\prime \prime}\left(x_{m-1}\right)+\frac{5}{6} g^{\prime \prime}\left(x_{m}\right)-\frac{2}{h}\left[g^{\prime}\left(x_{m}\right)-\frac{g\left(x_{m}\right)-g\left(x_{m-1}\right)}{h}\right] \\ =& \frac{h}{6} g^{\prime \prime \prime}\left(x_{m}\right)-\frac{h^{3}}{90} g^{(5)}\left(x_{m}\right)+\frac{h^{4}}{180} \int_{0}^{1} g^{(6)}\left(x_{m}-s h\right) \alpha(s) d s \\ =& \frac{h}{6} g^{\prime \prime \prime}\left(x_{m}\right)-\frac{h^{3}}{90} g^{(5)}\left(x_{m}\right)+\frac{h^{4}}{240} g^{(6)}\left(x_{m}-\theta_{m} h\right), \quad \theta_{m} \in(0,1) \end{aligned}$

$\begin{aligned} & \frac{1}{12}\left[g^{\prime \prime}\left(x_{i-1}\right)+10 g^{\prime \prime}\left(x_{i}\right)+g^{\prime \prime}\left(x_{i+1}\right)\right]-\frac{1}{h^{2}}\left[g\left(x_{i+1}\right)-2 g\left(x_{i}\right)+g\left(x_{i-1}\right)\right] \\ =& \frac{h^{4}}{360} \int_{0}^{1}\left[g^{(6)}\left(x_{i}+s h\right)+g^{(6)}\left(x_{i}-s h\right)\right] \alpha(s) d s \\ =& \frac{h^{4}}{240} g^{(6)}\left(x_{i}+\theta_{i} h\right), \quad \theta_{i} \in(-1,1) \end{aligned}$

設 $v \in \mathcal{U}_{h} .$ 定義平均值算子 $\mathcal{A}$ :

$\mathcal{A} v_{i}= \begin{cases}\frac{5}{6} v_{0}+\frac{1}{6} v_{1}, & i=0, \\ \frac{1}{12}\left(v_{i-1}+10 v_{i}+v_{i+1}\right), & 1 \leqslant i \leqslant m-1 \\ \frac{5}{6} v_{m}+\frac{1}{6} v_{m-1}, & i=m .\end{cases}$

設 $(1.67)-(1.70)$ 存在解 $w(x, t) \in C^{6,4}([0, L] \times[0, T])$ . 定義網格函數

$U_{i}^{k}=u\left(x_{i}, t_{k}\right), \quad w_{i}^{k}=W\left(x_{i}, t_{k}\right), \quad 0 \leqslant i \leqslant m, \quad 0 \leqslant k \leqslant n$

由 (1.67) 和 (1.69) 可得

$\begin{aligned} &w_{x x x}(0, t)=0, \quad w_{x x x}(L, t)=0, \quad 0 \leqslant t \leqslant T, \\ &w_{x x x x x}(0, t)=0, \quad w_{x x x x x}(L, t)=0, \quad 0 \leqslant t \leqslant T \end{aligned}\tag{1.71-1.72}$

在點 $\left(x_{i}, t_{k+\frac{1}{2}}\right)$ 處考慮方程 (1.67), 有

$w_{t}\left(x_{i}, t_{k+\frac{1}{2}}\right)-\nu w_{x x}\left(x_{i}, t_{k+\frac{1}{2}}\right)=0, \quad 0 \leqslant i \leqslant m, 0 \leqslant k \leqslant n-1$

應用引理 $1.2$ , 可得

$w_{t}\left(x_{i}, t_{k+\frac{1}{2}}\right)-\frac{\nu}{2}\left[w_{x x}\left(x_{i}, t_{k}\right)+w_{x x}\left(x_{i}, t_{k+1}\right)\right]=\left(R_{x t} w\right)_{i}^{k+\frac{1}{2}}$

$0 \leqslant i \leqslant m, 0 \leqslant k \leqslant n-1$

其中

$\left(R_{x t} w\right)_{i}^{k+\frac{1}{2}}=-\frac{\nu \tau^{2}}{8} \frac{\partial^{4} w}{\partial x^{2} \partial t^{2}}\left(x_{i}, t_{k}+\eta_{i}^{k} \tau\right), \quad \eta_{i}^{k} \in(0,1), 0 \leqslant i \leqslant m$

用算子 $\mathcal{A}$ 作用上述等式的兩邊, 得

$\begin{gathered} \mathcal{A} w_{t}\left(x_{i}, t_{k+\frac{1}{2}}\right)-\frac{\nu}{2}\left[\mathcal{A} w_{x x}\left(x_{i}, t_{k}\right)+\mathcal{A} w_{x x}\left(x_{i}, t_{k+1}\right)\right]=\mathcal{A}\left(R_{x t} w\right)_{i}^{k+\frac{1}{2}} \\ 1 \leqslant i \leqslant m-1,0 \leqslant k \leqslant n-1 \end{gathered}$

應用引理 $1.4$ 并注意到 (1.71)-(1.72), 得到

$\begin{aligned} &\mathcal{A} \delta_{t} W_{0}^{k+\frac{1}{2}}-\nu\left(\frac{2}{h} \delta_{x} W_{\frac{1}{2}}^{k+\frac{1}{2}}\right)=R_{0}^{k+\frac{1}{2}}, \quad 0 \leqslant k \leqslant n-1 \\ &\mathcal{A} \delta_{t} W_{i}^{k+\frac{1}{2}}-\nu \delta_{x}^{2} W_{i}^{k+\frac{1}{2}}=R_{i}^{k+\frac{1}{2}}, \quad 1 \leqslant i \leqslant m-1,0 \leqslant k \leqslant n-1 \\ &\mathcal{A} \delta_{t} W_{m}^{k+\frac{1}{2}}-\nu\left(-\frac{2}{h} \delta_{x} W_{m-\frac{1}{2}}^{k+\frac{1}{2}}\right)=R_{m}^{k+\frac{1}{2}}, \quad 0 \leqslant k \leqslant n-1 \end{aligned}\tag{1.73-1.75}$

存在常數 $c_{8}$ 使得

$\left|R_{i}^{k+\frac{1}{2}}\right| \leqslant c_{8}\left(\tau^{2}+h^{4}\right), \quad 0 \leqslant i \leqslant m, 0 \leqslant k \leqslant n-1 \tag{1.76}$

注意到初值條件

$W_{i}^{0}=\tilde{\varphi}\left(x_{i}\right), \quad 0 \leqslant i \leqslant m. \tag{1.77}$

在 (1.73)-(1.75) 中略去小量項, 對 (1.67)-(1.69) 建立如下差分格式

$\begin{aligned} &\mathcal{A} \delta_{t} w_{0}^{k+\frac{1}{2}}-\nu\left(\frac{2}{h} \delta_{x} w_{\frac{1}{2}}^{k+\frac{1}{2}}\right)=0, \quad 0 \leqslant k \leqslant n-1 \\ &\mathcal{A} \delta_{t} w_{i}^{k+\frac{1}{2}}-\nu \delta_{x}^{2} w_{i}^{k+\frac{1}{2}}=0, \quad 1 \leqslant i \leqslant m-1,0 \leqslant k \leqslant n-1, \\ &\mathcal{A} \delta_{t} w_{m}^{k+\frac{1}{2}}-\nu\left(-\frac{2}{h} \delta_{x} w_{m-\frac{1}{2}}^{k+\frac{1}{2}}\right)=0, \quad 0 \leqslant k \leqslant n-1 \\ &w_{i}^{0}=\widetilde{\varphi}\left(x_{i}\right), \quad 0 \leqslant i \leqslant m \end{aligned}\tag{1.78-1.81}$

$1.4.3$ 差分格式解的存在性和唯一性}

引理 $1.5$ 設 $v=\left(v_{0}, v_{1}, \cdots, v_{m}\right) \in \mathcal{U}_{h}$ , 則有

$\begin{aligned} &(\mathcal{A} v, v)=\|v\|^{2}-\frac{h^{2}}{12}|v|_{1}^{2} \\ &\frac{2}{3}\|v\|^{2} \leqslant(\mathcal{A} v, v) \leqslant\|v\|^{2} \end{aligned}$

證明由

$(\mathcal{A} v)_{i}= \begin{cases}v_{0}+\frac{h}{6} \delta_{x} v_{\frac{1}{2}}, & i=0 \\ v_{i}+\frac{h^{2}}{12} \delta_{x}^{2} v_{i}, & 1 \leqslant i \leqslant m-1 \\ v_{m}-\frac{h}{6} \delta_{x} v_{m-\frac{1}{2}}, & i=m\end{cases}$

得到

$\begin{aligned} (\mathcal{A} v, v) &=h\left[\frac{1}{2}\left(\mathcal{A} v_{0}\right) v_{0}+\sum_{i=1}^{m-1}\left(\mathcal{A} v_{i}\right) v_{i}+\frac{1}{2}\left(\mathcal{A} v_{m}\right) v_{m}\right] \\ &=h\left[\frac{1}{2}\left(v_{0}+\frac{h}{6} \delta_{x} v_{\frac{1}{2}}\right) v_{0}+\sum_{i=1}^{m-1}\left(v_{i}+\frac{h^{2}}{12} \delta_{x}^{2} v_{i}\right) v_{i}+\frac{1}{2}\left(v_{m}-\frac{h}{6} \delta_{x} v_{m-\frac{1}{2}}\right) v_{m}\right] \\ &=h\left(\frac{1}{2} v_{0}^{2}+\sum_{i=1}^{m-1} v_{i}^{2}+\frac{1}{2} v_{m}^{2}\right) \end{aligned}$

$\begin{aligned} &+\frac{h^{2}}{12}\left[\left(\delta_{x} v_{\frac{1}{2}}\right) v_{0}+h \sum_{i=1}^{m-1}\left(\delta_{x}^{2} v_{i}\right) v_{i}-\left(\delta_{x} v_{m-\frac{1}{2}}\right) v_{m}\right] \\ =&\|v\|^{2}-\frac{h^{2}}{12} \cdot h \sum_{i=0}^{m-1}\left(\delta_{x} v_{i+\frac{1}{2}}\right)^{2} \\ =&\|v\|^{2}-\frac{h^{2}}{12}|v|_{1}^{2} . \end{aligned}$

易知

$(\mathcal{A} v, v) \leqslant\|v\|^{2}$

由

$|v|_{1}^{2} \leqslant \frac{4}{h^{2}}\|v\|^{2}$

易得

$(\mathcal{A} v, v) \geqslant\|v\|^{2}-\frac{h^{2}}{12} \cdot \frac{4}{h^{2}}\|v\|=\frac{2}{3}\|v\|^{2}$

定義 $\mathcal{U}_{h}$ 上的范數

$\|v\|_{\mathcal{A}}=\sqrt{(\mathcal{A} v, v)}$

由引理 $1.5$ 知 $\Vert v\Vert _{\mathcal{A}}$ 和 $\Vert v\Vert$ 等價.

定理 $1.11$ 差分格式 (1.78)-(1.81) 的解是存在唯一的.

證明第 0 層的值 $w^{0}$ 是由（1.81）給定. 設已得到第 $k$ 層的值 $w^{k}$ , 則由 (1.78)-(1.80) 可得關于第 $k+1$ 層值 $w^{k+1}$ 的線性方程組. 考慮其齊次方程組

$\begin{aligned} &\frac{1}{\tau} \mathcal{A} w_{0}^{k+1}-\nu \frac{1}{h} \delta_{x} w_{\frac{1}{2}}^{k+1}=0 \\ &\frac{1}{\tau} \mathcal{A} w_{i}^{k+1}-\nu \frac{1}{2} \delta_{x}^{2} w_{i}^{k+1}=0, \quad 1 \leqslant i \leqslant m-1 \\ &\frac{1}{\tau} \mathcal{A} w_{m}^{k+\frac{1}{2}}-\nu\left(-\frac{1}{h} \delta_{x} w_{m-\frac{1}{2}}^{k+1}\right)=0 \end{aligned}\tag{1.82-1.84}$

用 $\frac{1}{2} h w_{0}^{k+1}$ 與 $(1.82)$ 相乘, 用 $h w_{i}^{k+1}$ 與 $(1.83)$ 相乘, 用 $\frac{1}{2} h w_{m}^{k+\frac{1}{2}}$ 與 $(1.84)$ 相乘, 并將所得結果相加, 得

$\frac{1}{\tau}\left(\mathcal{A} w^{k+1}, w^{k+1}\right)-\frac{1}{2} \nu\left[w_{0}^{k+\frac{1}{2}} \delta_{x} w_{\frac{1}{2}}^{k+1}+h \sum_{i=1}^{m-1} w_{i}^{k+1} \delta_{x}^{2} w_{i}^{k+1}-w_{m}^{k+1} \delta_{x} w_{m-\frac{1}{2}}^{k+1}\right]=0$

由上式易得

$\frac{1}{\tau}\left\|w^{k+1}\right\|_{\mathcal{A}}^{2}+\frac{1}{2} \nu\left|w^{k+1}\right|_{1}^{2}=0$

因而

$w^{k+1}=0$

即 $(1.82)-(1.84)$ 只有零解. 于是 $(1.78)-(1.80)$ 唯一確定 $w^{k+1}$ .

$1.4.4$ 差分格式解的收斂性

定理 $1.12$ 設 $\left\{W_{i}^{k} \mid 0 \leqslant i \leqslant m, 0 \leqslant k \leqslant n\right\}$ 是問題 (1.67)-(1.69) 的解, $\left\{w_{i}^{k} \mid 0 \leqslant i \leqslant m, 0 \leqslant k \leqslant n\right\}$ 是差分格式 $(1.78)-(1.81)$ 的解. 令

$e_{i}^{k}=W_{i}^{k}-w_{i}^{k}, \quad 0 \leqslant i \leqslant m, 0 \leqslant k \leqslant n$

則存在常數 $c_{9}, c_{10}$ 使得

$\begin{aligned} &\left\|e^{k}\right\| \leqslant c_{9}\left(\tau^{2}+h^{4}\right), \quad 0 \leqslant k \leqslant n \\ &\left|e^{k}\right|_{1} \leqslant c_{10}\left(\tau^{2}+h^{4}\right), \quad 0 \leqslant k \leqslant n \\ &\left\|e^{k}\right\|_{\infty} \leqslant \frac{\sqrt{L}}{2} c_{10}\left(\tau^{2}+h^{4}\right), \quad 0 \leqslant k \leqslant n \end{aligned}\tag{1.85-1.87}$

證明將 (1.73)-(1.75), (1.77) 和 (1.78)-(1.81) 相減, 得到誤差方程組

$\begin{aligned} &\mathcal{A} \delta_{t} e_{0}^{k+\frac{1}{2}}-\nu \cdot \frac{2}{h} \delta_{x} e_{\frac{1}{2}}^{k+\frac{1}{2}}=R_{0}^{k+\frac{1}{2}}, \quad 0 \leqslant k \leqslant n-1 \\ &\mathcal{A} \delta_{t} e_{i}^{k+\frac{1}{2}}-\nu \delta_{x}^{2} e_{i}^{k+\frac{1}{2}}=R_{i}^{k+\frac{1}{2}}, \quad 1 \leqslant i \leqslant m-1,0 \leqslant k \leqslant n-1, \\ &\mathcal{A} \delta_{t} e_{m}^{k+\frac{1}{2}}-\nu\left(-\frac{2}{h} \delta_{x} e_{m-\frac{1}{2}}^{k+\frac{1}{2}}\right)=R_{m}^{k+\frac{1}{2}}, \quad 0 \leqslant k \leqslant n-1 \\ &e_{i}^{0}=0, \quad 0 \leqslant i \leqslant m . \\ &\text { (I) 用 } \frac{1}{2} h e_{0}^{k+\frac{1}{2}} \text { 與 }(1.88) \text { 相乘, 用 } h e_{i}^{k+\frac{1}{2}} \text { 與 }(1.89) \text { 相乘, 用 } \frac{1}{2} h e_{m}^{k+\frac{1}{2}} \text { 與 } \end{aligned}\tag{1.88-1.91}$

相乘, 并將結果相加, 得到

$\begin{aligned} & \frac{1}{2 \tau}\left(\left\|e^{k+1}\right\|_{\mathcal{A}}^{2}-\left\|e^{k}\right\|_{\mathcal{A}}^{2}\right)-\nu\left(e_{0}^{k+\frac{1}{2}} \delta_{x} e_{\frac{1}{2}}^{k+\frac{1}{2}}+h \sum_{i=1}^{m-1} e_{i}^{k+\frac{1}{2}} \delta_{x}^{2} e_{i}^{k+\frac{1}{2}}-e_{m}^{k+\frac{1}{2}} \delta_{x} e_{m-\frac{1}{2}}^{k+\frac{1}{2}}\right) \\ =&\left(R^{k+\frac{1}{2}}, e^{k+\frac{1}{2}}\right), \quad 0 \leqslant k \leqslant n-1 \end{aligned}$

即

$\frac{1}{2 \tau}\left(\left\|e^{k+1}\right\|_{\mathcal{A}}^{2}-\left\|e^{k}\right\|_{\mathcal{A}}^{2}\right)+\nu\left|e^{k+\frac{1}{2}}\right|_{1}^{2}=\left(R^{k+\frac{1}{2}}, e^{k+\frac{1}{2}}\right), \quad 0 \leqslant k \leqslant n-1$

對上式右端用 Cauchy-Schwarz 不等式, 并應用引理 $1.5$ , 得

$\begin{aligned} & \frac{1}{2 \tau}\left(\left\|e^{k+1}\right\|_{\mathcal{A}}^{2}-\left\|e^{k}\right\|_{\mathcal{A}}^{2}\right) \\ \leqslant &\left\|R^{k+\frac{1}{2}}\right\| \cdot\left\|e^{k+\frac{1}{2}}\right\| \\ \leqslant & \sqrt{\frac{3}{2}}\left\|R^{k+\frac{1}{2}}\right\| \cdot\left\|e^{k+\frac{1}{2}}\right\|_{\mathcal{A}} \\ \leqslant & \sqrt{\frac{3}{2}}\left\|R^{k+\frac{1}{2}}\right\| \cdot \frac{\left\|e^{k+1}\right\|_{\mathcal{A}}+\left\|e^{k}\right\|_{\mathcal{A}}}{2} \end{aligned}$

兩邊約去 $\frac{1}{2}\left(\left\|e^{k+1}\right\|_{\mathcal{A}}+\left\|e^{k}\right\|_{\mathcal{A}}\right),$ 得到

$\frac{1}{\tau}\left(\left\|e^{k+1}\right\|_{\mathcal{A}}-\left\|e^{k}\right\|_{\mathcal{A}}\right) \leqslant \sqrt{\frac{3}{2}}\left\|R^{k+\frac{1}{2}}\right\|, \quad 0 \leqslant k \leqslant n-1$

因而

$\left\|e^{k+1}\right\|_{\mathcal{A}} \leqslant\left\|e^{0}\right\|_{\mathcal{A}}+\sqrt{\frac{3}{2}} \tau \sum_{l=0}^{k}\left\|R^{l+\frac{1}{2}}\right\|, \quad 0 \leqslant k \leqslant n-1$

由 $(1.76)$ 和 $(1.91)$ , 得到

$\left\|e^{k+1}\right\|_{A} \leqslant \sqrt{\frac{3}{2}}(k+1) \tau \sqrt{L} c_{8}\left(\tau^{2}+h^{4}\right) \leqslant \sqrt{\frac{3}{2}} T \sqrt{L} c_{8}\left(\tau^{2}+h^{4}\right), \quad 0 \leqslant k \leqslant n-1$

再次應用引理 $1.5$ , 得

$\left\|e^{k}\right\| \leqslant \sqrt{\frac{3}{2}}\left\|e^{k}\right\|_{\mathcal{A}} \leqslant \frac{3}{2} T \sqrt{L} c_{8}\left(\tau^{2}+h^{2}\right), \quad 1 \leqslant k \leqslant n$

(II) 用 $\frac{1}{2} h \delta_{t} e_{0}^{k+\frac{1}{2}}$ 與 $(1.88)$ 相乘, 用 $h \delta_{t} e_{i}^{k+\frac{1}{2}}$ 與 (1.89) 相乘, 用 $\frac{1}{2} h \delta_{t} e_{m}^{k+\frac{1}{2}}$ 與 (1.90) 相乘, 并將結果相加, 得

$\left\|\delta_{t} e^{k+\frac{1}{2}}\right\|_{\mathcal{A}}^{2}-\nu\left[\left(\delta_{x} e_{\frac{1}{2}}^{k+\frac{1}{2}}\right) \delta_{t} e_{0}^{k+\frac{1}{2}}+h \sum_{i=1}^{m-1}\left(\delta_{x}^{2} e_{i}^{k+\frac{1}{2}}\right) \delta_{t} e_{i}^{k+\frac{1}{2}}-\left(\delta_{x} e_{m-\frac{1}{2}}^{k+\frac{1}{2}}\right) \delta_{t} e_{m}^{k+\frac{1}{2}}\right]$

$=\left(R^{k+\frac{1}{2}}, \delta_{t} e^{k+\frac{1}{2}}\right), \quad 0 \leqslant k \leqslant n-1$

即

$\begin{aligned} &\left\|\delta_{t} e^{k+\frac{1}{2}}\right\|_{\mathcal{A}}^{2}+\nu \cdot \frac{1}{2 \tau}\left(\left|e^{k+1}\right|_{1}^{2}-\left|e^{k}\right|_{1}^{2}\right) \\ =&\left(R^{k+\frac{1}{2}}, \delta_{t} e^{k+\frac{1}{2}}\right) \\ \leqslant &\left\|R^{k+\frac{1}{2}}\right\| \cdot\left\|\delta_{t} e^{k+\frac{1}{2}}\right\| \\ \leqslant & \frac{2}{3}\left\|\delta_{t} e^{k+\frac{1}{2}}\right\|^{2}+\frac{3}{8}\left\|R^{k+\frac{1}{2}}\right\|^{2} \\ \leqslant &\left\|\delta_{t} e^{k+\frac{1}{2}}\right\|_{\mathcal{A}}^{2}+\frac{3}{8}\left\|R^{k+\frac{1}{2}}\right\|^{2}, \quad 0 \leqslant k \leqslant n-1 \end{aligned}$

因而

$\frac{1}{2 \tau}\left(\left|e^{k+1}\right|_{1}^{2}-\left|e^{k}\right|_{1}^{2}\right) \leqslant \frac{3}{8 \nu}\left\|R^{k+\frac{1}{2}}\right\|^{2}, \quad 0 \leqslant k \leqslant n-1$

遞推可得

$\left|e^{k+1}\right|_{1}^{2} \leqslant\left|e^{0}\right|_{1}^{2}+\frac{3}{4 \nu} \tau \sum_{l=0}^{k}\left\|R^{l+\frac{1}{2}}\right\|^{2} \leqslant \frac{3}{4 \nu}(k+1) \tau L c_{8}^{2}\left(\tau^{2}+h^{4}\right)^{2}, \quad 0 \leqslant k \leqslant n-1$

于是

$\left|e^{k}\right|_{1} \leqslant \frac{1}{2} \sqrt{\frac{3}{\nu} T L} c_{8}\left(\tau^{2}+h^{4}\right), \quad 1 \leqslant k \leqslant n$

(III) 由 $(1.86)$ 和引理 1.1(b) 得

$\left\|e^{k}\right\|_{\infty} \leqslant \frac{\sqrt{L}}{2}\left|e^{k}\right|_{1} \leqslant \frac{\sqrt{L}}{2} c_{10}\left(\tau^{2}+h^{4}\right), \quad 1 \leqslant k \leqslant n$

$1.4.5$ 原問題解的計算

設 $g(x) \in C^{5}\left[x_{0}, x_{m}\right]$ , 則存在常數 $\hat{c}_{1}, \hat{c}_{2}, \cdots, \hat{c}_{6}$ 使得

$g^{\prime}\left(x_{1}\right)=\hat{c}_{1} \delta_{x} g_{\frac{1}{2}}+\hat{c}_{2} \delta_{x} g_{\frac{3}{2}}+\hat{c}_{3} \delta_{x} g_{\frac{5}{2}}+\hat{c}_{4} \delta_{x} g_{\frac{7}{2}}+O\left(h^{4}\right)\tag{1.92}$

$g^{\prime}\left(x_{i}\right)=\hat{c}_{5} \delta_{x} g_{i-\frac{3}{2}}+\hat{c}_{6} \delta_{x} g_{i-\frac{1}{2}}+\hat{c}_{6} \delta_{x} g_{i+\frac{1}{2}}+\hat{c}_{5} \delta_{x} g_{i+\frac{3}{2}}+O\left(h^{4}\right)$

$2 \leqslant i \leqslant m-2, \tag{1.93}$

$g^{\prime}\left(x_{m-1}\right)=\hat{c}_{1} \delta_{x} g_{m-\frac{1}{2}}+\hat{c}_{2} \delta_{x} g_{m-\frac{3}{2}}+\hat{c}_{3} \delta_{x} g_{m-\frac{5}{2}}+\hat{c}_{4} \delta_{x} g_{m-\frac{7}{2}}+O\left(h^{4}\right)\tag{1.94}$

由變換 (1.66) 有

$u\left(x_{i}, t_{k}\right)=-2 \nu \frac{w_{x}\left(x_{i}, t_{k}\right)}{w\left(x_{i}, t_{k}\right)}$

利用 $(1.92)-(1.94)$ 可得

$\begin{gathered} U_{1}^{k}=-\frac{2 \nu}{W_{1}^{k}}\left(\hat{c}_{1} \delta_{x} W_{\frac{1}{2}}^{k}+\hat{c}_{2} W_{\frac{3}{2}}^{k}+\hat{c}_{3} \delta_{x} W_{\frac{5}{2}}^{k}+\hat{c}_{4} \delta_{x} W_{\frac{7}{2}}^{k}\right)+\hat{R}_{1}^{k}, \quad 1 \leqslant k \leqslant n \\ U_{i}^{k}=-\frac{2 \nu}{W_{i}^{k}}\left(\hat{c}_{5} \delta_{x} W_{i-\frac{3}{2}}^{k}+\hat{c}_{6} \delta_{x} W_{i-\frac{1}{2}}^{k}+\hat{c}_{6} \delta_{x} \hat{W}_{i+\frac{1}{2}}^{k}+\hat{c}_{5} \delta_{x} W_{i+\frac{3}{2}}^{k}\right)+\hat{R}_{i}^{k} \\ 2 \leqslant i \leqslant m-2,1 \leqslant k \leqslant n . \\ U_{m-1}^{k}=-\frac{2 \nu}{W_{m-1}^{k}}\left(\hat{c}_{4} \delta_{x} W_{m-\frac{7}{2}}^{k}+\hat{c}_{3} \delta_{x} W_{m-\frac{5}{2}}^{k}+\hat{c}_{2} \delta_{x} W_{m-\frac{3}{2}}^{k}+\hat{c}_{1} \delta_{x} W_{m-\frac{1}{2}}^{k}\right)+\hat{R}_{m-1}^{k} \\ 1 \leqslant k \leqslant n \end{gathered}\tag{1.95-1.97}$

存在常數 $c_{11}$ 使得

$\left|\hat{R}_{i}^{k}\right| \leqslant c_{11} h^{4}, \quad 1 \leqslant i \leqslant m-1,1 \leqslant k \leqslant n. \tag{1.98}$

在 (1.95)-(1.97) 中略去小量項, 得到如下計算格式

$\begin{gathered} u_{1}^{k}=-\frac{2 \nu}{w_{1}^{k}}\left(\hat{c}_{1} \delta_{x} w_{\frac{1}{2}}^{k}+\hat{c}_{2} \delta_{x} w_{\frac{3}{2}}^{k}+\hat{c}_{3} \delta_{x} w_{\frac{5}{2}}^{k}+\hat{c}_{4} \delta_{x} w_{\frac{7}{2}}^{k}\right), \quad 1 \leqslant k \leqslant n \\ u_{i}^{k}=-\frac{2 \nu}{w_{i}^{k}}\left(\hat{c}_{5} \delta_{x} w_{i-\frac{3}{2}}^{k}+\hat{c}_{6} \delta_{x} w_{i-\frac{1}{2}}^{k}+\hat{c}_{6} \delta_{x} w_{i+\frac{1}{2}}^{k}+\hat{c}_{5} \delta_{x} w_{i+\frac{3}{2}}^{k}\right) \\ 2 \leqslant i \leqslant m-2,1 \leqslant k \leqslant n, \\ u_{m-1}^{k}=-\frac{2 \nu}{w_{m-1}^{k}}\left(\hat{c}_{4} \delta_{x} w_{m-\frac{7}{2}}^{k}+\hat{c}_{3} \delta_{x} w_{m-\frac{5}{2}}^{k}+\hat{c}_{2} \delta_{x} w_{m-\frac{3}{2}}^{k}+\hat{c}_{1} \delta_{x} w_{m-\frac{1}{2}}^{k}\right) \\ 1 \leqslant k \leqslant n \end{gathered}\tag{1.99-1.101}$

利用定理 $1.12$ 的結果可以證明

$\sqrt{h \sum_{i=1}^{m-1}\left(U_{i}^{k}-u_{i}^{k}\right)^{2}} \leqslant c_{12}\left(\tau^{2}+h^{4}\right), \quad 1 \leqslant k \leqslant n$

$1.5$ 小結與延拓

本章討論了 Burgers 方程的差分方法. 首先證明了問題 (1.1)-(1.3) 的解滿足能量守恒性. 接著在 $1.2$ 節(jié)和 $1.3$ 節(jié)分別介紹了二層非線性差分格式和三層線性化差分格式. 證明了差分格式解的存在性、唯一性抚芦、有界性和收斂性. 三層線性化差分格式的有關結果主要取材于 [30].

對于問題 (1.1)-(1.3) 可建立如下二層線性化差分格式

$\begin{aligned} &\delta_{t} u_{i}^{k+\frac{1}{2}}+\frac{1}{2}\left(u_{i}^{k} \Delta_{x} u_{i}^{k+1}+u_{i}^{k+1} \Delta_{x} u_{i}^{k}\right)=\nu \delta_{x}^{2} u_{i}^{k+\frac{1}{2}}, \\ &\quad 1 \leqslant i \leqslant m-1,0 \leqslant k \leqslant n-1 \\ &u_{i}^{0}=\varphi\left(x_{i}\right), \quad 1 \leqslant i \leqslant m-1 \\ &u_{0}^{k}=0, \quad u_{m}^{k}=0, \quad 0 \leqslant k \leqslant n . \end{aligned}\tag{1.102-1.104}$

可以證明差分格式 (1.102)-(1.104) 是唯一可解的, 在無窮范數下關于時間步長和空間步長均是二階收玫的.

文 [39] 研究了二維 Burgers 方程的二階差分方法.

應用 Hopf-Cole 變換可將 Burgers 方程的初邊值問題 (1.1)-(1.3) 變?yōu)榫€性的熱傳導方程的初邊值問題 (1.67)-(1.69). 對 $(1.67)-(1.69)$ , 我們建立了緊致差分格式 $(1.78)-(1.81) .$ 證明了 (1.78)-(1.81) 解的存在性和唯一性以及解在無窮范數下關于時間步長 2 階、空間步長 4 階的收斂性. 如果在 (1.78)-(1.81) 中用單位算子 $\mathcal{I}$ 代替平均值算子 $\mathcal{A}$ , 得到如下格式

$\begin{aligned} &\delta_{t} w_{0}^{k+\frac{1}{2}}-\nu\left(\frac{2}{h} \delta_{x} w_{\frac{1}{2}}^{k+\frac{1}{2}}\right)=0, \quad 0 \leqslant k \leqslant n-1 \\ &\delta_{t} w_{i}^{k+\frac{1}{2}}-\nu \delta_{x}^{2} w_{i}^{k+\frac{1}{2}}=0, \quad 1 \leqslant i \leqslant m-1,0 \leqslant k \leqslant n-1 \\ &\delta_{t} w_{m}^{k+\frac{1}{2}}-\left(-\frac{2}{h} \delta_{x} w_{m-\frac{1}{2}}^{k+\frac{1}{2}}\right)=0, \quad 0 \leqslant k \leqslant n-1 \\ &w_{i}^{0}=\widetilde{\varphi}\left(x_{i}\right), \quad 0 \leqslant i \leqslant m \end{aligned}$

該差分格式是唯一可解的, 在無窮范數下關于時間步長和空間步長均是二階收斂的.

我們借助于 Browder 定理證明了非線性方程組 (1.19)-(1.20) 解的存在性. 與 Browder 定理相伴的還有一個 Leray-Schauder 定理 [43]. 設 $H$ 是一個有限維內積空間, $\Vert\cdot\Vert$ 是導出范數. 考慮 $H \rightarrow H$ 的算子 $T_{\lambda}(w)$ , 其中 $\lambda \in[0,1]$ 為參數. 如果 $T_{\lambda}(w)$ 滿足如下條件:

(a) $T_{\lambda}(w)$ 是 $H$ 上的連續(xù)算子;

(b) $T_{0}(w)=0$ 有唯一解;

則對任意 $\lambda \in[0,1], T_{\lambda}(w)=0$ 存在解. 特別地, $T_{1}(w)=0$ 存在解.

現在用上述結論來證明定理 $1.4$ , 即證明 $(1.19)-(1.20)$ 存在解.

令 $H=\dot{\mathcal{U}}_{h} .$ 對任意的 $w \in \stackrel{\circ}{\mathcal{U}}_{h}$ , 定義

$\begin{aligned} &T_{\lambda}(w)_{i}=\frac{2}{\tau}\left(w_{i}-u_{i}^{k}\right)+\lambda \psi(w, w)_{i}-\nu \delta_{x}^{2} w_{i}, \quad 1 \leqslant i \leqslant m-1 \\ &T_{\lambda}(w)_{0}=0 \\ &T_{\lambda}(w)_{m}=0 \end{aligned}$

易知 (a) $T_{\lambda}(w)$ 是連續(xù)的; (b) $T_{0}(w)_{i}=0, i=1,2, \cdots, m-1$ 是一個嚴格對角占優(yōu)的三對角線性方程組, 故有唯一解. 現在來檢驗 (c). 設 $w$ 是 $T_{\lambda}(w)=0$ 可能的解. 用 $w$ 和 $T_{\lambda}(w)=0$ 作內積, 得

$\frac{2}{\tau}\left((w, w)-\left(u^{k}, w\right)\right)+\lambda(\psi(w, w), w)-\nu\left(\delta_{x}^{2} w, w\right)=0$

利用

$(\psi(w, w), w)=0, \quad-\left(\delta_{x}^{2} w, w\right)=|w|_{1}^{2}$

得

$\frac{2}{\tau}\left(\|w\|^{2}-\left(u^{k}, w\right)\right)+\nu|w|_{1}^{2}=0$

于是

$\|w\|^{2} \leqslant\left(u^{k}, w\right) \leqslant\left\|u^{k}\right\| \cdot\|w\|$

易知

$\|w\| \leqslant\left\|u^{k}\right\|$

條件 (c) 滿足. 由 Leray-Schauder 定理. (1.19)-(1.20) 存在解.

?著作權歸作者所有,轉載或內容合作請聯系作者

人面猴
序言：七十年代末亦歉，一起剝皮案震驚了整個濱河市，隨后出現的幾起案子，更是在濱河造成了極大的恐慌截碴，老刑警劉巖溶诞，帶你破解...
沈念sama閱讀 207,113評論 6贊 481
死咒
序言：濱河連續(xù)發(fā)生了三起死亡事件，死亡現場離奇詭異拧篮，居然都是意外死亡词渤，警方通過查閱死者的電腦和手機，發(fā)現死者居然都...
沈念sama閱讀 88,644評論 2贊 381
救了他兩次的神仙讓他今天三更去死
文/潘曉璐我一進店門串绩，熙熙樓的掌柜王于貴愁眉苦臉地迎上來缺虐，“玉大人，你說我怎么就攤上這事礁凡〕锰遥” “怎么了皮获？”我有些...
開封第一講書人閱讀 153,340評論 0贊 344
道士緝兇錄：失蹤的賣姜人
文/不壞的土叔我叫張陵，是天一觀的道長。經常有香客問我撩穿，道長，這世上最難降的妖魔是什么摩瞎？我笑而不...
開封第一講書人閱讀 55,449評論 1贊 279
?港島之戀（遺憾婚禮）
正文為了忘掉前任篮绿，我火速辦了婚禮，結果婚禮上运挫，老公的妹妹穿的比我還像新娘状共。我一直安慰自己，他們只是感情好谁帕，可當我...
茶點故事閱讀 64,445評論 5贊 374
惡毒庶女頂嫁案：這布局不是一般人想出來的
文/花漫我一把揭開白布峡继。她就那樣靜靜地躺著，像睡著了一般匈挖。火紅的嫁衣襯著肌膚如雪碾牌。梳的紋絲不亂的頭發(fā)上，一...
開封第一講書人閱讀 49,166評論 1贊 284
城市分裂傳說
那天儡循，我揣著相機與錄音舶吗，去河邊找鬼。笑死贮折，一個胖子當著我的面吹牛裤翩，可吹牛的內容都是我干的。我是一名探鬼主播，決...
沈念sama閱讀 38,442評論 3贊 401
雙鴛鴦連環(huán)套：你想象不到人心有多黑
文/蒼蘭香墨我猛地睜開眼踊赠，長吁一口氣：“原來是場噩夢啊……” “哼呵扛！你這毒婦竟也來了？” 一聲冷哼從身側響起筐带，我...
開封第一講書人閱讀 37,105評論 0贊 261
萬榮殺人案實錄
序言：老撾萬榮一對情侶失蹤今穿，失蹤者是張志新（化名）和其女友劉穎，沒想到半個月后伦籍，有當地人在樹林里發(fā)現了一具尸體蓝晒，經...
沈念sama閱讀 43,601評論 1贊 300
?護林員之死
正文獨居荒郊野嶺守林人離奇死亡，尸身上長有42處帶血的膿包…… 初始之章·張勛以下內容為張勛視角年9月15日...
茶點故事閱讀 36,066評論 2贊 325
?白月光啟示錄
正文我和宋清朗相戀三年帖鸦，在試婚紗的時候發(fā)現自己被綠了芝薇。大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
茶點故事閱讀 38,161評論 1贊 334
活死人
序言：一個原本活蹦亂跳的男人離奇死亡作儿，死狀恐怖洛二，靈堂內的尸體忽然破棺而出，到底是詐尸還是另有隱情攻锰，我是刑警寧澤晾嘶，帶...
沈念sama閱讀 33,792評論 4贊 323
?日本核電站爆炸內幕
正文年R本政府宣布，位于F島的核電站娶吞，受9級特大地震影響垒迂，放射性物質發(fā)生泄漏。R本人自食惡果不足惜妒蛇，卻給世界環(huán)境...
茶點故事閱讀 39,351評論 3贊 307
男人毒藥：我在死后第九天來索命
文/蒙蒙一机断、第九天我趴在偏房一處隱蔽的房頂上張望。院中可真熱鬧材部，春花似錦毫缆、人聲如沸。這莊子的主人今日做“春日...
開封第一講書人閱讀 30,352評論 0贊 19
一樁弒父案，背后竟有這般陰謀
文/蒼蘭香墨我抬頭看了看天上的太陽浸颓。三九已至物臂，卻和暖如春，著一層夾襖步出監(jiān)牢的瞬間产上，已是汗流浹背棵磷。一陣腳步聲響...
開封第一講書人閱讀 31,584評論 1贊 261
情欲美人皮
我被黑心中介騙來泰國打工，沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留晋涣，地道東北人仪媒。一個月前我還...
沈念sama閱讀 45,618評論 2贊 355
代替公主和親
正文我出身青樓，卻偏偏與公主長得像谢鹊，于是被迫代替她去往敵國和親算吩。傳聞我的和親對象是個殘疾皇子留凭，可洞房花燭夜當晚...
茶點故事閱讀 42,916評論 2贊 344

Hopf-Cole 變換與高階差分格式

非線性發(fā)展方程的有限差分方法