• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.741-750
• /
• 1995

The turbulent buoyancy-driven flow in 2-dimensional enclosed cavities heated from the vertical side is numerically calculated for both cases of a Rayleigh number of 5*10$^{10}$ for air and 2.5*10$^{10}$ for water. Three different turbulence models are considered : standard k-.epsilon. model of Ozoe and low-Reynolds-number model of Lam and Bremhorst, and another low-Reynolds-number model of Davidson. The results indicate that the use of low-Reynolds number models is recommended for the indoor airflow computation, and the results from Davidson model are reasonably close to the reported experimental data. A sensitivity study shows that the amounts of wall-heat transfer and the velocity profiles with the Lam and Bremhorst model largely depend on the choice of the wall function for .epsilon..

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.911-918
• /
• 1997

The turbulent viscous wake flows behind a single airfoil, two-dimensional stationary blade row and three-dimensional rotating blade row were calculated, and the numerical results were compared with experimental ones. The numerical technique was based on the SIMPLE algorithm using three turbulent closure models, standard k-.epsilon. model(WFM), low Reynolds number k-.epsilon. model(LRN) and Reynolds stress model (RSM). In the case of a single airfoil, WFM, LRN and RSM presented fairly good velocity distributions in the wake compared with experimental data. In the case of the stationary blade row, LRN and RSM presented better results than WFM for wake velocity distribution, and especially LRN showed best results among these three turbulent models. In the case of the rotating blade row, WFM and LRN showed fairly good agreement with experimental data of the three-dimensional velocity component distributions in the range from hub to mid span region. LRN was also superior to WFM in accuracy of prediction for the wake velocity distribution as same with the cases of a airfoil and the stationary blade row.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.6
• /
• pp.733-739
• /
• 1999

A modified low-Re $k-\varepsilon$ model is used for the calculation of drag-reducing turbulent flow by polymer injection in a pipe. With the viscoelastic model, molecular viscosity in the definition of turbulent viscosity is related to elongations viscosity of the solution to account for the effects of drag reduction. Finite volume method is used for the discretization, and power-law scheme is used as a numerical scheme. Computed dimensionless velocity profiles are in good agreements with the experimental data in case of low drag reductions. However, in case of high drag reductions, they deviate largely from the measurements in the central zone of the flow field.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.12
• /
• pp.1681-1690
• /
• 2003

The effect of roughness is a change in the velocity and turbulence distributions near the surface. Turbulence models with surface roughness effect are applied to the fully developed flow in a two-dimensional, rough wall channel. Modified wall function model, low-Reynolds number k-$\varepsilon$ model, and k-$\omega$ model are selected for comparison. In order to make a fair comparison, the calculation results are compared with the experimental data. The modified wall function model and the low-Reynolds number k-$\varepsilon$ model require further refinement, while the k-$\omega$ model of Wilcox performs remarkably well over a wide range of roughness values.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.9
• /
• pp.1106-1112
• /
• 1999

Substantial losses behind axial flow rotor are generated by the wake, various vortices in the hub region and the tip leakage vortex in the tip region. Particularly, the leakage vortex formed near blade tip is one of the main causes of the reduction of performance, generation of noise and aerodynamic vibration in downstream. In this study, the three-dimensional flow fields in an axial flow rotor were calculated with varying tip clearance under various flow rates, and the numerical results were compared with experimental ones. The numerical technique was based on SIMPLE algorithm using standard $k-{\varepsilon}$ model(WFM) and Launder & Sharma's Low Reynolds Number $k-{\varepsilon}$ model(LRN). Through calculations, the effects of tip clearance and attack angle on the 3-dimensional flow fileds behind a rotor and leakage flow/vortex were investigated. The presence of tip leakage vortex, loci of vortex center and its behavior behind the rotor for various tip clearances and attack angles was described well by calculation.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.428-433
• /
• 2000

The characteristics of unsteady heat transfer and boundary layer flow in the SSME turbine rotor passage are investigated with LRN $k-{\varepsilon}$ turbulence model. The unsteady flow and heat transfer in a rotor blade passage as a result of wake/blade interaction is modeled by the inviscid/boundary-layer flow approach. The relevant governing equations are discretized to a system of finite different equations by means of a BTBCS implicit method. These equations have been solved numerically, for the velocity and temperature fields using TDMA method. Heat flux on the blade surface and flow parameters in the rotor passage are calculated with wake interaction. Numerical results show that velocity, pressure, turbulent kinetic energy and heat flux on the blade surface are varied periodically by wake passing.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.1
• /
• pp.61-69
• /
• 2013

Numerical simulations have been performed for assessment of compressibility correction of two-equation turbulent models suitable for large scale separation flows perturbed by a pintle strokes. Two-equation turbulence models, the low Reynolds k-${\varepsilon}$ and the k-${\omega}$ SST models with or without compressibility correction proposed by Wilcox and Sarkar are evaluated. The detail flow structures are observed and static pressures along nozzle wall are compared with experimental results. Mach disk location and pressure recovery profiles in flow separation region are noticeably distinct between turbulent models of k-${\varepsilon}$ and k-${\omega}$ SST. The compressible effect corrections to those models improve resolving of separation flow behaviors. The compressibility corrections to k-${\varepsilon}$ model have provided very comparable results with test data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1971-1978
• /
• 1992

A turbulent flow in a conical diffuser with total divergence angle of 8.deg. was numerically studied. The low Reynolds number k-.epsilon. model(Launder-Sharma model) was adopted to simulate the turbulence. The continuity and time averaged Navier-Stokes equations in a nonorthogonal coordinate system were solved by a finite volume method based on the fully elliptic formulation. The low Reynolds number k-.epsilon. model reasonably simulates the pressure recovery and the mean velocity components. However, there are also considerable discrepancies between predicted and measured shear stress distribution on the wall and turbulent kinetic energy distributions. It is necessary to investigate the flow structure at the entry of the diffuser, numerically as well as experimentally.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.12 no.2
• /
• pp.150-160
• /
• 2000

Turbulent natural convective flow and heat transfer in a square enclosure with horizontal partition are investigated numerically. The enclosure is composed of a lower hot and a upper cold horizontal walls and adiabatic vertical walls. Partitions carried with the upward, downward, and both control plates are attached perpendicularly to the one of the vertical insulated walls, respectively. The low Reynolds number $k-\varepsilon$ model is adopted to calculate the turbulent thermal convection. The governing equations are solved by using the finite element method with Galerkin method. The computations have been carried out by varying the length of partition, the position of control plates, and the Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height for water(Pr=4.95). When the control plates are attached at the edge of partition, the stability of oscillating flow grows wrose with the increase of Rayleigh number and the partition length. The heat transfer rate has been reducer than that of no control plate due to the restraint of control plates with the increase of Rayleigh number.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.20 no.4
• /
• pp.59-67
• /
• 2016

Numerical simulations have been performed to investigate thrust characteristics of a through-type pintle nozzle with or without flow separation at various operating altitudes. The low Reynolds number $k-{\varepsilon}$ with compressibility correction proposed by Sarkar are applied. The detail flow structures are observed and static pressures along nozzle wall are compared with experimental results. The flow separation in the pintle nozzle disappears and jet plume strongly expands as its operating altitude increases. To evaluate the thrust characteristics, the momentum term and pressure term of thrust are analyzed. Thrust and thrust coefficient at altitude 20 km are about 10% more than them at the ground 0km.