• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.3
• /
• pp.158-165
• /
• 2005

In order to analyze incompressible two-phase flow, Level Set method was applied on pseudo-compressibility formulation. Level Set function is defined as a signed distance function from the phase interface, and gives the information of the each phase location and the geometric data to the flow. In this study, Level Set function transport equation was coupled with flow conservation equations, and owing to pseudo-compressibility technique we could solve the resultant vector equation iteratively. Two-phase flow analysis code was developed on general curvilinear coordinate, and numerical tests of bubble dynamics and surging wave problems demonstrate its capability successfully.

• Bulletin of the Korean Mathematical Society
• /
• v.55 no.3
• /
• pp.881-898
• /
• 2018

This article deals with implementation of a high-order finite difference scheme for numerical solution of the incompressible Navier-Stokes equations on curvilinear grids. The numerical scheme is based on pseudo-compressibility approach. A fifth-order upwind compact scheme is used to approximate the inviscid fluxes while the discretization of metric and viscous terms is accomplished using sixth-order central compact scheme. An implicit Euler method is used for discretization of the pseudo-time derivative to obtain the steady-state solution. The resulting block tridiagonal matrix system is solved by approximate factorization based alternating direction implicit scheme (AF-ADI) which consists of an alternate sweep in each direction for every pseudo-time step. The convergence and efficiency of the method are evaluated by solving some 2D benchmark problems. Finally, computed results are compared with numerical results in the literature and a good agreement is observed.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.147-157
• /
• 2006

This paper represents a numerical study of the flow field due to the interactions between a pair of vortices produced by vortex generators in a rectangular channel flow. In order to analyze longitudinal vortices induced by the vortex generators, the pseudo-compressibility method is introduced into the Reynolds-averaged Navier-Strokes equations of a 3-dimensional unsteady, incompressible viscous flow. A two-layer $k-{\epsilon}$ turbulence model is applied to a flat plate 3-dimensional turbulence boundary to predict the flow structure and turbulence characteristics of the vortices. The computational results predict accurately the vortex characteristics related to the flow field, the Reynolds shear stresses and turbulent kinetic energy. Also, in the prediction of skin friction characteristics the computational results are reasonably close to those of the experiment obtained from other researchers.

• The KSFM Journal of Fluid Machinery
• /
• v.6 no.1 s.18
• /
• pp.44-50
• /
• 2003

Flow through a turbine flow meter is simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo-compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved steadily in rotating reference frames, and the centrifugal force and the Coriolis force are added to the equation of motion. The standard $k-{\epsilon}$model is employed to evaluate turbulent viscosity. Computational results yield quantitative as well as qualitative information on the design of turbine flow meters by showing the distributions of pressure and velocity around the turbine blades.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1624-1629
• /
• 2008

The governing equations in generalized curvilinear coordinates for a 3D pulsatile flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. The computational technique implements the implicit approximate factorization method of the Beam and Warming method (1978), which is the extension of the Alternate Direction Implicit (ADI) method. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Rindt & Steenhoven, 1991).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.802-813
• /
• 2000

This paper represents numerical computations of the interaction between the longitudinal vortex and a flat plate 3-D turbulent boundary layer. In the present study, the main interest is in the behavior of longitudinal vortices introduced in turbulent boundary layers. The flow field behind vortex generator is modeled by the information that is available from studies on the delta winglet. Also, the Reynolds-averaged Navier-Stoke equations for three-dimensional turbulent flows, together with a two-layer turbulence model to resolve the near-wall flow, is solved by the method of pseudo compressibility. The present results show that the boundary layer is thinned in the regions where the secondary flow is directed toward the wall and thickened where it is directed away from the wall, and have a good agreement with the experimental data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1392-1396
• /
• 2007

The governing equations in generalized curvilinear coordinates for 3D laminar flow are the Incompressible Navier-Stokes (INS) equations with the artificial dissipative terms. and continuity equation discretized using a second-order accurate, finite volume method on the nonstaggered computational grid. This method adopts a dual or pseudo time-stepping Artificial Compressibility (AC) method integrated in pseudo-time. Multigrid methods are also applied because solving the equations on the coarse grids requires much less computational effort per iteration than on the fine grid. The algorithm yields practically identical velocity profiles and secondary flows that are in excellent overall agreement with an experimental measurement (Humphrey et al., 1977).

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.17 no.9
• /
• pp.799-807
• /
• 2005

The flow characteristics and the heat transfer rate on a surface by the interaction of a pair of vortices are studied numerically. To analyze the common flow up produced by vortex generators in a rectangular channel flow, the pseudo-compressibility viscous method is introduced into the Reynolds-averaged Navier-Stokes equation for 3-dimensional unsteady, incompressible viscous flows. To predict turbulence characteristics, a two-layer $k-\varepsilon$ turbulence model is used on the flat plate 3-dimensional turbulence boundary The computational results predict accurately Reynolds stress, turbulent kinetic energy and flow field generated by the vortex generators. The numerical results, such as thermal boundary layers, skin friction characteristics and heat transfers, are also reasonably close to the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.83-84
• /
• 2006

• Journal of Advanced Marine Engineering and Technology
• /
• v.24 no.1
• /
• pp.31-40
• /
• 2000

This paper describes the numerical computations of the interaction between the streamwise vortex and a flat plate 3-D turbulent boundary layer. In the present study, the main interest is in the behavior of the streamwise vortices introduced in turbulent boundary layers. The flow behind a vortex generator is modeled by the information that is avilable from studies on the dalta winglet. An algorithm of the solution of the incompressible Navier-Stokes equations for three-dimensional turbulent flows, together with a two layer turbulent model to resolve the near-wall flow, is based on the method of artificial compressibility. The present results show boundary layer distortion due to vortices, such as strong spanwise flow divergence and boundary thinning, and have a good agreement with the experimental data.