logo

非線形1階偏微分方程式の特性方程式 📂偏微分方程式

非線形1階偏微分方程式の特性方程式

  • xとpに関して、偏微分方程式の変数を強調する時は普通のフォントでx,pRnx,p \in \mathbb{R}^{n}と表記し、ssの関数であることを強調する時は太字でx,pRn\mathbf{x}, \mathbf{p} \in \mathbb{R}^{n}と表記する。

特性メソッド1

開集合 ΩRn\Omega \subset \mathbb{R}^{n}が与えられたとしよう。uC2(Ω)u\in C^{2}(\Omega)が以下の非線形1階偏微分方程式の解だとする。

F(Du, u, x)=0 F(Du,\ u,\ x)=0

そしてx,z,p\mathbf{x}, z, \mathbf{p}を以下のように置こう。

x(s)=(x1(s),,xn(s))C1(I;Ω)(sIR)z(s)=u(x(s))p(s)=Du(x(s)) \begin{align*} \mathbf{x}(s) &=(x^{1}(s), \dots, x^{n}(s)) \in C^1(I;\Omega)\quad (s\in I \subset \mathbb{R}) \\ z(s) &= u(\mathbf{x}(s)) \\ \mathbf{p}(s) &= Du(\mathbf{x}(s)) \end{align*}

ここでx\mathbf{x}が以下の式を満たすと仮定しよう。

x˙(s)=DpF(p(s),z(s),x(s)) \dot{\mathbf{x}}(s) = D_{p}F\big(\mathbf{p}(s), z(s), \mathbf{x}(s) \big)

するとp(s)\mathbf{p}(s)z(s)z(s)はそれぞれ以下の常微分方程式の解となる。

{p˙(s)=DxF(p(s), z(s), x(s))DzF(p(s), z(s), x(s))p(s)z˙(s)=DpF(p(s), z(s), x(s))p(s)x˙(s)=DpF(p(s), z(s), x(s)) \begin{cases} \dot{\mathbf{p}} (s) = -D_{x}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big)-D_{z}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big)\mathbf{p}(s) \\ \dot{z}(s) = D_{p}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big) \cdot \mathbf{p}(s) \\ \dot{\mathbf{x}}(s) = D_{p}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big) \end{cases}

説明

特性メソッドは非線形1階偏微分方程式を解く方法の一つで、一つの偏微分方程式を連立常微分方程式に表して解くことである。

導出

以下のような非線形1階偏微分方程式が与えられたとする。

F(Du,u,x)=0 \begin{equation} F(Du, u, x) = 0 \label{eq1} \end{equation}

このとき、ΩRn\Omega \subset \mathbb{R}^{n}開集合であり、FC(Rn×R×Ω)F\in C^{\infty}(\mathbb{R}^n \times \mathbb{R} \times \Omega)と仮定する。そして、以下のような境界条件が与えられたとする。

u=g on Γ \begin{equation} u=g \quad \text{ on } \Gamma \label{eq2} \end{equation}

特性メソッドのアイデアは、与えられた境界条件を利用して固定されたxΩx\in\Omegax0Γx^0 \in \Gammaを結ぶΩ\Omega内の線を見つけることである。言い換えれば、uC2u \in C^{2}(eq1),(eq2)\eqref{eq1}, \eqref{eq2}を満たす解であるとき、次の図のような線に沿ってuuの値を得ることが目標である。

1.JPG

図の線が以下の関数で表されるとしよう。

x(s)=(x1(s), , xn(s)),sIR \mathbf{x}(s)=\big( x^1(s),\ \cdots,\ x^n(s) \big), \quad s\in I\subset \mathbb{R}

x(0)=x0\mathbf{x}(0)=x^{0}として、x\mathbf{x}ssが増加するにつれてΩ\Omega内の線に沿って動くと考えよう。ssによるuuの値をzzという関数で表現しよう。

z(s):=u(x(s)) \begin{equation} z(s):= u(\mathbf{x}(s)) \label{eq3} \end{equation}

同様に、ssによるDuDuの値をp\mathbf{p}としよう。

p(s):=Du(x(s)) \mathbf{p}(s) := Du(\mathbf{x}(s))

するとp(s)=(p1(s),,pn(s))\mathbf{p}(s)=(p^{1}(s), \dots, p^{n}(s))のとき、各成分は以下のようになる。

pi(s)=uxi(x(s)) \begin{equation} p^i(s)=u_{x_{i}}( \mathbf{x}(s)) \label{eq4} \end{equation}

ssに対する微分を簡単にdfds=f˙\dfrac{d f}{ds} = \dot{f}と表記しよう。dpi(s)=uxix1dx1+uxixndxndp^i(s)=u_{x_{i}x_{1}}dx^1+\cdots u_{x_{i}x_{n}}dx^nなのでp˙i(s)\dot{p}^i(s)は以下のようになる。

