• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.316-323
• /
• 2000

A nonlinear low-Reynolds-number heat transfer model is developed to predict turbulent flow and heat transfer in separated and reattaching flows. The $k-{\varepsilon}-f_{\mu}$ model of Park and Sung (1997) is extended to a nonlinear formulation, based on the nonlinear model of Gatski and Speziale (1993). The limiting near-wall behavior is resolved by solving the $f_{\mu}$ elliptic relaxation equation. An improved explicit algebraic heat transfer model is proposed, which is achieved by applying a matrix inversion. The scalar heat fluxes are not aligned with the mean temperature gradients in separated and reattaching flows; a full diffusivity tensor model is required. The near-wall asymptotic behavior is incorporated into the $f_{\lambda}$ function in conjunction with the $f_{\mu}$ elliptic relaxation equation. Predictions of the present model are cross-checked with existing measurements and DNS data. The model preformance is shown to be satisfactory.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.190-193
• /
• 2010

The turbulent mixing and combustion of a shear coaxial injector under supercritical pressures have been theoretically/numerically investigated. Turbulent numerical model is based on large eddy simulation with real-fluid transport and thermodynamics over the entire pressure; Soave modification of Redlich-Kwong equation of state, Chung's model for viscosity/conductivity, and Fuller's theorem for diffusivity to take account Takahashi's compressible effect. The results are compared with previous researcher's. The large-scale vortices shedding from the outer rim into the recirculation region to react with gaseous oxygen was investigated.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.04a
• /
• pp.103-108
• /
• 2011

The turbulent mixing of a kerosene/liquid oxygen coaxial swirl injector under supercritical pressures have been numerically investigated. Kerosene surrogate models are proposed for the kerosene thermodynamic properties. Turbulent numerical model is based on LES(Large Eddy Simulation) with real-fluid transport and thermodynamics over the entire pressure range; Soave modification of Redlich-Kwong equation of state, Chung's model for viscosity/conductivity, and Fuller's theorem for diffusivity to take account Takahashi's compressible effect. The effect of operating pressure on thermodynamic properties and mixing dynamics inside an injector and a combustion chamber are investigated. Power spectral densities of pressure fluctuations in the injector under various chamber pressure are analyzed.

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

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.12 no.4
• /
• pp.390-399
• /
• 2007

In order to investigate the mixing of sea waters on the continental shelf in the northern East China Sea, Korea Ocean Research and Development Institute conducted hydrographic surveys including turbulence measurements using the R/V Eardo in August 2005 and August 2006. The turbulent kinetic energy dissipation rates based on velocity shear measurements are estimated to be $10^{-7}{\sim}10^{4}$, $10^{-7}{\sim}10^{-6}$, and $10^{-7}$ W/kg in the surface layer, bottom layer, and lower thermocline, respectively. The data sets suggest that surface layer water is being constantly mixed by winds. High dissipation rate in the lower thermocline seems to be caused by internal waves. The bottom layer with high dissipation rate also shows high turbidity, indicating the effect of tidal stirring turbulence. The vertical eddy diffusivities are $10^{-3}{\sim}10^{-2}m^2/s$ near the bottom, and these high values appear to arise from both the low stability and high turbulent mixing.

• Journal of Mechanical Science and Technology
• /
• v.19 no.3
• /
• pp.913-923
• /
• 2005

The dynamic mixed model (DMM) combined with a box filter of Zang et. al. (1993) has been generalized for passive scalar transport and applied to large eddy simulation of turbulent channel flows with Prandtl number up to 10. Results from a priori test showed that DMM is capable of predicting both subgrid-scale (SGS) scalar flux and dissipation rather accurately for the Prandtl numbers considered. This would suggest that the favorable feature of DMM, originally developed for the velocity field, works equally well for scalar transport problem. The validity of the DMM has also been tested a posteriori. The results of the large eddy simulation showed that DMM is superior to the dynamic Smagorinsky model in the prediction of scalar field and the model performance of DMM depends to a lesser degree on the ratio of test to grid filter widths, unlike in the a priori test.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.9
• /
• pp.2443-2453
• /
• 1994

