• 한국전산유체공학회지
• /
• 제4권3호
• /
• pp.12-20
• /
• 1999

The computations of the flowfield and pollutant dispersion over a flat plate and the Russian hills of various slopes are described. The Gaussian plume and the puff model have been used to calculate concentration of pollutant. The Reynolds-averaged unsteady incompressible Navier-Stokes equation with low Reynolds κ-ε model has been used to calculate the flowfield. The flow data of a flat plate and the Russian hills from Navier-Stokes equation solutions has been used as the input data for the puff model. The computational results of flowfield agree well with experimental results of both a flat plate and Russian hills. The concentration prediction by the Gaussian plume model and the Gaussian puff model also agrees flirty well with experiments.

• 한국전산유체공학회지
• /
• 제13권1호
• /
• pp.28-34
• /
• 2008

A particle method recognized as one of the gridless methods has been developed to investigate the nonlinear free-surface motions interacting to the structures. The method is more feasible and effective than convectional grid methods for solving the non-linear free-surface motion with complicated boundary shapes. The right-handed side of the governing equations for incompressible fluid, which includes gradient, viscous and external force terms, can be replaced by the particle interaction models. In the present study, the developed method is applied to the dam-broken problem on dried- and wet-floor and its adequacy will be discussed by the comparison with the experimental results.

• 대한기계학회논문집A
• /
• 제22권2호
• /
• pp.365-371
• /
• 1998

In the present study, we attempted to add the incompatible functions as additional variable terms to the conventional u-p formulation. It is derived from the four-field generalized variational principle that encompasses velocity, pressure, velocity strains and stress fields as independent interpolated variables. As a severe test of the present formulation, we have investigated the driven cavity with the corner velocity singularity like leaky lid. Through the test, the present element performs very well without unstable oscillation of pressure distribution.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제15권3호
• /
• pp.233-251
• /
• 2011

The Gauge-Uzawa method [GUM] in [9] which is a projection type algorithm to solve evolution Navier-Stokes equations has many advantages and superior performance. But this method has been studied for backward Euler time discrete scheme which is the first order technique, because the classical second order GUM requests rather strong stability condition. Recently, the second order time discrete GUM was modified to be unconditionally stable and estimated errors in [12]. In this paper, we contemplate several GUMs which can be derived by the same manner within [12], and we dig out properties of them for both stability and accuracy. In addition, we evaluate an stability condition for the classical GUM to construct an adaptive GUM for time to make free from strong stability condition of the classical GUM.

• 한국유체기계학회 논문집
• /
• 제2권2호
• /
• pp.57-63
• /
• 1999

The flow through multistage turbomachinery is affected by the interaction between a rotor and a stator. The interaction is due to the inviscid potential effect and viscous effect between closely spaced rotor and stator airfoils. Three-dimensional, unsteady, incompressible Navier-Stokes equations with a standard $k-{\epsilon}$ model are solved using a non-staggered grid system. This method is applied to the flow through a multistage compressor measured by Stauter et al. The results have shown strong interaction between the rotating and stationary flow field. The decay of rotor wake and the pressure profiles agree very well with experimental data. The wake produced by rotor causes unsteady pressure on the surface of a stator. The rotor/stator interaction produces the unsteady pressure force on the rotor and stator blades.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제14권1호
• /
• pp.43-56
• /
• 2010

The representative numerical algorithms to solve the time dependent Navier-Stokes equations are projection type methods. Lots of projection schemes have been developed to find more accurate solutions. But most of projection methods [4, 11] suffer from inconsistency and requesting unknown datum. E and Liu in [5] constructed the gauge method which splits the velocity $u=a+{\nabla}{\phi}$ to make consistent and to replace requesting of the unknown values to known datum of non-physical variables a and ${\phi}$. The errors are evaluated in [9]. But gauge method is not still obvious to find out suitable combination of discrete finite element spaces and to compute boundary derivative of the gauge variable ${\phi}$. In this paper, we define 4 gauge algorithms via combining both 2 decomposition operators and 2 boundary conditions. And we derive variational derivative on boundary and analyze numerical results of 4 gauge algorithms in various discrete spaces combinations to search right discrete space relation.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• pp.71-76
• /
• 2004

An extension of our recently developed locally linear reconstruction scheme to 2 dimensional incompressible flow solver is presented. The solver is based on a semi-implicit fractional step method in which the convective term is discretized by Adams-Bashforth method and the diffusion term by Crank-Nicolson method. Several numerical examples are tested to demonstrate the mesh type independent accuracy of the solver, which include decaying vortex flow, square cavity flow, and flow around a circular cylinder. The above examples are solved on quadrilateral or hybrid meshes. For all numerical examples, we obtained reasonable results.

• 한국전산유체공학회지
• /
• 제9권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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 추계 학술대회논문집
• /
• pp.169-177
• /
• 2005

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 추계 학술대회논문집
• /
• pp.179-182
• /
• 2005

We have developed a numerical scheme to reproduce the unsteady flows with cavitation by the finite-difference method. The evolution of cavitation is represented by the source/sink of vapor phase in the incompressible liquid flow. The pressure-velocity coupling is based on the fractional-step method for incompressible fluid flows, in which the compressibility is taken into account through the low Mach number assumption. We applied our method for the cavitating flows in a two-dimensional cascade, which approximates the portion near the tip of inducer in liquid-fuel engine. Particular attention was focused on the influence of turbulence model in this report. Using an eddy viscosity model, although it was not an optimized one for our purpose, the agreement with the experimental observation was improved.