• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1596-1608
• /
• 1993

Low Reynolds number second moment turbulence model which be applicable to the fine gird near the wall region was developed. In this model, turbulence model coefficients in the pressure strain model of the Reynolds stress equation was expressed as functions of turbulence Reynolds number $R_{t}\equivk^{2}/(\nu\varepsilon)).$ In the derivation procedure of the present low Reynolds number algebraic stress model, Laufer's near wall experimental data on Reynolds stresses were curve fitted as functions of R$_{t}$ and the resulting simultaneous equations of the model coefficients were solved by using the boundary conditions at wall and high Reynolds number limiting conditions. Predicted Reynolds stresses and dissipation rate of turbulent kinetic energy etc. in the 2 dimensional parallel, plane channel flow and pipe flow were compared with the preditions obtained by employing the Launder-Shima model, standard algebraic stress model and several experimental data. Results show that all the Reynolds stresses and dissipation rate of turbulent kinetic energy predicted by the present low Reynolds number algebraic stress model agree better with the experimental data than those predicted by other algebraic stress models.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.5
• /
• pp.1203-1217
• /
• 1994

In the present study, in order to analyze a turbulent flow in a square sectiond 180.deg. bend, Kim's low Reynolds number second moment turbulence closure is adopted. In this model, turbulence model constants in the wall region are modified as functions of turbulent Reynolds number by use of near wall turbulent universal properties based on Laufer's experimental results of Reynolds stress distriburions. Algebraic stress model and Reynolds stress equation model are used to verify the low Reynolds number second moment closure. The application of the present low Reynolds number algebraic stress model to the prediction of a square sectioned 180.deg. bend flow gives improved velocities and Reynolds stresses profiles compared with those obtained by using the van Driest mixing length model and present low Reynolds number Reynolds stress equation model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2572-2592
• /
• 1996

A low-Reynolds-number second moment turbulence closure was developed with the aid of DNS data. Model coefficients of nonlinear return to isotropy term were derived by use of Cayley-Hamilton theorem and two component turbulence limit condition as the functions of invariances of anisotropy and turbulent Reynolds number. Launder and Tselepidakis' cubic mean pressure strain model was modified to fit the predicted pressure-strain components to the DNS data. Two component turbulence limit condition was the precondition to be satisfied in developing the second moment turbulence closure for the realizable Reynolds stress prediction. But the satisfactions of Reynolds stress level and pressure-strain level of each component were compromised because the satisfaction of both levels was impossible.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.5
• /
• pp.1661-1670
• /
• 1996

Incompressible flow over a backward-facing step is computed by low Reynolds number turbulence models in order to compare with direct simulation results. In this study, selected low Reynolds number 1st and 2nd (Algebraic Stress Model : ASM) moment closure turbulence models are adopted and compared with each other. Each turbulence model predicts different flow characteristics, different re-attachment point, velocity profiles and Reynolds stress distribution etc. Results by .kappa.-.epsilon. turbulence models indicate that predicted re-attachment lengths are shorter than those by standard model. Turbulent intensity and eddy viscosity by low Reynolds number .kappa.-.epsilon. models are still greater than DNS results. The results by algebraic stress model (ASM) are more reasonable than those by .kappa.-.epsilon. models. The convective scheme is QUICK (Quadratic Upstream Interpolation for Convective Kinematics) and SIMPLE algorithm is adopted. Reynolds number based on step height and inlet free stream velocity is 5100.

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

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• 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.10
• /
• pp.1229-1239
• /
• 1999

A computational study on unsteady compressible flows has been performed by adopting a low Reynolds number $k-{\omega}$ turbulence model in conjunction with dual time stepping scheme. An explicit four-stage Runge-Kutta scheme for the Navier-Stokes equations and an approximate factorization scheme for the $k-{\omega}$ turbulence model equations are used. Computational results obtained for blade surface pressure distributions in the process of rotor-stator interaction in a turbine stage are in good agreement with extant experimental data. The effects of the wake from the stator on the boundary-layer transition over the rotor blade surface are discussed by showing that high intensity turbulence of the stator wake induces an early transition.

• Journal of Korean Society for Atmospheric Environment
• /
• v.13 no.1
• /
• pp.65-77
• /
• 1997

A theoretical and numerical investigation on the boundary-layer flow over a two- or three-dimensional hill is presented. The numerical model is based on the finite volume method with boundary-fitted coordinates. The k-$\varepsilon$ turbulence model with modified wall function and the low-Reynolds-number model are employed. The hypothesis of Reynolds number independency for the atmospheric boundary-layer flow over aerodynamically rough terrain is confirmed by the numerical simulation. Comparisons of the mean velocity profiles and surface pressure distributions between the numerical predictions and the wind-tunnel experiments on the flow over a hill show good agreement. The linear theory provides generally good prediction of speed-up characteristics for the gentle-sloped hills. The flow separation occurs in the hill slope of 0.5 and the measured reattachment points are compared with the numerical prediction. It is found that the k- $\varepsilon$ turbulence model is reasonably accurate in predicting the attached flow, while the low- Reynolds-number model is more suitable to simulate the separated flows.ows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.821-827
• /
• 2003

Fully developed turbulent flow in a square duct is numerically predicted with two nonlinear low-Reynolds-number ${\kappa}-{\varepsilon}$ models. Typical predicted quantities such as axial and secondary velocities, turbulent kinetic energy and Reynolds stresses are compared in detail with each other. It is found that the nonlinear low-Reynolds-number ${\kappa}-{\varepsilon}$ model adopted in a commercial code is unable to predict accurately duct flows involving turbulence-driven secondary motion with the prediction level of secondary flows one order less than that of the experiment.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.172-177
• /
• 2004

In the present work, we have interests on the modification of parallel implemented with MPI(Message Passing Interface) programming method, 3-Dimensional, unsteady, incompressible Navier-Stokes equation solver to analyze the low-Reynolds number flow In order to accurate calculation aerodynamic coefficients in low-Reynolds number flow field, we modified the two-equation turbulence model. This paper describes the development and validation of a new two-equation model for the prediction of flow transition. It is based on Mentor's low Reynolds $\kappa-\omega$ model with modifications to include Total Stresses Limitation (TSL) and Separation Transition Trigger (STT)