• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• Journal of the Korean Society of Safety
• /
• v.10 no.3
• /
• pp.68-80
• /
• 1995

A series of parametric calculations have been performed in order to investigate the short-term and short-range plume and puff behavior of toxic gaseous and solid pollutant dispersion in an open atmosphere. The simulation is made by the use of the computer program developed by this laboratory, in which a control-volume based finite-difference method is used together with the SIMPLEC algorithm for the resolution of the pressure-velocity coupling appeared In Wavier-Stokes equation. The Reynolds stresses are solved by the standard two-equation k-$\varepsilon$ model modified for buoyancy together with the RNG(Renormalization Group) k-$\varepsilon$ model. The major parameters considered in this calculation are pollutant gas density and temperature, the relative velocity of pollutants to that of the surrounding atmospheric air, and particulate size and density together with the height released. The flow field is typically characterized by the formation of a strong recirculation region for the case of the low density gases such as $CH_4$ and air due to the strong buoyancy, while the flow is simply declining pattern toward the downstream ground for the case of heavy molecule like the $CH_2C1_2$and $CCl_4$, even for the high temperature, $200^{\circ}C$. The effect of gas temperature and velocity on the flow field together with the particle trajectory are presented and discussed in detail. In general, the results are physically acceptable and consistent.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.2 no.3
• /
• pp.20-35
• /
• 1998

Characteristics of turbulent flow with wall transpiration is analyzed. The wall transpiration includes both of suction and injection and extends their range to 0~160 of absolute magnitude of Re$_{w}$ . Reynolds number based on inlet velocity also covers wide range of 3${\times}$$10^3$~8${\times}$$10^4$. The turbulent flow with wall transpiration induces change of wall boundary layer and rapid change of turbulent field. This, in turn, leads the change of whole flow field. For predicting this complicated flow field properly, newly modified $\kappa$-$\varepsilon$ model is utilized, which is formed by modifying dissipation rate equation. The modified $\kappa$-$\varepsilon$ model of Chien is also adopted for the comparison of model performance. Analysis shows the newly modified $\kappa$-$\varepsilon$ model is successfully able to reflect the characteristics of turbulent flow field with wall transpiration.ion.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.12 no.9
• /
• pp.3808-3814
• /
• 2011

The prediction of heat transfer and pressure drops in the exchanger passages is a clue to the problem of heat exchanger design. In order to make such predictions for non-Newtonian fluids, it is necessary to know the relation between the viscous properties of the fluid and the wall shear rate in the duct. This study deals with the limits of validity of the power law equation. The useful methodology of the present research involves a consideration of a more general equation which has power law and Newtonian behavior as asymptotes. It isconcluded that use of the power law equation outside of its applicability range can lead to serious errors inpredicting the heat transfer and pressure drops. The present computational results of the friction factors times Reynolds number for shear-thinning fluid flows in a triangular duct are compared with previous published results, showing agreement with 0.13 % in Newtonian region and 2.85 % in power law region. These shear-thinning fluid results also showed the 12% increase of convective heat transfer enhancement compared with Newtonian heat transfer.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.9
• /
• pp.907-915
• /
• 2008

The computational fluid dynamic analysis has been conducted for the thermo-chemical flow field in an arcjet thruster with mono-propellant Hydrazine (N2H4) as a working fluid. The Reynolds Averaged Navier-Stokes (RANS) equations are modified to analyze compressible flows with the thermal radiation and electric field. the Maxwell equation, which is loosely coupled with the fluid dynamic equations through the Ohm heating and Lorentz forces, is adopted to analyze the electric field induced by the electric arc. The chemical reactions of Hydrazine were assumed to be infinitely fast due to the high temperature field inside the arcjet thruster. The chemical and the thermal radiation models for the nitrogen-hydrogen mixture and optical thick media respectively, were incorporated with the fluid dynamic equations. The results show that performance indices of the arcjet thruster with 1kW arc heating are improved by amount of 180% in thrust and 200% in specific impulse more than frozen flow. In addition thermo-physical process inside the arcjet thruster is understood from the flow field results.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• Journal of Mechanical Science and Technology
• /
• v.17 no.3
• /
• pp.440-448
• /
• 2003

In the present paper, numerical simulations of buoyant plume dispersion in a neutral and stable atmospheric boundary layer have been carride out. A Lagrangian Stochastic Model (LSM) with a Non-Linear Eddy Viscosity Model (NLEVM) for turbulence is used to generate a Reynolds stress field as an input condition of dispersion simulation. A modified plume-rise equation is included in dispersion simulation in order to consider momentum effect in an initial stage of plume rise resulting in an improved prediction by comparing with the experimental data. The LSM is validated by comparing with the prediction of an Eulerian Dispersion Model (EDM) and by the measured results of vertical profiles of mean concentration in the downstream of an elevated source in an atmospheric boundary layer. The LSM predicts accurate results especially in the vicinity of the source where the EDM underestimates the peak concentration by 40% due to inherent limitations of gradient diffusion theory. As a verification study, the LSM simulation of buoyant plume dispersions under a neutral and stable atmospheric condition is compared with a wind-tunnel experiment, which shows good qualitative agreements.

• Journal of Mechanical Science and Technology
• /
• v.14 no.3
• /
• pp.350-357
• /
• 2000

Experiments were performed by using PF-5060 and water to investigate the influence of an aspect ratio of a horizontal rectangular channel on the cooling characteristics from an in-line $6{\times}1$ array of discrete heat sources which were flush mounted on the top wall of the channel. The experimental parameters were aspect ratio of rectangular channel, heat flux of simulated VLSI chip, and channel Reynolds number. The chip surface temperatures decreased with the aspect ratio at the first and sixth rows, and decreased more rapidly at a high heat flux than at a low heat flux. The measured friction factors at each aspect ratio for both water and PF-5060 gave a good agreement with the values predicted by the modified Blasius equation within ${\pm}7%$. The Nusselt number increased as the aspect ratio decreased, but the increasing rate of Nusselt number reduced as the aspect ratio decreased. A 5:1 rectangular channel yields the most efficient cooling performance when the heat transfer and pressure drop in the test section were considered simultaneously.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1487-1494
• /
• 1990

In this paper the dynamic characteristics of self acting air lubricated slider bearing of hard disk/head system are investigated. The dynamic equations of magnetic head mechanism considering both parallel and pitch motion and the time dependent modified Reynolds equation are analyzed and the dynamic pressure distribution of air film is numerically calculated in frequency domain by small perturbation method and finite difference scheme with variable grid. The dynamic response of the slider spacing is obtained accordingly as the moving recording surface vibrates in parallel mode.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.