• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.291-293
• /
• 2011

In this paper, the problem of transonic aerodynamic characteristics of a NACA0012 airfoil is numerically investigated in the inviscid gas-droplet two-phase flow with the compressible two-fluid model. In the present study, the airfoil flight in the cloud is simulated by taking account of the viscous drag of the droplets, the heat transfer, the phase change, and the droplet fragmentation The two-fluid equation system is solved by the fractional-step method and the WAF-HIL scheme. The effects of size and volume fraction of the droplets on the flow characteristics of the airfoil in the cloud are elaborated and discussed.

• Tribology and Lubricants
• /
• v.13 no.3
• /
• pp.56-62
• /
• 1997

Sealing an oil-air mixture plays important roles to have an enhanced lubrication for high speed spindle. High speed spindles require non-contact type sealing mechanism. In this study, an optimum seal design to minimize leakage is concerned in the aspect of flow control. This paper categorizes geometries of mostly used non-contact type seals and analyzes each leakage characteristics to minimize a leakage on sealing area. Effect of minimum clearance and its position are considered according to variation of detail geometry. The estimation of non-leaking property is determined by amount of pressure drop in the leakage path assuming constant leakage flow. To simulate an oil jet or oil mist type high speed spindle lubrication, the working fluid is regarded as two phases that are mixed flow of oil phase and air phase. Both of the turbulence and the compressible flow model were introduced in CFD(Computational Fluid Dynamics) analysis. Design parameters has been induced to minimize leakage in limited space, and a methodological study on geometrical optimization has been conducted.

• Journal of computational fluids engineering
• /
• v.18 no.4
• /
• pp.9-16
• /
• 2013

Cavitation on an axi-symmetric hemispherical head-form body was studied using an Reynolds-averaged Navier-Stokes equations solver based on a cell-centered finite volume method. To consider compressibility effects on the vapor phase and cavity interface, a pressure-based compressible flow CFD code was developed. To validate the developed CFD code, cavitating flow around the hemispherical head-form body was simulated using pressure-based incompressible and compressible CFD codes and validated against existing experimental data in the three-way comparison. The cavity shedding behavior, length of re-entrant jet, drag history, and Strouhal number of the hemispherical head-form body were compared between two CFD codes. The results, in this paper, suggested that the computations of cavitating flow with compressibility effects improve the description of cavity dynamics.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• Journal of KSNVE
• /
• v.6 no.2
• /
• pp.163-177
• /
• 1996

An investigation on thermohydraulic stability of flow oscillations in the CANada Deuterium Uranium-600(CANDU-6) heat transport system has been conducted. Flow oscillations in reactor coolant loops, comprising two heat sources and two heat sinks in series, are possibly caused by the response of the pressure to extraction of fluid in two-phase region. This response consists of two contributions, one arising from mass and another from enthalpy change in the two-phase region. The system computer code used in the investigation os SOPHT, which is capable of simulating steady states as well as transients with varying boundary conditions. The model was derived by linearizing and solving one-dimensional, homogeneous single- and two-phase flow conservation equations. The mass, energy and momentum equations with boundary conditions are set up throughout the system in matrix form based on a node-link structure. Loop stability was studied under full power conditions with interconnecting the two compressible two phase regions in the figure-of-eight circuit. The dominant function of the interconnecting pipe is the transfer of mass between the two-phase regions. Parametric survey of loop stability characteristics, i. e., damping ratio and period, has been made as a function of geometrical parameters of the interconnection line such as diameter, length, height and orifice flow coefficient. The stability characteristics with interconnection line has been clarified to provide a simple criterion to be used as a guide in scaling of the pipe.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.115-119
• /
• 2009

We numerically investigated propagation of various waves in the two-phase flows such as sound wave, shock wave, rarefaction wave, and contact discontinuity in terms of pressure, void fraction, velocity and density of the two phases. The waves have been generated by a hydrodynamic shock tube, a pair of symmetric impulsive expansion, impulsive pressure and impulsive void waves. The six compressible two-fluid two-phase conservation laws with interfacial friction terms have been solved in two fractional steps. The first PDE Operator is solved by the HLL scheme and the second Source Operator by the semi-implicit stiff ODE solver. In the HLL scheme, the fastest wave speeds were estimated by the analytic eigenvalues of an approximate Jacobian matrix. We have discussed how the interfacial friction terms affect the wave structures in the numerical solution.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.653-662
• /
• 2007

Due to the difficulty raised from the coupling of cavitation modeling with turbulent flow, numerical simulation for two phase flow remains one of the challenging issues in the society. This research focuses on the development of numerical code to deal with incompressible two phase flow around 2D hydrofoil by combing the cavitation model suggested by Kunz et al. with $k-{\varepsilon}$ turbulent model. The simulation results are compared to experimental data to verify the validity of the developed code. Also, the comparison of the calculation results is made with LES results to evaluate the capability of $k-{\varepsilon}$ turbulence model. The calculation results show very good agreement with experimental observations even though this code can not grasp the small scaled bubbles in the calculation wheres LES can hold the real physics. This code will be extended to 3D compressible two phase flow for the study on the fluid dynamics in the inducers and impellers.

• Journal of the Korean Mathematical Society
• /
• v.51 no.1
• /
• pp.163-187
• /
• 2014

We present a multi-stage Roe-type numerical scheme for a model of two-phase flows arisen from the modeling of deflagration-to-detonation transition in granular materials. The first stage in the construction of the scheme computes the volume fraction at every time step. The second stage deals with the nonconservative terms in the governing equations which produces states on both side of the contact wave at each node. In the third stage, a Roe matrix for the two-phase is used to apply on the states obtained from the second stage. This scheme is shown to capture stationary waves and preserves the positivity of the volume fractions. Finally, we present numerical tests which all indicate that the proposed scheme can give very good approximations to the exact solution.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.362-364
• /
• 2010

Shock wave interaction with droplet-gas medium is investigated in this paper. In the present computation, the shock wave is initially started in a pure gas and reflected from the wedge to interact with the droplet-ridden gas flows. We used the compressible two-fluid two-phase model that is solved by the two-fluid version of the HLL scheme. The interfacial drag force and heat transfer were included to model the interaction between continuous and dispersed phases. The parametric effect of void fraction on the shock wave reflection in the two-phase media was investigated.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.264-267
• /
• 2010

This research focuses on the development of numerical code to deal with compressible two phase flow around three dimensional objects combined with cavitation model suggested by Weishyy et al. with k-e turbulent model. The cryogenic cavitation is carried out by considering the thermodynamic effect on physical properties of cryogenic fluids in physical point of view and implementing the temperature sensitivity in the energy equation of the government equations in numerical point of view, respectively. The formulation has been extensively validated for both liquid nitrogen and liquid hydrogen by simulating the experiments of Hord on hydrofoils. Then, simulations of cavitating turbopump inducers at their design flow rate are presented. Results over a broad range of Nss numbers extending from single-phase flow conditions through the critical head break down point are discussed. In particular, thermal depression effects arising from cavitation in cryogenic fluids are identified and their impact on the suction performance of the inducer quantified.