• 방사성폐기물학회지
• /
• 제10권4호
• /
• pp.281-294
• /
• 2012

현재까지 개발된 고준위폐기물 심지층처분장의 열-수리-역학적 복합거동 해석을 위한 전산 코드의 현황을 조사하고, 문헌에 보고된 각 코드에 의한 계산치와 현장실험 측정치의 비교 결과를 이용하여, 기존 전산 코드들의 신뢰도를 분석하였다. 개발된 전산코드들은 완충재가 없는 처분장에서는 붕괴열에 따른 암반의 열-수리-역학적 거동을 비교적 잘 모사하였으나, 포화 경암층에 위치한 완충재가 존재하는 처분장의 공학적방벽시스템 내에서 일어나는 열-수리-역학적 복합거동의 예측은 만족스럽지 못하였다. 현재 제안된 열-수리-역학적 복합거동 해석모델을 고준위폐기물 처분장 공학적방벽시스템의 거동 해석에 적용하기 위해서는 완충재 내의 수분함량 및 전 압력 분포를 보다 정교하게 모사할 수 있도록 수학적 모델의 개선이 필요하다.

• 터널과지하공간
• /
• 제32권2호
• /
• pp.144-159
• /
• 2022

본 연구에서는 지반의 열-수리-역학적 복합거동을 모델링하기 위한 순차적 접근법 기반의 시뮬레이터를 개발하고 적용된 연계해석 알고리즘의 계산성능을 분석하였다. 본 연구의 순차적 연계해석에서는 다공성 매질의 열 및 유체거동 분석을 위한 오픈소스 기반의 OpenGeoSys 수치코드와 역학해석을 위한 상용 소프트웨어 FLAC3D가 연동되었다. 해석해가 주어진 열-수리-역학적 복합거동 문제를 토대로 개발된 시뮬레이터에 대한 벤치마크 테스트가 수행되었다. 적용된 벤치마크 문제는 완전포화된 지반 내 점열원 작용 시 지반거동(시간에 따른 온도, 간극수압, 응력, 변형 변화)과 관계된다. 해석해와 수치해석 시뮬레이션 결과를 비교 분석하고 연계해석 시뮬레이터의 적정성을 조사하였다.

• 터널과지하공간
• /
• 제25권2호
• /
• pp.155-167
• /
• 2015

본 연구에서는 TOUGH2-FLAC3D 연계해석기법을 이용하여 암반공동에 고온의 열에너지를 30년간 저장하는 경우 주변 암반에 야기되는 열-수리-역학적 연계거동을 살펴보았다. 열에너지저장에 따른 암반의 거동 특성 및 환경 영향을 예측하고 이에 대한 제어기준을 수립하기 위한 기초 연구로서, 저장소 주변 암반에서 발생하는 열-수리 흐름과 역학적 거동의 상호작용에 대하여 검토하였다. 기본해석으로서 결정질 암반 내 원통형 공동에$350^{\circ}C$의 대용량 열에너지를 저장하는 경우를 모델링하였으며, 열에너지저장소의 단열성능은 고려하지 않았다. 암반 내 열전달의 주요 메카니즘은 암반의 전도에 의한 것으로 판단되며, 암반의 역학적 거동은 수리적 요소보다는 열적 요소에 지배적인 영향을 받는 것으로 나타났다. 암반과 지하수 가열에 따른 유효응력 재분포 양상과 열팽창으로 인한 암반 변위 및 지표 융기를 검토하였으며, 주변 암반에서의 전단파괴 위험도를 정량적인 수치를 통해 제시하였다. 암반 가열에 따른 열팽창으로 인하여 지표면에서 수 cm의 융기가 발생하였으며, 저장공동 상부에 인장응력이 크게 발달하면서 전단파괴의 위험도가 증가하는 것으로 나타났다.

• 터널과지하공간
• /
• 제25권2호
• /
• pp.168-185
• /
• 2015

