Reynoldsの輸送定理
空間表示した函数$ \bm\theta(\bm x,t)の体積分の時間変化に関する変換則 $ \frac{\partial}{\partial t}\int_{B(t)}\bm\theta(\bm x,t)\mathrm dv=\int_{B(t)}\left(\frac{\mathrm D\bm\theta}{\mathrm Dt}+\bm \theta\bm\nabla\cdot\bm v\right)\mathrm dv
解釈
$ \frac{\partial}{\partial t}\int_{B(t)}\bm\theta(\bm x,t)\mathrm dv=\int_{B(t)}\left(\frac{\mathrm D\bm\theta}{\mathrm Dt}+\bm \theta\bm\nabla\cdot\bm v\right)\mathrm dv
$ =\int_{B(t)}\left(\frac{\partial\bm\theta}{\partial t}+\bm v\cdot\bm\nabla\bm\theta+\bm \theta\bm\nabla\cdot\bm v\right)\mathrm dv
$ =\int_{B(t)}\left(\frac{\partial\bm\theta}{\partial t}+\bm\nabla\cdot(\bm v\bm\theta)\right)\mathrm dv
$ =\underbrace{\int_{B(t)}\frac{\partial\bm\theta}{\partial t}\mathrm dv}_{領域内の時間変化}+\underbrace{\int_{\partial B(t)}\bm\theta\bm v\cdot\mathrm d\bm s}_{領域外からの流入}
運動する連続体$ B(t)上の量$ \bm\thetaの体積分の時間変化は、
時刻$ tにおける領域内の$ \bm\thetaの時間変化の合計
の和に等しい
別の言い方をすると、「時刻tにおけるB(t)内の物理量に、dt間の領域内の生成分と流入分を加えると、t+dtにおける領域の変形も加味した物理量になる」
$ \frac{\partial}{\partial t}\int_{B(t)}\bm\theta(\bm x,t)\mathrm dvは物質微分を用いて$ \frac{\mathrm D}{\mathrm D t}\int_{B(t)}\bm\theta(\bm x,t)\mathrm dvと書かれることがほとんどだが、ここでは時間の偏微分$ \frac{\partial}{\partial t}を用いる $ \int_{B(t)}\bm\theta(\bm x,t)\mathrm dvは空間位置$ \bm xでの体積分しているので、$ \bm xは束縛変数扱いになる
そのため、$ \int_{B(t)}\bm\theta(\bm x,t)\mathrm dvは$ tのみの函数となる
導出
ある時刻$ tにおける連続体の現位置の領域を$ B(t)とする $ B(t)=\{\bm x\in\R^3|\exist\bm X\in B_0;\bm x=\bm\phi(\bm X,t)\}という関係が成り立つ
ある空間表示された物質量$ \bm\theta(\bm x,t)を、$ B(t)で体積分した量は$ \int_{B(t)}\bm\theta(\bm x,t)\mathrm dvである
$ \mathrm dv(\bm x,t)：微小体素の体積の空間表示
これの時間変化率$ \frac{\partial}{\partial t}\int_{B(t)}\bm\theta(\bm x,t)\mathrm dvを考えたい
このままだと、積分範囲$ B(t)が$ tに依存しているので、$ \frac{\partial}{\partial t}を積分の中に入れられない
そこで、$ \bm xを$ \bm Xに置換積分する
$ \frac{\partial}{\partial t}\int_{B(t)}\bm\theta(\bm x,t)\mathrm dv
$ = \frac{\partial}{\partial t}\int_{\bm\phi(\bm X,t)\in B(t)}\bm\theta(\bm\phi(\bm X,t)\,t)(\det\bm\nabla\bm\phi)\mathrm dV
$ \underbrace{\mathrm dv(\bm\phi(\bm X,t),t)}_{時刻tでの微小体素の体積}=\det\bm\nabla\bm\phi(\bm X,t)\underbrace{\mathrm dV(\bm X)}_{初期時刻における微小体素の体積}
2024-04-04今なら$ \underbrace{\left.\mathrm dv\right|_{\bm x=\bm\phi(X,t)}}_{現配置の体積要素}=J\underbrace{\mathrm dV}_{基準配置の体積要素}と説明するかなtakker.icon
$ \mathrm dv(\bm X,0)=\mathrm dV(\bm X)という関係にある
$ = \frac{\partial}{\partial t}\int_{\bm X\in B_0}\bm\theta(\bm\phi(\bm X,t)\,t)J\mathrm dV
$ J:=\det\bm\nabla\bm\phi
$ = \int_{\bm X\in B_0}\frac{\partial}{\partial t}(\bm\theta(\bm\phi(\bm X,t)\,t)J)\mathrm dV
$ = \int_{\bm X\in B_0}\left(\left.\frac{\mathrm D\bm\theta}{\mathrm D t}\right|_{\bm x=\bm\phi(\bm X,t)}J+\bm\theta(\bm\phi(\bm X,t)\,t)\frac{\partial J}{\partial t}\right)\mathrm dV
$ = \int_{\bm X\in B_0}\left(\left.\frac{\mathrm D\bm\theta}{\mathrm D t}\right|_{\bm x=\bm\phi(\bm X,t)}+\bm\theta(\bm\phi(\bm X,t)\,t)\left.\bm\nabla\cdot\bm v\right|_{\bm x=\bm\phi(\bm X,t)}\right)J\mathrm dV
体積変化率の時間導函数$ \frac{\partial J}{\partial t}=J\left.\bm\nabla\cdot\bm v\right|_{\bm x=\bm\phi(\bm X,t)}を使った Jacobianの時間変化が、領域内への流入を担っているみたい？
$ = \int_{\bm\phi(\bm X,t)\in B(t)}\left.\left(\frac{\mathrm D\bm\theta}{\mathrm D t}+\bm\theta\bm\nabla\cdot\bm v\right)\right|_{\bm x=\bm\phi(\bm X,t)}\mathrm dv(\bm\phi(\bm X,t)\,t)
置換積分して、空間表示に直した
$ \underline{= \int_{B(t)}\left(\frac{\mathrm D\bm\theta}{\mathrm D t}+\bm\theta\bm\nabla\cdot\bm v\right)\mathrm dv\quad}_\blacksquare
$ \bm\theta=\rho\bm\psiのときに成り立つ式
$ \frac{\partial}{\partial t}\int_{B(t)}\rho\bm\psi\mathrm dv=\int_{B(t)}\left(\frac{\mathrm D\rho\bm\psi}{\mathrm Dt}+\rho\bm\psi\bm\nabla\cdot\bm v\right)\mathrm dv
$ =\int_{B(t)}\left(\rho\frac{\mathrm D\bm\psi}{\mathrm Dt}+\bm\psi\frac{\mathrm D\rho}{\mathrm Dt}+\rho\bm\psi\bm\nabla\cdot\bm v\right)\mathrm dv
$ =\int_{B(t)}\left(\rho\frac{\mathrm D\bm\psi}{\mathrm Dt}+\bm\psi\left(\frac{\mathrm D\rho}{\mathrm Dt}+\rho\bm\nabla\cdot\bm v\right)\right)\mathrm dv
$ =\int_{B(t)}\left(\rho\frac{\mathrm D\bm\psi}{\mathrm Dt}+\bm0\right)\mathrm dv
$ \underline{\therefore\frac{\partial}{\partial t}\int_{B(t)}\rho\bm\psi\mathrm dv=\int_{B(t)}\rho\frac{\mathrm D\bm\psi}{\mathrm Dt}\mathrm dv\quad}_\blacksquare
一般に、体積変化に不変($ \frac{\partial}{\partial t}\int_{B(t)}\xi\mathrm dv=0)な量$ \xiについて
$ \frac{\partial}{\partial t}\int_{B(t)}\xi\bm\psi\mathrm dv=\int_{B(t)}\xi\frac{\mathrm D\bm\psi}{\mathrm Dt}\mathrm dv
が成り立つ
$ \because\left(\forall B.\frac{\partial}{\partial t}\int_{B(t)}\xi\mathrm dv=0\right)\implies\frac{\mathrm D\xi}{\mathrm Dt}+\xi\bm\nabla\cdot\bm v=0