• 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 computational fluids engineering
• /
• v.5 no.1
• /
• pp.33-42
• /
• 2000

The LiBr-H₂O absorption process on a horizontal tube has been analyzed numerically. The flow field, which was calculated in the authors' previous study by solving the fully elliptic Navier-Stokes equations with accurate free-surface-tracking method, is used to solve the temperature and concentration distributions in the absorption film. With the assumption that the absorbent is linear, calculations have been made for various inlet temperature and flow-rate conditions. For low inlet temperature, the absorption rate is large in the upstream region but the mean temperature also increases and as a result the absorption decreases as the film flows to downstream while high-inlet-temperature case does the opposite. The difference in the absorption rate due to the inlet temperature change becomes smaller in the downstream than that in the upstream. For large flow rate, the heat transfer to the wall becomes poor due to the thick film and so does the absorption rate. The analyses have also been carried out for multiple tube arrangement and the results show that the absorption rate converges after a few tube rows.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.27-33
• /
• 2009

In this paper, the flow past a rotary oscillating circular cylinder is simulated. The high-order and high-resolution numerical schemes with the characteristic boundary conditions are used for the compressible Navier-Stokes equation. The frequencies of rotating oscillation are $0.19\;{\leq}\;S_f\;{\leq}\;0.25$ for the maximum angular $\theta_{max}=10^{\circ}$ and $17^{\circ}$. The flow conditions are Mach number of 0.3 and Reynolds number of 1000. At Lock-on and Non-lock-on region which are defined by the relation between the vortex shedding frequency and the oscillating frequency, the drag and lift coefficient are analyzed.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.90-95
• /
• 2011

The Navier-Stokes equations are solved to study the flow characteristics of a micro viscous pump. The viscous micropump is consisted of a stationary disk with a spiral shaped channel and a rotating disk. A simple geometrical model for the tip clearance is proposed and validated by comparing computed flow rate with corresponding experimental data. Present numerical solutions show satisfactory agreement with the corresponding experimental data. The tip clearance effect is found to become significant as the rotational speed increases. As the pressure load increases, a reversed flow region is seen to form near the stationary disk. The height of the channel is shown to be optimized in terms of the flow rate for a given rotational speed and pressure load. The optimal height of the channel becomes small as the rotational speed decreases or the pressure load increases. The flow rate of the pump is found to be in proportion to the width of channel.

• Journal of the Korean Mathematical Society
• /
• v.38 no.1
• /
• pp.1-23
• /
• 2001

This paper is concerned with an optimal shape control problem for the stationary Navier-Stokes system. A two-dimensional channel flow of an incompressible, viscous fluid is examined to determine the shape of a bump on a part of the boundary that minimizes the viscous drag. by introducing an artificial compressibility term to relax the incompressibility constraints, we take the penalty method. The existence of optima solutions for the penalized problem will be shown. Next, by employing Lagrange multipliers method and the material derivatives, we derive the shape gradient for the minimization problem of the shape functional which represents the viscous drag.

• 유체기계공업학회:학술대회논문집
• /
• 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 computational fluids engineering
• /
• v.1 no.1
• /
• pp.112-122
• /
• 1996

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.35-43
• /
• 1997

An efficient inverse method for 1.he supersonic/hypersonic axisymmetric body design is developed for the parabolized Navier-Stokes equations. The developed method is examined numerically for three extreme testcases in the supersonic(M/sub ∞/=3.0) and hypersonic(M/sub ∞/=6.28) speeds. The first one is a negative pressure distribution near a vacuum pressure and the second one is a positive pressure distribution over the whole region of the body. The last one is the case of abrupt change of pressure distribution to zero in the forward region of the body. These testcases show the robustness of the method. By introducing a regular-falsi method and by using a not-fully converged inverse solution, the convergence behavior was greatly improved.