logo

非線形一次微分方程式の境界の線形化 📂偏微分方程式

非線形一次微分方程式の境界の線形化

ビルドアップ

非線形1次微分方程式特性方程式を簡単に解く方法の一つに、定義域Ω\Omegaの境界であるΩ\partial \Omegaの小さい部分Γ\Gammaを直線にすることがある。これは常に可能であるので、境界上の点x0x^{0}の近くでは、最初から境界が直線だと仮定して問題にアプローチすることができる。これを境界の直線化という。

1.JPG

ΩRn\Omega \subset \mathbb{R}^{n}開集合とし、Ω\partial \OmegaC2C^{2}とする。そして、偏微分方程式FC1(Rn×R×Ωˉ)F \in C^{1}(\mathbb{R}^{n} \times \mathbb{R} \times \bar \Omega)が与えられたとする。また、次のような境界条件が与えられたとする。

{F(Du, u, x)=0in Ωu=gon Γ \begin{equation} \left\{ \begin{aligned} F(Du,\ u,\ x)&=0 && \text{in } \Omega \\ u&=g && \text{on } \Gamma \end{aligned} \right. \end{equation}

このときΓΩ\Gamma \subset \partial \Omegaであり、g:ΓRg : \Gamma \to \mathbb{R}である。

説明

境界上に固定された点x0Ωx^{0} \in \partial \Omegaがある。そして変換Φ : RnRn\Phi\ :\ \mathbb{R}^{n} \to \mathbb{R}^{n}を次のように定義する。

{Φ(Ω):=VΦ(x):=(Φ1(x), ,Φn(x))=(y1, , yn)=y,yVyi=Φi(x):=xi,xRn (i=1,,n1)yn=Φn(x):=xnγ(x1,,xn1),xRn \left\{ \begin{align*} \Phi (\Omega) &:= V \\ \Phi (x) &:= \left( \Phi^{1}(x),\ \cdots, \Phi^{n}(x) \right) = (y_{1},\ \cdots,\ y_{n})=y, \quad y \in V \\ y_{i}=\Phi^{i}(x) &:= x_{i}, \quad x \in \mathbb{R}^{n}\ (i=1,\cdots, n-1) \\ y_{n}=\Phi^{n}(x) &:= x_{n}-\gamma (x_{1},\cdots,x_{n-1}), \quad x \in \mathbb{R}^{n} \end{align*} \right.

すなわち、Φ\Phiは境界の特定の部分のnn番目の座標を00にする変換である。これは定義により全単射であることは明白である。従って、逆変換が存在し、これをΨ\Psiとする。

{lΨ(y):=Φ1(y)Ψi(y)=xi=yi,yRn (i=1,,n1)Ψn(y)=xn=yn+γ(x1,,xn1),xRn \left\{ \begin{align*} {l}\Psi (y) &:= \Phi^{{}-1}(y) \\ \Psi^{i}(y) &= x_{i}=y_{i}, \quad y \in \mathbb{R}^{n}\ (i=1,\cdots, n-1) \\ \Psi^{n}(y) &= x_{n}=y_{n}+\gamma (x_{1},\cdots,x_{n-1}), \quad x \in \mathbb{R}^{n} \end{align*} \right.

図で表すと、以下のようになる。

2.JPG

さて、ΓΩ\Gamma \subset \partial \Omegaが開集合であり、gC(Γ)g\in C(\Gamma)とする。そして固定された点x0Γx^{0} \in \Gammaが与えられたとする。そしてuC1(Ω)C(Ωˉ)u \in C^{1} (\Omega)\cap C(\bar \Omega)が境界条件(1)(1)を解く解であると仮定する。そしてvvを以下のように定義する。

v(y):=u(Ψ(y)) yV v(y) := u(\Psi (y)) \quad \forall\ y\in V

つまり、VVuuと同じ関数値を持つようにvvを定義したわけである。すると、次が成立する。

u(x)=v(Φ(x)) xΩ u(x)=v(\Phi (x)) \quad \forall\ x\in \Omega

では、Du,u,xDu, u, xVVでどのようになるか見てみよう。まずuxiu_{x_{i}}から計算してみると、次のようになる。

uxi(x)=k=1nvyk(Φ(x))Φxik(x) u_{x_{i}}(x)=\sum \limits _{k=1}^{n} v_{y_{k}} \left( \Phi (x) \right) \Phi^{k}_{x_{i}}(x)

従って、

Du(x)=(k=1nvyk(Φ(x))Φx1k(x), , k=1nvyk(Φ(x))Φxnk(x))=(vy1Φx11++vynΦx1n, , vy1Φxn1++vynΦxnn)=(vy1vy2vyn)(Φx11Φx21Φxn1Φx12Φx22Φxn2Φx1nΦx2nΦxnn)=Dv(Φ(x))DΦ(x) \begin{align*} Du(x) &= \left( \sum \limits _{k=1}^{n} v_{y_{k}} \left( \Phi (x) \right) \Phi^{k}_{x_{1}}(x),\ \cdots,\ \sum \limits _{k=1}^{n} v_{y_{k}} \left( \Phi (x) \right) \Phi^{k}_{x_{n}}(x) \right) \\ &= (v_{y_{1}}\Phi^{1}_{x_{1}}+\cdots+v_{y_{n}}\Phi^{n}_{x_{1}},\ \cdots ,\ v_{y_{1}}\Phi^{1}_{x_{n}}+\cdots+v_{y_{n}}\Phi^{n}_{x_{n}} ) \\ &= \begin{pmatrix} v_{y_{1}} & v_{y_{2}} & \cdots & v_{y_{n}} \end{pmatrix} \begin{pmatrix} \Phi^{1}_{x_{1}} & \Phi^{1}_{x_{2}} & \cdots &\Phi^{1}_{x_{n}} \\ \Phi^{2}_{x_{1}} & \Phi^{2}_{x_{2}} & \cdots & \Phi^{2}_{x_{n}} \\ \vdots & \vdots & \ddots & \vdots \\ \Phi^{n}_{x_{1}} & \Phi^{n}_{x_{2}} & \cdots & \Phi^{n}_{x_{n}} \end{pmatrix} \\ &= Dv\left( \Phi (x) \right) D\Phi (x) \end{align*}

あるいは、

Du(Ψ(y))=Dv(y)DΦ(Ψ(y)) Du(\Psi (y)) = Dv(y) D\Phi (\Psi (y))

すると、(1)(1)の式は次のようになる。

F(Du(Ψ(y)),u(Ψ(y)),Ψ(y))=F(Dv(y)DΦ(Ψ(y)),v(y),Ψ(y))=0 F\Big( Du(\Psi (y) ), u( \Psi (y) ), \Psi (y) \Big) = F\Big( Dv(y)D\Phi (\Psi (y)), v(y), \Psi (y) \Big)=0

さて、次のように非線形1次偏微分方程式を定義しよう。

G(q,w,y):=F(qDΦ(Ψ(y)),w,Ψ(y)),(q,w,y)Rn×R×Vˉ G(q, w, y):=F \Big (q D\Phi (\Psi (y) ), w, \Psi (y) \Big), \quad \forall (q, w, y)\in\mathbb{R}^{n}\times\mathbb{R}\times\bar V

するとGC1(Rn×R×Vˉ)G\in C^{1}(\mathbb{R}^{n} \times \mathbb{R} \times \bar V)であり、上記の結果から次の式が成立する。

G(Dv(y), v(y), y)=0,yV G \Big( Dv(y),\ v(y),\ y \Big)=0, \quad \forall y\in V

そして、Δ:=Φ(Γ)\Delta:=\Phi (\Gamma)であり、h(y):=g(Φ(y))yΔh(y):=g(\Phi (y)) \quad y \in \Deltaと定義する。するとΔ\Deltaは開集合であり、ΔV\Delta \subset \partial Vである。そして、Δ\Deltay0y^{0}の近くで平らである。要するに、ここで定義したvC1(V)C(Vˉ)v \in C^{1}(V) \cap C(\bar V)は下記の境界条件を満たす解になる。

{G(Dv, v, y)=0in Vv=hon ΔV \left\{ \begin{align*} G(Dv,\ v,\ y) &= 0 && \text{in } V \\ v &= h && \text{on } \Delta \subset \partial V \end{align*} \right.

これは、(1)(1)と境界の任意に選ばれた部分が平らであること以外は、すべて同じである。そして、常にこのように境界を平らにすることができるので、最初から与えられた問題がこれであると仮定して解くことができる。