デュアメルの原理

数学の、特に偏微分方程式の分野で用いられるデュアメルの原理（デュアメルのげんり、英: Duhamel's principle）とは、熱方程式や波動方程式や振動板方程式（英語版）などの非同次（英語版）線形発展方程式の解を得るための一般的な手法である。薄い板を底から温める際の熱の分布のモデルとしての非同次熱方程式に対して初めてこの原理を利用した、ジャン＝マリー・デュアメル（英語版）の名にちなむ。デュアメルの原理は、調和振動子のような空間依存性を持たない線型発展方程式に対しては、線型同次常微分方程式を解く際に用いられる定数変化法に帰着される。 デュアメルの原理の根本となるアイデアは、コーシー問題（あるいは初期値問題）の解から非同次問題の解を得ることが可能、というものである。例えば、Rn 内の熱エネルギー u の分布をモデル化する熱方程式の例を考える。このときの初期値問題は { u t ( x , t ) − Δ u ( x , t ) = 0 ( x , t ) ∈ R n × ( 0 , ∞ ) u ( x , 0 ) = g ( x ) x ∈ R n {\displaystyle {\begin{cases}u_{t}(x,t)-\Delta u(x,t)=0&(x,t)\in \mathbf {R} ^{n}\times (0,\infty )\\u(x,0)=g(x)&x\in \mathbf {R} ^{n}\end{cases}}} となる。ただし g は初期の熱分布である。この熱方程式に対応する非同次の問題は { u t ( x , t ) − Δ u ( x , t ) = f ( x , t ) ( x , t ) ∈ R n × ( 0 , ∞ ) u ( x , 0 ) = 0 x ∈ R n {\displaystyle {\begin{cases}u_{t}(x,t)-\Delta u(x,t)=f(x,t)&(x,t)\in \mathbf {R} ^{n}\times (0,\infty )\\u(x,0)=0&x\in \mathbf {R} ^{n}\end{cases}}} のように表される。ここで ƒ(x,t)dt は各点に加えられる外的な熱エネルギーである。直感的に、この非同次問題は、各時間 t = t0 ごとに考えられる同次問題を集めたものと考えられるであろう。線型性により、その同次問題の解を時間 t0 毎に足し上げる（積分する）ことで、求めたい非同次問題の解を得ることが出来る。この考えがデュアメルの原理の本質である。 形式的に、ある函数 u : D × ( 0 , ∞ ) → R {\displaystyle u:D\times (0,\infty )\to \mathbf {R} } と Rn 内の空間領域 D に関する線型の非同次発展方程式として、次のような形状のものを考える。 { u t ( x , t ) − L u ( x , t ) = f ( x , t ) ( x , t ) ∈ D × ( 0 , ∞ ) u | ∂ D = 0 u ( x , 0 ) = 0 x ∈ D . {\displaystyle {\begin{cases}u_{t}(x,t)-Lu(x,t)=f(x,t)&(x,t)\in D\times (0,\infty )\\u|_{\partial D}=0&\\u(x,0)=0&x\in D.\end{cases}}} ここで L は、時間に関する微分は含まない線型の微分作用素である。 デュアメルの原理は、形式的に、この問題の解は u ( x , t ) = ∫ 0 t ( P s f ) ( x , t ) d s {\displaystyle u(x,t)=\int _{0}^{t}(P^{s}f)(x,t)\,ds} で与えられるということを述べたものである。ただし Psƒ は次の問題の解である。 { u t − L u = 0 ( x , t ) ∈ D × ( s , ∞ ) u | ∂ D = 0 u ( x , s ) = f ( x , s ) x ∈ D . {\displaystyle {\begin{cases}u_{t}-Lu=0&(x,t)\in D\times (s,\infty )\\u|_{\partial D}=0&\\u(x,s)=f(x,s)&x\in D.\end{cases}}} デュアメルの原理は（ベクトル値函数 u を伴う）線型システムに対しても同様に適用される。またデュアメルの原理は、波動方程式に現れるような高次の t 微分に対する一般化も与える（下記参照）。この原理の正当性は、同次問題を適切な函数空間の上で解くことが出来るか、またその解は積分が well-defined となるような良いパラメータ依存性を示すものであるか、という点に依存する。u および f に関する正しい解析的条件は、その適用される場面に依存するものである。 次の非同次の波動方程式を考える。 u t t − c 2 u x x = f ( x , t ) {\displaystyle u_{tt}-c^{2}u_{xx}=f(x,t)\,} ただし初期条件は u ( x , 0 ) = u t ( x , 0 ) = 0. {\displaystyle u(x,0)=u_{t}(x,0)=0.\,} とする。この解は次のように書ける。 u ( x , t ) = 1 2 c ∫ 0 t ∫ x − c ( t − s ) x + c ( t − s ) f ( ξ , s ) d ξ d s . {\displaystyle u(x,t)={\frac {1}{2c}}\int _{0}^{t}\int _{x-c(t-s)}^{x+c(t-s)}f(\xi ,s)\,d\xi \,ds.\,} デュアメルの原理とは、非同次の線型偏微分方程式の解は、初めにある入力ステップに対する解を見つけて、それをデュアメルの積分（英語版）を使って積み重ねることによって求めることが出来る、という結果である。次の定数係数 m 階非同次常微分方程式を考える。 P ( ∂ t ) u ( t ) = F ( t ) {\displaystyle P(\partial _{t})u(t)=F(t)\,} ∂ t j u ( 0 ) = 0 , 0 ≤ j ≤ m − 1 {\displaystyle \partial _{t}^{j}u(0)=0,\;0\leq j\leq m-1} ただし P ( ∂ t ) := a m ∂ t m + ⋯ + a 1 ∂ t + a 0 , a m ≠ 0. {\displaystyle P(\partial _{t}):=a_{m}\partial _{t}^{m}+\cdots +a_{1}\partial _{t}+a_{0},\;a_{m}\neq 0.} である。この非同次方程式の解を、以下に示す方法によってある同次方程式の解に書き下すことが出来る。以下の手順はすべて形式的に行われるが、解が well-defined となるために必要ないくつかの設定は無視している。 初めに、次の方程式の解 G を求める。 P ( ∂ t ) G = 0 , ∂ t j G ( 0 ) = 0 , 0 ≤ j ≤ m − 2 , ∂ t m − 1 G ( 0 ) = 1 / a m . {\displaystyle P(\partial _{t})G=0,\;\partial _{t}^{j}G(0)=0,\quad 0\leq j\leq m-2,\;\partial _{t}^{m-1}G(0)=1/a_{m}.} ここで H = G χ [ 0 , ∞ ) {\displaystyle H=G\chi _{[0,\infty )}} 　定める。ただし χ [ 0 , ∞ ) {\displaystyle \chi _{[0,\infty )}} は区間 [ 0 , ∞ ) {\displaystyle [0,\infty )} に対する指示函数である。このとき、 P ( ∂ t ) H = δ {\displaystyle P(\partial _{t})H=\delta } が超函数の意味で成立する。したがって u ( t ) = ( H ∗ F ) ( t ) {\displaystyle u(t)=(H\ast F)(t)} = ∫ 0 ∞ G ( τ ) F ( t − τ ) d τ {\displaystyle =\int _{0}^{\infty }G(\tau )F(t-\tau )\,d\tau } = ∫ − ∞ t G ( t − τ ) F ( τ ) d τ {\displaystyle =\int _{-\infty }^{t}G(t-\tau )F(\tau )\,d\tau } が、元の常微分方程式の解となる。 より一般に、次の非同次の定数係数偏微分方程式を考える。 P ( ∂ t , D x ) u ( t , x ) = F ( t , x ) {\displaystyle P(\partial _{t},D_{x})u(t,x)=F(t,x)\,} ただし D x = 1 i ∂ ∂ x {\displaystyle D_{x}={\frac {1}{i}}{\frac {\partial }{\partial x}}\,} とする。以下に示す方法で、この非同次方程式の解はある同次方程式の解へと書き下すことが出来る。すべての手順は形式的に行われるが、解が well-defined となるための必要な設定は無視している。 はじめに、x についてフーリエ変換を行うことで、 P ( ∂ t , ξ ) u