[論文レビュー] Model and Solution of Darcy's Law: Homogeneous and Inhomogeneous Media

V. Vázquez-Báez, Ana Beatriz Rubio-Arellano+2 more|arXiv (Cornell University)|Feb 3, 2018

Groundwater flow and contamination studies被引用数 1

ひとこと要約

本論文は、定常状態下での2次元の水中層における地下水の流れをシミュレートする計算モデルを、ダルシーの法則に基づいて開発した。均質的等方的および非等方的媒体を比較し、ラプラス方程式および非等方的流れの一般化された形式を解き、トゥースの解析的解と照合して高い適合精度を確認した。

ABSTRACT

Using a mathematical analytic scheme we present a model of groundwater dynamics, given as usual in terms of partial differential equations whose solution requires the implementation of boundary conditions under physical basis. The final aim is to model the water flow, thus allowing us to fix a suitable domain to predict sustainable explotaiton of a regional aquifer. Under conditions of stationary flow the most fitting model to describe the properties of water motion through porous media is Darcy's Law, this is done for two cases of a 2D aquifer with water table of constant slope assuming first we have an homogeneous isotropic media and after within a more realistic anisotropic soil. Taking into account some geophysical parameters we develop a computational routine that solves the resulting Laplace equation and the generalized equation corresponding to flow in the anisotropic medium, both in the stationary flow regime. Under certain conditions we calibrate the numerical model comparing with the analytic solution proposed by Toth.

研究の動機と目的

• 定常な流れ条件下における2次元の地域的水中層における地下水の流れのダイナミクスをダルシーの法則を用いてモデル化すること。
• 均質的等方的媒体とより現実的な非等方的多孔質媒体との間での流れの挙動を比較すること。
• 非等方的流れのための一般化されたラプラス方程式およびその一般化形式を解く計算ルーチンを開発すること。
• 精度の検証のため、数値モデルをトゥースによる解析的解と照合すること。
• 正確な流れ予測を通じて、地域的水中層の持続可能な利用を支援すること。

提案手法

• 定常状態の地下水の流れとしてダルシーの法則を偏微分方程式（PDE）として定式化する。
• 物理的制約に基づく境界条件を適用して水中層の定義域を設定する。
• 数値的手法を用いて均質的等方的媒体におけるラプラス方程式を解く。
• 一般化されたダルシーの法則方程式を解くことで、非等方的媒体へのモデルの拡張を実施する。
• 均質的および非均質的媒体における両方の流れをシミュレートする計算ルーチンを実装する。
• 検証のため、数値解をトゥースの解析的解と照合する。

実験結果

リサーチクエスチョン

• RQ1 均質的等方的媒体における2次元水中層の定常状態の地下水の流れを、ダルシーの法則に基づく数値モデルがどの程度正確に再現できるか。
• RQ2 2次元水中層において非等方的多孔質媒体と等方的媒体をモデル化した場合、流れのパターンにどのような差が生じるか。
• RQ3 一般化されたダルシーの法則方程式の数値解が、トゥースによる解析的解とどの程度一致するか。
• RQ4 地球物理学的パラメータの組み込みが、不均一な水中層における流れ予測の精度にどのように影響するか。
• RQ5 適合された数値モデルは、持続可能な地域的水中層管理を信頼性を持って支援できるか。

主な発見

• 数値モデルはトゥースが提唱した解析的解を正確に再現しており、高い適合適合性が確認された。
• 非等方的媒体における流れのパターンは、等方的仮定とは顕著に異なる傾向を示し、特に流れの方向および速度分布に顕著な差が生じた。
• 計算ルーチンは、2次元領域における非等方的条件下の一般化されたダルシーの法則方程式を効果的に解くことに成功した。
• さまざまな地球物理学的パラメータの入力に対して、モデルは一貫した収束性と安定性を示した。
• 均質的媒体におけるラプラス方程式の解は、理論的期待と正確に一致した。
• 本研究は、地域的水中層における現実的な地下水の流れの予測には、非等方性のモデル化が不可欠であることを確認した。