• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

• 유체기계공업학회:학술대회논문집
• /
• 2001.11a
• /
• pp.334-338
• /
• 2001

In order to investigate the effects of various injection hole shapes on the film cooling of turbine blade, three test models having cylindrical and shaped holes were used. A three-dimensional Navier-Stokes code with standard k-$\epsilon$ model was used to compute the film cooling coefficient on the film cooled turbine blade. Over 330,000 grids were used to compute the flow over the blade. Mainstream Reynolds number based on the cylinder diameter was $7.1{\times}10^4$. The turbulence intensity kept at $5.0\%$ for all inlets. The effect of coolant blowing ratio was studied for various blowing ratios. For each blowing ratios, wall temperatures around the surface of test model were calculated. Temperature was visualized using cartesian cut-cell method to obtain traces of the injected secondary air on the test surface, so we could interpret the film effectiveness as temperature distributions.

• 유체기계공업학회:학술대회논문집
• /
• 2001.11a
• /
• pp.454-459
• /
• 2001

The present study is concerned with the flow patterns induced by other impellers in a rectangular tank Impellers are FBT(Flat blade turbine), PBT(Pitched blade turbine), Shroud turbine, Rushton Turbine, and Helical ribbon turbine. The solution of flows in moving reference frames requires the use of 'moving' cell zone. The moving zone approaches are MRF(Multiple reference frame), which is a steady-state approximation and Sliding method, which is a unsteady-state approximation. Numerical results using two moving zone approaches are compared with experiments by Ranade & Joshi, which have done extensive LDA measurements of the flow generated by a standard six-bladed Rushton turbine in a cylindrical baffled vessel. In this paper we simulated the flow patterns with above mentioned moving zone approaches and impellers. Turbulence model is RNG k-$\epsilon$ model.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.1
• /
• pp.43-53
• /
• 1996

Numerical simulation was performed for the 3-dimensional flow filed of gas and particle phase for cyclone dust collector. FVM(Finite Volume Method) was employed for gas phase. The flow was solved suing the k-.varepsilon. epsilon turbulence model. The particle exit at the bottom of the cone was treated as a solid wall in this model because the gas flow through the effective dust exit is usually insignificant. The major parameters considered in this study was vortex finder diameter, effective dust exit diameterm vortex finder length, inlet type for dimension performance. Particle trajectory calculations were made for three different, particle sizes of 1, 25 and 50 .mu.m. The results obtained from this study give some physical insight of dust particle collection mechanism together with the indication of the collection efficiency. The simulation results were in generally good agreement with empirical knowledge. The application of this kind of computer program looks promising as a potential tool for the design of cyclone and determination of optimum operating condition.

