■FEMプログラム■
ここでは、運動方程式とHeat equation を解析するシステムとして、以下の3つのソフトを用意しました。

以下にプログラム名と役割や特徴を示します。作図に関するプログラムは、前に紹介しておりますので、詳細については、そちらを参考にしてください。

解析用プログラムについて、追加説明をしておきます。水の密度は、subroutine BODYF の中で定義してあります。定義方法は、文関数を用いています。重力加速度もsubroutine BODYF で定義してあります。単位は、cm/sec²で値は-981です。したがって、座標値の単位もcmになりますので注意して下さい。

プログラムNSTEQ8EE.FORでは、t時の温度分布を参考に、要素毎に平均温度を計算し、その温度をベースに粘性係数と熱伝導係数を計算しています。そして、t+Δt時の U, V, T にその粘性係数と熱伝導係数が使われます。係数と U, V, T の間にΔtの時間差がありますが、Δt間のU, Vの変化が小さいため、解析結果に大きな誤差を招くことはないと思います。