特性曲線法のソースを表示

[[数学]]において'''特性曲線法'''（とくせいきょくせんほう、{{Lang-en-short|method of characteristics}}）とは、[[偏微分方程式]]に対する一つの解法である。一般には[[一階偏微分方程式]]に対して適用されるが、任意の[[双曲型偏微分方程式]]に対するより一般の特性曲線法も存在する。この方法では偏微分方程式を、常微分方程式の族に書き下し、適切な[[超曲面]]上で与えられたいくつかの初期データより積分されることによってその線に沿った解が得られる。

== 一階偏微分方程式の特性曲線 ==

一階の[[偏微分方程式]]（PDE）に対する特性曲線法では、それが[[常微分方程式]]（ODE）となるようなある曲線（特性曲線あるいは単に特性線と呼ばれる）を探すことになる。そのようなODEが見つかれば、特性曲線に沿って解いた後に元のPDEに対して解を変換すれば良いことになる。

ここで、二つの独立変数 ''x'' と ''y'' の函数のケースを取り上げる。次の形の[[準線型]]{{要曖昧さ回避|date=2018-8-17}}偏微分方程式を考える：

{{NumBlk|:|<math>a(x,y,z) \frac{\partial z}{\partial x}+b(x,y,z) \frac{\partial z}{\partial y}=c(x,y,z).</math>|{{EquationRef|1}}}}

ここで、解 ''z'' が得られたとして、'''R'''<sup>3</sup> 内の曲面のグラフ ''z''&nbsp;=&nbsp;''z''(''x'',''y'') を考える。この曲面に対する[[法線ベクトル]]は次で与えられる。

:<math>\left(\frac{\partial z}{\partial x}(x,y),\frac{\partial z}{\partial y}(x,y),-1\right).\,</math>

これは次のようにして分かる。x,y方向の接ベクトルをそれぞれ<math>n_1</math>,<math>n_2</math>とすると、これらは例えば<math>n_1=(1,0,\partial z /\partial x ) dx</math>,<math>n_2 = (0,1,\partial z / \partial y ) dy</math>となる。これらの外積から上述の法線ベクトル（と平行なベクトル）が得られる。

したがって<ref>{{harvnb|John|1991}}</ref>式 ({{EquationNote|1}}) は、ベクトル場

:<math>(a(x,y,z),b(x,y,z),c(x,y,z))\,</math>

が全ての点において曲面 ''z''&nbsp;=&nbsp;''z''(''x'', ''y'') に接するという幾何学的な内容を意味する。言い換えると、解はこのベクトル場の[[積分曲線]]の合併となる。これらの積分曲線は、元の偏微分方程式の'''特性曲線'''と呼ばれる。

特性曲線の方程式は、'''ラグランジュ＝シャルピ方程式'''によって次のように不変な形で表すことが出来る<ref>{{harvnb|Delgado|1997}}</ref>：

:<math>\frac{dx}{a(x,y,z)} = \frac{dy}{b(x,y,z)} = \frac{dz}{c(x,y,z)}.</math>

また、この曲線のパラメータ化 ''t'' が固定された場合、これらの方程式は ''x''(''t''), ''y''(''t''), ''z''(''t'') に対する次の連立常微分方程式として書くことが出来る。

:<math>
\begin{align}
\frac{dx}{dt}&=a(x,y,z),\\
\frac{dy}{dt}&=b(x,y,z),\\
\frac{dz}{dt}&=c(x,y,z).
\end{align}
</math>

これらを元の偏微分方程式の'''特性方程式''' (characteristic equation) という<!-- [[Characteristic equations]] redirects to this article -->。

=== 線型と準線型の場合 ===

次の形式のPDEを考える。

:<math>\sum_{i=1}^n a_i(x_1,\dots,x_n,u) \frac{\partial u}{\partial x_i}=c(x_1,\dots,x_n,u).</math>

このPDEを[[線型性|線型]]とするためには、係数 ''a''<sub>''i''</sub> は空間変数のみに依存し、''u'' には独立とすればよい。準線型とするためには、''a''<sub>''i''</sub> はその函数の値にも依存するが、導函数には依存しないものとすればよい。これら二つのケースの区別は、ここでの議論では本質的ではない。

線型あるいは準線型のPDEに対し、特性曲線はパラメータ的に次で与えられる。

