• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3033-3038
• /
• 2007

In this research, Rayleigh instability of gas-liquid flow in annular pipe is studied in film boiling using viscous potential flow. Viscous potential flow is a kind of approximation of gas-liquid interface considering velocity field as potential including viscosity. A dispersion relation is obtained including the effect of heat and mass transfer and viscosity. New expression for dispersion relation in film boiling and critical wave number is obtained. Viscosity and heat and mass transfer have a stabilizing effect on instability and its effect appears in maximum growth rate and critical wave number. And the existence of marginal stability region is shown.

• Ocean Systems Engineering
• /
• v.1 no.1
• /
• pp.31-57
• /
• 2011

Impact pressure due to sloshing is of great concern for the ship owners, designers and builders of the LNG carriers regarding the safety of LNG containment system and hull structure. Sloshing of LNG in partially filled tank has been an active area of research with numerous experimental and numerical investigations over the past decade. In order to accurately predict the sloshing impact load, a new numerical method was developed for accurate resolution of violent sloshing flow inside a three-dimensional LNG tank including wave breaking, jet formation, gas entrapping and liquid-gas interaction. The sloshing flow inside a membrane-type LNG tank is simulated numerically using the Finite-Analytic Navier-Stokes (FANS) method. The governing equations for two-phase air and water flows are formulated in curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered grid. Simulations were performed for LNG tank in transverse and longitudinal motions including horizontal, vertical, and rotational motions. The predicted impact pressures were compared with the corresponding experimental data. The validation results clearly illustrate the capability of the present two-phase FANS method for accurate prediction of impact pressure in sloshing LNG tank including violent free surface motion, three-dimensional instability and air trapping effects.

• Nuclear Engineering and Technology
• /
• v.56 no.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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.47 no.3
• /
• pp.204-211
• /
• 2019

The analysis on two-phase flow in a Lean Direct Injection(LDI) combustor has been investigated. Linearized Instability Sheet Atomization(LISA) and Aerodynamically Progressed Taylor Analogy Breakup(APTAB) breakup models are applied to simulate the droplet breakup process in hollow-cone spray. Breakup model is validated by comparing penetration length and Sauter Mean Diameter(SMD) of the experiment and simulation. In the LDI combustor, Precessing Vortex Core(PVC) is developed by swirling flow and most droplets are atomized along the PVC. It has been confirmed that all droplets have Stokes number less than 1.0.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• 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.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• 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
• /
• v.54 no.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.03b
• /
• 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
• /
• v.53 no.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
• /
• v.81 no.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.