• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.777-784
• /
• 2000

A Navier-Stokes code based on an unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number ${\kappa}-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for the whole experimental domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. The first harmonics of the velocity in the boundary layer shows local peak value inside the boundary layer and also local minimum near the edge of boundary layer. It is intensified as it develops along the blade surface. This is shown to be caused as the unsteadiness inside the boundary layer is being convected at a speed less than the free stream value. It is also shown that there is negligible mixing of the unsteadiness between the boundary layer and the free stream.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.1
• /
• pp.87-95
• /
• 2001

A numerical study was made of heat transfer in locally-forced turbulent separated and reattaching flow over a backward-facing step. The local forcing was given to the flow by means of sinusoidally oscillating jet from a separation line. A Rhee and Sung version of the unsteady $\kappa$-$\varepsilon$-f(sub)u model and the diffusivity tensor heat transfer model were employed. The Reynolds number was fixed at Re(sub)H=33,000 and the forcing frequency was varied in the range 0$\leq$fH/U(sub)$\infty$$\leq$2. The condition of constant heat flux was imposed at the bottom wall. The predicted results were compared and validated with the experimental data of Chun and Sung and Vogel and Eaton. The enhancement of heat transfer in turbulent separated and reattaching flow by local forcing was evaluated and analyzed.

• Environmental Engineering Research
• /
• v.11 no.3
• /
• pp.164-171
• /
• 2006

Numerical modeling of the flow velocity fields for the near corona wire electrohydrodynamic (EHD) flow was conducted. The steady, two-dimensional momentum equations have been computed for a wire-plate type electrostatic precipitator (ESP). The equations were solved in the conservative finite-difference form on a fine uniform rectilinear grid of sufficient resolution to accurately capture the momentum boundary layers. The numerical procedure for the differential equations was used by SIMPLEST algorithm. The Phoenics (Version 3.5.1) CFD code, coupled with Poisson's electric field, ion transport equations and the momentum equation with electric body force were used for the numerical simulation and the Chen-Kim ${\kappa}-{\varepsilon}$ turbulent model numerical results that an EHD secondary flow was clearly visible in the downstream regions of the corona wire despite the low Reynolds number for the electrode ($Re_{cw}=12.4$). Secondary flow vortices caused by the EHD increases with increasing discharge current or EHD number, hence pressure drop of ESP increases.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.6 no.3
• /
• pp.267-281
• /
• 1994

The Combustion gas behavior induced by fire in a building is numerically investigated. The typical building for this analysis is partially divided by a vertical baffle projecting from the ceiling. The solution procedure includes the low Reynolds number ${\kappa}-{\varepsilon}$ model for the turbulent flow and the discrete ordinates method is used for the calculation of radiative heat transfer equation. The effects of the location and size of fire source and baffle length on velocity and temperature distributions, species mass fraction and flame location are analyzed. As the results of this study, it is found that the case when the fire source is located at the vertical wall is more dangerous than at the bottom wall in view of the combustion products and flame location. It is also found that the radiation effect cannot be neglected in analyzing the building in fire.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.3
• /
• pp.31-39
• /
• 1999

In this paper, numerical experiment is attempted to analyze and compare the combustion efficiency of the burning sprays due to OFO, FOF triplet / FOOF split doublet injectors. Preconditioned Wavier-Stokes equation system with low Reynolds number $\kappa$-$\varepsilon$ model for turbulence closure, is LU-SGS time-integrated. Spray processes are modeled by DSF analysis with experimentally determined injection characteristics. n-heptane/air global reaction model approximates the combustion for simplicity, and the influence of turbulence on the chemical reaction is included using eddy dissipation model. The results showed the FOF triplet injector of highest combustion efficiency, whereas the OFO type of poet performance. It was also observed that the droplet mean diameter and the average gas temperature due to the mixing efficiency, are the representative parameters for the performance design of combustion.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.1
• /
• pp.21-29
• /
• 2002

Hybrid propulsion systems provide many advantages in terms of stable operation and safety. However, classical hybrid rocket motors have lower fuel regression rate and combustion efficiency compared to solid propellant rocket motor. The recent research efforts are focused on the improvement of volume limitation and regression rate in the hybrid rocket engine. The present study has numerically investigated the combustion processes in the hybrid rocket engine. The turbulent combustion is represented by the eddy breakup model and Hiroyasu and Nagle and Strickland-Constable model are used for soot formation and soot oxidation. Radiative heat transfer is modeled by finite volume method. To reduce the uncertainties for convective heat transfer near solid fuel surface having strong blowing effect, the Low Reynolds number $\kappa-\varepsilon$ turbulent model is employed. Based on numerical results, the detailed discussion has been made for the turbulent combustion processes in the vortex hybrid rocket engine.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.2
• /
• pp.1-9
• /
• 1999

Applicability of the pressure gradient method which is formulated based on pressure gradient is verified against turbulent flow analysis. In the pressure gradient method, pressure gradient instead of pressure itself is obtained using continuity constraint. Since correct pressure gradient is found only when mass conservation is satisfied, pressure gradient method can reflect physics of flow field properly The pressure gradient method is formulated with semi-staggered grid system which locates each primitive variables on the same grid point but evaluates pressure gradient in-between. This grid system ensures easy programming and reflection of correct physics in analysis. For verifying applicability of this method, the pressure gradient method is applied to turbulent flow analysis with low Reynolds number $\kappa$-$\varepsilon$ model. Turbulent flows include fully developed channel flow, backward-facing step flow, and conical diffuser flow. Prediction results show that the pressure gradient method can be applied to turbulent flow analysis. However, the pressure gradient method requires somewhat long computation time. Proper way to find optimum under-relaxation factor, $\gamma$, is also need to be developed.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.6 no.1
• /
• pp.1-10
• /
• 1996

A numerical study was conducted on the effects of the cusp magnetic field in 8" silicon single crystal grwoth by Czochralski method. For a damping effects simulation by magnetic field, low reynolds number ${\kappa} - {\varepsilon}$ model was adopted. Symmetrci cusp magnetic field has a effect of damping streamline crystal, is lowerd with the increasing cusp magnetic field intensity. The uniformity of the oxygen concentration was improved. The asymmetirc cusp magnetic field increased the oxygen concentration however, oxygen concentration distribution in the radial direction was remained uniform. Suitable combination of symmetric and asymmetric cusp magnetic fields could give uniform and low oxygen concentration in the axial direction.tion.