:<math>(x_1,\dots,x_n,u) = (x_1(s),\dots,x_n(s),u(s))</math>

但し次の常微分方程式系が満たされるものとする。

{{NumBlk|:|<math>\frac{dx_i}{ds} = a_i(x_1,\dots,x_n,u)</math>|{{EquationRef|2}}}}
{{NumBlk|:|<math>\frac{du}{ds} = c(x_1,\dots,x_n,u).</math>|{{EquationRef|3}}}}

式 ({{EquationNote|2}}) と ({{EquationNote|3}}) が、元のPDEの特性曲線である。

=== 完全に非線型の場合 ===

次の偏微分方程式を考える。

{{NumBlk|:|<math>F(x_1,\dots,x_n,u,p_1,\dots,p_n)=0</math>|{{EquationRef|4}}}}

ここで変数 ''p''<sub>i</sub> は次の偏微分を略記したものである。

:<math>p_i = \frac{\partial u}{\partial x_i}.</math>
'''R'''<sup>''n''+1</sup> 内の超曲面 (''x<sub>i</sub>'', ''u'') が偏微分方程式の解であるとする。解の超曲面の上にある任意の滑らかな(微分可能な)曲線を特性曲線と言い、''s'' を曲線長さに沿うパラメータとして、曲線上の各点は次のように表されるものとする。
:<math>u(s) = u(x_1(s),\dots,x_n(s)).</math>
また、解の曲面の方向は、この特性曲線の各点での接線の傾き <math>p_i = \frac{\partial u}{\partial x_i}</math> により指定されているとする。解に沿って ({{EquationNote|4}})を ''s'' に関して微分すると、次が得られる。
{{NumBlk|:|<math>\sum_i(F_{x_i} + F_u p_i)\dot{x}_i + \sum_i F_{p_i}\dot{p}_i = 0</math>|{{EquationRef|5}}}}
{{NumBlk|:|<math>\dot{u} - \sum_i p_i \dot{x}_i = 0</math>|{{EquationRef|6}}}}
{{NumBlk|:|<math>\sum_i (\dot{x}_i dp_i - \dot{p}_i dx_i)= 0.</math>|{{EquationRef|7}}}}
式(6)は、解 ''u'' に対して連鎖律を適用することで得られる。また式(7)の括弧内は
:<math>\dot{x_i}dp_i-\dot{p_i}dx_i=\frac{dx_i}{ds}\frac{dp_i}{ds}ds-\frac{dp_i}{ds}\frac{dx_i}{ds}ds=0.</math>
だから式(7)が成立する。

''λ'' をある定数として ''λ''×(5)+(7)/''ds''を 作ると次の式(8)となる。
{{NumBlk|:|<math>\sum_i (\lambda (F_{x_i}+F_u p_i) + \dot{p}_i)(dx_i/ds)+ \sum_i(\lambda F_{p_i}-\dot{x}_i)(dp_i/ds)=0</math>|{{EquationRef|8}}}}
式(8)はその点を通る任意の特性曲線 ''x''<sub>''i''</sub> (''s'')に対して成り立つ。特性曲線 ''x''<sub>''i''</sub>(''s'')が任意に変わると、''dx''<sub>''i''</sub>/''ds''および''dp''<sub>''i''</sub>/''ds''はそれに応じて変わってしまう変数である。それでも式(8)が成り立つためには、d''x''<sub>''i''</sub>/''ds''およびd''p''<sub>''i''</sub>/''ds''の係数は 0 でなければならない。よって
:<math>\quad\lambda(F_{x_i}+F_up_i)+\dot{p}_i=0,\quad\lambda F_{p_i}-\dot{x}_i=0,.</math>

:<math>\dot{x}_i=\lambda F_{p_i},\quad  \dot{p}_i=-\lambda(F_{x_i}+F_up_i).</math>

この<math>\dot{x}_i=\lambda F_{p_i}</math>を式(6)に入れると
:<math>\quad \dot{u}-\lambda\sum_i p_iF_{p_i}=0,\quad \dot{u}=\lambda\sum_i p_iF_{p_i}.</math>

ここに ''λ'' はある定数である。これらの式をより対称的に書くと、特性曲線に対する次のラグランジュ＝シャルピ方程式が得られる。

