• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.131-133
• /
• 2003

The streamfunction-vorticity equations for two-dimentional cavity flow are solved by a new finite element method which uses finite spectral basis functions as interpolation functions for rectangular elements. Results for several cases with different Renold's number are compared with benchmark solutions and found to be in well agreement.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.803-806
• /
• 2002

The dispersion of Lagrangian fluid particles in a turbulent channel flow is studied by a direct numerical simulation. Four points Hermite interpolation in the homogeneous direction and Chebyshev polynomials in the inhomogeneous direction is adopted by assesing the acceleration of fluid particles. In order to characterize the inhomogeneous Lagrangian statistics, accurate single particle Lagrangian statistics are obtained along the wall normal direction. Integral time scales of Lagrangian velocity can be normalized by Eulerian mean shear stresses.

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

In this study, it is performed numerical simulation on multi-phase flow by means of CIP-CSI2 scheme. It is applied In a two-phase free surface flow problem at a high density ratio equivalent to that of an air-water system, for examining the computational capability. The method that is being developed and improved is a CIP(Constrained Interpolation Profile) and CSL2(Conservative Semi-Lagrangian) based Cartesian Grid Method.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.1-10
• /
• 2011

A new computational program, which is based on the CIP/CCUP(Constraint Interpolation Profile/CIP Combined Unified Procedure) method, has been developed to numerically analyse sloshing phenomena dealt as multiphase-flow problems. For the convection terms of Navier-Stokes equations, the RCIP(Rational function CIP) method was adopted and the THINC-WLIC(Tangent of Hyperbola for Interface Capturing-Weighted Line Interface Calculation) method was used to capture the air/water interface. To validate the present numerical method, two-dimensional dam-breaking and sloshing problems in a rectangular tank were solved by the developed method in a stationary Cartesian grid system. In the case of sloshing problems, simulations by using a improved MPS(Moving Particle Simulation) method, which is named as PNU-MPS(Pusan National University-MPS), were also carried out. The computational results are compared with those of experiments and most of the comparisons are reasonably good.

• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.21-27
• /
• 2015

Flow field around the KRISO 3600TEU container ship is computed using a surface generated based on interpolations of station lines, which are given in a body plan of the ship, without using any CAD program. An interpolation method is suggested based on inscribed circles to generate curves between two neighboring station lines. The interpolated surface is saved in a STL format to use the snappyHexMesh utility of the openfoam. Computed resistance of the ship is compared with experimental and other computational results and the effects of the interpolation of neighboring station lines on the computed resistance are investigated. The suggested method is applied to calculate the flow field around a submarine with appendages. The surface triangulations for the hull and the appendages are generated without consideration of each other, then those surface triangulations are simply combined to provide a grid generator with the body boundary. The junctures of the hull and the appendages are identified automatically during the grid generation procedure. Tip vortex is captured, which travels downstream from the tip of the appendages.

• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.603-607
• /
• 2007

Recently, the frequency of unexpecting heavy rains has been increased due to abnormal climate and extreme rainfall. There was a limit to analyze one dimension or two dimension stream flow of domestic rivers that was applied simple momentum equation and fixed energy conservation. Therefore, hydrodynamics flow analysis in rivers has been needed three dimensional numerical analysis for correct stream flow interpolation. In this study, CFD model on FLOW-3D was applied to stream flow analysis, which solves three dimension RANS(Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behavior and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as $k-{\backepsilon}$, RNG $k-{\backepsilon}$ and LES. Those numerical analysis results have been illustrated to bends and junctions by the turbulence energy effects, velocity of flow distributions, water level pressure distributions and eddy flows.

• Journal of computational fluids engineering
• /
• v.6 no.3
• /
• pp.19-26
• /
• 2001

An Eulerian-Lagrangian method, so called immersed boundary method, is used for analysing viscous flow around arbitrary bodies, where governing equations are discretized on a regular grid by using a finite volume method. To improve the accuracy of flow near body boundaries, a second-order accurate interpolation scheme is used and a level-set based grid deformation method is presented to construct the adaptive grids around body boundaries. The present scheme is used to simulate steady flow around a semicircular cylinder mounted on the bottom of flow domain and calculated results are validated by results of a body fitted grid method. Finally, present method is applied to a complex flow around multi body and the usefulness is checked by investigating calculated results.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.105-110
• /
• 2001

An Eulerian-Lagrangian method, so called immersed boundary method, is used for analysing viscous flow around arbitrary bodies, where governing equations are discretized on a regular grid by using a finite volume method. To improve the accuracy of flow near body boundaries, a second-order accurate interpolation scheme is used and a level-set based grid deformation method is presented to construct the adaptive grids around body boundaries. The present scheme is used to simulate steady flow around a semicircular cylinder mounted on the bottom of flow domain and calculated results are validated by results of a body fitted grid method. Finally, present method is applied to a complex flow around multi body and the usefulness is checked by investigating calculated results.

• Nuclear Engineering and Technology
• /
• v.44 no.1
• /
• pp.71-78
• /
• 2012

In many semi-empirical analyses of flow boiling heat transfer, an annular flow is often assumed as a model flow and the local liquid film thickness is a key parameter in the analysis. This work considers a simple electrical conductance technique to estimate the local liquid film thickness in two-phase annular flows. In this approach, many electrodes are mounted flush with the inner wall of the pipe. Voltage differences between two neighboring electrodes for concentric annular flows with various liquid film thicknesses are obtained before the main experiments and logged in a look-up table. For an actual application in the annual flow, voltage differences of neighboring electrodes are measured and then corresponding local film thicknesses are determined by the interpolation of the look-up table. Even though the proposed technique is quite simple and straightforward, the numerical and static phantom experiments support its usefulness.