logo

非粘性バーガース方程式における質量保存の法則 📂偏微分方程式

非粘性バーガース方程式における質量保存の法則

定理

{ut+uux=0,t>0u(t,x)=f(x),t=0 \begin{cases} u_{t} + u u_{x} = 0 & , t>0 \\ u(t,x) = f(x) & , t=0 \end{cases}

非粘性ブルガース方程式の解 uuについて、区間 [a,b][a,b]までの質量 MMを次のように定義する。

Ma,b(t):=abu(t,x)dx M_{a,b}(t) := \int_{a}^{b} u(t,x) dx

そして、破裂時間をtt_{\ast}とすると、t(0,t)t \in ( 0 , t_{\ast})に対して次が成立する。

ddtMa,b(t)=(12u2(t,b)12u2(t,a)) {{d} \over {dt}} M_{a,b}(t) = - \left( {{1} \over {2}} u^2 (t,b) - {{1} \over {2}} u^2 (t,a) \right)

説明

破裂時間は数学的には関数ではなく、物理的には同時に多数の状態が重なる瞬間を意味する。

式を簡単に言うと、質量の変化量は入ってくるものと出ていくものの純合計と同じになるという意味だ。

右辺の 12u2(t,b)\displaystyle {{1} \over {2}} u^2 (t,b)流出と呼び、12u2(t,a)\displaystyle {{1} \over {2}} u^2 (t,a)流入と呼ぶ。時間による Ma,b(t)M_{a,b}(t) の変化量が負であるということは、 [a,b][a,b] 上で uu が減少しているということで、つまり流出が流入よりも大きいということだ。逆に正であれば質量が増加しているということで、流入が流出よりも大きいということになり、表現がうまく一致する。

これを総括すると、F(u):=12u2\displaystyle F(u) := {{1} \over {2}} u^2流量関数と呼ぶことができる。

ut+uxF(u)=0 {{\partial u} \over {\partial t}} + {{\partial u} \over {\partial x}} F(u) =0

この感覚でブルガース方程式を**11次元における質量保存の法則**とも呼ぶ。

導出

uu が連続関数であるため、次が成立する。

ddtabu(t,x)dx=abtu(t,x)dx {{d } \over {dt} } \int_{a}^{b} u(t,x) dx = \int_{a}^{b} {{\partial } \over { \partial t} } u(t,x) dx

ut=uuxu_{t} = - u u_{x}であるため、次が成立する。

abtu(t,x)dx=abuuxdx \int_{a}^{b} {{\partial } \over { \partial t} } u(t,x) dx = - \int_{a}^{b} u u_{x} dx

uux=x(12u2)\displaystyle u u_{x} = {{\partial} \over {\partial x}} \left( {{1 } \over {2}} u^2 \right)であり、かつut=uuxu_{t} = - u u_{x}であるため、次が成立する。

abuuxdx=abx(12u2)dx=(12u2(t,b)12u2(t,a)) - \int_{a}^{b} u u_{x} dx = - \int_{a}^{b} {{\partial} \over {\partial x}} \left( {{1} \over {2}} u^2 \right) dx = - \left( {{1} \over {2}} u^2 (t,b) - {{1} \over {2}} u^2 (t,a) \right)


特にf(x)dx<\displaystyle \int_{-\infty}^{\infty} f(x) dx < \inftyの場合、t[0,t)t \in [0,t_{\ast})に対して次が成立する。

u(t,x)dx=f(x)dx \int_{-\infty}^{\infty} u(t,x) dx = \int_{-\infty}^{\infty} f(x) dx