• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.20-27
• /
• 2007

Two-dimensional stagnation flow toward a plane wall coated with magnetic fluid of uniform thickness is investigated. The flow field is represented as a similarity solution of the Navier-Stokes equation for this incompressible laminar flow. The resulting third order ordinary differential equation is solved numerically by using the shooting method and by determining two shooting parameters so as to satisfy the boundary and interface conditions. Features of the flow including streamline patterns are investigated for the varying values of density ratio, viscosity ratio, and Reynolds number. An adverse flow with double eddy pair in magnetic fluid region is found to emerge as the Reynolds number becomes higher than a threshold value. The results for the interface velocity, interface and wall shear stress, and boundary layer and displacement thickness are also presented.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.216-221
• /
• 2004

In two dimensional incompressible flaws, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint and hence achieve an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flaw analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flaw showed the present procedure to be stable, very efficient and useful in flaw analyses in conjunction with hierarchical elements.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.5 s.54
• /
• pp.11-18
• /
• 2003

In two-dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure (HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint hence, achieving an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared (CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper and efficient iterative procedure for higher order finite element formulation has not been available, thus far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient, and useful in flow analyses, in conjunction with hierarchical elements.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Journal of Ocean Engineering and Technology
• /
• v.14 no.3
• /
• pp.35-40
• /
• 2000

SMAC method is, one of the numerical simulation techniques, modified from the original MAC for the time-dependent variation of fluid flows. The Navier-Stokes equation for incompressible time-dependent viscous flow is applied and, also marker particles which move with the fluid are used. Two-dimensional numerical computations of fluid flow are carried out in a virtual water tank. This paper has shown very well the movements of marker particles using SMAC method.

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

Simulations of atmospheric diffusion process under stable and unstable conditions were carried out using both numerical and experimental methods. Results from the previous study show that numerical simulation using 3-dimensional incompressible Navier-Stokes equation and density deviation are in good agreement with typical plume pattern. In this study, we use experimental data of temperature and wind profile obtained from a thermally stratified wind tunnel as initial conditions for numerical simulation and compare the results.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.139-143
• /
• 2002

The purpose of this paper is to develop the CFD code to consider the viscous flow features of the marine propeller. The flow of the marine propeller has been numerically analyzed by using three dimensional viscous incompressible Navier-Stokes equation. The model used in this study is Screw B with 4 blades whose pitch ratio is 1 in Ka-4-55 screw series. The result of the analysis was compared with panel method.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.410-415
• /
• 2001

An analytical and numerical examination of second-order fractional-step methods and boundary condition for the incompressible Navier-Stokes equations is presented. In this study, the compatibility condition for pressure Poisson equation and its boundary conditions, stability, and numerical accuracy of canonical fractional-step methods has been investigated. It has been found that satisfaction of compatibility condition depends on tentative velocity and pressure boundary condition, and that the compatible boundary conditions for type D method and approximately compatible boundary conditions for type P method are proper for divergence-free velocity for type D and approximately divergence-free for type P method. Instability of canonical fractional-step methods is induced by approximation of implicit viscous term with explicit terms, and the stability criteria have been founded with simple model problems and numerical experiments of cavity flow and Taylor vortex flow. The numerical accuracy of canonical fractional-step methods with its consistent boundary conditions shows second-order accuracy except $D_{MM}$ condition, which make approximately first-order accuracy due to weak coupling of boundary conditions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.44-52
• /
• 2002

Aerodynamic characteristics of a flapping airfoil in low Reynolds number flows are numerically studied using the unsteady, incompressible Navier-Stokes flow solver with a two-equation turbulence model. For more efficient computation of unsteady flows over flapping airfoil, the flow solver is parallel-implemented by MPI programming method Unsteady computations are performed for low Reynolds number flows over a NACA four-digit series airfoils. Effects of pitching, plunging, and flapping motion with different reduced frequency, amplitude, thickness and camber on aerodynamic characteristics are investigated. Present computational results yield a better agreement in thrust at various reduced frequency with experimental data.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.