HydroGeoSphere/多孔介質

預設多孔介質傳輸值

預設情況下，域中所有多孔介質區域（和單元）將被分配預設的多孔介質傳輸屬性，這些屬性列在表 5.21中。這裡包括用於修改多孔介質的引數，使其充當模擬傳輸的雙重孔隙介質，如第 2.6.1.3 節所述，以及用於同位素分餾的引數，如第 2.6.1.4 節所述。

表 5.21：多孔介質傳輸屬性的預設值

引數	值	單位
縱向彌散度 ${\alpha }_{l}$	1.0	m
橫向彌散度的水平分量 ${\alpha }_{t}$	0.1	m
橫向彌散度的垂直分量 ${\alpha }_{t}$	0.1	m
堆積密度 ${\rho }$	2031.25	kg m−3
曲折度 ${\tau }$	0.1	-
非流動區孔隙度 ${\theta }_{Imm}$	0.0	-
非流動區質量傳遞係數 ${\alpha }$	0.0	s−1
逆向分餾速率 $k_{r}$	0.0	s−1
分餾因子 ${\alpha }_{r}$	0.0	-
固體相與水相的質量比 $x_{r}$	0.0	-
固體熱導率 $k_{s}$	3.0	W m−1 K−1
固體比熱容 $c_{s}$	738.0	J kg−1 K−1
固體密度 ${\rho }_{s}$	2650.0	kg m−3

以下說明可應用於多孔介質，如第 5.8.1 節所述，以修改預設的傳輸引數。對於每條說明，我們將指出其範圍（即.grok, .mprops）。請記住，如果在字首中使用說明.grok檔案，它將影響當前選擇的區域集，而在屬性（例如.mprops）檔案中，它將僅影響其中包含的命名材料。

縱向彌散度

範圍：.grok .mprops

1. val 縱向彌散度 [L]，${\alpha }_{l}$ 在公式 2.97 中。

• • •

橫向彌散度

範圍：.grok .mprops

1. val 橫向彌散度的水平分量 [L]，${\alpha }_{t}$ 在公式 2.97 中。

• • •

垂直橫向彌散度

範圍：.grok .mprops

1. val 橫向彌散度垂直分量 [L]，${\alpha }_{t}$ 在方程式 2.97 中。

• • •

曲折度

範圍：.grok .mprops

1. val 曲折度，${\tau }$ 在方程式 2.97 中。

• • •

各向異性曲折度比

範圍：.grok .mprops

1. y_tortratio y 方向上的曲折度比。預設值為 1。
2. z_tortratio z 方向上的曲折度比。預設值為 1。

預設情況下，曲折度是各向同性的，因為 y 和 z 方向的比率值都設定為 1。可以透過輸入大於 0 小於 1 的值使曲折度各向異性。這些值將用於分別乘以 y 和 z 方向的曲折度，${\tau }$，以獲得方向值。

• • •

堆積密度

範圍：.grok .mprops

1. val 堆積密度 [M L⁻³]，${\rho }_{b}$ 在方程式 2.96 中。如果使用此指令，則之前為該材料儲存的固體密度值將被從堆積密度、水密度和孔隙度計算出的固體密度覆蓋，$[{\rho }_{s}=({\rho }_{b}-{\theta }_{s}{\rho })/(1-{\theta }_{s})]$。

• • •

預設情況下，多孔介質充當單孔隙介質（即固定區不活躍），因為孔隙度和傳質係數設定為零。為了啟用雙孔隙特性，可以使用以下兩條指令輸入非零值。

固定區孔隙度

範圍：.grok .mprops

1. val 固定區孔隙度，${\theta }_{Imm}$ 在方程式 2.101、2.135 和 2.136 中。

• • •

固定區傳質係數

範圍：.grok .mprops

1. val 固定區傳質係數 [T⁻¹]，${\alpha }_{Imm}$ 在方程式 2.134 中。

• • •

同位素分餾資料...結束

範圍: .mprops

導致 grok 開始讀取一組同位素分餾指令，直到它遇到結束指令。

如果未發出進一步的指令，則將使用 表 5.21 中列出的默認同位素分餾引數值。

• • •

以下三個指令可用於更改這些值

反向速率

範圍: .mprops

1. val 反向分餾速率 [L⁻¹]，$k_{r}$ 在方程式 2.143 中。

• • •

分餾因子

範圍: .mprops

1. val 分餾因子，${\alpha }_{r}$ 在方程式 2.143 中。

• • •

岩石-水質量比

範圍: .mprops

1. val 同位素岩石-水質量比，$x_{r}$，見公式 2.102。

• • •

接下來的四條指令可用於更改多孔介質的熱學性質

固體熱導率

範圍：.grok .mprops

1. val 固體的溫度不變熱導率 [W L⁻¹ K⁻¹]。系統內部根據固相和液相的體積分數計算總體熱導率。

• • •

固體的溫度依賴熱導率

範圍：.grok .mprops

1. k_s1 溫度為 $t_{s1}$ 時的熱導率 [W L⁻¹ K⁻¹]。
2. t_s1 熱導率等於 $k_{s1}$ 時的溫度 [°C]。

如果指定了該指令，則固相的熱導率與溫度相關。總體熱導率也與溫度相關，並根據固相和液相的體積分數在系統內部計算得出。這裡假設固體的熱導率以每升高 10°C 降低 1% 的恆定速率下降，熱導率和溫度之間的關係由 $(k_{s1},t_{s1})$ 這對值定義。

• • •

固體的比熱容

範圍：.grok .mprops

1. val 固相的比熱容 [J kg⁻¹ K⁻¹]。預設值為 730.0 J kg⁻¹ K⁻¹。

• • •

請注意，現在從總體密度和孔隙率自動計算多孔介質固相的密度。