2人の大学院生による雑記ブログ
僕はあまり気にしたことがありませんが、女子にとって日焼けは天敵らしいです。 運動が好きな子は外で行うと日焼けで黒くなってしまうので、運動したいけど黒くなりたくないと板挟み状態です。 ふと思いました。肌を白くすればいいのではと。 もちろん顔は化粧という魔法がありますが、全身を白くできなければいけません。 そこで『みかんを食べると黄色くなる』という原理を使えば健康的に内側から白くなるのではないかと。 そもそもなんで『みかんを食べると黄色くなる』のか調べてみました。 みかんにはβカロチンなるオレンジ色の物質が含まれるそうです。このβカロチンは体内に摂取すると血液中に蓄積するため、皮膚にその色が現れる。（柑皮症（かんぴしょう）） という仕組みらしいです。 なんと血液.....
コードが長くなってきましたので、何をしているのか補足したいと思います。 前回の記事と合わせてご参照ください。 ①仮想セルについて（13~14,94~98行目） 有限体積法では、検査体積の中心を代表値として算出しています。つまり解析領域の端に代表点はなく、セル境界が解析領域の端になってしまいます。 そこで上の図のように一つ外側にセルを用意し境界のバランスを調整します。 左端には「初期値で一定」のディリクレ条件を与えます。隣のセルとの平均値を境界の値として考え、下記のように外側のセルの値を決めます。 $$\frac{Q_{仮想セル}+Q_1}{2}=Q_L ⇔ Q_{仮想セル}=2Q_L-Q_1 $$ 右端では「勾配0」のノイマン条件を与えます。上記と同様に考え.....
$$\frac{\partial Q}{\partial t}+ \frac{\partial F}{\partial x} =0$$ $$Q=\left\{ \begin{array}{c} \rho \\ \rho u \\ e \end{array}\right\} , F=\left\{ \begin{array}{c} \rho u \\ p+\rho u^2 \\ (e+p)u \end{array}\right\} $$ 今回は今までのまとめとして、Sod shock tubeと呼ばれる衝撃波管に対して数値解析を行います。この衝撃波管は厳密解が求められており、数値解の検証によく使われます。 pythonのコードは最後に書いておきます。 衝撃波管.....
前回は圧縮性流体の支配方程式としてオイラー方程式を取り上げました。この式は保存形で書かれており、有限体積法で解くことができます。 しかし、圧縮性流体において重要な問題があり、その一つがリーマン(Riemann)問題です。 次の図のように、ある領域において不連続な初期値をもつ場合の方程式をリーマン問題と言います。 離散化した際は当然ながら、上記のように不連続な値を取らざるを得ないです。 しかし、現実で衝撃波のように不連続な値が発生する場合、いままで通り風上差分を用いることができるのでしょうか。 オイラー方程式では3種の式が連動していますが、それぞれの伝播する方向は常に等しいわけではありません。 そこでFVS(Flux Vector Splitting)と呼ばれ.....
$$\frac{\partial Q}{\partial t}+ \frac{\partial F}{\partial x} =0$$ $$Q=\left\{ \begin{array}{c} \rho \\ \rho u \\ e \end{array}\right\} , F=\left\{ \begin{array}{c} \rho u \\ p+\rho u^2 \\ (e+p)u \end{array}\right\} $$ 第5回では圧縮性流体の代表である衝撃波管(Sod shock tube)を解きたいと思います． この Sod shock tube は厳密解が求められており，数値解との比較ができるため圧縮性流体コードの妥当性検証が行えます． .....