• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.999-1011
• /
• 1989

A consideration on the trubulence models for describing the redeveloping turbulent boundary layer beyond separation-reattachment in the flow over a backward-facing step is given through experimental and numerical studies. By considering the blance among the measured values of respective terms in the transport equations for the turbulent kinetic energy and the turbulent shear stress, the recovering process of the redeveloping boundary layer from non-equilibrium to equilibrium has been investigated, which takes place slowly over a substantial distance in the downstream direction. In the numerical study, the standard K-.epsilon. model and the Reynolds stress model have been applied to two kinds of flow regions, one for the entire downstream region after the backward-facing step and another for the downstream region after reattachment. Then the results are compared to a meaningful extent, with the experimental values of the turbulent kinetic energy k, the turbulent energy production term P, the dissipation term K-.epsilon. model, a necessity for a new modelling has been brought forward, which can be also applied to the case of the nonequlibrium turbulent flow.

• Journal of the Korean Society for Precision Engineering
• /
• v.6 no.4
• /
• pp.43-51
• /
• 1989

Free jet was investigated experimentally and numerically in range of Reynolds number from 9900 to 21000. The working fluid was air; the mean velocity components and turbulent quantities were measured by a hot-wire anemometer. In numerical computations, the governing partial differential equations of elliptic type were solved with conventional k- ${\epsilon}$ turbulence model. The measurements show that the jet increased linearly in flow direction, and that similarity for each turbulent quantity such as Reynolds shear stress, or turbulent kinetic energy was revealed in the fully developed region. The computational results show good agreements with experiments.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.1
• /
• pp.184-192
• /
• 1994

The parameters such as Richardson numbers or stability parameters are widely used to account for the extra straining effects due to three-dimensionality, curvature, rotation, swirl and others arising in paractical complex flows. Existing expressions for the extra strain in turbulence models such as $k-{\epsilon}$ models, however, do not satisfy the tensor invariant condition representing the coordinate indifference. In the present paper, considering the characteristics of both the mean strain rate and the mean vorticity, a new parameter to deal with the extra straining effects is proposed. The new parameter has a simple form and satisfies the tensor invariant condition. A semi-quantitative analysis between the present and previous parameters for several typical complex flows suggests that the newly proposed parameter is more general and adequate in representing the extra straining effects than the previous ad-hoc parameters.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2315-2328
• /
• 1993

The turbulence effect of particle deposition on a horizontal free-standing wafer in a vertical flow has been studied numerically by using the low-Reynolds-number k-.epsilon. turbulence model. For both the upper and lower surfaces of the wafer, predictions are made of the averaged particle deposition velocity and its radial distribution. Thus, it is now possible to obtain local information about the particle deposition on a free-standing wafer. The present result indicates that the particle deposition velocity on the lower surface of wafer is comparable to that on the upper one in the diffusion controlled deposition region in which the particle sizes are smaller than $0.1{\mu}m$. And it is found in this region that, compared to the laminar flow case, the averaged deposition velocity under the turbulent flow is about two times higher, and also that the local deposition velocity at the center of wafer is high equivalent to that the wafer edge.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.9
• /
• pp.1744-1759
• /
• 1992

A multidimensional simulation of turbulent flow and combustion with swirl in the cylinder of SI engine is implemented to clarify the effects of swirl. present simulation employs the ICED-ALE numerical technique, the skew-upwind difference scheme, a modified k_.epsilon. turbulence model, a combustion model of Arrhenius type and turbulence-mixing-control type. First, the calculations for swirling flow in an axisymmetric cylinder are carried out. The results are compared with the experimental data to validate the numerical analysis. Second, the calculations for intake, compression and combustion processes in an axisymmetric cylinder are performed. The effects of swirl on turbulent flow and combustion are examined through the parametric study of swirl number 0.0, 0.6, 1.2 and 2.4. As a result, it is numerically shown that the turbulent kinetic energy and the swirl velocity, which are produced during the intake process, affect the combustion process.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.2
• /
• pp.287-298
• /
• 1989

