[논문 리뷰] Modelling and Numerical Simulation of Gas Migration in a Nuclear Waste Repository

Alain Bourgeat, Mladen Jurak외 1명|arXiv (Cornell University)|2010. 06. 15.

Hydraulic Fracturing and Reservoir Analysis참고 문헌 31인용 수 5

한 줄 요약

이 논문은 깊은 지하 핵 폐기물 저장소에서 수소 이동을 위한 새로운 비열화 없는 복합 두상 유동 모델을 제시한다. 주로 총 수소 질량 농도와 액체 압력을 주요 변수로 사용한다. 이 구성은 포화 및 비포화 조건에서 상의 출현과 소멸을 성공적으로 시뮬레이션하며, 수치 시험을 통해 날카로운 포화 전면과 동적 평형 전이를 다룰 때의 강건성을 확인한다.

ABSTRACT

We present a compositional compressible two-phase, liquid and gas, flow model for numerical simulations of hydrogen migration in deep geological radioactive waste repository. This model includes capillary effects and the gas diffusivity. The choice of the main variables in this model, Total or Dissolved Hydrogen Mass Concentration and Liquid Pressure, leads to a unique and consistent formulation of the gas phase appearance and disappearance. After introducing this model, we show computational evidences of its adequacy to simulate gas phase appearance and disappearance in different situations typical of underground radioactive waste repository.

연구 동기 및 목표

• 핵 폐기물 저장소 시뮬레이션 중 다공성 매체에서 기체 상의 출현과 소멸을 모델링하는 데 도전하는 문제를 해결한다.
• 상 포화도 방정식에서 발생하는 수치적 열화 문제를 해결하기 위해 비열화 변수를 사용해 시스템을 재구성한다.
• 포화 및 비포화 유동 영역을 아우르는 수소 이동을 일관되고 안정적으로 시뮬레이션한다.
• 모델이 캡일러리 효과, 기체 확산, 압력 불연속성과 같은 복잡한 물리 현상을 다룰 수 있는지를 검증한다.
• 모델이 지하 고형 방사성 폐기물 저장소에서 흔히 발생하는 실제적이고 도전적인 시나리오에 적용 가능함을 보여준다.

제안 방법

• 액체 상과 기체 상을 포함한 복합적이고 압축 가능한 두상 유동 모델을 수립하며, 구성 성분으로는 물과 수소(H₂)를 중심으로 한다.
• 포화도 방정식에서의 열화를 피하기 위해 총 수소 질량 농도와 액체 압력을 주요 변수로 사용한다.
• 캡일러리 압력 관계와 기체 법칙을 기반으로 기체 상을 모델링하고, 액체 상에서 용해된 수소의 확산을 포함한다.
• 다르시-마스크트 법칙을 사용해 상 유량을 기술하며, 상대 투과도 함수는 상의 포화도에 따라 결정된다.
• 포화 전면에서의 불연속성과 급격한 기울기를 다룰 수 있도록 제어 방정식의 약한 형태를 적용한다.
• 적응형 시간 스텝과 공간 이산화를 사용한 수치 시뮬레이션을 통해 포화 및 비포화 상태 간의 동적 전이를 해상한다.

실험 결과

연구 질문

• RQ1 기체 상의 출현과 소멸이 발생하는 동안 수치적으로 안정적인 두상 유동 모델을 어떻게 구성할 수 있는가?
• RQ2 포화 및 비포화 조건에서 포화도 기반 방정식의 일관성과 열화를 방지하기 위해 어떤 변수 선택이 필요한가?
• RQ3 캡일러리 효과와 용해된 수소의 확산은 기체 상의 포화도와 압력 프로파일의 진화에 어떻게 영향을 미치는가?
• RQ4 이동하는 포화 전면을 동반한 포화 상태에서 비포화 상태로의 전이를 모델이 정확하게 시뮬레이션할 수 있는가?
• RQ5 대류 및 확산 유량은 포화 전면 근처에서 수소 농도 기울기를 어떻게 형성하는가?

주요 결과

• 모델은 1.92 × 10⁵ s 동안 x = 0.5 m 에서 출발한 기체 상 포화 전면이 왼쪽 경계(Γ_in)로 이동하는 것을 성공적으로 시뮬레이션하였으며, 전면 높이는 약 0.16에서 약 0.03으로 감소하였다.
• 포화 영역에서는 압력이 급격히 증가하여 t = 10³ s 에 약 1 MPa 에서 약 1.6 MPa 로 상승한 후 서서히 감소하는 경향을 보였으며, 이는 상부 기체 상에 의한 동적 압축을 나타낸다.
• 포화 전면에서 용해된 수소 농도에 강한 기울기가 발생하였으며, 이는 좌측에서 우측으로의 대류 운반과 반대 방향의 확산 유량 간의 경쟁으로 인해 t = 1.92 × 10⁵ s 까지 지속되었다.
• t = 1.92 × 10⁵ s 이후 전체 도메인이 비포화 상태(S_g > 0)가 되었고, 모든 변수(압력, 수소 밀도, 포화도)가 공간적으로 균일한 정적 상태로 단조롭게 수렴하였다.
• 계는 압력과 농도 프로파일이 연속적이며, 고립된 도메인 내 질량 보존 법칙과 일치하는 평형 상태에 도달하였다.
• 모델은 압력 연속성의 즉각적인 복귀와 농도 기울기의 점진적 평탄화를 정확히 포착하였으며, 상 전이 역학을 다룰 때의 강건성을 검증하였다.