• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.44-53
• /
• 2012

Immersed boundary method(IBM) is a numerical scheme proposed to simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies, the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. The weight coefficients of the bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies. For an analysis of moving boundary, we smulated an oscillating circular cylinder with Re=100 and KC(Keulegan-Carpenter) number of 5. The predicted flow fields were compared with experimental data and they also showed good agreements.

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

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

The three-dimensional unsteady compressible Full Navier-Stokes equation solver with sliding multi-block method has been applied to analyze three dimensional characteristics of the flow field and compression wave around the high speed train which Is entering into a tunnel. The numerical scheme of AF + ADI was used to efficiently solve Navier-Stokes equations in the curvilinear coordinate system. The vortex formation around the nose region was found and the generation of compression wave due to the blockage effects was observed ahead of the train in the form of plane wave. The three dimensional characteristics of the flow field compared to the analytic results were discussed in detail. The variation of pressure of tunnel wall surface and velocity profile of the train are identified as the train enters into a tunnel. The changes in aerodynamic forces and streamlines of each specific sections are also discussed and presented.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.423-426
• /
• 2006

An analysis is made of flow and rocket motion during a supersonic separation stage of air-launching rocket(ALR) from the mother plane. Three-dimensional compressible Navier-Stokes equations is numerically solved to analyze the steady/unsteady flow field around the rocket which is being separated from the mother plane configuration(F-4E Phantom). The simulation results clearly demonstrate the effect of shock-expansion wave interaction between the rocket and the mother plane. To predict the behavior of the ALR according to the change of the C.G., three cases of numerical analysis are performed. As a result, a design-guideline of supersonic air-launching rocket for the safe separation is proposed.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.4
• /
• pp.28-36
• /
• 1999

An efficient inverse design method 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 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 propeller. It has been found that they are well converged to their targeting shapes.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.117-124
• /
• 2005

Splitting algorithms of the incompressible Navier-Stokes equations using P1P1/P2P1 finite element formulation are newly proposed. P1P1 formulation allocates velocity and pressure at the same nodes, while P2P1 formulation allocates pressure only at the vertex nodes and velocity at both the vertex and mid nodes. For comparison of the elapsed time and accuracy of the two methods, they have been applied to the well-known benchmark problems. The three cases chosen are the two-dimensional steady and unsteady flows around a fixed cylinder, decaying vortex, and impinging slot jet. It is shown that the proposed P2P1 semi-splitting method performs better than the conventional P1P1 splitting method in terms of both accuracy and computation time.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.13-21
• /
• 2013

Cavitating flow is usually formed on the surface of a high speed underwater object. When a object moves near a free surface at very high speed, the cavity signature becomes one of the major factors to be overcome by sensors of military satellite. The present work was to study the free surface effect on the ventilated cavitation. The governing equations were Navier-Stokes equations based on a homogeneous mixture model. The multiphase flow solver used an implicit preconditioning method in the curvilinear coordinate system. The cavitation model used here was the one first presented by Merkle et al.(2006) and redeveloped by Park & Ha(2009). Computations considered the free surface effects were carried out with a NACA0012 hydrofoil and the corresponding results were compared with the experimental data to have a good agreement. Calculations were then performed considering the ventilated cavitation, including the effect of non-condensable gas under the free surface effects.

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

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.60-66
• /
• 2011

Engineering interests of submerged bodies and turbomachinery has led researchers to study various cavitation models for decades. The governing equations used for the present work are the two-phase Navier-Stokes equations with homogeneous mixture model. The solver employed on implicit dual time preconditioning algorithm in curvilinear coordinates. Three different cavitation models were applied to two axisymmetric cylinders and compared with experiments. It is concluded that the Merkle's new cavitation model has successfully accounted for cavitating flows and well captured the re-entrant jet phenomenon over the 0-caliber cylinder.

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

The numerical procedure has been developed for simulating incompressible viscous flow around a turbine stage with rotor-stator interaction. This study solves 2-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system. The Marker-and-Cell concept is applied to efficiently solve continuity equation. To impose an accurate boundary condition, O-H multiblocked grid system is generated. O-type grid and H-type grid is generated near and outer rotor-stator The cubic-spline interpolation is applied to handle a relative motion of a rotor to the stator. Turbulent flows have been modeled by the Baldwin- Lomax turbulent model. To validate present procedure, the time averaged pressure coefficients around the rotor and stator are compared with experiment and a good agreement obtained.