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

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

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

Bubble rising phenomenon is widely founded in many industrial applications such as a stream generator in power plant. Many experimental and numerical researches have been already performed to predict dynamic behavior of the bubble rising process. Recently numerical approaches are getting popular since it can offer much detailed information which is almost impossible to obtain from the experiments. Rising bubble could penetrate through the top free surface which makes the problem much more complicate in addition to the phase changing effect even with latest numerical techniques. In this paper, the top free surface effect on rising bubble has been investigated. The gas-liquid interface was explicitly tracked using high-order Level Contour Reconstruction Method(LCRM) which is a hybridization of Front-Tracking and Level-Set method. Break-up behavior of rising bubble at free surface showed different characteristics with initial diameter of bubble.

• Journal of Environmental Health Sciences
• /
• v.31 no.6
• /
• pp.489-497
• /
• 2005

A numerical method for the mathematical water modeling in 2-dimensional flow has been developed. The model based on a split operator technique, in which, the advection term is calculated using the upwind scheme. The diffusion term is one- dimensionalized and calculated using Crank-Nicholson's implicit finite difference scheme to reduce the numerical errors from large time steps and variable spacings. It also provides a relatively simple and economic method for more accurate simulation of pollutant dispersion. Water depths and flow velocities in the Boreyong reservoir during the normal water periods were predicted by numerical experiments with a 2-dimensional flow model so as to provide current field data for the study of advection and diffusion of pollutants. Developed 2-dimensional water quality model is applied to Boreyong reservoir to simulate a spatial and periodical changes of water quality.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.152-157
• /
• 1999

For the simulation of the blood flow in coronary artery, the system modeling of coronary hemodynamics is combined with CFD technique. The blood flow in coronary artery interacts with the global coronary circulation. Especially in case of the coronary artery with stenosis, the interaction plays an important role in the hemodynamics of the circulation. In this study we present a combined numerical approach using both the CFD technique for flow simulation and the global system model of coronary circulation. We use a lumped parameter model for the global simulation of coronary circulation whereas the finite element method is employed to compute the viscous flow field in stenosed coronary artery, The time variation of the pressure drop due to stenosis is obtained from the proposed numerical method. Numerical results shows that the flow resistance and pressure drop due to stenosis has a relatively large value in systole.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.156-161
• /
• 1998

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.3
• /
• pp.130-138
• /
• 1996

The transient incompressible flow behind the widely-spaced co-axial jet is numerically simulated using the random vortex method(RVM). This numerical approach is based on the Lagrangian approach for the vorticity formulation of the unsteady Navier-Stokes equations, utilizing vortex elements to account for the convection and diffusion processes. The effects of the mass flow rate of an annular air jet and a central fuel jet on the co-axial jet flow dynamics is investigated. To validate the present procedure, the numerical results are compared with the available experimental data the present procedure, the numerical results are compared with the available experimental data in terms of the centerline and off-centerline profiles of the mean axial velocity. Discrepancies between the RVM results and the measurements are discussed in detail.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.2 no.1
• /
• pp.27-39
• /
• 1998

We describe a fast numerical method for non-separable elliptic equations in self-adjoin form on irregular adaptive domains. One of the most successful results in numerical PDE is developing rapid elliptic solvers for separable EPDEs, for example, Fourier transformation methods for Poisson problem on a square, however, it is known that there is no rapid elliptic solvers capable of solving a general nonseparable problems. It is the purpose of this paper to present an iterative solver for linear EPDEs in self-adjoint form. The scheme discussed in this paper solves a given non-separable equation using a sequence of solutions of Poisson equations, therefore, the most important key for such a method is having a good Poison solver. High performance is achieved by using a fast high-order adaptive Poisson solver which requires only about 500 floating point operations per gridpoint in order to obtain machine precision for both the computed solution and its partial derivatives. A few numerical examples have been presented.

• Journal of Aerospace System Engineering
• /
• v.2 no.2
• /
• pp.20-27
• /
• 2008

The characteristics of compressible flow are different from those of incompressible flow from the mathematical and physical point of view. Therefore, the way to solve the flowfield is different between compressible flow and incompressible flow. In general, density-based numerical algorithm is mainly used for compressible flow solver development. On the other hand, incompressible flow solver prefers to use pressure-based numerical algorithm. In this research, a compressible Navier-Stokes flow solver is developed by means of extending from pressure-based incompressible numerical algorithm to handle both compressible and incompressible flows using the same flow solver. The present flow solver is tested at various speed ranges and compared with the solutions of density-based compressible flow solver. Numerical results show a good agreement between two flow solvers.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.6
• /
• pp.427-435
• /
• 2015

This paper provides the numerical results of a simulation of the flow around a propeller working beneath the free surface. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes (URANS) equations, where the wave-making problem is solved using a volume-of-fluid (VOF) method. The numerical analysis focuses on the propeller wake structure affected by the free surface, where we consider another free surface boundary condition that treats the free surface as a rigid wall surface. The propeller wake under the effect of these two free surface conditions shows a reduction in the magnitude of the longitudinal and vertical flow velocities, and its vortical structures strongly interact with the free surface. The thrust and torque coefficient under the free surface effect decrease about 3.7% and 3.1%, respectively. Finally, the present numerical results show a reasonable agreement with the available experimental data.

• Advances in aircraft and spacecraft science
• /
• v.1 no.2
• /
• pp.145-175
• /
• 2014

The results of a series of numerical experiments are presented to verify some of the important developments made in the first part of this paper. Firstly, the static solution of an algebraic system obtained through Strong Formulation Finite Element Method (SFEM) is presented. Secondly, the stress and strain recovery procedure is descripted for the present technique. It will be clear that the present approach is suitable for any strong formulation finite element methodology, due to the presented general approach based on the unknown displacements and on the elasticity equations. Thirdly, the numerical solutions for some classical and other numerical results found in literature are exposed. Finally, an arbitrarily shaped composite plate is solved and good agreement is observed for all the presented cases.