JAEA-Research-2010-034.pdf:16.23MB - JAEAの研究開発成果 ...
JAEA-Research-2010-034.pdf:16.23MB - JAEAの研究開発成果 ...
JAEA-Research-2010-034.pdf:16.23MB - JAEAの研究開発成果 ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
�<br />
�t<br />
� � � �<br />
� � �cG<br />
n �<br />
ac � �aDGij<br />
G n<br />
�x<br />
�<br />
i �<br />
�x<br />
j<br />
�<br />
<strong>JAEA</strong>-<strong>Research</strong> <strong>2010</strong>-034<br />
ここに、a は気相率、cGは気相化学種 n の濃度、DG は気相化学種の拡散係数である。<br />
- 31 -<br />
(3.5-12)<br />
また THMC モデルでは、定置後の緩衝材が閉鎖環境に置かれる場合のことを想定して、不飽和<br />
から飽和に至る過程で要素内の飽和度の値が 0.9999 を超えた場合には飽和状態であると判断し、<br />
要素内の気相化学種の質量を強制的に液相中に溶解させ、気相化学種を解析条件から排除するモ<br />
デルとしている。<br />
(2)拡散項の取扱い<br />
THMC モデルでは既存の独立した解析コードを連成させながら質量保存を満足させるために、<br />
地下水中の物質移行問題に関する移流項の質量粒子法と同様の方法を採用しており、分散項と同<br />
様に有限体積法によって離散化される(日本原子力研究開発機構、2009)。<br />
解析コード間で受け渡す変数は次式によって濃度から質量へ変換することで時間ステップ間で<br />
の質量を保持される。<br />
t<br />
gi<br />
t<br />
gi<br />
t<br />
i<br />
t<br />
i<br />
M � C � a �Vol<br />
(3.5-13)<br />
ここに、a は気相率、�は体積含水率、Cg は気相化学種の濃度、Mg は気相化学種の質量、Vol<br />
は要素体積である。<br />
式(3.5-13)の質量は次の時間ステップにおける要素内水分量に応じた気相の体積変化を考慮し<br />
て、次式で濃度へ変換される。<br />
M<br />
t�1<br />
gi<br />
C gi �<br />
(3.5-14)<br />
a �Vol<br />
t�1<br />
gi<br />
t<br />
t�1<br />
i<br />
3.5.3 安定基準について<br />
数値的安定を確保するための基準として、移流に対するクーラン基準もしくは分散に対するフ<br />
ォン・ノイマン基準を満足するための時間刻みが自動設定される。クーラン基準は双曲形の移流<br />
方程式に対して陽解法による差分法を用いる場合の安定基準であり、Δt の間に許容される移流距<br />
離が格子 1 個分未満であることを意味し、フォン・ノイマン基準は放物形の非定常拡散方程式に<br />
対して陽解法による差分法を用いる場合の安定基準である(越塚、1997)。それぞれの時間刻み Δt<br />
に対する条件は次式で表わされる。<br />
① クーラン基準<br />
C x � C y � C z �1.<br />
0<br />
(3.5-15)<br />
�t<br />
�t<br />
�t<br />
C x � Vx<br />
、 C y � V y 、 C z � Vz<br />
(3.5-16)<br />
�x<br />
�y<br />
�x<br />
ここに、C は移流方程式の差分計算上現れる無次元数であるクーラン数、Δt は時間刻み、Δx、<br />
Δy、Δz は格子幅、V は実流速を表す。