• Fire Science and Engineering
• /
• v.13 no.4
• /
• pp.20-29
• /
• 1999

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP(Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k- epsilon turbulence model with buoyancy term. Compressibility is assumed and the perfect gas law is used. The results of the calculated upper-layer average temperature and smoke layer interface height has shown reasonable agreement compared with the zone models. The zone models used are the CFAST developed at the Building and Fire Research Laboratory NIST U.S.A. and the NBTC one-room of FIRECALC developed at CSIRO, Australia. The smoke layer interface heights that are important in fire safety were not as sensitive as the smoke layer temperature to the nature of ceiling heat flux condition.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.1-11
• /
• 2006

The finite volume based multi-block RANS code, WAVIS developed at MOERI is applied to the numerical self-propulsion test. WAVIS uses the cell-centered finite volume method for discretization of the governing equations. The realizable $k-{\epsilon}$ turbulence model with a wall function is employed for the turbulence closure. The free surface is captured with the two-phase level set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. The propeller forces are obtained using an unsteady lifting surface method based on potential flow theory. The numerical procedure followed the self-propulsion model experiment based on the 1978 ITTC performance prediction method. The self-propulsion point is obtained iteratively through balancing the propeller thrust, the ship hull resistance and towing force that is correction for Reynolds number difference between the model and full scale. The unsteady lifting surface code is also iterated until the propeller induced velocity is converged in order to obtain the propeller force. The self-propulsion characteristics such as thrust deduction, wake fraction, propeller efficiency, and hull efficiency are compared with the experimental data of the practical container ship. The present paper shows that hybrid RANS and potential flow based numerical method is promising to predict the self-propulsion parameters of practical ships as a useful tool for the hull form and propeller design.

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

Turbulent flows around two-dimensional vehicle-like bodies in ground proximity are numerically simulated. The Reynolds averaged Navier-Stokes equations with a k-.epsilon. turbulence model are numercally solved, and a body-fitted coordinate system is used. It is shown that the simulation is acceptable in comparison with limitted data measured in the wind-tunnel. According to numerical simulations, drag coefficients are under-estimated and lift coefficients are over-estimated during the model test in the wind-tunnel if the ground is fixed. Such ground effects are reduced as Reynolds number is increased. Reducing the gap between the vehicle and the ground make drag coefficients smaller and lift coefficients larger. The changes in static pressure distributions on the bottom and the rear surface play dominent roles in determination of the drag and the lift of the body in ground proximity. Drag component less than 10% of the total amount is contributed by skin-frictions. When the slant-angle of the body is reduced, the drag shows its minimum value and the lift shows its maximum value at about 22 degree.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.5
• /
• pp.1319-1332
• /
• 1995

Numerical analysis was performed to characterize the particle deposition behavior on a horizontal free-standing wafer with thermophoretic effect under the turbulent flow field. A low Reynolds number k-.epsilon. turbulence model was used to analyze the turbulent flow field around the wafer, and the temperature field for the calculation of the thermophoretic effect was predicted from the energy equation introducing the eddy diffusivity concept. The deposition mechanisms considered were convection, diffusion, sedimentation, turbulence and thermophoresis. For both the upper and lower surfaces of the wafer, the averaged particle deposition velocities and their radial distributions were calculated and compared with the laminar flow results and available experimental data. It was shown by the calculated averaged particle deposition velocities on the upper surface of the wafer that the deposition-free zone, where the deposition velocite is lower than 10$^{-5}$ cm/s, exists between 0.096 .mu.m and 1.6 .mu.m through the influence of thermophoresis with positive temperature difference of 10 K between the wafer and the ambient air. As for the calsulated local deposition velocities, for small particle sizes d$_{p}$<0.05 .mu.m, the deposition velocity is higher at the center of the wafer than at the wafer edge, whereas for particle size of d$_{p}$ = 2.0 .mu.m the deposition takes place mainly on the inside area of the wafer. Finally, an approximate model for calculating the deposition velocities was recommended and the calculated deposition velocity results were compared with the present numerical solutions, those of Schmidt et al.'s model and the experimental data of Opiolka et al.. It is shown by the comparison that the results of the recommended model agree better with the numerical solutions and Opiolka et al.'s data than those of Schmidt's simple model.

• 유체기계공업학회:학술대회논문집
• /
• 2001.11a
• /
• pp.63-68
• /
• 2001

A numerical analysis to predict the flow characteristics in turbo fan for vacuum cleaner has been performed by using CFD. The rotating reference frame method is applied to compute the impeller-diffuser interaction and the characteristics of two models, 460 and 380, are calculated for various rotating speeds and flow rates. The flow in impeller, diffuser and return channel is assumed as steady and compressible. STAR-CD with k- $\epsilon$ turbulence model is used to solve the Navier-Stokes equations. Computed relative velocity, absolute pressure and flow angles are shown and compared with measurement results. The good agreement between the predictions and measurement results confirms the validity of this study.

• Journal of the Korean Solar Energy Society
• /
• v.27 no.4
• /
• pp.1-9
• /
• 2007

Sensitivity analysis of wind resource micrositing has been performed through the application case at the Antarctic King Sejong station with the most representative micrositing softwares: WAsP, WindSim and Meteodyn WT. The wind data obtained from two met-masts separated 625m were applied as a climatology input condition of micro-scale wind mapping. A tower shading effect on the met-mast installed 20m apart from the warehouse has been assessed by the CFD software Fluent and confirmed a negligible influence on wind speed measurement. Theoretically, micro-scale wind maps generated by the two met-data located within the same wind system and strongly correlated meteor-statistically should be identical if nothing influenced on wind prediction but orography. They, however, show discrepancies due to nonlinear effects induced by surrounding complex terrain. From the comparison of sensitivity analysis, Meteodyn WT employing 1-equation turbulence model showed 68% higher RMSE error of wind speed prediction than that of WindSim using the ${\kappa}-{\epsilon}$ turbulence model, while a linear-theoretical model WAsP showed 21% higher error. Consequently, the CFD model WindSim would predict wind field over complex terrain more reliable and less sensitive to climatology input data than other micrositing models. The auto-validation method proposed in this paper and the evaluation result of the micrositing softwares would be anticipated a good reference of wind resource assessments in complex terrain.