Geologic disposal at deep depth is an acceptable way to dispose of high-level radioactive waste and isolate it from the biosphere. The geological repository system comprises an engineered barrier system (EBS) and the host rock. The system aims to delay radionuclide migration through groundwater flow, and also, the flow affects the saturation of the bentonite in the EBS. The thermal conductivity of bentonite is a function of saturation, so the temperature in the EBS is directly related to the flow system. High-temperature results in the two-phase flow, and the two-phase flow system also affects the flow system. Therefore, comprehending the influencing parameters on the flow system is critical to ensure the safety of the disposal system. Various studies have been performed to figure out the complex two-phase flow characteristics, and numerical simulation is considered an effective way to predict the coupled behavior. DECOVALEX (DEvelopment of COupled models and their VALidation against EXperiments) is one of the most famous international cooperating projects to develop numerical methods for thermo-hydro-mechanicalchemical interaction, and Task C in the DECOVALEX-2023 has the purpose of simulating the Fullscale Emplacement (FE) experiment at the Mont-Terri underground research laboratory. We used OGS-FLAC, a self-developed numerical simulator combining OpenGeoSys and FLAC3D, for the simulation and targeted to analyze the effecting parameters on the two-phase flow system. We focused on the parameters of bentonite, a key component of the disposal system, and analyzed the effect of compressibility and air entry pressure on the flow system. Compressibility is a parameter included in the storage term, defining the fluid storage capacity of the medium. While air entry pressure is a crucial value of the water retention curve, defining the relation between saturation and capillary pressure. From a series of sensitivity analyses, low compressibility resulted in faster flow due to low storage term, while low air entry pressure slowed flow inflow into the bentonite. Low air entry pressure means the air easily enters the medium; hence the flow rate becomes lower based on the relativity permeability definition. Based on the sensitivity analysis, we further investigate the effect of shotcrete around the tunnel and excavation damaged zone. Also, long-term analysis considering heat decay of the radioactive waste will be considered in future studies.

In high-level radioactive waste disposal, a high temperature is generated from the canister containing the waste in the engineered barrier, while groundwater flows into the buffer system from the host rock. The temperature increase and groundwater inflow result in the water phase change and saturation variation. Saturation change is related to the thermal conductivity of buffer material; hence the phase change and saturation strongly interact with the temperature evolution. The complex coupled behavior affects the stability of the whole disposal system, and the security of the repository is critical to human-being life. However, it is difficult to predict the long-term coupled behavior in the disposal system due to the considerable field-test scale, and therefore a numerical simulation is a suitable method having repeatability and cost-effectiveness. DECOVALEX is an international cooperating project for developing numerical methods and models for thermo-hydro-mechanical-chemical (THMC) interaction. DECOVALEX has a four-year cycle with various topics. At the current phase, Task C aims to simulate the full-scale emplacement (FE) experiment performed at Mont Terri underground rock laboratory. Nine research groups are participating in the task, and among them, KAERI simulates the experiment using OGS-FLAC. The simulator combines OpenGeoSys for TH simulation and FLAC3D for M simulation. Through the benchmark simulation, we verified OGS-FLAC for the two-phase flow analysis in the disposal system and finally modeled the FE experiment with a three-dimensional grid. We performed a simple sensitivity analysis to investigate the effect of input parameters on the two-phase flow system and confirmed that the compressibility and permeability affected the flow behavior. We also compared the simulation results to the field data and obtained well-matched results from a series of simulation.

본 연구에서는 이수의 유동 특성을 분석하기 위한 기초 연구로서 상용 코드인 ANSYS CFX 14.5를 이용하여 고체-액체 2상 유동 에 대한 수치해석적 연구를 수행하였다. 고체-액체 2상 유동 현상을 모사하기 위해서 균질류 모델과 분리류 모델을 사용하였다. 분리류 모델에서는 Gidaspow의 항력모델을 적용하였으며, 고체 입자에 운동 이론 모델을 적용하였다. 기존의 실험 결과를 기반으로 본 연구에서 사용한 수치해석 모델의 유효성을 검토하였으며, 수치해석은 직경 54.9 mm, 길이 3 m의 수평관에서 체적 분율 0.1~0.5, 속도 1~5 m/s 범위에서 수행되었다. 그리고 압력강하와 고체 입자의 체적 분율 분포를 확인하였으며, 압력강하는 균질류 모델과 분리류 모델이 각각 MAE 17.04%, 8.98 % 이내에서 실험결과를 잘 예측하였다. 관의 하부에서 높은 체적 분율이 나타나며, 상부로 갈수록 체적 분율은 감소하였다. 그리고 속도가 증가할수록 높이 변화에 따른 체적 분율 분포의 변화는 감소하였으며, 수치해석 결과는 이러한 유동 특성을 잘 예측하였다.

As the fabrication technology used in FPDs(flat-panel displays) advances, the size of these panels is increasing and the pattern size is decreasing to the um range. Accordingly, a cleaning process during the FPD fabrication process is becoming more important to prevent yield reductions. The purpose of this study is to develop a FPD cleaning system and a cleaning process using a two-phase flow. The FPD cleaning system consists of two parts, one being a cleaning part which includes a two-phase flow nozzle, and the other being a drying part which includes an air-knife and a halogen lamp. To evaluate the particle removal efficiency by means of two-phase flow cleaning, silica particles 1.5μm in size were contaminated onto a six-inch silicon wafer and a four-inch glass wafer. We conducted cleaning processes under various conditions, i.e., DI water and nitrogen gas at different pressures, using a two-phase-flow nozzle with a gap distance between the nozzle and the substrate. The drying efficiency was also tested using the air-knife with a change in the gap distance between the air-knife and the substrate to remove the DI water which remained on the substrate after the two-phase-flow cleaning process. We obtained high efficiency in terms of particle removal as well as good drying efficiency through the optimized conditions of the two-phase-flow cleaning and air-knife processes.

