• Journal of Theoretical and Applied Mechanics
• /
• v.2 no.1
• /
• pp.15-29
• /
• 1996

A fair portion of the dream to acquire the solutions to the Navier-Stokes equations has come true through the remarkable development of computers and solution algorithms in recent years. However, it is also true that there still remain serious hurdles in simulating general fluid flows. A few numerical trials to overcome the existing difficulties are introduced. The issues in numerical simulations of high-Reynolds-number flows, flows characterized by complex body geometry, and multi-phase flows, are scrutinized. The future of computational fluid dynamics as a promising tool for flow analyses is illuminated by this review.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.19-26
• /
• 1991

The complete, fully-elliptic Reynolds-averaged Navier-Stokes equations have been solved using a two-layer model, in the $\kappa-\varepsilon$ turbulence model, for the axisymmetric body. Numerically generated boundary-fitted coordinate system and the finite analytic methods are used to solve the governing equations. Calculations are started after the middle body with given inlet conditions. The velocities and the turbulent quantities at the inlet section are specified by solving the boundary layer equations or by standard flat-plate boundary profiles. The effects of the inlet conditions on the solution are investigated.

• Journal of Korean Society for Atmospheric Environment
• /
• v.9 no.E
• /
• pp.347-357
• /
• 1993

The Lagrangian grid-free numerical method, the discrete vortex method, was applied to solve the Navier-Stokes euqations. This method avoids the introduction of numerical viscosity swamping the real physical viscosity at high Reynolds number, unlike Eulerian method, e.g. finite difference and element methods. The boundary integral equation method for the potential flow solution was included to make the discrete vortex method more feasible for complex geometries. The fast adaptive multipole expansion method was incorporated to reduce the computational time from $O(N^2)$ to O(N) for the computations of vortex-vortex interactions. The test problems were air flow around one circular cylinder and two circular cylinders in tandem with various gaps. The numerical results were in excellent gareement with the experimental and other computational results. The applicabilty of the method was discussed with the indoor and the outdoor air pollution problems, especially the contaminant transport in the recirculation regions.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.1
• /
• pp.63-72
• /
• 1995

Convective heat transfer in a two-dimensional horizontal and vertical channel with isothermal rectangular beams attached to one adiabatic wall is investigated from the numerical solution of Navier-Stokes and energy equations. The solutions have been obtained for dimensionless beam spacings, S/L=1~4, aspect ratios of beam, H/B=0.25~4, Reynolds numbers, Re=50~1000 and Grashof numbers, $Gr=0{\sim}5{\times}10^4$. The total mean Nusselt number, Nu_T for horizontal and vertical channels shows same value at Gr=0. As Gr increases, Nu_T for horizontal channel increases, but Nu_T for vertical channel shows similar value at S/L=2, H/B=0.25, Re=100. The total mean Nusselt number for horizontal channel is higher than that for vertical channel. As H/B increases, $Nu_T$ for both channel decrease at $Gr=10^4$, Re=100.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.1
• /
• pp.23-30
• /
• 2004

In the Lagrangian vortex particle method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations, a numerical scheme for calculating pressure fields is presented. Implementation of the numerical method is directly connected with the well-established surface panel methods, just by dealing with the dynamic coupling among vorticity field. Assuming the vorticity and the velocity fields are to be calculated in time domain analysis, the pressure calculation for a complete set of solution at present time step is performed in a similar way to the one used in the Eulerian description. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsive started circular cylinder for Reynolds number 550. The comparative study with the Eulerian finite Volume method provides an extensive understanding and application of the mesh-free Lagrangian vortex methods for numerical simulation of viscous flows around arbitrary bodies of general shape.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.66-72
• /
• 1999

A numerical simulation method has been under development for solving turbulent flows around a ship model in maneuvering motion using the Reynolds Averaged Navier-Stokes equations. The method used second-order finite differences, collocated grids, pressure-Poisson equation and four-stage Runge-Kutta scheme as key components of the solution method. A modified Baldwin-Lomax model is used for the turbulence closure. This paper presents a preliminary result of the computational study on turbulent flows past a ship model in drift motion. Calculations are carried out for a Series 60 $C_B=0.6$ ship model, for which detailed experimental data are available. The results of the present calculations are compared with the experimental data for hydrodynamic forces acting on the model as well as velocity distributions at longitudinal sections. Only fair agreements has been achieved. The computational results show the complex asymmetrical shear flow patterns including three-dimensional separations followed by formation of bilge vortices both in bow and stern regions.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.125-130
• /
• 1999

A flowfields around a NACA0012 airfoil pitching about a 1/4 chord and plunging in vertical displacement are analyzed by solving two-dimensional compressible Navier-Stokes equations. A steady solution was solved first as a validation of the code used and the results were compared with experimental data. Then as a unsteady case, the oscillatory airfoil was solved to compare the results with experimental data. Oscillating rate of pitching and plunging motion was set to have analogy and the magnitude of plunging was set using the magnitude of pitching angle of attack. Finally combined pitching and plunging motion was solved to show the effect of 2 different types of oscillating motion of the airfoil.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.74-84
• /
• 2007

The flowfield characteristics in a straight rectangular channel have been investigated through a numerical model to analyze the regenerative cooling system that is used in rocket engine cooling. The supercritical hydrogen coolant introduces strong property variations that have a major influence on the developing flow and heat transfer characteristics. Of particular interest is the improved understanding of the physical characteristics of such flows through parametric studies. The approach used is a numerical solution of the full Navier-Stokes equations in the three dimensional form including the arbitrary equation of state and property variations. The present study compares constant and variable property solutions for both laminar and turbulent flow. For laminar flow, the variation of aspect ratio is examined, while for turbulent flow, the effects of variation of channel length and Reynolds number are discussed.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.336-340
• /
• 2017

When supersonic flow is through an internal duct, there forms a flow structure called pseudo-shock. Pseudo-shock is a result of shockwave-boundary layer interaction(SBLI) and to simulate pseudo-shock correctly, one needs to correctly anticipate not only the strength of the shock but also the boundary layer behavior as well. In this study, pseud-shockwave structure at a rectangular duct will be numerically simulated using dedicated inlet boundary conditions to obtain accurate solution in terms of its structure and pressure rise pattern.