• Nuclear Engineering and Technology
• /
• 제56권3호
• /
• pp.1002-1012
• /
• 2024

To remove insoluble fission products, which could possibly cause reactor instability and significantly reduce heat transfer efficiency from primary system of molten salt reactor, a helium bubbling method is employed into a passive molten salt fast reactor. In this regard, two-phase flow behavior of molten salt and helium bubbles was investigated experimentally because the helium bubbles highly affect the circulation performance of working fluid owing to an additional drag force. As the helium flow rate is controlled, the change of key thermal-hydraulic parameters was analyzed through a two-phase experiment. Simultaneously, to assess the applicability of numerical model for the analysis of two-phase flow behavior, the numerical calculation was performed using the OpenFOAM 9.0 code. The accuracy of the numerical analysis code was evaluated by comparing it with the experimental data. Generally, numerical results showed a good agreement with the experiment. However, at the high helium injection rates, the prediction capability for void fraction of helium bubbles was relatively low. This study suggests that the multiphaseEulerFoam solver in OpenFOAM code is effective for predicting the helium bubbling but there exists a room for further improvement by incorporating the appropriate drag flux model and the population balance equation.

• 한국항공우주학회지
• /
• 제47권3호
• /
• pp.204-211
• /
• 2019

희박연료 직접분사(Lean Direct injection(LDI)) 가스터빈 연소기에 대한 이상유동 특성을 해석하였다. LDI 연소기에 적용된 환형분사기(hollow-cone spray injector)의 분열을 모사하기 위해 분열모델(Linearized Instability Sheet Atomization(LISA), Aerodynamically Progressed Taylor Analysis Breakup(APTAB)을 적용하였다. 침투깊이와 평균입도(Sauter Mean Diameter(SMD))를 통해 분열모델을 검증하였으며, LDI 연소기에 적용하여 이상유동특성을 분석하였다. 스월인젝터로 인해 Precessing Vortex Core(PVC)가 발생하였으며, 액적들이 PVC를 따라 미립화되는 것을 확인하였다. SMD 결과를 통해 PVC가 회전하는 영역의 외곽으로 즉, 빠른 속도 영역에 액적들이 분포하며, 스톡스수(Stokes number)는 1보다 작다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2005년도 제24회 춘계학술대회논문집
• /
• pp.334-342
• /
• 2005

A numerical analysis is performed using two dimensional axisymmetric RANS (Reynolds Averaged Navier-Stokes) equations system on the transition sequence of the Integrated Rocket Ramjet and the unsteady reacting flow-field in a ramjet combustor during unstable combustion. The mode transition of an axisymmetric ramjet is numerically simulated starting from the initial condition of the boost end phase of the entire ramjet. The unsteady reacting flow-field within combustor is computed for varying injection area. In calculation results of the transition, the terminal normal shock is occurred at the downstream of diffuser throat section and no notable combustor pressure oscillation is observed after certain time of the inlet port cover open. For the case of a small injection area at the same equivalence ratio, periodic pressure oscillation in the combustor leads to the terminal shock expulsion from the inlet and hence the buzz instability occurred.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회B
• /
• pp.2524-2529
• /
• 2008

The thermal system like a combustion chamber is believed to experience a significant instability problem with vibration in case that the thermal energy or the acoustic energy are transformed into a different form through a relevant path. This study deals with a numerically- predicted, Thermoacoustic instability in a Rijke tube by using a non-linear model for a heat source. The heating part where the energy transformation occurs actively is modeled after simulating two-dimensional cylinder case with constant surface temperature, and a nonlinear model that accounts for the transfer function of magnitude- and phase-characteristics is properly implemented so as to be dependent on the pulsation strength in the tube. The heat source model is observed to result in equivalent Thermoacoustic instabilities in the Rijke tube except low flow-rate cases in which the natural convection is dominant.

• Nuclear Engineering and Technology
• /
• 제54권5호
• /
• pp.1851-1859
• /
• 2022

Experimental research on flow characteristics of open natural circulation system was performed, to figure out the mechanism of the open natural circulation behaviors. The influence factors, such as the heating power, the inlet subcooled and the level of cooling tank on the flow characteristics of the system were examined. It was shown that within the scope of the experimental conditions, there are five flow types: single-phase stable flow, flash and geyser coexisting unstable flow, flash stable flow, flash unstable flow, and flash and boiling coexisting unstable flow. The geyser flow in flash and geyser coexisting unstable flow is different from classic geysers flow. The flow oscillation period and amplitude of the former are more regular, is a newly discovered flow pattern. By drawing the flow instability boundary diagram and sorting out the flow types, it is found that the two-phase unstable flow is mainly characterized by boiling and flash, which determine the behavior of open natural circulation respectively or jointly. Moreover, compared with full liquid level system, non-full liquid level system is more prone to boiling phenomenon, and the range of heat flux density and undercooling degree corresponding to unstable flow is larger.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 춘계학술대회논문집
• /
• pp.572-575
• /
• 2008

Two-dimensional multiphase flows due to density difference such as the Rayleigh-Taylor instability problem and the droplet splash are simulated by an in-house solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present code simulates complex free surface flows such as multiphase flows due to density difference efficiently and accurately.

• Nuclear Engineering and Technology
• /
• 제53권1호
• /
• pp.61-68
• /
• 2021

Experiments of vertically upward steam-water two-phase flow have been carried out in single-side heated narrow rectangular channel with a gap of 3 mm. Flow patterns were identified and classified through visualization directly. Slug flow was only observed at 0.2 MPa but replaced by block-bubble flow at 1.0 MPa. Flow pattern maps at the pressure of 0.2 MPa and 1.0 MPa were plotted and the difference was analyzed. The experimental data has been compared with other flow pattern maps and transition criteria. The results show reasonable agreement with Hosler's, while a wide discrepancy is observed when compared with air-water two-phase experimental data. Current criteria developed based on air-water experiments poorly predict bubble-slug flow transition due to the different formation and growth of bubbles. This work is significant for researches on heat transfer, bubble dynamics and flow instability.

• Structural Engineering and Mechanics
• /
• 제81권4호
• /
• pp.523-527
• /
• 2022

The objective of this study is to simulate the flow in pipes with various boundary conditions. Free-pressure fluid model, is used in the pipe based on Navier-Stokes equation. The models are solved by using the numerical method. A problem called "stability of pipes" is used in order to compare frequency and critical fluid velocity. When the initial conditions of problem satisfied the instability conditions, the free-pressure model could accurately predict discontinuities in the solution field. Employing nonlinear strains-displacements, stress-strain energy method the governing equations were derived using Hamilton's principal. Differential quadrature method (DQM) is used for obtaining the frequency and critical fluid velocity. The results of this paper are analyzed by hyperbolic numerical method. Results show that the level of numerical diffusion in the solution field and the range of well-posedness are two important criteria for selecting the two-fluid models. The solutions for predicting the flow variables is approximately equal to the two-pressure model 2. Therefore, the predicted pressure changes profile in the two-pressure model is more consistent with actual physics. Therefore, in numerical modeling of gas-liquid two-phase flows in the vertical pipe, the present model can be applied.

• Nuclear Engineering and Technology
• /
• 제56권5호
• /
• pp.1902-1912
• /
• 2024

This paper presents results from system-level modeling of a water-based reactor cavity cooling system using RELAP5-3D. The computational model is benchmarked with experimental data from a half-scale RCCS test facility at Argonne National Laboratory. The model prediction is first compared with a two-phase oscillatory baseline experimental case where mixed accuracy is obtained. The model shows reasonable prediction of mass flow rate, pressure, and temperature but significant overprediction of void fraction. The model prediction is then compared with a fault case where the inlet of the risers is gradually reduced using a throttling valve. As the valve is closed, the model is able to predict some major flow phenomena observed in the experiment such as the dampening of oscillations, the reintroduction of oscillations, as well as boiling, flashing, and geysering in the risers. However, the timeline of these events are not well captured by the model. The model is also used to investigate the evolution of flow regime in the chimney. This work highlights that the semi-empirical constitutive relations used in RELAP-3D could have a strong influence on the accuracy of the model in two-phase oscillatory flows.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.118-120
• /
• 2003

The behavior of viscous fingerings caused by an external force is investigated by using a two­phase lattice Boltzmann method. The effects of the modified capillary number, the viscosity contrast, and the modified Darcy-Rayleigh number on the instability of interfaces are found. The calculated wave numbers are in good agreement with the theoretical ones in the range of wave numbers smaller than 10, but the calculated ones tend to become smaller than the theoretical ones in higher wave numbers.