• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.126-134
• /
• 2008

We developed a three dimensional Navier-Stokes code based on the level set method to simulate two phase flows with high density ratio. The Navier-Stokes equations with consideration of the surface tension effects are solved by using SIMPLE algorithm on a non-staggered grid. The present code is validated by simulating two test problems. First one is to simulate a rising bubble inside a cube. The thickness of the interface of the bubble is shown to affect the pressure distribution around the interface. As the thickness decreases, the pressure field around the interface becomes more oscillatory. As the bubble rises, a ring vortex is shown to form around the interface and the bubble eventually develops into an ellipsoidal shape. Merge of two bubbles inside a container is secondly tested to show the robustness of the present code for two phase flow simulation. Numerical results show stable and reliable behavior during the process of merging of two bubbles. The velocity and pressure fields around the interface of bubbles are shown oscillation free during the merging of two bubbles.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.5
• /
• pp.19-29
• /
• 2006

In this study, numerical simulation of airframes separating from a missile system has been preformed. For the time-accurate trajectory simulation, six D.O.F equations of motion of multiply connected bodies were derived and these equations have been coupled with the unstructured overset mesh technique for the treatment of independent mesh blocks moving with each body component. Applications were made for the simulation of the airframe separation at missile angles of attack of 0 and 5 degrees. It was demonstrated that the present method is efficient and robust for the prediction of unsteady time-accurate flow fields involving multiple bodies in relative motion.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1348-1358
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation fur flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among seven filling patterns at each tetrahedral control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. The collapse of a water dam and the filling of a fluidity spiral have been analyzed. The numerical results have been in good agreement with the experimental results and the efficiency of the adaptive grid refinement and filling pattern techniques have been verified.

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

This paper presents three-dimensional flow analysis for a mixed-flow pump which consists of a rotor and a stator. Reynolds-averaged Navier-Stokes equations with shear stress transport turbulence model are discretized by finite volume approximations and solved by the commercial CFD code CFX 11.0. Structured grid system is constructed in the computational domain, which has O-type grids near the blade surfaces and H-type grids in other regions. Validation of the numerical results was performed with experimental data for head coefficients and hydraulic efficiencies at different flow coefficients. This paper shows that the pump characteristics can be predicted effectively by numerical analysis.

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

An analysis is made of flow and rocket motion during a supersonic separation stage of air-launching rocket from the mother plane. Three-dimensional Euler and Navier-Stokes equations are numerically solved to analyze the steady/unsteady flow field around the rocket which is being separated from two cases of mother plane configuration: one is an idealized ogive-cylinder body and the other is a real F-4E Phantom. The simulation results clearly demonstrate the effect of shock-expansion wave interaction between the rocket and the mother plane. As a result, a design-guideline of supersonic air-launching rocket for the safe separation is proposed.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.175-181
• /
• 1995

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.45-48
• /
• 2007

In the present work, flow analysis has been performed for side channel type double stage ring blower by solving three-dimensional Reynolds-averaged Navier-Stokes equation. Shear stress transport model is used as turbulent closure. The commercial CFD code CFX 11.0 is used for the calculations. Each of two stage is calculated separately and the second stage inlet flow is same as the first stage outlet flow so that consecutive calculation is possible. Velocity and pressure fields have been analyzed at the midplane between blades. The numerical results are validated with experimental data for head coefficients at different flow coefficients.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.79-82
• /
• 2002

Three-dimensional pipe flows with elbows and tees for few different pipe fittings are calculated to estimate the effect of swirling flow on measuring accuracy of orifice flow meter. It is evaluated how the pressure difference across the orifice is dependent on the length of upstream straight pipe in a branch and how swirl intensity, swirl angel and axial velocity distribution affect the measuring error of orifice flowmeter. From the results, it is found that, regardless of flow rate specified in this calculation, the effect of the straight pipe length can be neglected for the lengths larger than thirty diameters although there still remain significant swirl at the orifice

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.95-98
• /
• 2002

Three-dimensional flow analyses for two different ratios of radius to height of sedimentation bed are implemented to evaluate the effect of blockage ratio of center feed wall and angle of distributor on sedimentation efficiency, and to find the optimal value of those parameters. Sedimentation efficiencies for three different shapes are compared with and without rotation speed. And then, five different combinations of blockage ratio of center feed wall and angle of distributor are compared It reveals that the effect of blockage ratio of center fled wall and angle of distributor is considerable to sedimentation efficiency while rotation effect can be neglected and $0.55 and 33^{\circ}$for blockage ratio of center food wall and angle of distributor, respectively, ive the best sedimentation efficiency.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.71-74
• /
• 2002

Fluid flow in a rectangular duct system are measured by W laser doppler velocity meter, and also computed by commercial software of STAR-CD for comparison between then First, for a rectangular duct with 90 degree metered elbow, the fluid flow with Reynolds numbs's of 1,508 is predicted by assumption of both laminar and turbulent models. But, even though the Reynolds number is less than 2,300-3,000, the computation by turbulent model is close to the experimental data. Moeover, the computation by turbulent model for Reynolds number of 11,751 also predicts the experimental data satisfactorily. Second, for a rectangular duct with two branch ducts, the ratios between flow rates in the two branches are invariant to Reynolds number according to both of numerical and experimental results.