A numerical work was performed to study the flow behaviors of the open channel type flow with its geometric boundary conditions being similar to that of the Multi-Stage-Flash evaporator with and without a baffle. For the analysis, two-dimensional steady turbulent flow was assumed and the widely known k-.epsilon. turbulence model was usded. SIMPLE algorithm and the power difference scheme were used for the numerical approach. Numerical results generally agree with the previous experimental results though there are some uncertainties at far downstream and near the free surface due to the three dimensionality of the flow and surface waves. Without a baffle, the flow has basically the shape of the submerged plane wall jet with its upper boundary at downstream being sharply curved toward the free surface. For the flow with a baffle, recirculation flow patterns are observed at the upper inlet portion and at the backside of the baffle. For the case without a baffle, it was also confirmed that the ratio between the liquid level and the gate opening height is the most important parameter to determine the flow behavior.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.2 s.32
• /
• pp.138-146
• /
• 1999

The viscous flow around a ship hull is calculated by the use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stresses are midelled by using the k-${epsilon}$ turbulence model and the law is applied near the body. Body fitted corrdinates are introduced for the treatment of the complex boundary of the ship hull form and the governing equations in the physical domain transformed into ones in the computational domain. The transformed equations are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implicit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the sidcretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). To assure the proprietty of this computing method, HSVA tanker and Dyne hull are calculated ar both model and ship scale Reynolds number. Their reaults of pressure distributions on fore and aft body, axial velocity contours and transverse velocity velocity vectors and viscous resistance coefficients are compared with other's experiments and calculations.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.2
• /
• pp.467-474
• /
• 1993

Viscous flow and heat transfer phenomena in an axisymmetric cylinder which models a diesel engine have been numerically studied. In order to search for a way to minimize numerical diffusion, the effectiveness and the appropriateness of two selected numerical schemes for convective terms in the governing equations have been tested. They are Linear Upwind Difference Scheme and Hybrid Scheme. Using a standard k-.epsilon. turbulence model, the calculation has been carried out basically up to 180.deg. of crank angle. As a result, it was shown from comparison with previous experimental data that Linear Upwind Difference Scheme is less influenced than Hybrid Scheme by the numerical diffusion and it was suggested that these effects of numerical diffusion can be more significant than those due to turbulence modeling.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.81-86
• /
• 2012

The objective of the present study is to investigate the effects of various turbulence models on the aerodynamic noise of an air-conditioner (AC) indoor unit. The results from URANS (unsteady Reynolds-averaged Navier-Stokes) simulations with the standard k-$\varepsilon$, k-$\omega$ shear stress transport (SST) and Spalart-Allmaras (S-A) turbulence models were analyzed and compared with the noise data from the experiments. The frequency spectra of the far-field acoustic pressure were computed using the Farrasat equation derived from the Ffowcs Williams-Hawkings (FW-H) equation based on the acoustic analogy model. Two fixed fan casings and the rotating cross-flow fan were used as the source surfaces of the dipole noise in the Farrasat equation. The result with the standard k-$\epsilon$ model showed a much better agreement with the experimental data compared to the k-w SST and S-A models. The differences in the pressure spectra from the different turbulence models were discussed based on the instantaneous vorticity fields. It was found that the over-estimated power spectra with the k-w SST and S-A models are related to the emphasized small-scale vortices produced with these models.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.28 no.1
• /
• pp.1-6
• /
• 2016

This paper studies the wind pressure distribution over the Commonwealth Advisory Aeronautical Council building model (CAARC model) using CFD. We also considered the interaction between the CAARC model and other buildings. The Reynolds number based on the building height was 380,000. The number of sells for the simulation was about 500,000. The wind pressure was lowest when the wind direction was blowing at an angle 45 degrees of the CAARC model. When the gap between the two buildings in front of the CAARC was over 1/2 the horizontal length of the CAARC model, the wind pressure was higher than the pressure without the two buildings. When the distance between the two front buildings and the CAARC was less than 1.5 times the vertical length of the CAARC model, the wind pressure increased. Accordingly, the relative distance between two buildings or the distance from the CAARC model should be considered when extra wind exists due to other buildings.