• Journal of Mechanical Science and Technology
• /
• v.14 no.4
• /
• pp.448-455
• /
• 2000

The bow wave breaking and the viscous interaction of stern wave are studied by simulating the free-surface flows. The Navier-Stokes equation is solved by a finite difference method in which the body-fitted coordinate system, the wall function and the triple-grid system are invoked. After validation, the calculations are extended to turbulent flows. The wave elevation at the Reynolds number of $10^4$ is much less than that at $10^6$ although the Froude number is the same. The numerical appearance of the sub-breaking waves is qualitatively supported by experimental observation. They are also applied to study the stern flow of S-103 for which extensive experimental data are available. Although the interaction between separation and the stern wave generation are not yet clear, the effects of the bow wave on the development of the boundary layer flows are concluded to be significant.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.27-32
• /
• 1996

A previously developed pressure based calculation procedure for incompressible flows was modified and applied to transonic and supersonic flows. It uses pressure as a primary variable in preference to density and body fitted coordinate and non-staggered grid system. The discretized momentum equations were rearranged as a system of equations with respect to covariant velocity components. Three different discretization schemes, QUICK hybrid and first order upwind, were used to approximate the convective terms and compared. Present approach was tested far two transonic flow and one supersonic flow problems. Comparison with previous results show that present approach can be used as a solver for compressible flows.

• Journal of Advanced Marine Engineering and Technology
• /
• v.33 no.6
• /
• pp.873-879
• /
• 2009

The steady breakers at an infant stage are investigated through the numerical simulation. The appearing condition and characteristics of the sub-breaking waves are reviewed by analysing bow waves. The instability analysis is possibly done through the relationship between the free-surface curvature and circumferential force, which is obtained from the momentum equations. Navier-Stokes equations are solved by a finite difference method where the body-fitted coordinate system, the wall function and the advanced mesh system are invoked. The numerical result shows that the gradient of M/$U_s$ is greatly influenced by the Froude number and the decrease of M/$U_s$ indicates that the flows are unstable. Additionally flows with plunging or spilling are simulated successfully, but the application of breakers to the severely broken wave still remains to be settled in the future.

• Journal of Mechanical Science and Technology
• /
• v.18 no.7
• /
• pp.1203-1212
• /
• 2004

The purpose of this study is to predict and understand the hydrodynamic forces and their nonlinear behaviors of fluid motion around the submerged plate oscillating near a free surface. To achieve this objective, we have developed a composite grid method for the solution of a radiation problem. The domain is divided into two different grids; one is a moving grid system and the other is a fixed grid system. The moving grid is employed for the body fitted coordinate system and moves with the body. This numerical method is applied to calculation of radiation forces generated by the submerged plate oscillating near a free surface. In order to investigate the characteristics of the radiation forces, the forced heaving tests have been performed with several amplitudes and different submergences near a free surface. These experimental results are compared with the numerical ones obtained by the present method and a linear potential theory. As a result, we can confirm the accuracy of the present method. Finally, the effect of nonlinear and viscous damping has been evaluated on the hydrodynamic forces acting on the submerged plate.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.474-481
• /
• 2001

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

• Journal of Navigation and Port Research
• /
• v.28 no.8
• /
• pp.753-757
• /
• 2004

The free-surface flow is simulated to make clear the viscous interaction of stem waves and the sub-breaking phenomena around a high speed vehicle. The Navier-Stokes equation is solved by a finite difference method where the body-fitted coordinate system, the wall function and the triple-grid system are invoked They are applied to study precisely on the stem flow of S-103 as to which extensive experimental data are available. Computations are extended to the submerged revolutional body. The numerical result shows that the gradient of M/Us is greatly influenced by the submerged depth And the stem wave is influenced by the separation due to the bow wave.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.9-25
• /
• 1994

In the present study, the Reynolds-averaged Navier-Stokes equations, together with the equations of the $k-{\varepsilon}$ model of turbulence, were solved numerically in a general body-fitted coordinate system for three-dimensional turbulent flows around the six basic shapes of the magnetically levitated train(MAGLEV). The numerical computations were conducted on the MAGLEV model configurations to provide information on shapes of this type very near the elevated track at a constant Reynolds number of $1.48{\times}10^{6}$ based on the body length. The coordinate system was generated by numerically solving a set of Poisson equations. The convective transport equations were discretized using the finite-analytic scheme which employed analytic solutions of the locally-linearized equations. A time marching algorithm was employed to enable future extensions to be made to handle unsteady and fully-elliptic problems. The pressure-velocity coupling was treated with the SIMPLER-algorithm. Of particular interests were wall effect by the elevated track on the aerodynamic forces and flow characteristics of the six models calculated. The results indicated that the half-circle configuration with extended sides and with smooth curvature of sides was desirable because of the low aerodynamic forces and pitching moment. And it was found that the separation bubble was occured at wake region in near the elevated track.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.49-56
• /
• 2003

Three-dimensional endwall flow within a linear cascade passage of high performance turbine blade is simulated with a 3-D Navier-Stokes CFD code (MOSA3D), which is based on body-fitted coordinate system, pressure-correction and finite volume method. The endwall flow characteristics, including the development and generation of horseshoe vortex, passage vortex, etc. are clearly simulated, consistent with the generally known tendency. The effects of both turbulence model and convective differencing scheme on the prediction performance of endwall flow are systematically analyzed in the present paper. The convective scheme is found to have stronger effect than the turbulence model on the prediction performance of endwall flow. The present simulation result also indicates that the suction leg of the horseshoe vortex continues on the suction side until it reaches the trailing edge.

• The KSFM Journal of Fluid Machinery
• /
• v.1 no.1 s.1
• /
• pp.72-80
• /
• 1998

Computational analysis of incompressible flows by numerically solving Navier-Stokes equations using multi-block finite volume method is conducted on a parallel computing system. Numerical algorithms adopted in this study $include^{(1)}$ QUICK upwinding scheme for convective $terms,^{(2)}$ central differencing for other terms $and^{(3)}$ the second-order Euler differencing for time-marching procedure. Structured grids are used on the body-fitted coordinate with multi-block concept which uses overlaid grids on the block-interfacing boundaries. Computational code is parallelized on the MPI environment. Numerical accuracy of the computational method is verified by solving a benchmark test case of the flow inside two-dimensional rectangular cavity. Computation in the axial compressor cascade is conducted by using 4 PE's md, as results, no numerical instabilities are observed and it is expected that the present computational method can be applied to the turbomachinery flow problems without major difficulties.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.1-6
• /
• 2003

Three-dimensional endwall flow within a linear cascade passage of high performance turbine blade is simulated with a 3-D Wavier-Stokes CFD code (MOSA3D), which is based on body-fitted coordinate system, pressure-correction and finite volume method. the endwall flow characteristics, including the development and generation of horseshoe vortex, passage vortex, etc. are clearly simulated, consistent with the generally known tendency The effects of both turbulence model and convective differencing scheme on the Prediction performance of endwall flow are systematically analyzed in the present paper. The convective scheme is found to have stronger effect than the turbulence modei on the prediction performance of endwall flow. The present simulation result also indicates that the suction leg of the horseshoe vortex continues on the suction side until it reaches the trailing edge.