• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.777-785
• /
• 2003

A numerical model for the solution of two-dimensional free surface flow is developed on unstructured grid. By using fractional step method, the two-dimensional shallow water equations (SWE) are treated as two one-dimensional problems. Thus, it is possible to simulate computational hydraulic problems with higher computational efficiency. The one-dimensional problems are solved using upwind TVD version of second-order Weighted Averaged Flux (WAF) scheme with HLLC approximate Riemann solver. The numerical oscillations which are common with second-order numerical scheme are controlled by exploiting WAF flux limiter, Some idealized test problems are solved using this model and very accurate and stable solutions are obtained. It can be concluded as an efficient implement for the computation of SWE including dam break problems that concerning discontinuities, subcritical and supercritical flows and complex domain.

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

The tidal bores of the Qiantang River on the East coast of China are simulated numerically based on the shallow water theory. The governing equations, which were traditionally formulated using water depth, are formulated in terms of water surface level, and the fractional-step method is applied in conjunction with a Godunov-type scheme. In addition, the source terms due to bottom gradient are discretized centrally to exactly balance the flux terms. Our numerical simulation produces tidal bores in excellent agreement with field measurements.

• Proceedings of the Korea Contents Association Conference
• /
• 2012.05a
• /
• pp.143-144
• /
• 2012

전통적인 VOF 기법을 이용한 내부 조파 방법을 레블셋 기법에 적용하였다. 기하학적으로 유리한 유한요소법을 이용하여, Navier-Stokes 방정식의 공간차분에는 Characteristic Galerkin 기법을, 시간차분에는 Fractional Four-step 기법을 적용하였다. 중심(x=0)에서 전파하는 경우, 외부조파에 의한 영역내 재반사 문제가 해결되어 선형파를 의도한 바대로 잘 조파할 수 있었다.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.704-708
• /
• 2000

A new efficient numerical method for computing unsteady, incompressible flows, DRANS (Decoupled Reynolds-Averaged Navier-Stokes), is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used fer both the diffusion and convection terms. is adopted. Based on decomposition method, the velocity-turbulent quantity decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully-implicit time advancement scheme. Since the iterative procedures for the momentum, ${\kappa}\;and\;{\varepsilon}$ equations are not required, the components decouplings bring fourth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to turbulent boundary layer with local forcing.

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

The effect of the Dirichlet boundary condition for the redistance equation of level set method on the solutionof sloshing problem is investigated by adopting four Dirichlet boundary conditions. For the solution of the incompressible Navier-Stokes equations, P1P1 four-step fractional finite element method is employed and a least-square finite element method is used for the solutions of the two hyperbolic type equations of level set method; advection and redistance equation. ALE (Arbitrary Lagrangian Eulerian) method is used to deal with a moving computational domain. It has been shown that the free surface motion in a sloshing tank is strongly dependent on the type of the Dirichlet boundary condition and the results of broken dam and sloshing problems using various Dirichlet boundary conditions are discussed and compared with the existing experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.15-18
• /
• 2004

A numerical method for simulating two-phase flows including surface force is presented. The method is based on fractional step method of finite volume formulation and the interface is tracked with PLIC VOF method. In the CSF model, as color function, f, representing the location of interface varies steeply in the interface region, we need to use smoothed function f to get accurate unit normal and the curvature. Peskin kernel is used to get smoothed function f. A spherical drop in static equilibrium and three-dimensional merging of gas bubble are tested, resulting in the validation of this method

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.910-915
• /
• 2003

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

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

We present an immersed boundary (IB) method for 3D simulation of flappingflags in a uniform flow. The proposed formulation is manipulated on the basis of an efficient Navier-Stokes solver adopting the fractional step method and a staggered Cartesian grid system. A direct numerical method is developed to calculate the flag motion, with the elastic force treated implicitly. The fluid motion defined on an Eulerian grid and the flag motion defined on a Lagrangian grid are independently solved and the mass of flag is handled in a natural way. An additional momentum forcing is formulated from the flag motion equation in a way similar with the direct-forcing IB formulation and acts as the interaction force between the flag and ambient fluid. A series of numerical tests are performed and the present results are compared qualitatively and quantitatively with previous studies. The instantaneous flag motion is analyzed under different conditions and surrounding vortical structures are identified. The effects of physical parameters on the flapping frequency are studied.

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

The Immersed boundary method(IBM) is one of CFD techniques which can simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In these cases, pressure buildup near IB was found to occur when linear interpolation and stadard mass conservation is used and the interpolation scheme became complicated when higher order of interpolation is adopted. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. Bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2005-2011
• /
• 2003

Introducing the interfacial pressure jump terms based on the surface tension into the momentum equations of two-phase two-fluid model, the system of governing equations is turned mathematically into the hyperbolic system. The eigenvalues of the equation system become always real representing the void wave and the pressure wave propagation speeds as shown in the previous manuscript. To solve the interfacial pressure jump terms with void fraction gradients implicitly, the conventional semi-implicit method should be modified as an intermediate iteration method for void fraction at fractional time step. This advanced semi-implicit method (ASIM) then becomes stable without conventional additive terms. As a consequence, including the interfacial pressure jump terms with the advanced semi-implicit method, the numerical solutions of typical two-phase problems can be more stable and sound than those calculated exclusively by using any other terms like virtual mass, or artificial viscosity.