The structure of turbulence of fully developed flow through four concentric annuli with the rough outer wall was investigated experimentally for a Reynolds number range Re=15, 000-93, 000. Turbulence intensities were measured in three(u, v, w) directions, and turbulence shear stresses in annuli of radius=0.13, 0.26, 0.4 and 0.56, respectively. Due to the square roughness element attached periodically along the axial direction, the radial velocity fluctuations show similar distribution regardless of the different .alpha.cases. However, the axial and circumferential velocity fluctuation profiles demonstrate the longitudinal turbulence structures are strongly influenced by the .alpha. values. The turbulent eddy viscosity deduced form mean velocity distributions and the measured Reynolds shear stresses are also presented and discussed.

• Membrane and Water Treatment
• /
• v.3 no.2
• /
• pp.77-98
• /
• 2012

A two-dimensional (2D) steady state numerical model of concentration polarisation (CP) phenomena in a membrane channel has been developed using the commercially available computational fluid dynamics (CFD) package CFX (Ansys, Inc., USA). The model incorporates the transmembrane pressure (TMP), axially variable permeate flux, variable diffusivity and viscosity, and osmotic pressure effects. The model has been verified against several benchmark analytical and empirical solutions from the membrane literature. Additionally, the model is able to predict the rejection of an arbitrary solute by the membrane using a pore model, given some basic knowledge of the geometry of the solute molecule or particle, and the membrane pore geometry. This allows for predictive design of membrane systems without experimental determination of the membrane rejection for the specified operating conditions. A demonstration of the model is presented against experimental results for two uncharged test compounds (sucrose and PEG1000) from the literature. The model will be extended to incorporate charge effects, transient simulations, three-dimensional (3D) geometry and turbulent effects in future work.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.1
• /
• pp.214-224
• /
• 1990

A "two-fluid" model using thermal eddy diffusivity concept and Lumley's drag reduction theory, is proposed to analyze heat transfer of the turbulent dilute gas-particle flow in a vertical pipe with constant wall heat flux. The thermal eddy diffusivity is derived to be a function of the ratio of the heat capacity-density products .rho. over bar $C_{p}$ of the gaseous phase and the particulate phase and also of the ratio of thermal relaxation time scale to that of turbulence. The Lumley's theory dictates the variation of the viscous sublayer thickness depending on the particle loading ratio Z and the relative particle size $d_{p}$/D. At low loading ratio, the size of viscous sublayer thickness is important for suspension heat transfer, while at higher loading, the effect of the ratio .rho. $_{p}$ over bar $C_{p}$$_{p}$/ .rho. $_{f}$ over bar $C_{p}$$_{f}$ is dominant. The major cause of decrease in the suspension Nusselt number at lower loading ratio is found to be due to the increase of the viscous sublayer thickness caused by the suppression of turbulence near the wall by the presence of solid particles. Predicted Nusselt numbers using the present model are in satisfactory agreements with available experimental data both in pipe entrance and the fully developed regions.

• Journal of Korean Society for Atmospheric Environment
• /
• v.13 no.5
• /
• pp.383-396
• /
• 1997

Numerical prediction of the pollutant dispersion over a two-dimensional hilly terrain is presented. The dispersion model used in the present work is based on the gradient diffusion theory and the finite-volume method on a non-orthogonal boundary-fitted grid system. The numerical model is validated by comparing the results with the available experimental data for the flat-floor dispersion within a turbulent boundary-layer. The numerical error analysis is performed based on the guideline of Kasibhatla et al.(1988) for the elevated-source dispersion in the flat-floor boundary layer having a power-law velocity and linear eddy-diffusivity profile. The influences of the two-dimensional hilly terrain on the dispersion from a continuously released source are numerically investigated by changing the emission locations and heights. It is found that the distributions of ground-level concentration are strongly influenced by the source location and the emission height. Hence, the terrain amplification factor is greatly enhanced when the pollutant source is located within a flow separation region. Dispersion from a source of short duration is also simulated and the duration time of the pollutant is compared at several downstream locations on a hilly terrain. The results of the numerical prediction are applied to the evaluation of environmental impacts due to the automobile exhausts at the seashore highway with a parallel mountain range.