역학손상모델을 이용한 1차원 기체 주입 시험 모델링: 국제공동연구 DECOVALEX-2019 Task A Stage 1A Numerical Modelling of One Dimensional Gas Injection Experiment using Mechanical Damage Model: DECOVALEX-2019 Task A Stage 1A원문보기
고준위방사성폐기물처분장의 공학적 방벽에서는 다양한 원인으로 인해 기체가 발생한다. 만약 기체 생성 속도가 기체 확산 속도보다 빠를 경우 기체의 압력이 증가하게 되고, 기체 유입 압력(gas entry pressure)을 넘어서게 되면 기체가 급격히 벤토나이트 완충재를 통과하는 기체 이동 현상(gas migration)이 발생하게 되며 이는 사람과 주변 환경을 방사능에 노출시킬 수 있기 때문에, 공학적 방벽의 장기 건전성 확보 측면에서 기체 이동 현상을 명확히 규명하는 것이 매우 중요하다. 특히 벤토나이트 완충재와 같이 점토 물질을 다량 함유한 매질에서만 나타나는 매우 중요한 기체 흐름 현상인 팽창 흐름에 대한 수리-역학적 메커니즘을 규명하고, 기체 이동 현상의 정량적 평가를 위한 새로운 수치 해석 기법 개발 및 검증이 필수적이다. 따라서 본 연구에서는 공학적 방벽에서의 기체 이동 현상을 모사하고자 역학 손상 모델 및 손상도를 고려한 2상 유동 모델을 개발하였으며, 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험 모사를 통해 개발된 모델의 적용성을 검토하였다. 수치 해석 결과 공극 수압 및 응력, 기체 유출량이 팽창 흐름 발생 시 급격히 증가하는 현상을 모사할 수 있었다.
고준위방사성폐기물처분장의 공학적 방벽에서는 다양한 원인으로 인해 기체가 발생한다. 만약 기체 생성 속도가 기체 확산 속도보다 빠를 경우 기체의 압력이 증가하게 되고, 기체 유입 압력(gas entry pressure)을 넘어서게 되면 기체가 급격히 벤토나이트 완충재를 통과하는 기체 이동 현상(gas migration)이 발생하게 되며 이는 사람과 주변 환경을 방사능에 노출시킬 수 있기 때문에, 공학적 방벽의 장기 건전성 확보 측면에서 기체 이동 현상을 명확히 규명하는 것이 매우 중요하다. 특히 벤토나이트 완충재와 같이 점토 물질을 다량 함유한 매질에서만 나타나는 매우 중요한 기체 흐름 현상인 팽창 흐름에 대한 수리-역학적 메커니즘을 규명하고, 기체 이동 현상의 정량적 평가를 위한 새로운 수치 해석 기법 개발 및 검증이 필수적이다. 따라서 본 연구에서는 공학적 방벽에서의 기체 이동 현상을 모사하고자 역학 손상 모델 및 손상도를 고려한 2상 유동 모델을 개발하였으며, 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험 모사를 통해 개발된 모델의 적용성을 검토하였다. 수치 해석 결과 공극 수압 및 응력, 기체 유출량이 팽창 흐름 발생 시 급격히 증가하는 현상을 모사할 수 있었다.
In the engineering barriers of high-level radioactive waste disposal, gases could be generated through a number of processes. If the gas production rate exceeds the gas diffusion rate, the pressure of the gas increases and gases could migrate through the bentonite buffer. Because people and the envi...
In the engineering barriers of high-level radioactive waste disposal, gases could be generated through a number of processes. If the gas production rate exceeds the gas diffusion rate, the pressure of the gas increases and gases could migrate through the bentonite buffer. Because people and the environment can be exposed to radioactivity, it is very important to clarify gas migration in terms of long-term integrity of the engineered barrier system. In particular, it is necessary to identify the hydro-mechanical mechanism for the dilation flow, which is a very important gas flow phenomenon only in medium containing large amounts of clay materials such as bentonite buffer, and to develop and validate new numerical approach for the quantitative evaluation of the gas migration phenomenon. Therefore, in this study, we developed a two-phase flow model considering the mechanical damage model in order to simulate the gas migration in the engineered barrier system, and validated with 1D gas flow modelling through saturated bentonite under constant volume boundary conditions. As a result of numerical analysis, the rapid increase in pore water pressure, stress, and gas outflow could be simulated when the dilation flow was occurred.
In the engineering barriers of high-level radioactive waste disposal, gases could be generated through a number of processes. If the gas production rate exceeds the gas diffusion rate, the pressure of the gas increases and gases could migrate through the bentonite buffer. Because people and the environment can be exposed to radioactivity, it is very important to clarify gas migration in terms of long-term integrity of the engineered barrier system. In particular, it is necessary to identify the hydro-mechanical mechanism for the dilation flow, which is a very important gas flow phenomenon only in medium containing large amounts of clay materials such as bentonite buffer, and to develop and validate new numerical approach for the quantitative evaluation of the gas migration phenomenon. Therefore, in this study, we developed a two-phase flow model considering the mechanical damage model in order to simulate the gas migration in the engineered barrier system, and validated with 1D gas flow modelling through saturated bentonite under constant volume boundary conditions. As a result of numerical analysis, the rapid increase in pore water pressure, stress, and gas outflow could be simulated when the dilation flow was occurred.
* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.
문제 정의
The authors appreciate and thank theDECOVALEX-2019 Funding Organisations Andra, BGR/UFZ, CNSC, US DOE, ENSI, JAEA, IRSN, KAERI, NWMO, RWM, SÚRAO, SSM and Taipower for their financial and technical support of the work described in this paper.
특히 벤토나이트 완충재와 같이 점토 물질을 다량 함유한 매질에서만 나타나는 매우 중요한 기체 흐름 현상인 팽창 흐름에 대한 수리-역학적 메커니즘을 규명하고, 기체 이동 현상의 정량적 평가를 위한 새로운 수치 해석 기법 개발 및 검증이 필수적이다. 따라서 본논문에서는 공학적 방벽에서의 기체 이동 현상을 모사하고자역학 손상 모델 및 손상도를 고려한 2상 유동 모델을 개발하였으며, 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험 모사를 통해 개발된 모델의 적용성을 검토하였다.
본 연구는 DECOVALEX-2019 Task A의 일환으로 수행되었으며, 그 일부인 Stage 1A에 해당하는 1차원 기체 주입 시험 모델링에 관한 연구이다. 본 논문에서는 공학적 방벽에서의 기체 이동 현상을 모사하고자 역학 손상 모델 및 손상도를 고려한 2상 유동 모델을 개발하였으며, 일정 체적 경계 조건(constant volume boundary condition) 하에서의 1차원 기체 주입 시험 모사를 통해 개발된 모델의 적용성을 검토하였다.
벤토나이트 완충재와 같은 저투수율 매질에서의 기체 이동과 관련한 국제적인 관심사를 토대로 국제 공동 연구 프로젝트인DEvelopment of COupled models and their VALidation against EXperiments-2019 (DECOVALEX-2019) Task A가 수립되었으며, 기체 압력에 의해 매질 입자 사이가 벌어지면서 발생하는 팽창 흐름을 적절하게 모델링하는 데 필요한 수리-역학적(HM)메커니즘을 규명하고 정량적인 기체 이동 현상 모델 개발 및 검증이 수행되고 있다. 본 연구는 DECOVALEX-2019 Task A의 일환으로 수행되었으며, 그 일부인 Stage 1A에 해당하는 1차원 기체 주입 시험 모델링에 관한 연구이다. 본 논문에서는 공학적 방벽에서의 기체 이동 현상을 모사하고자 역학 손상 모델 및 손상도를 고려한 2상 유동 모델을 개발하였으며, 일정 체적 경계 조건(constant volume boundary condition) 하에서의 1차원 기체 주입 시험 모사를 통해 개발된 모델의 적용성을 검토하였다.
본 연구에서는 벤토나이트 완충재에서의 기체 이동 메커니즘 중 기체 압력에 의해 매질 입자 사이가 벌어지면서 발생하는 팽창 흐름을 모사하기 위해 수리-역학적 상호작용을 고려한 모델을 개발하였으며, 주요 지배방정식 및 구성방정식은 다음과 같다.
가설 설정
3 단계의 경우 54일까지는 실험실 시험 결과에 비해 낮은 기체 주입량을 모사하였으며, 그 이후에는 시험 결과와 유사한 주입량을 모사하였다. 3 단계 중 46 일부터54일까지 낮은 주입량으로 모델링을 수행한 이유는 실험실 시험 결과와 동일한 주입량으로 모델링을 하였을 경우 3 단계 초기에 기체 압력이 급격하게증가하는 현상이 발생하였기 때문에, 실제 시험 결과에 비해 낮은 값으로 가정 후 수치해석을 수행하였다.
, 2006), 벤토나이트 완충재와 같이 점토 물질을 다량 함유한 매질에 대한 역학 손상 모델은 아직 제안된 바가 없다. 따라서 본 연구에서는 벤토나이트 완충재의 역학 손상 거동을 탄성 손상 거동으로 가정하였으며, 매게 변수 불확실성(parameter uncertainty)을 줄이기 위해 가장 단순한 모델 중 하나인 Tang et al. (2002)이 제안한 탄성 손상 모델을 적용하였으며이는 다음과 같다.
5 MPa의 응력이 걸리는 것으로 확인되었지만, 수치 모델에서는 포화도의 변화가 발생하지 않기 때문에 추가적인 팽창압이 발생하지 않는다. 따라서 이를 보완하기 위해 초기 응력을 5.5 MPa로 가정하였다.
0 MPa의 수압 및 98 % 포화도를 갖는다. 또한 5.5 MPa의 팽윤압(swelling pressure)이 존재한다고 가정하였다. 초기 응력 측정 결과 기체 주입을 실시하기 전 벤토나이트 시료의 팽윤(swelling)으로 인해 5.
본 해석에 사용된 물성은 Table 3에 나타나 있다. 탄성 계수, 포아송비, 공극률, 고유 투수율은 Task 제안 기관인 BGS에서 제공한 물성이고, 압축 강도, 잔류 압축 강도, Biot 상수, 입자 밀도, 잔류 포화도, 비손상 투수율 계수, 최대 손상도는 참고문헌에서 얻은값이며, 그외 나머지 물성은 가정한 값이다.
제안 방법
BGS에서 수행한 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험을 모사하기 위해 총 4 단계로 나누어 해석을 진행하였다(Table 4). 1 단계에서는 벤토나이트 시료를 포화시키기 위해 기체 주입 필터(injection filter) 및 배압 필터(backpressure filter)를1 MPa의 일정 수압 조건으로 유지시켰으며, 2 단계에서는 기체 주입을 실시하기 위해 기체 주입 필터에서 3 MPa의 일정 압력 조건으로 기체를 주입하였다. 3 단계에서는 선형으로 증가하는 기체 주입량 조건을 이용해 본격적으로 기체 주입을 실시하였으며, 4 단계에서는 통로 밀봉 현상(pathway sealing)을 확인하기 위해 기체 주입을 실시하지 않은 상태를 유지하였다.
모델 적용성 검토를 위해 3차원 모델을 이용하여일정 체적 경계 조건 하에서의 1차원 기체 주입 시험 모사를 수행하였다. 4 단계로 나누어 기체 주입을 실시하였으며, 3 단계 초기의 경우 과도한 기체 주입 압력 증가 현상을 방지하기 위해 실험실 시험 결과보다낮은 주입량으로 가정 후 수치해석을 진행하였다. 수치해석 결과 기체 주입 필터에서 공극 수압 결과는 실험실 시험 결과와 유사한경향을 나타냈으며, 특히 최대 공극 수압 및 피크 이후 거동이 거의 일치하는 것으로 나타났다.
BGS에서 수행한 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험을 모사하기 위해 총 4 단계로 나누어 해석을 진행하였다(Table 4). 1 단계에서는 벤토나이트 시료를 포화시키기 위해 기체 주입 필터(injection filter) 및 배압 필터(backpressure filter)를1 MPa의 일정 수압 조건으로 유지시켰으며, 2 단계에서는 기체 주입을 실시하기 위해 기체 주입 필터에서 3 MPa의 일정 압력 조건으로 기체를 주입하였다.
먼저 TOUGH2에서는 전통적인 다공성 매질에서의 2상 유동 해석을 수행하였으며, 해석 결과인 공극 수압을 이용해 FLAC3D에서 유효응력을 변화시킨 후 역학 손상 모델을 이용하여 역학 해석을 수행하였다. FLAC3D 해석 결과인 평균 응력과 매질의 손상도를 바탕으로 수리 물성인 공극률과 투수율을 변화시켜 TOUGH2 해석을 다시 수행하였으며, 이 과정을 반복적으로 수행해 전체적인 해석을 진행하였다.
마지막으로 손상도에 따른 수리 물성 변화를 고려함으로써 매질의 변형에 의한 기체의 팽창 흐름을 모사하였다. 개발된 모델의 적용성 검토를 위해 지반공학 및 암반공학 분야의 범용 수치해석 코드인 TOUGH-FLAC을 이용하였으며(Rutqvist et al., 2002), 일부 코드 수정을 통해 본 연구에서 개발된 모델을 적용하였다(Fig. 3).
또한 기체 압력이 기준 응력 수준을 넘어설 때 매질 입자 사이가 기체 압력에 의해 벌어지면서 미세 균열의 형성 및 매질 팽창 현상이 발생하며, 이러한 미세 균열을 따라 기체가 갑작스럽게 이동하는 기체 이동 현상을 모사하기 위해 손상도 개념을 도입하였다. 기체 압력으로 인해 미세 균열의 형성 및 매질 팽창 정도를 손상도로 나타내었으며, 공극 수압의 변화로 인한 매질의 손상도(damage)를 판단하기 위해 유효응력 개념을 고려한 역학 손상 모델을 적용하였다. 마지막으로 손상도에 따른 수리 물성 변화를 고려함으로써 매질의 변형에 의한 기체의 팽창 흐름을 모사하였다.
9는 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험에 대한 기체 주입량과 유출량을 나타낸다. 기체 주입량의 경우 4 단계로 나누어 주입을 실시하였으며, 3 단계 초기의 경우 과도한 기체 주입 압력 증가 현상을 방지하기 위해 실험실 시험 결과보다 낮은 주입량으로 가정 후 수치해석을 진행하였다 (Fig. 6). 배압 필터에서의 기체 유출량이 팽창 흐름 발생 시 급격히 증가하는 현상을모사할 수 있었지만, 최대 기체 유출량 값이 실험실 시험 결과에 비해 약 1.
벤토나이트 완충재의 공극에서의 기체상 유동 및 액체상에 녹아있는 기체의 유동을 고려하기 위해 다공성 매질에서의 2상 유동을 고려하였다. 또한 기체 압력이 기준 응력 수준을 넘어설 때 매질 입자 사이가 기체 압력에 의해 벌어지면서 미세 균열의 형성 및 매질 팽창 현상이 발생하며, 이러한 미세 균열을 따라 기체가 갑작스럽게 이동하는 기체 이동 현상을 모사하기 위해 손상도 개념을 도입하였다. 기체 압력으로 인해 미세 균열의 형성 및 매질 팽창 정도를 손상도로 나타내었으며, 공극 수압의 변화로 인한 매질의 손상도(damage)를 판단하기 위해 유효응력 개념을 고려한 역학 손상 모델을 적용하였다.
모델 개발의 경우 벤토나이트 완충재의 공극에서의 기체상 유동 및 액체상에 녹아있는 기체의 유동을 고려하기 위해 다공성 매질에서의 2상 유동을 고려하였으며, 공극 수압의 변화로 인한 매질의 손상도(damage)를 판단하기 위해 유효응력 개념을 고려한 역학손상 모델을 적용하였다. 또한 손상도에 따른 수리 물성 변화를 고려함으로써 매질의 변형에 의한 기체의 팽창 흐름을 모사하였다.
기체 압력으로 인해 미세 균열의 형성 및 매질 팽창 정도를 손상도로 나타내었으며, 공극 수압의 변화로 인한 매질의 손상도(damage)를 판단하기 위해 유효응력 개념을 고려한 역학 손상 모델을 적용하였다. 마지막으로 손상도에 따른 수리 물성 변화를 고려함으로써 매질의 변형에 의한 기체의 팽창 흐름을 모사하였다. 개발된 모델의 적용성 검토를 위해 지반공학 및 암반공학 분야의 범용 수치해석 코드인 TOUGH-FLAC을 이용하였으며(Rutqvist et al.
먼저 TOUGH2에서는 전통적인 다공성 매질에서의 2상 유동 해석을 수행하였으며, 해석 결과인 공극 수압을 이용해 FLAC3D에서 유효응력을 변화시킨 후 역학 손상 모델을 이용하여 역학 해석을 수행하였다. FLAC3D 해석 결과인 평균 응력과 매질의 손상도를 바탕으로 수리 물성인 공극률과 투수율을 변화시켜 TOUGH2 해석을 다시 수행하였으며, 이 과정을 반복적으로 수행해 전체적인 해석을 진행하였다.
모델 개발의 경우 벤토나이트 완충재의 공극에서의 기체상 유동 및 액체상에 녹아있는 기체의 유동을 고려하기 위해 다공성 매질에서의 2상 유동을 고려하였으며, 공극 수압의 변화로 인한 매질의 손상도(damage)를 판단하기 위해 유효응력 개념을 고려한 역학손상 모델을 적용하였다. 또한 손상도에 따른 수리 물성 변화를 고려함으로써 매질의 변형에 의한 기체의 팽창 흐름을 모사하였다.
모델 적용성 검토를 위해 3차원 모델을 이용하여일정 체적 경계 조건 하에서의 1차원 기체 주입 시험 모사를 수행하였다. 4 단계로 나누어 기체 주입을 실시하였으며, 3 단계 초기의 경우 과도한 기체 주입 압력 증가 현상을 방지하기 위해 실험실 시험 결과보다낮은 주입량으로 가정 후 수치해석을 진행하였다.
벤토나이트 완충재의 공극에서의 기체상 유동 및 액체상에 녹아있는 기체의 유동을 고려하기 위해 다공성 매질에서의 2상 유동을 고려하였다. 또한 기체 압력이 기준 응력 수준을 넘어설 때 매질 입자 사이가 기체 압력에 의해 벌어지면서 미세 균열의 형성 및 매질 팽창 현상이 발생하며, 이러한 미세 균열을 따라 기체가 갑작스럽게 이동하는 기체 이동 현상을 모사하기 위해 손상도 개념을 도입하였다.
본 연구에서는 기체 압력 증가에 따른 미세 균열 발생 및 이에 따른 역학적 거동 변화를 모사하기 위해 손상도 개념을 적용하였다. 미세 균열 발생에 의한 매질의 손상도는 0과 1 사이의 값을 갖으며, 미세 균열이 많이 발생할수록1에 가까운 값을 갖는다.
본 연구에서는 매질의 팽창에 의한 우선 경로를 직접적으로 고려하지 않고 손상도에 따른 등가 투수율 개념으로 고려하였다. 만약 매질의 손상 혹은 파괴가 발생한다면 공극 구조의 변화를 유발하기 때문에 유체 흐름 현상이 변화하게 된다.
본 연구에서는 수분 보유 곡선을 고려하지 않았는데, 1차원 기체 주입 시험에서 벤토나이트 시료가 초기부터 시험 종료 시점까지 포화도가 98 %로 완전 포화 상태에 가깝게 유지되었기 때문에 모세관압에 의한 영향은 무시하였다.
역학 경계 조건은 일정 체적 경계 조건을 모사하기 위해 고정 경계 조건(Fixed boundary condition)을 이용하였다. 수리 경계 조건은 Fig. 5에서와 같이 스테인리스 강 튜브가 없는 면은 기체 주입을 위해 일정 유량 경계 조건(constant flux boundary condition)을 적용하였으며, 반대면 및 반경 필터(radial filter)에서는 일정 수압 경계 조건(constant pressure boundary condition)을, 반경 필터(radial filter)를 제외한 나머지 옆면은 no flow 경계 조건을 적용하였다.
시험에 사용된 벤토나이트 시료는 MX-80을이용하여 직경 60 mm 및 길이 120 mm 크기로 제작되었으며, 약 98 % 포화된 상태로 시험을 수행하였다. 약 120 일간 한쪽 면에는 헬륨 기체를 주입하고 반대면에는 일정한 수압을 걸어줌으로써 1차원 기체 주입 시험을 수행하였으며, 공극 수압을 비롯하여 축방향 응력 및 기체 유출양을 측정함으로써 기체 이동 메커니즘을 규명하고자 하였다.
DECOVALEX-2019 Task A: modElliNg Gas INjection ExpERiments (ENGINEER)는 벤토나이트 및 점성토에서의 기체이동 현상 모델링을 위해 영국 지질조사소(British Geological Survey(BGS))가 제안한 국제 공동 프로젝트이며, 벤토나이트 완충재를 포함한 공학적 방벽에서의 기체 이동 현상의 정량적 평가를 위한 새로운 수치 해석 기법 개발 및 검증이 주목적이다. 이 프로젝트에는 총 8개국의 9개 기관이 참여하고 있으며(Table 1), 1년에 2회 개최되는 워크숍을 통해 서로의 연구결과를 공유하고 의견을 교류함으로써 벤토나이트 및 점성토에서의 기체 이동 현상에 대한 수치해석 방법론 및 모델을 개발하고 있다. DECOVALEX-2019 Task A는 총 4개의 단계로 구성되어 있으며, 각 단계 특성에 따라 필수적으로 수행해야 하는 과제와 선택적으로 수행하는 과제로 구성되어 있다(Table 2).
대상 데이터
2와 같다(Daniels and Harrington, 2017). 시험에 사용된 벤토나이트 시료는 MX-80을이용하여 직경 60 mm 및 길이 120 mm 크기로 제작되었으며, 약 98 % 포화된 상태로 시험을 수행하였다. 약 120 일간 한쪽 면에는 헬륨 기체를 주입하고 반대면에는 일정한 수압을 걸어줌으로써 1차원 기체 주입 시험을 수행하였으며, 공극 수압을 비롯하여 축방향 응력 및 기체 유출양을 측정함으로써 기체 이동 메커니즘을 규명하고자 하였다.
일정 체적 경계 조건 하에서의 1차원 기체 주입 시험 모사를 위해 3차원 모델을 이용하였으며, 벤토나이트 시료의 크기는 직경 60mm, 길이 120 mm 이고, 스테인리스 강 튜브의 크기는 직경 6.35 mm, 길이 64.8 mm 이다(Fig. 5). 스테인리스 강 튜브는 벤토나이트 시료 중앙에 기체를 주입하기 위해 고안되었으나, 본 시험에서는 벤토나이트 시료의 측면에서 기체 주입을 실시하였기 때문에경계 조건의 역할만을 수행한다.
스테인리스 강 튜브는 벤토나이트 시료 중앙에 기체를 주입하기 위해 고안되었으나, 본 시험에서는 벤토나이트 시료의 측면에서 기체 주입을 실시하였기 때문에경계 조건의 역할만을 수행한다. 해석 모델의 총 요소의 개수는 7560 개 이며, 6면체 요소를 사용하였다.
이론/모형
이 때, 유체의 총 부피는 매질의 공극률과 같으며 포화도와 질량비(mass fraction)에 따라 각 상의 부피가 나뉘게 되고 각각의 밀도를 이용해 전체 유체의 질량을 계산한다. 또한 단위 면적 당 질량 유속은 다상 Darcy의 법칙 및 Fick의 법칙을 통해 계산되는데, 이를 식 (1)에 적용하면 다음과 같은 각 상에 대한 질량 보존 방정식을 얻을 수 있다.
마지막으로 공극수에 녹아있는 기체의 양을 계산하기 위해 Henry의 법칙을 적용하였으며, 이는 다음과 같다.
미세 균열 발생에 의한 매질의 손상도는 0과 1 사이의 값을 갖으며, 미세 균열이 많이 발생할수록1에 가까운 값을 갖는다. 만약 매질의 손상이 발생하였을 경우 매질의 역학적 거동 변화가 발생하는데, 이를 모사하기 위해 손상도에 따른 역학적 구성 모델(constitutive model)을 적용하였다. 암석과 같은 지반 재료 및 콘크리트에서의 손상에 따른 역학적 거동을 모사하기 위해 다양한 역학 손상 모델(mechanical damage model)이 개발되었으나(Mazars, 1986, Simo and Ju, 1987, Lee and Fenves, 1998, Jirásek,2004, Jason et al.
는 점성도(viscosity in phase β) [Pa·s]이다. 상대 투수율은 Fatt and Klikoff(1959) 식을 이용하여 계산하였으며, 이는 다음과 같다.
압축 혹은 전단파괴를 모사하기 위해, Mohr-Coulomb 파괴 기준식을 두 번째 파괴 조건으로 고려하였으며, 이때의 손상도는 다음과 같다.
역학 경계 조건은 일정 체적 경계 조건을 모사하기 위해 고정 경계 조건(Fixed boundary condition)을 이용하였다. 수리 경계 조건은 Fig.
첫 번째 메커니즘인 공극수에 용해된 기체의 이류-확산 이동은, 생성된 기체 분자가 공극수에 용해되어 액체의 흐름에 따라 이동하거나(이류) 용해된 기체 분자의 농도차로 인해 이동하는 현상(확산)을 의미한다. 이 중, 이류는 압력 구배에 의한 유체 유동을 나타내는 Darcy의 법칙을 따르며, 확산은 공극수의 기체 농도 구배에 의한 유체 유동을 나타내는 Fick의 법칙 및 공극수에 용해되는 정도를 나타내는 Henry의 법칙을 따른다.
따라서 손상도에 따른 투수율의 변화를 필수적으로 고려하여야만 한다. 이를 위해 손상도와 투수율의 관계를 설명한 다양한 연구들이 진행되어 왔지만(Gawin et al., 2001, souley et al., 2001, Tang et al., 2002), 본 연구에서는 Meschke and Grasberger(2003)이 제안한 방법론에따라 투수율을 계산하였다. 만약 특정 요소에서 매질의 손상이 발생한 경우, 이 요소의 등가 투수율은 손상되지 않은 매질의 투수율(undamaged permeability)과 손상된 매질의 투수율(damaged permeability)의 합으로 계산되며(Fall et al.
성능/효과
1차원 기체 주입 시험 결과, 기체 주입으로 인해 공극 수압이 증가하였으며 특정 기체 주입 압력(약 8.5 MPa)을 넘었을 때 기체가급격히 벤토나이트 시료로 들어가게 되는 팽창 흐름 현상이 관찰되었다. 벤토나이트 시료 부피를 일정하게 유지시키기 위한 구속경계 조건으로 인해 벤토나이트에서의 미세 균열 생성 및 전파로 인한 요소 부피 증가 현상 대신 응력이 증가하였으며, 팽창 흐름으로 인해 배압 필터에서의 기체 유출량의 급격히 발생하는 현상이 관찰되었다.
따라서 반경 필터 주변의 요소들이 인장 파괴가 먼저 발생하게 되고, 인장 파괴에 의해 투수율이 급격히 증가함으로써 반경 필터 방향으로의 유체 흐름이 강화되는 현상이 나타났다. 결과적으로 스테인리스 강 튜브 주변부에서의 공극 수압은 파괴를 일으킬 만큼 증가하지 않아, 스테인리스 강 튜브 주변부는 상대적으로 손상이 거의 발생하지 않은 것으로 판단된다.
7은 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험에 대한 기체 주입 필터, 배압 필터, 그리고 반경 필터에서의 공극 수압 변화를 나타낸그림이다. 기체 주입 필터에서 공극 수압 결과는 실험실 시험 결과와 유사한 경향을 나타냈으며, 특히 최대 공극수압 및 피크 이후(post-peak) 거동이 거의 일치하는 것으로 나타났다. 하지만 최대 공극 수압 이전의 거동의 경우, 수치 모델링은초기에 급격한 압력 증가 현상을 보이다가 점차 기울기가 감소하였는데 실내 시험 결과는 기울기가 점차 증가하는 형태로 나타났다.
수치 해석 결과 반경 필터에서의 갑작스러운 공극 수압 증가 현상을 확인할 수 있었으며, 최대 공극 수압 및피크 이후 거동도 실험실 시험 결과와 거의 일치하는 것으로 나타났다. 다만 실험실 시험 결과는 필터의 위치에 상관없이 거의 동시에 압력이 급격히 증가하였으나, 수치 해석 결과는 필터의 위치에 따라 압력 증가 현상 발생 시점이 다르게 나타났다. 수치 해석의경우 손상도가 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하기 때문에 반경 필터의 위치에 따라 압력 증가 발생 시점이다르게 나타났으며, 이는 본 연구에서 개발된 모델의 한계점을 나타낸다.
반경 필터는 시험 기간 동안 1 MPa의 압력으로 유지되기 때문에 완전한 1차원 유체 흐름이 아닌 반경 필터 방향으로의 유체 흐름이 나타났으며, 이에 따라 스테인리스 강 튜브 주변부 보다 반경 필터에 가까운 요소들의 유체 압력이 높게 나타났다. 따라서 반경 필터 주변의 요소들이 인장 파괴가 먼저 발생하게 되고, 인장 파괴에 의해 투수율이 급격히 증가함으로써 반경 필터 방향으로의 유체 흐름이 강화되는 현상이 나타났다. 결과적으로 스테인리스 강 튜브 주변부에서의 공극 수압은 파괴를 일으킬 만큼 증가하지 않아, 스테인리스 강 튜브 주변부는 상대적으로 손상이 거의 발생하지 않은 것으로 판단된다.
수치해석 결과 기체 주입 필터에서 공극 수압 결과는 실험실 시험 결과와 유사한경향을 나타냈으며, 특히 최대 공극 수압 및 피크 이후 거동이 거의 일치하는 것으로 나타났다. 또한 수치 해석 결과 응력 증가 현상및 최대 응력 이후 감소하는 경향은 모사할 수 있었으며, 배압 필터에서의 기체 유출량이 팽창 흐름 발생 시 급격히 증가하는 현상을 모사할 수 있었다. 하지만 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하는 손상도로 인해 공극 수압 증가 시점에서 차이가 발생하였으며, 응력 증가 기울기나 최대 응력 이후 응력 감소 기울기가 실험실 결과에 비해 다소 완만하게 나타나는 등 현재 개발된 해석 모델이 수리-역학적 상호작용의 영향을 과소평가 하는 한계를 확인할 수 있었다.
6). 배압 필터에서의 기체 유출량이 팽창 흐름 발생 시 급격히 증가하는 현상을모사할 수 있었지만, 최대 기체 유출량 값이 실험실 시험 결과에 비해 약 1.0E-8 m3/s 정도 낮게 나타남을 확인할 수 있었다. 이는 실험실 시험의 경우 특정 통로를 통해서만 기체 흐름이 발생하지만 수치 해석의 경우 손상도가 벤토나이트 시료 전 영역에 걸쳐 발생하기 때문에, 수치 해석에서 벤토나이트 시료가 저장하고 있는 기체의 양이 실험실 시험에 비해 많기 때문인 것으로 판단된다.
5 MPa)을 넘었을 때 기체가급격히 벤토나이트 시료로 들어가게 되는 팽창 흐름 현상이 관찰되었다. 벤토나이트 시료 부피를 일정하게 유지시키기 위한 구속경계 조건으로 인해 벤토나이트에서의 미세 균열 생성 및 전파로 인한 요소 부피 증가 현상 대신 응력이 증가하였으며, 팽창 흐름으로 인해 배압 필터에서의 기체 유출량의 급격히 발생하는 현상이 관찰되었다.
반경 필터의 경우 특정 시점에서 공극 수압이 갑작스럽게 증가하는 현상을 확인할 수 있는데, 이는 기체 압력에 따른 팽창 흐름이발생하였기 때문이다. 수치 해석 결과 반경 필터에서의 갑작스러운 공극 수압 증가 현상을 확인할 수 있었으며, 최대 공극 수압 및피크 이후 거동도 실험실 시험 결과와 거의 일치하는 것으로 나타났다. 다만 실험실 시험 결과는 필터의 위치에 상관없이 거의 동시에 압력이 급격히 증가하였으나, 수치 해석 결과는 필터의 위치에 따라 압력 증가 현상 발생 시점이 다르게 나타났다.
팽창 흐름 발생 이후 응력 또한 급격히 증가하는 현상이 나타났으며, 최대 응력 이후 점차 응력이 감소하였다. 수치 해석 결과 응력 증가 현상 및 최대 응력 이후 감소하는 경향은 모사할 수 있었으나, 최대 응력 값에서 차이가 발생하였고 응력 증가 기울기나 최대 응력 이후 응력 감소 기울기가 실험실 결과에 비해 다소 완만하게나타났다. 최대 응력 값의 경우 기체 주입 필터의 로드셀에서는 수치 해석 결과가 실험실 시험 결과에 비해 약 1 MPa 정도 높게 나타났으며, 배압 필터의 로드셀에서는 약 2 MPa 정도 낮게, 반경 로드셀의 경우 약 3 MPa 정도 높은 결과를 나타냈다.
최대 응력 값의 경우 기체 주입 필터의 로드셀에서는 수치 해석 결과가 실험실 시험 결과에 비해 약 1 MPa 정도 높게 나타났으며, 배압 필터의 로드셀에서는 약 2 MPa 정도 낮게, 반경 로드셀의 경우 약 3 MPa 정도 높은 결과를 나타냈다. 수치 해석의경우 손상도가 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하기 때문에 기체 주입 필터에서 응력이 가장 높게 나타났으며, 배압 필터쪽으로 갈수록 응력이 낮게 나타남을 확인할 수 있었다. 수치 해석 결과 응력 증가 및 응력 감소 기울기가 실험실 시험결과에 비해 완만하게 나타났는데, 이는 본 연구에서 개발된 해석 모델에서 수리-역학적 상호작용의 영향이 과소평가됨을 의미하기에, 모델 개선에 관한 추가적인 연구가 필요하다.
4 단계로 나누어 기체 주입을 실시하였으며, 3 단계 초기의 경우 과도한 기체 주입 압력 증가 현상을 방지하기 위해 실험실 시험 결과보다낮은 주입량으로 가정 후 수치해석을 진행하였다. 수치해석 결과 기체 주입 필터에서 공극 수압 결과는 실험실 시험 결과와 유사한경향을 나타냈으며, 특히 최대 공극 수압 및 피크 이후 거동이 거의 일치하는 것으로 나타났다. 또한 수치 해석 결과 응력 증가 현상및 최대 응력 이후 감소하는 경향은 모사할 수 있었으며, 배압 필터에서의 기체 유출량이 팽창 흐름 발생 시 급격히 증가하는 현상을 모사할 수 있었다.
8은 일정 체적 경계 조건 하에서의 1차원 기체 주입 시험에 대한 기체 주입 필터, 배압 필터, 그리고 반경 필터에서의 응력 변화를 나타낸 그림이다. 실험실 시험의 경우 시간이 지남에 따라 팽윤압이 계속 발생하면서 팽창 흐름이 발생하기 이전 최대 약 10MPa에 달하는 응력이 측정되었으나, 앞서 설명한 대로 수치 해석에서는 포화도의 변화가 발생하지 않기에 추가적인 팽윤압이 발생하지 않아 초기 조건으로 가정한 5.5 MPa의 응력을 유지하였다. 팽창 흐름 발생 이후 응력 또한 급격히 증가하는 현상이 나타났으며, 최대 응력 이후 점차 응력이 감소하였다.
수치 해석 결과 응력 증가 현상 및 최대 응력 이후 감소하는 경향은 모사할 수 있었으나, 최대 응력 값에서 차이가 발생하였고 응력 증가 기울기나 최대 응력 이후 응력 감소 기울기가 실험실 결과에 비해 다소 완만하게나타났다. 최대 응력 값의 경우 기체 주입 필터의 로드셀에서는 수치 해석 결과가 실험실 시험 결과에 비해 약 1 MPa 정도 높게 나타났으며, 배압 필터의 로드셀에서는 약 2 MPa 정도 낮게, 반경 로드셀의 경우 약 3 MPa 정도 높은 결과를 나타냈다. 수치 해석의경우 손상도가 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하기 때문에 기체 주입 필터에서 응력이 가장 높게 나타났으며, 배압 필터쪽으로 갈수록 응력이 낮게 나타남을 확인할 수 있었다.
또한 수치 해석 결과 응력 증가 현상및 최대 응력 이후 감소하는 경향은 모사할 수 있었으며, 배압 필터에서의 기체 유출량이 팽창 흐름 발생 시 급격히 증가하는 현상을 모사할 수 있었다. 하지만 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하는 손상도로 인해 공극 수압 증가 시점에서 차이가 발생하였으며, 응력 증가 기울기나 최대 응력 이후 응력 감소 기울기가 실험실 결과에 비해 다소 완만하게 나타나는 등 현재 개발된 해석 모델이 수리-역학적 상호작용의 영향을 과소평가 하는 한계를 확인할 수 있었다. 추후 DECOVALEX-2019 Task A 워크숍을 통해 국외 연구팀들과 서로의 연구결과를 공유하고 의견을 교류함으로써 본 연구에서 개발된 수치 해석 모델을 지속적으로 개선해 나가는 한편, 추가적인 해석을 통해 모델 적용성을 검토할 예정이다.
기체 주입 필터에서 공극 수압 결과는 실험실 시험 결과와 유사한 경향을 나타냈으며, 특히 최대 공극수압 및 피크 이후(post-peak) 거동이 거의 일치하는 것으로 나타났다. 하지만 최대 공극 수압 이전의 거동의 경우, 수치 모델링은초기에 급격한 압력 증가 현상을 보이다가 점차 기울기가 감소하였는데 실내 시험 결과는 기울기가 점차 증가하는 형태로 나타났다. 수치 모델링의 경우 벤토나이트 시료 옆면에 직접적으로 기체를 주입하는 반면, 실험실 시험의 경우 기체 주입 펌프를 이용해 기체 주입 필터를 거쳐 기체 주입이 이루어지기 때문에 압력 증가 반응에 차이가 발생한 것으로 판단된다.
후속연구
수치 해석의경우 손상도가 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하기 때문에 기체 주입 필터에서 응력이 가장 높게 나타났으며, 배압 필터쪽으로 갈수록 응력이 낮게 나타남을 확인할 수 있었다. 수치 해석 결과 응력 증가 및 응력 감소 기울기가 실험실 시험결과에 비해 완만하게 나타났는데, 이는 본 연구에서 개발된 해석 모델에서 수리-역학적 상호작용의 영향이 과소평가됨을 의미하기에, 모델 개선에 관한 추가적인 연구가 필요하다.
다만 실험실 시험 결과는 필터의 위치에 상관없이 거의 동시에 압력이 급격히 증가하였으나, 수치 해석 결과는 필터의 위치에 따라 압력 증가 현상 발생 시점이 다르게 나타났다. 수치 해석의경우 손상도가 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하기 때문에 반경 필터의 위치에 따라 압력 증가 발생 시점이다르게 나타났으며, 이는 본 연구에서 개발된 모델의 한계점을 나타낸다.
하지만 기체 주입 필터에서부터 배압 필터로 점진적으로 발생하는 손상도로 인해 공극 수압 증가 시점에서 차이가 발생하였으며, 응력 증가 기울기나 최대 응력 이후 응력 감소 기울기가 실험실 결과에 비해 다소 완만하게 나타나는 등 현재 개발된 해석 모델이 수리-역학적 상호작용의 영향을 과소평가 하는 한계를 확인할 수 있었다. 추후 DECOVALEX-2019 Task A 워크숍을 통해 국외 연구팀들과 서로의 연구결과를 공유하고 의견을 교류함으로써 본 연구에서 개발된 수치 해석 모델을 지속적으로 개선해 나가는 한편, 추가적인 해석을 통해 모델 적용성을 검토할 예정이다.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.