본 연구에서는 천부의 암반 공동에 대용량 고온의 열에너지를 저장하는 경우 주변 암반에 야기되는 열-수리-역학적 연계거동을 살펴보고, 이에 지하수위와 암반 투수계수 등 수리적 조건이 미치는 영향을 검토하였다. 해석대상을 투수계수가 비교적 낮은 수준($10^{-17}m^2$)인 결정질 암반으로 가정할 때 열에너지 저장으로 인한 암반 거동에 지하수가 미치는 영향은 크지 않을 것으로 예측되었다. 저장 공동이 지하수위 하부에 위치하는 경우의 온도, 주응력, 변위 분포 등은 저장공동이 불포화대에 위치하는 경우와 거의 동일하게 나타났다. 암반내 열전달 특성은 암반의 투수계수에 매우 큰 영향을 받는 것으로 나타났다. 암반의 투수계수를 $10^{-15}m^2$ 이하로 가정한 경우 열전달은 주로 암반에 전도에 의한 것으로 판단할 수 있었으나, 투수계수를 $10^{-12}m^2$으로 가정하는 경우 지하수 대류에 의한 상향 열유동이 뚜렷이 관찰되었다. 암반 투수계수의 크기에 따라 열수의 대류나 비등으로 인한 상변화 등 복합적인 유동 특성을 나타났으며, 온도, 압력, 포화도 분포가 상이하게 발달하였다.

• 터널과지하공간
• /
• 제32권6호
• /
• pp.518-529
• /
• 2022

DECOVALEX-2023 Task C에서는 6개국 9개 참여 기관들이 스위스 Mont Terri 지하처분연구시설에 서 수행된 FE 실험을 대상으로 열-수리-역학 복합거동 모사를 위한 해석코드 개발을 수행하고 있다. 현재 현장시험 결과와 비교 분석을 위한 Step 1이 진행되고 있으며, 본 연구진은 OGS-FLAC 해석 시뮬레이터를 활용하여 일련의 해석을 진행하였다. 해석 결과 히터 가열에 따른 온도 상승이 잘 구현되었고, 상 변화에 따른 완충재 내 포화도 변화를 관측할 수 있었다. 반면 완충재 흡입력의 과대평가로 완충재 내 상대습도, 온도 변화 및 Opalinus 점토암 내 압력 변화가 현장 결과와 다소간 차이를 나타내는 것을 확인하였다. 이를 통해 완충재 흡입력이 처분시스템 해석 시 유동 해석 결과에 지배적인 영향을 미침을 확인할수 있었으며, 향후 지보재 및 초기 수압 모사 개선을 통해 향상된 결과를 도출하고자 한다. 또한, Opalinus 점토암의 열, 수리, 역학적 이방성이 잘 구현되었으며 해석 결과를 통해 OGS-FLAC 시뮬레이터의 처분시스템 해석 적용성을 확인하였다.

• 방사성폐기물학회지
• /
• 제9권1호
• /
• pp.23-32
• /
• 2011

한국형 기준 처분시스템의 공학적 방벽에서의 열-수리-역학 복합 현상을 실증하기 위한 공학적 규모 실증실험 장치인 KENTEX에서 얻은 열, 수리, 역학적 실험 데이터를 이용하여 벤토나이트의 포화공정을 해석하였다. ABAQUS를 사용한 모델계산의 함수율과 실험 결과의 비교에서 불포화 영역에서는 온도상승으로 인해 초기 수분이 감소하는 수분 재분포 공정을 모델에 포함시키지 않아 함수율의 차가 컸다. 포화 영역에서는 실험에서 초기 수분보다 낮은 함수율에서부터 지하수로 포화가 진행되지만 모델과 실험에서 얻은 함수율 값의 차이가 점점 감소해 완전포화에 도달할 때에는 두 함수율 값이 거의 비슷한 결과를 보여주였다. 포화도 약 95%에 이르는 시간은 실험결과와 계산 결과가 서로 비슷한 약 500일 정도로 예측할 수 있었다. 그리고 불포화 영역의 수분 재분포가 벤토나이트의 완전포화에 도달하는 시간에는 큰 영향을 미치지 않는 것으로 분석되었다. 따라서 본 해석기법을 사용하면 지하처분연구시설의 완충재인 벤토나이트의 포화시간을 예측할 수 있을 것으로 판단된다.