• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.169-177
• /
• 2005

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

• Journal of the Korean Applied Science and Technology
• /
• v.27 no.4
• /
• pp.415-420
• /
• 2010

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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• pp.671-686
• /
• 2005

The transient nature and complex geometries of two-phase gas-liquid flows cause fundamental difficulties when measuring flow velocity using an electromagnetic flowmeter. Recently, a current-sensing flowmeter was introduced to obtain measurements with high temporal resolution (Ahn et al.). In this study, current-sensing flowmeter theory was applied to measure the fast velocity transients in slug flows. The velocity fields of axisymmetric gas-liquid slug flow in a vertical pipe were obtained using Volume-of-Fluid (VOF) method, and the virtual potential distributions for the electrodes of finite size were also computed using the finite volume method for simulating slug flow. The output signal prediction for slug flow was carried out from the velocity and virtual potential (or weight function) fields. The flowmeter was numerically calibrated to obtain the cross-sectional liquid mean velocity at an electrode plane from the predicted output signal. Two calibration parameters are proposed for this procedure: a flow pattern coefficient and a localization parameter. The flow pattern coefficient was defined by the ratio of the liquid resistance between the electrodes for two-phase flow with respect to that for single-phase flow, and the localization parameter was introduced to avoid errors in the flowmeter readings caused by liquid acceleration or deceleration around the electrodes. These parameters were also calculated from the computed velocity and virtual potential fields. The results can be used to obtain the liquid mean velocity from the slug flow signal measured by a current-sensing flowmeter.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1951-1956
• /
• 2004

The transient nature and complex flow geometries of two-phase gas-liquid flows cause fundamental difficulties when measuring flow velocity using an electromagnetic flowmeter. Recently, a current-sensing flowmeter was introduced to obtain measurements with high temporal resolution (Ahn et $al.^{(1)}$). In this study, current-sensing flowmeter theory was applied to measure the fast velocity transients in slug flows. To do this, the velocity fields of axisymmetric gas-liquid slug flow in a vertical pipe were obtained using Volume-of-Fluid (VOF) method and the virtual potential distributions for the electrodes of finite size were also computed using the finite volume method for the simulated slug flow. The output signal prediction for slug flow was carried out from the velocity and virtual potential (or weight function) fields. The flowmeter was numerically calibrated to obtain the cross-sectional liquid mean velocity at an electrode plane from the predicted output signal. Two calibration parameters are required for this procedure: a flow pattern coefficient and a localization parameter. The flow pattern coefficient was defined by the ratio of the liquid resistance between the electrodes for two-phase flow with respect to that for single-phase flow, and the localization parameter was introduced to avoid errors in the flowmeter readings caused by liquid acceleration or deceleration around the electrodes. These parameters were also calculated from the computed velocity and virtual potential fields. The results can be used to obtain the liquid mean velocity from the slug flow signal measured by a current-sensing flowmeter.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.6
• /
• pp.714-724
• /
• 1985

Two-phase(air-solid, air-liquid droplet) turbulent round jet has been analyzed numerically using two equation turbulence model. The mean motion of suspending particles in air has been treated as the secondary fluid with virtual density and eddy viscosity. In this paper, the local mean velocity of secondary fluid is not assumed to be the same as that of the primary one. Dissipation rate of turbulent kinetic energy which arises because the particles can not catch up with the turbulent fluctuations of the primary fluid has been modelled by using the concept of Kolmogorov's spectral energy transfer. Numerical computations were performed for flows with different volume fraction of the dispersed phase and the diameter of particle. Results show that the total rate of turbulent energy dissipation, turbulent intensities and spreading rate of jets are reduced by the increase of volume fraction of dispersed phase. However it does not show consistent tendency with increasing the particle diameter. This investigation also shows that presence of particles in the fluid modifies the structure of the primary fluid flow significantly. Predicted velocity profiles and turbulence properties qualitatively agree with available data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.6
• /
• pp.446-454
• /
• 2008

The purpose of this paper is to analyze two-phase flows of the hydrogen recirculation system. Two-phase flow modeling is one of the great challenges in the classical sciences. As with most problems in engineering, the interest in two-phase flow is due to its extreme importance in various industrial applications. In hydrogen recirculation systems of fuel cell, the changes in pressure and temperature affect the phase change of mixture. Therefore, two-phase flow analysis of the hydrogen recirculation system is very important. Two-phase computation fluid dynamics (CFD) calculations, using a commercial CFD package FLUENT 6.2, were employed to calculate the gas-liquid flow. A two-phase flow calculation was conducted to solve continuity, momentum, energy equation for each phase. Then, the mass transfer between water vapor and liquid water was calculated. Through an experiment to measure production of liquid water with change of pressure, the analysis model was verified. The predictions of rate of condensed liquid water with change of pressure were within an average error of about 5%. A comparison of experimental and computed data was found to be in good agreement. The variations of performance, properties, mass fraction and two-phase flow characteristic of mixture with resepct to the fuel cell power were investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.

• Journal of the Korean Geotechnical Society
• /
• v.38 no.4
• /
• pp.5-20
• /
• 2022

To analyze the flow and density variations in debris flows, a two-phase finite volume model simplified with momentum equations was constructed in this study. The Hershel-Buckley rheology model was employed in this model to account for the internal and basal friction of debris flows and was utilized to analyze complex topography and entrainments of basal soil beds. In order to numerically solve the debris flow analysis model, a finite volume model with the Harten-Lax-van Leer-Contact method was used to solve the conservation equation for the debris flow interface. Case studies of circular dam failure, non-Newtonian fluid dam failure, and multiple debris flows were analyzed using the proposed model to evaluate shock absorption capacity, numerical isotropy, model accuracy, and mass conservation. The numerical stability and correctness of the debris flow analysis of this analysis model were proven by the analysis results. Additionally, the rate of debris flow with various rheological properties was systematically simulated, and the effect of debris flow rheological properties on behavior was analyzed.

• Journal of Mechanical Science and Technology
• /
• v.15 no.7
• /
• pp.941-950
• /
• 2001

A study of counter-current two-phase flow in narrow rectangular channels has been performed. Two-phase flow regimes were experimentally investigated in a 760mm long and 100mm wide test section with 2.0 and 5.0mm gap widths. The resulting flow regime maps were compared with the existing transition criteria. The experimental data and the transition criteria of the models showed relatively good agreement. However, the discrepancies between the experimental data and the model predictions of the flow regime transition become pronounced as the gap width increased. As the gap width increased the transition gas superficial velocities increased. The critical void fraction for the bubbly-to-slug transition was observed to be about 0.25. The two-phase distribution parameter for the slug flow was larger for the narrower channel. The uncertainties in the distribution parameter could lead to a disagreement in slug-to-churn transition between the experimental findings and the transition criteria. For the transition from churn to annular flow the effect of liquid superficial velocity was found to be insignificant.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2037-2042
• /
• 2003

The purpose of this research is to develop the model of turbulence modulation due to the presence of particles in various types of particle-laden flows. Available experimental data were surveyed and the dependence of turbulence modulation of carrier-phase on particle size, concentration and particle Reynolds number were examined. This study takes into account the effect of wake produced by particle, the drag between phases and the velocity gradient in the wake to estimate the production of turbulence. The model of turbulence modulation using the mixing length theory under the assumption of equilibrium flow is proposed. Numerical results show that the model is successful in predicting the characteristics of the particle-laden flow in various conditions both qualitatively and quantitatively.