■2次元 Horizontal confined aquifer の例題■
ここでは、時間依存(非定常)の例題を紹介します。定常ですと、あまりのも簡単すぎますので。なにせLaplace equation を解くだけですから。

計算に用いたbasic paramater は、下表の通りです。計算例を示すことが目的ですから、transmissivity や storage coefficient は、実測の数値ではありません。

Hydraulic Head の初期条件は、h(x,y,0) = 0. としました。境界条件と要素分割、およびMXNTIME=100 (t=10)でのHydraulic Head を表す等ポテンシャル線は、下図の様になります。

Boundary conditions, Element discretization, and Equi-potential lines

この計算と作図には、下のプログラムを使いました。解析ソフトのIMP4Q.FOR を改良すると、様々なタイプの confined aquifer の解析に使えます。