• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.100-109
• /
• 2011

Large eddy simulation(LES) of fully developed turbulent pipe flow has been performed to investigate the effect of Reynolds number on the flow field at $Re_{\tau}$=180, 395, 590 based on friction velocity and pipe radius. A dynamic subgrid-scale model for the turbulent subgrid-scale stresses was employed to close the governing equations. The mean flow properties, mean velocity profiles and turbulent intensities obtained from the present LES are in good agreement with the previous numerical and experimental results currently available. The Reynolds number effects were observed in the higher-order statistics(Skewness and Flatness factor). Furthermore, the budgets of the Reynolds stresses and turbulent kinetic energy were computed and analyzed to elucidate the effect of Reynolds number on the turbulent structures.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.254-263
• /
• 2015

In this study, numerical simulations are carried out to investigate the dynamic SGS model effects in a Kerosene-LOx coaxial swirl injector under high pressure conditions. The turbulent model is based on large-eddy simulation (LES) with real-fluid transport and thermodynamics. To assess the effect of the dynamic subgrid-scale (SGS) model, the dynamic SGS model is compared with that of the algebraic SGS model. In a swirl injector under supercritical pressure, the characteristics of temporal pressure fluctuation and power spectral density (PSD) present comparable discrepancies dependant on the SGS models, which affect the mixing characteristics. Mixing efficiency and the probability density (PDF) function are conducted for a statistical description of the turbulent flow fields according to the SGS models. The back-scattering of turbulent kinetic energy is estimated in terms of the film thickness of the swirl injector.

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

Local mean fluctuating velocity components were measured in the separating and reattaching axisymmetrc region of turbulent boundary layer over the wall of convex cylinders placed in a water tunnel by using 2-color 4-beam fiber optics laser Doppler velocimetry. Measurements were made with three different diameters of cylinders with four different diameters of cylinders with four different diameter of the obstructions. The range of Reynolds number based on step height was between 5,000 to 25,200. The study demonstrates that the reattachment length decreases with decreasing cylinder radius and is always shorter than that for the two-dimensional backward-facing step flow at the condition of the same step height. It was also observed that the turbulent kinetic energy in the recirculating region increases with an increases in the radius of convex curvature. The measured velocity field suggests that the transverse curvature can effect definitely the formation of corner eddy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.12
• /
• pp.1014-1021
• /
• 2009

A new eddy viscosity equation was formulated from assumption of turbulence length scale equation and specific dissipation ratio equation. Then, a set of turbulence model equations for the turbulent kinetic energy ${\kappa}$, the viscosity ${\nu}_t$, and the intermittency factor ${\gamma}$ is proposed by considering the entrainment effect. Closure coefficients are determined by experimental data and resorting to numerical optimization. Present model has been applied to compute four representative cases of free shear flows and successfully compared with experimental data. In particular, the spreading rate, the centreline mean velocity and the profiles of intermittency are calculated with improved accuracy. Also, the proposed ${\nu}_t-{\kappa}-{\gamma}$ model was applied to channel flow by considering the wall effect and the results show good agreements with the Direct Numerical Simulation data.

• Ocean and Polar Research
• /
• v.23 no.2
• /
• pp.133-139
• /
• 2001

Ocean currents estimated from sea height anomalies derived from inter-calibrated TP/ERS are compared with daily mean currents measured with satellite-tracked drifters. The correlation coefficient between the geostrophic current from TP/ERS and surface current at 15 m depth from drifter tracks was found to be about 0.5. Due to the limitation of satellite ground tracks, small scale eddies less than 80 km are poorly resolved from TP/ERS. One of the interesting results of this study is that coastal currents along the eastern coast of Korea were well reproduced from sea height anomalies when the coastal currents were developed in association with eddies near the South Korean coast. The eddy kinetic energy (EKE) estimated from drifters, TP/ERS, and a numerical model are also compared. The EKE estimated from drifters was about 22 % higher than EKE calculated from TP/ERS. The pattern of low EKE level in the northern basin and high EKE level in the southern East Sea is shown in the EKE estimates derived from both the drifters and TP/ERS.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.32-38
• /
• 2006

The objective of this study is to evaluate the prediction accuracy of development large eddy simulation(LES) program for turbulent flow behind a bluff-body. The LES solver was implemented on parallel computer consisting 16 processors. To verify the capability of LES code, the results were compared with those of Reynolds Averaged Navier-Stokes(RANS) using standard ${\kappa}-{\varepsilon}$ model as well as experimental data. The results showed that the LES and RANS qualitatively well predicted the experimental results, such as mean axial, radial velocities and turbulent kinetic energy. In the quantitative analysis, however, the LES showed a better prediction performance than RANS. Specially, the LES well described characteristics of the recirculation zones, such as air stagnation point and jet stagnation point. Finally, the unsteady phenomena on the Bluff-body, such as the transition of recirculation region and vorticity, was examined with LES methodology.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.22 no.6
• /
• pp.773-783
• /
• 2019

A numerical analysis of a turbulent flow with a 'large-eddy breakup device(LEBU)' was performed to investigate the influence of the device on the drag of underwater vehicle using commercial CFD code, FLUENT. In the present study, the vehicle drag was decomposed to skin-friction coefficient(C_f) and pressure coefficient(C_p). The variation of the vehicle C_f and C_p were observed with changing location of the device and Reynolds number. As a result, the device decreased the vehicle C_f because it suppressed the turbulent characteristics behind the device. The larger Reynolds number, the higher reduction effect when the device was placed in front part of, and near the vehicle. On the other hand, the device increased/decreased the vehicle C_p with increasing/decreasing turbulent kinetic energy at recirculating flow region behind the vehicle. The total drag change by the device was caused by C_p rather than C_f.

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

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.819-822
• /
• 2002

It is necessary for the numerical simulation of 3-dimensional incompressible isotropic decaying turbulence to construct 3-dimensional initial velocity field which resembles the fully developed turbulence. Although the previous velocity field generation method proposed by Rogallo(1981) satisfies continuity equation and 3-dimensional energy spectrum, it has limitation, as indicated in his paper, that it does not produce the higher velocity moments(e. g. velocity derivative skewness) characteristic of real turbulence. In this study, a new velocity field generation method which is able to control velocity derivative skewness of initial velocity field is proposed. Brief descriptions of the new method and a few parameters which is used to control velocity derivative skewness are given. A large eddy simulation(LES) of isotropic decaying turbulence using dynamic subgrid-scale model is carried out to evaluate the performance of the initial velocity field generated by the new method. It was shown that the resolved turbulent kinetic energy decay curve and the resolved enstrophy decay curve from the initial field of new method were more realistic than those from the initial field of Rogallo's method. It was found that the dynamic model coefficient from the former was initially half the stationary value and experienced relatively short transition period, though that from the latter was initially zero and experienced relatively longer transition period.