This study is focused on the channel design of bipolar plate in the electrode of hydrogen gas generator. The characteristics of hydrogen gas generation was studied in view of efficiency of hydrogen gas generation rate and a tendency of gas flow through the riv design of electrode. Since the flow rate and flow pattern of generated gas in the two phase flow system are the most crucial in determining the efficiency of hydrogen gas generator, we adopted the commercial analytical program of COMSOL MultiphysicsTM to calculate the theoretical flow rate of hydrogen gas from the outlet of gas generator and flow pattern of two phase fluid in the electrode. In this study, liquid electrolyte flows into the bipolar plate and decomposed into gas phase, two phase flow simulation is applied to measure the efficiency of hydrogen gas generation.

가스상의 체적분율과 압력강하는 기액이상류에 대한 이해와 예측에 있어서 매우 중요한 인자이다. 또한 그것들은 산업용 대용량의 열교환시스템 및 선박에 설치되는 보일러 및 냉동시스템의 설계에 있어서 필수적인 항목이다. 따라서 본 논문에서는 파이프의 모든 경사각도에서 기액이상류 가스상의 체적분율과 압력손실을 예측할 수 있는 이론적 해석 방법을 제시한다. 여기서의 이론적 해석은 2유체 층상류 모델을 기초로 하고 있다. 또한 이론적 해석결과와 기존의 실험결과와 비교한 결과에 대해서도 제시한다.

Our goal is to present a simple volume-of-fluid type interface-tracking algorithm to compressible two-phase flow in two space dimensions. The algorithm uses a uniform underlying Cartesian grid with some cells cut by the tracked interfaces into two subcells. A volume-moving procedure that consists of two basic steps: (1) the update of volume fractions in each grid cell at the end of the time step, and (2) the reconstruction of interfaces from discrete set of volume fractions, is employed to follow the dynamical behavior of the interface motion. As in the previous work with a surface-tracking procedure for general front tracking (LeVeque & Shyue 1995, 1996), a high resolution finite volume method is then applied on the resulting slightly nonuniform grid to update all the cell values, while the stability of the method is maintained by using a large time step wave propagation approach even in the presence of small cells and the use of a time step with respect to the uniform grid cells. A sample preliminary numerical result for an underwater explosion problem is shown to demonstrate the feasibility of the algorithm for practical problems.

원자로 내 사고발생 시 냉각수의 비등으로 기포가 발생하고, 기포율을 측정하기 위하여 열수력 안전 분 야에서는 주로 Optical Fiber Probe(OFP)나 광학 카메라를 이용하여 측정하지만 기하학적 구조의 한계로 인 해 17 × 17 배열의 봉 다발 내에 장비를 설치하는 것에는 어려움이 있다. 본 연구는 예비 연구로서 봉 다발 에 적용하기 전 X선 시스템과 다양한 모사 팬텀을 이용하여 연구 가능성 평가를 수행하였다. 라디오그라피 및 토모그라피 실험을 통해 X선 발생 장치의 관전압 130 kVp, 관전류 1 mA가 적합하였다. 또한, 기포 해상 도 팬텀을 통해 가시적으로 1 mm 크기의 구멍에 대해 측정이 가능하였으며 막대 팬텀을 이용한 대조도 평 가의 경우 프레온 내부에서 대조도가 상대적으로 떨어짐을 확인할 수 있었다. 그러나 영상 재구성 시 일그 러짐이 없는 좋은 영상을 획득할 수 있었다. 기포 발생 팬텀 실험을 통해 기포의 유동 방향 확인 및 단층 영상을 획득할 수 있었고, Image J 툴을 이용하여 하나의 단층영상에 대해 18 %의 기포율을 측정할 수 있 었다. 본 연구는 핵연료 주변 기포율 측정을 위한 선행 연구를 수행하였고 지속적인 연구를 위한 기초 연 구로서 활용할 수 있을 것이다.

단일절리에서 2상유체의 거동을 모의하기 위해 개발된 수치모형의 검증을 위해서 상대투과계수 특성식을 수치모형에 적용하여 가스와 물의 동시거동을 해석한 후, 수치모의 결과를 모형실험결과와 비교하였다. 절리면의 거칠기와 간극의 크기에 댸한 민감도 분석을 실시한 결과, 가스의 이동속도는 절리면의 거칠기와 반비례하였으며, 절리간극의 크기와는 상대투과계수 특성식의 영향으로 단상유체의 흐름에서와 같은 간극크기의 제곱에 비례하는 향상은 보이지 않았다. 수치모형의 현장

본 연구에서는 단일절리에서 2상유체 동시거동을 해석하기 위해서 2차원 유한차분 수치모형을 개발하였다. 개발된 모형은 압력에 따른 점성의 변화가 포화도에 따른 상대투과계수의 변화를 절리간극의 크기별로 고려할 수 있다. 수치기법으로는 IMPES해법을 적용하여 물과 가스의 압력변화량과 포화도를 차례로 구하였다. 개발된 수치 모형에 이용할 상대투과계수의 특성식 도출을 위해서 일곱가지 경우의 평판모형실험을 실시하였다. 실험으로부터 도출된 상대투과계수 특성곡선은