{{NumBlk|:|<math>\frac{\dot{x}_i}{F_{p_i}}=-\frac{\dot{p}_i}{F_{x_i}+F_up_i}=\frac{\dot{u}}{\sum p_iF_{p_i}}.</math>|{{EquationRef|9}}}}
幾何学的に、完全に非線型の場合の特性曲線法は、微分方程式の{{仮リンク|モンジュ錐|en|Monge cone}}が至る所で解のグラフに接することを要求するものとして解釈される。

== 例 ==
一例として、次の[[移流|移流方程式]]が挙げられる（この例ではPDEの記法や基本的なODEの解についてはよく知っているものと仮定する）。

:<math>a \frac{\partial u}{\partial x} + \frac{\partial u}{\partial t} = 0\,</math>

ここに <math>a\,</math> は定数で、<math>u\,</math> は <math>x\,</math> と <math>t\,</math> の函数である。この線型一階偏微分方程式を、適切な曲線に沿って常微分方程式に変換することを考える。すなわち、次の形状の常微分方程式を考える。

:<math> \frac{d}{ds}u(x(s), t(s)) = F(u, x(s), t(s)) </math>,

ここに <math>(x(s),t(s))\,</math> は特性曲線である。はじめに、連鎖律より次が得られる。

:<math>\frac{d}{ds}u(x(s), t(s)) = \frac{\partial u}{\partial x} \frac{dx}{ds} + \frac{\partial u}{\partial t} \frac{dt}{ds}</math>

今、<math> \frac{dx}{ds} = a</math> and <math>\frac{dt}{ds} = 1</math> とすると、次が得られる。

:<math> a \frac{\partial u}{\partial x} + \frac{\partial u}{\partial t}  \,</math>

これははじめのPDEの左辺である。したがって

:<math>\frac{d}{ds}u = a \frac{\partial u}{\partial x} + \frac{\partial u}{\partial t}  = 0</math>

が得られる。したがって、特性曲線 <math>(x(s), t(s))\,</math> に沿って、元の PDE は ODE <math>u_s = F(u, x(s), t(s)) = 0\,</math> になる。すなわち、特性曲線に沿って解は定数である。すると、<math>u(x_s, t_s) = u(x_0, 0)\,</math> となる。但し <math>(x_s, t_s)\,</math> と <math>(x_0, 0)\,</math> は同一の特性曲線上に存在する。したがって、一般解を決定する上で、次のような常微分方程式の特性システムを解けば十分である。

* <math>\frac{dt}{ds} = 1</math>。<math>t(0)=0\,</math> とすることで <math>t=s\,</math> が得られる。
* <math>\frac{dx}{ds} = a</math>。<math>x(0)=x_0\,</math> とすることで <math>x=as+x_0=at+x_0\,</math> が得られる。
* <math>\frac{du}{ds} = 0</math>。<math>u(0)=f(x_0)\,</math> とすることで <math>u(x(t), t)=f(x_0)=f(x-at)\,</math> が得られる。

この場合、特性曲線は傾きが <math>a\,</math> であるような直線で、任意の特性曲線に沿って <math>u\,</math> の値は定数となる。

== 線型微分作用素の特性曲線 ==

''X'' を[[可微分多様体]]とし、''P'' を次数 ''k'' の線型[[微分作用素]]

:<math>P : C^\infty(X) \to C^\infty(X)</math>

とする。局所座標系 ''x''<sup>''i''</sup> においては

:<math>P = \sum_{|\alpha|\le k} P^{\alpha}(x)\frac{\partial}{\partial x^\alpha}</math>

とする。ここで α は[[多重指数]]である。''P'' の[[微分作用素の表象|主表象]]は  σ<sub>''P''</sub> と表され、これはそれらの局所座標系において定義される[[余接束]] T<sup>∗</sup>''X'' に関する次の函数である。

:<math>\sigma_P(x,\xi) = \sum_{|\alpha|=k} P^\alpha(x)\xi_\alpha</math>

ここに ξ<sub>''i''</sub> は座標微分 d''x''<sup>''i''</sup> によって導かれる余接束上のファイバー座標である。これは特定の座標系を用いることで定義されるが、ξ<sub>''i''</sub> と ''x''<sup>''i''</sup> を関連づける変換則は σ<sub>''P''</sub> がその余接束上で well-defined な函数であることを保証するものである。

