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

To improve the efficiency of internal combustion engines, it is necessary to understand mixed air-fuel in-cylinder flow processes accurately at intake and compression strokes. There is experimental and numerical methods to analyse in-cylinder flow process. In numerical method, standard $k-{\varepsilon}$ model with wall function was mostly adopted in in-cylinder flow process. But this type model was not efficiently predicted in the near wall region. Therefore in the present study, low Reynolds number $k-{\varepsilon}$ model was adopted near the cylinder wall and standard $k-{\varepsilon}$ model in other region. Also QUICK scheme was used for convective difference scheme. This study takes axisymmetric reciprocating model engine motored at 200rpm with a centrally located valve, incorporated 60 degree seat angie, and flat piston surface excluding inlet port. Because in-cylinder flow processes are undergoing unsteady and compressible, averaged cylinder pressure and inlet velocity at arbitrary crank angle are determined from thermodynamic analytic method and incylinder states at that crank angle are iteratively determined from the numerical analytic method.

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

In this study, flow of a model intake port/valve system is analyzed by using low Reynolds number $k-{\varepsilon}$ model. Discharge coefficient was obtained from computational results for the various cases of valve lifts. Discharge coefficient becomes maximum when the valve lift is 20mm, and does not increase or decrease in proportional to valve lift. Most of pressure drop and production of turbulent kinetic energy occur at the edge points of the valve and the valve seat Thus, in order to improve discharge coefficient, rounding of edge points in valve and valve seat is recommended. As valve lift is increased, the velocity of the intake jet in the valve passage decreases, and the direction of the jet is more inclined toward the valve seat.

• Journal of Mechanical Science and Technology
• /
• v.17 no.7
• /
• pp.1054-1062
• /
• 2003

The numerical simulations of flow and pollutant particle dispersion are described for two-dimensional bell shaped hills with various aspect ratios. The Reynolds-averaged incompressible Navier-Stokes equations with low Reynolds number $\kappa$-$\varepsilon$ turbulent model are used to simulate the flowfield. The gradient diffusion equation is used to solve the pollutant dispersion field. The code was validated by comparison of velocity, turbulent kinetic energy, Reynolds shear stress, speed-up ratio, and ground level concentration with experimental and numerical data. Good agreement has been achieved and it has been found that the pollutant dispersion pattern and ground level concentration have been strongly influenced by the hill shape and aspect ratio, as well as the location and height of the source.

• 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.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.37-42
• /
• 1995

In the present study, an improved version of 4-equation low-Reynolds-number 4-equation model is proposed. The equations of the temperature variance ($k_{\theta}$) and its dissipation rate(${\varepsilon}_{\theta}$) are solved, in concert with the equations of the turbulent kinetic energy(k) and its dissipation rate(${\varepsilon}$). In the present model, the near-wall effect and the non-equilibrium effect are fully taken into consideration. The validation of the model is then applied to the turbulent flow behind a backward-facing step and the flow over a blunt body. The predicted results of the present model are compared and evaluated with the relevant experiments.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.2
• /
• pp.287-297
• /
• 1995

In this study, the geometry consists of a two-dimensional rectangular enclosure with localized heating from below. The size and the location of the heater on the floor has been varied, and one of the vertical walls remains at a low temperature simulating a cold window. The governing equations for momentum, energy and continuity, which are coupled with turbulent equations have been solved using a finite volume method. A low Reynolds number $k-{\varepsilon}$ model has been incorporated to solve the turbulent kinetic energy and the dissipation rate. The heat transfer characteristics and the thermal environmental characteristics of the room have been obtained for various system parameters in a room with a partially heated floor.

• 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.

• Nuclear Engineering and Technology
• /
• v.27 no.4
• /
• pp.491-502
• /
• 1995

In accurate prediction of cross-flow based on detailed knowledge of the velocity field in subchannels of a nuclear fuel assembly is of importance in nuclear fuel performance analysis. In this study, the low-Reynolds number k-$\varepsilon$ turbulence model has been adopted in too adjacent subchannels with cross-flow. The secondary flow is accurately estimated by the anisotropic algebraic Reynolds stress model. This model was numerically calculated by the finite element method and has been verified successfully through comparison with existing experimental data. Finally, with the numerical analysis of the velocity Held in such subchannel domain, an analytical correlation of the lateral loss coefficient is obtained to predict the cross-flow rate in subchannel analysis codes. The correlation is expressed as a function of the ratio of the lateral How velocity to the donor subchannel axial velocity, recipient channel Reynolds number and pitch-to-diameter.

• International Journal of Aeronautical and Space Sciences
• /
• v.1 no.1
• /
• pp.36-47
• /
• 2000

Three-dimensional compressible turbulent flow fields within the passage of a diffusing S-duct have been simulated by solving the Navier-Stokes equations with SIMPLE scheme. The average inlet Mach number is 0.6 and the Reynolds number based on the inlet diameter is $1.76{\times}10^6$ The extended $k-{\varepsilon}$ turbulence model is applied to modeling the Reynolds stresses. Computed results of the flow in a circular diffusing S-duct provide an understanding of the flow structure within a typical engine inlet system. These are compared with experimental wall static-pressure, total-pressure fields, and secondary velocity profiles. Additionally, boundary layer thickness, skin friction values, and streamlines in the symmetric plane are presented. The computed results depict the interaction between the low energy flow by the flow separation and the high energy flow by the reversed duct curvature. The computed results obtained using the extended $k-{\varepsilon}$ turbulence model.