• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.575-581
• /
• 1988

To integrate the two-dimensional Navier-Stokes equations numerically in multiply-connected flow regions, the vorticity-stream function formulation is used. The steady stream function value at the surface of the multibody, initially unknown, has been determined interactively by introducing a line integral which requires the single-valuedness of pressure at each interaction step. This procedure is relatively simpler and more efficient than the primitive variable formulation which requires much more computing time and shows poor convergence. Three doubly-connected flow problems are defined and numerically analyzed by the present method. The results have been compared either with earlier existing ones or with the experimental interferograms to demon-strate the validity of the presented method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.15 no.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.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.37-45
• /
• 1998

A new and efficient inverse design method based on the numerical optimization technique has been developed. The 2-D incompressible Navier-Stokes equations are solved for obtaining the objective functions and coupled with the optimization procedure to perform the inverse design. The steepest descent and the conjugate gradient method have been applied to find the searching direction. The golden section method was applied to compute the design variable intervals. It has been found that the airfoil and the pump impellers are well converged to their targeting shapes.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.200-207
• /
• 1999

An efficient inverse design technique based on the MGM (Modified Garabedian-McFadden) method has been developed. The 2-D Navier-Stokes equations are solved for obtaining the surface pressure distributions and coupled with the MGM method to perform the inverse design. The solver is parallelized by using the domain decomposition method and the standard MPI library for communications between the processors. The MGM method is a residual-correction technique, in which the residuals are the difference between the desired and the computed pressure distribution. The developed code was applied to several airfoil shapes and the axial blade. It has been found that they are well converged to their target pressure distribution.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.8-13
• /
• 2008

We present a numerical simulation technique and some preliminary results of the impact and spreading of a droplet containing particles on the solid substrate in 2D. We used the 2nd-order Adams-Bashforth / Crank-Nicholson method to solve the Navier-Stokes equation and employed the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension. The impact velocity has been generated by the instantaneous gravity. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles and the discontinuous Galerkin method has been used for the stabilization of the interface advection equation. We investigated the droplet spreading by the inertial force and discussed effects of the presence of particles on the spreading behavior using an example problem. We observed reduced oscillation and spread for the particulate droplet.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.3
• /
• pp.9-16
• /
• 1987

In this paper, the fluid flow in the two-dimensional tank is analyzed by the Finite Difference Method. The Navier-Stokes equation is modified for the tank fixed coordinate system. For the treatment of the free surface, the Volume of Fluid Method by Hirt and Nichols is adopted. The continuity equation and the Poisson equation which is derived from the Navier-Stokes equation to find the pressure are solved by the Successive-Line-Overrelaxation Method. The comparison of the calculated results with experimental data show a favorable agreement. The fluid flow in the two-dimensional tank can be predicted reasonably before the free surface reaches breaking by this numerical method.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.152-157
• /
• 2003

As computer capacity has been progressed continuously, the studies of the flow characteristics have been performing by the numerical methods actively. Recent numerical simulation has a tendency to require the higher-order accuracy in time, as well as in space. This tendency is more true in LES and acoustic noise simulation. In this study, 3-dimensional unsteady Incompressible Navier-Stokes equation was solved by numerical method using the fractional step method with the fourth order compact pade scheme to achieve high accuracy To validate the present code and algorithm, 3D flow-field around a cylinder was simulated. The drag coefficient and lift coefficient were computed and, then, compared with experiment. The present code will be tailored to LES simulation for more accurate turbulent flow analysis.

• International Journal of Ocean System Engineering
• /
• v.1 no.1
• /
• pp.37-45
• /
• 2011

The nonlinear free-surface motions interacting with a floating body were investigated using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka and Oka [6] for incompressible flow. In the numerical method, more realistic Lagrangian moving particles were used for solving the flow field instead of the Eulerian approach with a grid system. Therefore, the convection terms and time derivatives in the Navier-Stokes equation can be calculated more directly, without any numerical diffusion, instabilities, or topological failure. The MPS method was applied to a numerical simulation of predicting the efficiency of floating-type breakwater interacting with waves.

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.661-669
• /
• 2007

Eyer since the Prandtl's experiment in 1934 and X-21 airjet test in 1950 both attempting to reduce drag, it was found that controlling the velocities of surface for extremely fast-moving object in the air through suction or injection was highly effective and active method. To obtain the right amount of suction or injection, however, repetitive trial-and error parameter test has been still used up to now. This study started from an attempt to decide optimal amount of suction and injection of incompressible Navier-Stokes by employing optimization techniques. However, optimization with traditional methods are very limited, especially when Reynolds number gets high and many unexpected variables emerges. In earlier study, we have proposed an algorithm to solve this problem by using step by step method in analysis and introducing SQP method in optimization. In this study, we propose more effective and robust algorithm and techniques in solving flow optimization problem.