방사성 물질 임시 저장소의 이동 메커니즘을 정밀 시뮬레이션

초록

본 연구는 멕시코 페냐 블랑카 지역의 임시 방사성 폐기물 저장소를 대상으로, 포화 다공성 매질에서의 방사성 동위 원소 이동을 2차원 유한요소법(FEM)으로 시뮬레이션하였다. 다르시 법칙으로 유속장을 계산하고, 이를 이용해 대류‑확산‑붕괴 방정식을 풀어 238U, 234U, 230Th의 전파 경로와 농도 분포를 예측하였다. 실험 데이터와의 비교를 통해 모델을 검증했으며, 균열이 존재하는 매질에서의 비정상적인 플럼 형성을 확인하였다.

상세 분석

이 논문은 전통적인 연속체 가정에 기반한 Bear(1979) 모델을 현대적인 유한요소 해석 프레임워크와 결합한 점이 특징이다. 먼저, 포화 다공성 매질에 대해 다르시 법칙을 적용해 압력 및 중력 구배로부터 유속(u)을 도출하고, 연속 방정식에 삽입해 토양 포화도(θs)와 유체 밀도(ρf)를 고려한 흐름 해석을 수행한다. 두 번째 단계에서는 대류‑확산‑붕괴 방정식(식 2)을 풀어 방사성 물질의 농도 변화를 기술한다. 여기서 기계적 분산 텐서(Dm)는 유속에 비례하는 장축·횡축 분산 계수를 포함하고, 분자 확산(Dd)은 아이소템의 확산 계수와 결합해 전체 확산 텐서(DL)를 형성한다. 흡착‑탈착 평형은 분배계수(kd)를 통해 용액과 고체 상 사이의 농도 비를 고정시켜, 흡착된 농도(Cpi)와 용해 농도(Ci) 사이의 상호작용을 단순화한다. 방사성 붕괴는 반감기(λi)를 이용해 액체와 고체 상 모두에 적용되며, 238U→234U→230Th 순서의 체인 붕괴를 가정해 장기 거동을 평가한다. 모델 검증을 위해 Nopal I 우라늄 매장지에서 보고된 실험적 농도 프로파일과 비교했으며, 시뮬레이션 결과는 플럼의 우선 경로, 균열에 의한 비정상적 전파, 그리고 시간에 따른 농도 감소를 실측과 일치시켰다. 특히, 균열이 존재할 경우 전단 흐름이 강화되어 전통적인 균일 매질 가정으로는 설명되지 않는 고농도 영역이 형성되는 점을 강조한다. 이러한 결과는 현장 설계 시 균열 네트워크와 매질 이질성을 정량적으로 반영해야 함을 시사한다.