• 한국전산유체공학회:학술대회논문집
• /
• 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 computational fluids engineering
• /
• v.9 no.3
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have 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 and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.78-83
• /
• 2000

An iterative time marching procedure for solving incompressible internal flow has been applied to the flow around a rotor-stator. This procedure solves three-dimensional incompressible Reynolds-averaged Navier-Stokes equation on a moving, time-deforming, non-orthogonal body-fitted grid using second-order accurate schemes for the time derivatives and third/second-order accurate schemes for the spatial derivatives. To handle rotationg geometry, the multiblock technique is applied and the overall flow domain is subdivided into two blocks. In each block, a grid is generated and flowfield is solved independently of the other blocks. The boundary data for each block is provided by the neighboring blocks using bilinear interpolation technique.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.3 s.258
• /
• pp.300-305
• /
• 2007

Many practical flow problems are defined with the circular boundary. Fluid flows within a circular boundary are however susceptible to a singularity problem when the cylindrical coordinates are employed. To remove this singularity a method has been developed in this study which uses the local harmonic functions in discretization of derivatives as well as interpolation. This paper describes the basic reason for introducing the harmonic functions and the overall numerical methods. The numerical methods are evaluated in terms of the accuracy and the stability. The Lamb-dipole flow is selected as a test flow. We will see that the harmonic-function method indeed gives more accurate solutions than the conventional methods in which the polynomial functions are utilized.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.35-42
• /
• 2007

This paper describes a recent development on the divergence free basis function based on a hermite stream function and verifies its validity by comparing results with those from a modified residual method known as one of stabilized finite element methods. It can be shown that a proper choice of degrees of freedom at a node with a proper arrangement of the hermite interpolation functions can yield solenoidal or divergent free interpolation functions for the velocities. The well-known cavity problem has been chosen for validity of the present algorithm. The comparisons from numerical results between the present and the modified residual showed the present method yields better results in both the velocity and the pressure within modest Reynolds numbers(Re = 1,000).

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.41-49
• /
• 2003

The spatial accuracy becomes first-order when second-order conservation schemes including the non-conservative interpolation in general Chimera method are used. To ensure the solution accuracy, the discontinuities must be located away from the overlapped regions, and the length of overlapped region also must be proportional to the grid spacing. In this paper, a proposed method, cut-paste algorithm, is used to satisfy above constraints. The cut-paste algorithm can generate the optimal mesh inteface region automatically, To validate the spatial and temporal accuracy due to the non-conservative interpolation, inviscid and viscous problems are tested.

• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.41-48
• /
• 2008

The fluid-structure interaction analysis such as a static aeroelastic analysis requires the result of each analysis as an input to the other analysis. Usually the grids for the fluid analysis and the structural analysis are different, so the results should be transformed properly for each other. The Infinite Plate Spline(IPS) and the Thin Plate Spline(TPS) are used in interpolating the displacement and the pressure. In this study, such interpolation methods are compared with kriging which provides a precise response surface. The static aeroelastic analysis is performed for the supersonic flow field with shock waves and the pressure field is interpolated by the TPS and kriging. The TPS shows tendency to weaken the shock strength, whereas kriging preserves the shock strength.

• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.49-56
• /
• 2008

Chimera grid methods have been widely used in Computational Fluid Dynamics due to its simplicity in constructing grid systems over complex bodies, and suitability for unsteady flow computations with bodies in relative motion. However, the interpolation procedure for ensuring the continuity of the solution over overlapped regions fails when the so-called orphan cells are present. We have adopted the MLS(Moving Least Squares) method to replace commonly used linear interpolations in order to alleviate the difficulty associated with the orphan cells. MLS is one of the interpolation methods used in mesh-less methods. A number of examples with MLS are presented to show the validity and the accuracy of the method.