函数 σ<sub>''P''</sub> は変数 ξ について次数 ''k'' の[[斉次函数]]である。σ<sub>''P''</sub> のゼロ解は、T<sup>∗</sup>''X'' のゼロ切断とは離れた所にあり、''P'' の特性曲線である。式 ''F''(''x'')&nbsp;=&nbsp;''c'' によって定義される ''X'' の超曲面が、''x'' での特性超曲面であるとは、

:<math>\sigma_P(x,dF(x)) = 0</math>

が成立することを言う。特性超曲面は、{{仮リンク|法束|label=余法束|en|normal bundle}}が ''P'' の特性集合に属する超曲面である。

== 特性曲線の定性的解析 ==

特性曲線はまた、PDEへの定性的な洞察を得る上での強力な道具となる。

圧縮性流体におけるポテンシャルフローに対する[[衝撃波]]を見つけるために、特性曲線の交点を利用することが出来る。直感的に言うと、各特性曲線はそれ自身に沿った <math>u\,</math> の解を意味するものであると考えられる。したがって、二つの特性曲線が交わる場合には、函数は複数値ということになりこれは非物理的な解である。物理的にはこのような矛盾は、衝撃波の構成、接線の不連続性あるいは弱不連続性によって除外することが出来る。その結果、初めの仮定を満たさない非ポテンシャルフローが得られる。

特性曲線は、PDEの定義域の一部分をカバーしないこともある。この事実は{{仮リンク|希薄化|en|rarefaction}}と呼ばれ、弱い意味、すなわち[[積分方程式]]に対してのみ解が存在することを意味する。

特性曲線の方向は、上述の例で示したように、解に沿った値のフローを示すものである。この種の知識は、問題に対して[[有限差分]]が最適であるように示すものであるため、PDEを数値的に解く上で有用となる。

== 関連項目 ==
* {{仮リンク|量子特性曲線法|en|Method of quantum characteristics}}

== 脚注 ==
<references />

== 参考文献 ==
* {{citation|first1=Richard|last1=Courant|authorlink1=リヒャルト・クーラント|first2=David|last2=Hilbert|authorlink2=ダフィット・ヒルベルト|title=Methods of Mathematical Physics, Volume II|publisher=Wiley-Interscience|year=1962}}
*{{citation|jstor=2133111|last=Delgado|first=Manuel|title=The Lagrange-Charpit Method|journal=SIAM Review|volume=39|year=1997|pages=298–304|doi=10.1137/S0036144595293534|issue=2|bibcode = 1997SIAMR..39..298D }}
* {{citation|first=Lawrence C.|last=Evans|title=Partial Differential Equations|publisher=American Mathematical Society|publication-place=Providence|year=1998|isbn=0-8218-0772-2}}
* {{citation|first=Fritz|last=John|authorlink=:en:Fritz John|title=Partial differential equations|publisher=Springer|edition=4th |year=1991|isbn=978-0-387-90609-6}}
* {{citation|first1=A. D.|last1=Polyanin|first2=V. F.|last2=Zaitsev|first3=A.|last3=Moussiaux|title=Handbook of First Order Partial Differential Equations|publisher=Taylor & Francis|publication-place=London|year=2002|isbn=0-415-27267-X}}
* {{citation|first=A. D.|last=Polyanin|title=Handbook of Linear Partial Differential Equations for Engineers and Scientists|publisher=Chapman & Hall/CRC Press|publication-place=Boca Raton|year=2002|isbn=1-58488-299-9}}
* {{citation|last=Sarra|first=Scott|title=The Method of Characteristics with applications to Conservation Laws|journal=Journal of Online Mathematics and its Applications|year=2003}}.
*{{citation|last1=Streeter|first1=VL|last2=Wylie|first2=EB|title=Fluid mechanics|publisher=McGraw-Hill Higher Education|edition=International 9th Revised|year=1998}}

== 外部リンク ==
* [http://www.scottsarra.org/shock/shock.html Prof. Scott Sarra tutorial on Method of Characteristics]
* [http://www-solar.mcs.st-and.ac.uk/~alan/MT2003/PDE/node5.html Prof. Alan Hood tutorial on Method of Characteristics]
{{偏微分方程式の数値解法}}
{{DEFAULTSORT:とくせいきよくせんほう}}
[[Category:偏微分方程式]]
[[Category:数学に関する記事]]