流体の運動方程式(レベル1):ナビエ・ストークスの方程式のイントロダクション
https://www.youtube.com/watch?v=a87Yuc26q3E
内容
Navier-Stokes方程式にやられた人は水理学が苦手になり、逆にのめり込んだ人は水理学大好きになる
ε-δ論法で聞いたような話で草takker.icon
連続体の運動の記述方法
非圧縮性Newton流れにおけるNavier-Stokes方程式
非圧縮性流れかつNewton流体で成立する流体の運動方程式
粘性率が速度に依存しない粘性液体をNewton流体と読んだ。
↓を見る限りこう言ってしまっても問題なさそう
ニュートン流体 - Wikipedia
粘性とは何か - EMANの流体力学
$ \rho(\pmb{\nabla}\otimes\pmb{u})^\top\pmb{u}+\rho\frac{\partial \pmb{u}}{\partial t}=\pmb{f}-\pmb{\nabla}P+\rho\nu\pmb{\nabla}^2\pmb{u}
物質時間微分を用いた表記
$ \rho\frac{\mathrm{D} \pmb{u}}{\mathrm{D} t}=\pmb{f}-\pmb{\nabla}P+\rho\nu\pmb{\nabla}^2\pmb{u}
$ \nuを動粘性係数と呼ぶ
$ \mu:=\rho\nuを粘性係数と呼ぶ
ナビエ–ストークス方程式 - Wikipedia
この式は粘性が一定かつ非圧縮性流れの場合の式らしい
圧縮性流れにおけるNavier-Stokes方程式は↓となるようだ
$ \rho\frac{\mathrm{D} \pmb{u}}{\mathrm{D} t}=\pmb{f}-\pmb{\nabla}P+\rho\nu\pmb{\nabla}^2\pmb{u}+(\lambda+\rho\nu)\pmb{\nabla}(\pmb{\nabla}\cdot\pmb{u})
$ \lambdaは第二粘性率と呼ばれるもの
第二粘性率の意味 - EMANの流体力学
非圧縮性流れなら$ \pmb{\nabla}\cdot\pmb{u}=0で最後の項が消える
流速と圧力を求めることが目的
局所加速度と移流加速度
加速度項は慣性項とも呼ぶ
移流加速度(移流項、対流項)の必要性
定常流れだとしても、物質点は管路中の断面の広い箇所から狭い箇所に流れると加速される
この位置変化による加速が移流項として現れている
数式だと現位置を求める函数$ \pmb{\phi}を使えば一発で求まる
流速の物質表示と空間表示の微分
連続体力学 | 東北大学 材料力学研究室はこれを運動函数と呼んでいる
この呼称は便利そうだ。使わせていただこうtakker.icon
移流項の図解してみたいかもtakker.icon
流体粒子を追跡すれば解けそう
やってみた
https://kakeru.app/f20c2ae2b73b121bcfb5ef3efb189dd9 https://i.kakeru.app/f20c2ae2b73b121bcfb5ef3efb189dd9.svg
$ \pmb{x}にいる物質点$ \pmb{X}が$ \mathrm{d}t後に$ \pmb{x}+\mathrm{d}\pmb{x}に移動した場合を考える
仮に定常流れであったとしても、$ \mathrm{d}S\neq0であれば、連続の式より流速が位置によって変動する($ |\pmb{u}(\pmb{x})|\neq|\pmb{u}(\pmb{x}+\mathrm{d}\pmb{x})|)
よって定常流れ($ \frac{\partial \pmb{u}}{\partial t}=0)の場合、物質点$ \pmb{X}は$ \mathrm{d}t後に$ (\pmb{\nabla}\otimes\pmb{u})^\top\mathrm{d}\pmb{x}だけ加速される
これが移流項である
体積力$ \pmb{f}
重力や電磁気力など
よくわからない外力もとりあえずここに含めてもいいかも
面積力
圧力(垂直応力)項$ -\pmb{\nabla}Pと粘性(剪断応力)項$ \rho\nu\pmb{\nabla}^2\pmb{u}に分離して表記している
粘性項がNewtonの粘性法則の根拠となる
連続の式 (流体)
$ \pmb{\nabla}\cdot(\rho\pmb{u})+\frac{\partial \rho}{\partial t}=0
積分形:$ \int_{\partial V}\rho\pmb{u}\cdot\mathrm{d}\pmb{S}+\int_V\frac{\partial \rho}{\partial t}\mathrm{d}V=0
この動画では、非圧縮性流れにおける連続の式を連続の式と呼称している
粘性
Newtonの粘性法則の導出
Hagen-Poiseuille流れになるはず
定常流れ、$ \pmb{u}=u_x\pmb{e}_x、$ \pmb{f}=0とすると
$ \rho\frac{\partial u_x}{\partial x}u_x=-\frac{\partial P}{\partial x}+\rho\nu\pmb{\nabla}^2u_x
あとは『今日から使えるベクトル解析 (今日から使えるシリーズ)』を見てから書く
論理が逆。Newtonの粘性法則を満たす流体がNewton流体で、そこから非圧縮性Newton流れにおけるNavier-Stokes方程式の粘性項が導出される
導出は次回:流体の運動方程式(レベル2):ナビエ・ストークスの方程式の導出 part1
#2023-06-14 08:49:43
#2022-06-27 17:17:49
#2022-06-22 08:05:45