Wolfram言語™

エネルギー輸送

ナビエ・ストークス方程式と熱方程式を組み合せるマルチフィジックスの問題を解く．

ここで はベクトル値の速度場， は圧力， は単位行列， と はそれぞれ密度と動粘性係数である． は液体の熱容量であり は温度である．

パラメータと，切抜き部分のある矩形領域を指定する．

領域を可視化する．

プラントル(Prandtl)数とレイリー(Rayleigh)数を指定する．

ブシネスク(Boussinesq)近似を利用する熱方程式と連立した，粘性のあるナビエ・ストークス方程式を設定する．指定された材料パラメータを使う．

すべての境界壁の速度に対して，滑りなしの境界条件を設定する．

圧力の基準点を設定する．

左の壁と右の壁との温度差を指定する．

境界条件のパラメータを置き換える．

系が休止するように初期条件を設定する．

指定された間隔のメッシュを使い，粘性 と ，温度 を二次補間，圧力 を一次補間して，時間積分の進行状況と，偏微分方程式を解くのにかかる合計時間を監視する．

領域の境界を可視化する．

圧力の分布を可視化する．

温度の分布を可視化する．

速度場を可視化する．

温度の変化と速度の流線をアニメーションにする．

切抜きの位置を変えるだけで異なる速度が生成されることを調べる．