• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.200-202
• /
• 2003

In this paper, Gauss filter function is used to filter the N-S equation and the subgrid-scale Reynold stresses model is introduced to deduce the practical form of LES equation for 2-D case for flow calculation of hydraulic machine. Then the LES equation and its discrete form in computational field are obtained in the body-fitted coordinate system and the numerical calculation program is built. The incompressible turbulent flow in double-flow-conduits-sewage pump impeller is computed by using the abovementioned program, and then the distribution rules of velocity and pressure in flow field are obtained. Based on this, the designs of double-flow-conduits-sewage pump impeller are optimized.

• Journal of Mechanical Science and Technology
• /
• 제15권8호
• /
• pp.1188-1195
• /
• 2001

A numerical analysis has been carried out to investigate the influences of thermosolutal convection on the heat and mass transfer and solute segregation in crystals grown by the vertical Bridgman technique. The governing equations are solved by a finite-volume method using the power law scheme and the SIMPLE algorithm in which body-fitted coordinate system has been used. A primary convective cell driven by thermal gradients forms in the bulk of the domain, while a secondary convective cell driven by solutal gradients forms near interface. As the solutal Rayleigh number increases, secondary cell becomes to be stronger and has a great influence on the radial concentration along the interface.

• Journal of Mechanical Science and Technology
• /
• 제14권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년도 춘계 학술대회논문집
• /
• 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.

• 대한기계학회논문집B
• /
• 제22권12호
• /
• pp.1796-1804
• /
• 1998

A numerical investigation has been performed to discuss the radiation-affected steady-laminar natural convection induced by a hot inner cylinder under a large temperature difference in the annuli filled with a gray gas. To examine the effects of thermal radiation on thermo-fluid dynamic behaviors in the eccentric geometry, the generalized body-fitted coordinate system is introduced while the finite volume method (FVM) is used for solving the radiative transport equation. After validating the numerical results for the case without radiation, the detailed radiation effect has been discussed. Based on the results of this study, when there exists a large temperature difference between two cylinders, the existence of radiatively participating medium is found to incur a distinct difference in fluid dynamic as well as thermal behavior.

• 한국전산유체공학회지
• /
• 제8권1호
• /
• pp.8-15
• /
• 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 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.

• 한국연소학회지
• /
• 제5권1호
• /
• pp.55-66
• /
• 2000

A general purpose program for the analysis of flows in a gas turbine combustor is developed. The program uses non-staggered grids based on finite volume method and the cartesian velocities as primitive variables. A flow inside the C-type diffuser is simulated to check the boundary fitted coordinate. The velocity profiles at cross section agree well with experimental results. A turbulent diffusion flame behind a bluff body is simulated for the combustion simulation. Simulated results show good agreement with experimental data. Finally, a turbulent flow with swirl in a gas turbine combustor was simulated. The results show two recirculating region and simulated velocity fields agree well with experimental data. The distance between two recirculating regions becomes shorter as swirl angle increases. Swirl angle changes angular momentum and streamlines in flow fields.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2002년도 춘계 학술대회논문집
• /
• pp.102-108
• /
• 2002

Numerical analysis has been conducted in order to simulate the flow characteristics by tip clearance of turbine cascades. A 3-D Navier-Stokes CFD code based on body-fitted coordinate system, pressure-correction and finite volume method has been used along with a commercial CFD code. The present results have showed that the development and generation of leakage vortex, vortex within tip clearance, etc. are clearly simulated, consistent with the generally known tendency. The leakage vortex occurs mainly by a separation of leakage flow that arises due to a pressure difference between two surfaces of the blade at the tip.

• Journal of Advanced Marine Engineering and Technology
• /
• 제33권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.

• 한국자동차공학회논문집
• /
• 제2권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.