p˙i(s)=j=1nuxixj(x(s))x˙j(s) \begin{equation} \dot{p}^i(s)=\sum \limits_{j=1}^n u_{x_{i}x_{j}}( \mathbf{x}(s))\dot{x}^j(s) \label{eq5} \end{equation}

このときdF=Fp1dp1+Fpndpn+Fzdz+Fx1dx++FxndxndF=F_{p_{1}}dp_{1}+\cdots F_{p_{n}}dp_{n}+F_{z}dz+F_{x_{1}}dx_{} + \cdots + F_{x_{n}}dx_{n}なので(eq1)\eqref{eq1}xix_{i}に対して微分すると以下のようになる。

j=1nFpj(p(s),z(s),x(s))pxij(s)+Fz(p(s),z(s),x(s))zxi(s)+Fxi(p(s),z(s),x(s))=0 \sum \limits_{j=1}^n F_{p_{j}}\big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big)p^j_{x_{i}}(s)+F_{z} \big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big)z_{x_{i}}(s)+F_{x_{i}} \big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big) = 0

このとき(eq3),(eq4)\eqref{eq3}, \eqref{eq4}によってzxi=uxi=piz_{x_{i}}=u_{x_{i}}=p^{i}であり、(eq4)\eqref{eq4}によってpxij=uxjxip^{j}_{x_{i}}=u_{x_{j}x_{i}}なので、これを上記の式に代入すると以下を得る。

j=1nFpj(p(s),z(s),x(s))uxjxi(x(s))+Fz(p(s),z(s),x(s))pi(s)+Fxi(p(s),z(s),x(s))=0 \begin{equation} \sum \limits_{j=1}^{n} F_{p_{j}}\big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big)u_{x_{j}x_{i}}(\mathbf{x}(s))+F_{z}\big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big) p^{i}(s)+F_{x_{i}}\big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big) = 0 \label{eq6} \end{equation}

しかし、上記の式を見ると(eq5)\eqref{eq5}のためにuxixju_{x_{i}x_{j}}という2階微分項が生じたことが分かる。1階偏微分方程式を解くために2階微分を解くのは適切でないため、この項を取り除く必要がある。このために、x˙\dot{\mathbf{x}}の各成分が以下を満たすと仮定しよう。

x˙j(s)=Fpj(p(s),z(s),x(s)), sIR \begin{equation} \dot{x}^{j}(s)=F_{p_{j}}(\mathbf{p}(s), z(s), \mathbf{x}(s)), \quad \forall\ s\in I \subset \mathbb{R} \label{eq7} \end{equation}

(eq6)\eqref{eq6}(eq7)\eqref{eq7}を代入して(eq5)\eqref{eq5}を利用すると以下のような式を導き出せる。

p˙i(s)=Fz(p(s),z(s),x(s))pi(s)Fxi(p(s),z(s),x(s)), sI(i=1,,n) \begin{equation} \dot{p}^{i} (s) = -F_{z} \big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big)p^i(s)-F_{x_{i}} \big( \mathbf{p}(s), z(s), \mathbf{x}(s) \big), \quad \forall\ s\in I(i=1,\dots,n) \label{eq8} \end{equation}

また、(eq3)\eqref{eq3}ssに対して微分すると以下のようになる。

z˙(s)=j=1nuxj(x(s))x˙j(s) \dot{z}(s)=\sum \limits_{j=1}^n u_{x_{j}}(\mathbf{x}(s))\dot{x}^j(s)

ここに(eq4)pi=uxi\eqref{eq4} p^{i}=u_{x_{i}}(eq7)x˙j=Fpj\eqref{eq7} \dot{x}^{j} = F_{p_{j}}を代入すると以下を得る。

z˙(s)=j=1npj(s)Fpj(p(s), z(s), x(s)),sIR, \begin{equation} \dot{z}(s)=\sum \limits_{j=1}^n p^j(s)F_{p_{j}}\big( \mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big), \quad s\in I \subset \mathbb{R}, \label{eq9} \end{equation}

このような過程で得られた(eq7),(eq8),(eq9)\eqref{eq7}, \eqref{eq8}, \eqref{eq9}を以下のように整理して2n+12n+1の連立方程式として表現したものが特性方程式である。また、各未知数p(s),z(s),x(s)\mathbf{p}(s), z(s), \mathbf{x}(s)(eq1)\eqref{eq1}特性という。

{p˙(s)=DxF(p(s), z(s), x(s))DzF(p(s), z(s), x(s))p(s)z˙(s)=DpF(p(s), z(s), x(s))p(s)x˙(s)=DpF(p(s), z(s), x(s)) \begin{cases} \dot{\mathbf{p}} (s) = -D_{x}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big)-D_{z}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big)\mathbf{p}(s) \\ \dot{z}(s) = D_{p}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big) \cdot \mathbf{p}(s) \\ \dot{\mathbf{x}}(s) = D_{p}F\big(\mathbf{p}(s),\ z(s),\ \mathbf{x}(s) \big) \end{cases}


  1. Lawrence C. Evans, Partial Differential Equations (2nd Edition, 2010), p96-98 ↩︎