Time Derivative and Upwind Method
1 Dim. Convection-Diffusion Equation-3
■例題: No flux■
参考になりそうな、例題を計算してみましょう。この例題の諸条件は、次の様になっています。
| 上流の境界からの物質のfluxは、ゼロとする。 |
| Sediment の密度は、全てρ(x)=1.とする。 |
| V(Z)=0.2とする。 |
| 上流の境界から少し領域に入ったところに(z=x(5))、質量=1000を注入する。 |
| 領域の長さTL=10を100要素に分割する。 |
| Time step は、DT=0.05とする。そして、TMAX = 34.とする。 |
| 拡散係数は、D = 0.001とする。 |
| 計算結果のファイルへの書き出しは、40 x DT 毎とする。 |
上の事柄は、データを生成するプログラムSET1D1.FORにまとめてあります。このプログラムを実行すると、ONEDEE2.FOR用の入力データONEDEE.DATが作成されます。そして、ONEDEE2.FORを実行すると、SOLUTION.FEMが出力されます。さらに、この出力結果をベースに、アニメーションGIFを作成することが出来ます。ここをクリック
してみて下さい。アニメーションが表示されます。
してみて下さい。アニメーションが表示されます。