• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• 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 the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.104-114
• /
• 1993

In this paper, effects of free surface on viscous flow is investigated. Continuity equation coupled with Navier-Stokes equations are solved numerically by using an artificial compressibility method[1, 2]. The body-fitted generalized curvilinear coordinate system is employed to deal with arbitrary body shape. The IAF scheme with finite difference method is used to solve the equations, and a diagonal algorithm is applied to time-varying Jacobian matrices for the computational economics. Free surface shape is obtained by applying zero pressure condition to still water surface at each time step. A numerical test is made for larminar flow around a circular cylinder vertically piercing the free surface. Computed flow patterns are largely affected by the existance of free surface in low Reynolds number flows treated in this paper. Free surface causes viscous pressure drag to vary much in depth direction in accordance with the variations of flow pattern.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.149-154
• /
• 1997

Based on the full Navier-Stokes solutions, the thermohydrodynamic performance of a long journal bearing is investigated. A numerical method based on Galerkin's procedure and B-spline test functions has been presented for solving two-dimensional problems involving fluid flow and heat transfer. For numerical stability the artificial compressibility is employed to the conservation of mass. The discretized algebraic equations are solved by Newton's method. Effects of varying the speed of an inner cylinder to load carrying capacity are investigated. The results indicated that the increase of the speed of an inner cylinder has a significant effect on the temperature profile and ultimately on the performance.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.33-39
• /
• 1996

An incompressible flow stover based on the unstructured finite volume method has been developed. The flow domain is discretized by triangles in 2D or tetrahedra in 3D. The convective and viscous fluxes are obtained using edge connectivities of the unstructured meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. Laminar test flow problems are computed and presented with a comparison against other numerical solutions or experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.73-80
• /
• 1997

Objective of the present study is to examine the convergence characteristics of the various multi-stage Runge-Kutta methods in solving the incompressible Navier-Stokes equations of a time-marching from casted by the artificial compressibility method. Convergence characteristics are examined over 2-stage, 4-stage and hybrid type (using 4-, 3-, 2-stages sequentially) Runge-Kutta methods for a laminar lid-driven cavity flow, and also for a turbulent bump channel flow using Chien's low-Reynolds number turbulence model. Efforts are made to establish a stable and fast convergent multi-stage Runge-Kutta method with minimal artificial dissipations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.49 no.6
• /
• pp.449-456
• /
• 2021

In this paper, we present a strategy to extend solution capability of an existing low Mach number preconditioned compressible solver to incompressible flows with a little modification. To this end, the energy equation that is of the same form of the total energy equation of compressible flows is used. The energy equation is obtained by a linear combination of the thermal energy equation, the continuity equation and the mechanical energy equation. Subsequently, a modified artificial compressibility method in conjunction with a time marching technique is applied to these incompressible governing equations for steady flow solutions. It is found that the Roe average of the common governing equations is equally valid for both the compressible and incompressible flow conditions. The extension of an existing compressible solver to incompressible flows does not affect the original compressible flow analysis. Validity for incompressible flow analysis of the extended solver is examined for various inviscid, laminar and turbulent flows.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.48-54
• /
• 1998

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.