• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.745-750
• /
• 2000

A numerical code which can simulate unsteady, incompressible and 3-dimensional flows in an engine cylinder has been developed. The governing equations based on the cylindrical coordinate are discretized by the finite volume method with staggered variable arrangements. A geometric conservation rule is also incorporated into the simulation code in order to deal with a moving boundary problem. For the unsteady simulation, a fractional step method is adopted. The law of wall is applied to the wall boundaries and standard $k-\;{\varepsilon}$ model is used to describe the in-cylinder turbulent flow. The model cylinder has one eccentric port, flat piston and flat cylinder-head. The comparisons with experimental data show fairly well qualitative agreement.

• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.9-16
• /
• 1998

Predictive behaviors by the extended k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ turbulence model are compared. Grid dependency is tested with the H-type grid as well as the O-type grid. Computations have been performed on a circular-to-rectangular transition duct. The Reynolds number is 390,000 based on the bulk velocity at the inlet. The computed axial velocity contours, transverse velocity profiles, static pressure contours, peripheral skin friction coefficient, peripheral wall static pressure distributions and turbulence kinetic energy have been compared with experimental results. The computed results than those obtained with the standard k-${\varepsilon}$ turbulence model. Comparing to the computed results obtained with the H-type grid and O-type grid, those with H-type grid seem to agree well with experimental results.

• Journal of Mechanical Science and Technology
• /
• v.14 no.1
• /
• pp.93-102
• /
• 2000

The objective of the present study is to investigate the pressure-strain correlation terms of the Reynolds stress models for the three dimensional turbulent boundary layer in a $30^{\circ}$ bend tunnel. The numerical results obtained by models of Launder, Reece and Rodi (LRR) , Fu and Speziale, Sarkar and Gatski (SSG) for the pressure-strain correlation terms are compared against experimental data and the calculated results from the standard k-${\varepsilon}$ model. The governing equations are discretized by the finite volume method and SIMPLE algorithm is used to calculate the pressure field. The results show that the models of LRR and SSG predict the anisotropy of turbulent structure better than the standard k-${\varepsilon}$ model. Also, the results obtained from the LRR and SSG models are in better agreement with the experimental data than those of the Fu and standard k-${\varepsilon}$ models with regard to turbulent normal stresses. Nevertheless, LRR and SSG models do not effectively predict pressure-strain redistribution terms in the inner layer because the pressure-strain terms are based on the locally homogeneous approximation. Therefore, to give better predictions of the pressure-strain terms, non-local effects should be considered.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2593-2600
• /
• 1996

A modified k-$\varepsilon$ turbulence model having a generality is proposed in the present study, in which the constant $C_{\varepsilon2}$in the $\varepsilon$-equation is simply changed as a functional form of a new parameter both satisfying the tensor invariant condition and representing the extra straining effect on complex shear flows. With this model turbulent shear flows over two-dimensional obstacles placed in a channel are numerically studied for different blockage ratios and aspect ratios. Comparing with the available experimental data, the predicted results with the present model provide definite improvements over the standard model's results and work fairly well with the experimental data on the size of the recirculation zone, as well as mean velocity, wall static pressure, turbulent kinetic energy and Reynolds stresses.

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

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..

• Journal of the Korean Mathematical Society
• /
• v.47 no.2
• /
• pp.263-275
• /
• 2010

Let {$X_n;n\;\geq\;1$} be a strictly stationary sequence of negatively associated random variables with mean zero and finite variance. Set $S_n\;=\;{\sum}^n_{k=1}X_k$, $M_n\;=\;max_{k{\leq}n}|S_k|$, $n\;{\geq}\;1$. Suppose $\sigma^2\;=\;EX^2_1+2{\sum}^\infty_{k=2}EX_1X_k$ (0 < $\sigma$ < $\infty$). We prove that for any b > -1/2, if $E|X|^{2+\delta}$(0<$\delta$$\leq$1), then $$lim\limits_{\varepsilon\searrow0}\varepsilon^{2b+1}\sum^{\infty}_{n=1}\frac{(loglogn)^{b-1/2}}{n^{3/2}logn}E\{M_n-\sigma\varepsilon\sqrt{2nloglogn}\}_+=\frac{2^{-1/2-b}{\sigma}E|N|^{2(b+1)}}{(b+1)(2b+1)}\sum^{\infty}_{k=0}\frac{(-1)^k}{(2k+1)^{2(b+1)}}$$ and for any b > -1/2, $$lim\limits_{\varepsilon\nearrow\infty}\varepsilon^{-2(b+1)}\sum^{\infty}_{n=1}\frac{(loglogn)^b}{n^{3/2}logn}E\{\sigma\varepsilon\sqrt{\frac{\pi^2n}{8loglogn}}-M_n\}_+=\frac{\Gamma(b+1/2)}{\sqrt{2}(b+1)}\sum^{\infty}_{k=0}\frac{(-1)^k}{(2k+1)^{2b+2'}}$$, where $\Gamma(\cdot)$ is the Gamma function and N stands for the standard normal random variable.

• Wind and Structures
• /
• v.12 no.1
• /
• pp.21-47
• /
• 2009

The logarithmic form for turbulent flow analysis guarantees the positivity of the turbulence variables as ${\kappa}$ and ${\varepsilon}$ of the ${\kappa}-{\varepsilon}$ model by using the natural logarithm of these variables. In the present study, the logarithmic form is incorporated into the finite element solution procedure for the unsteady turbulent flow analysis. A backward facing step flow using the standard ${\kappa}-{\varepsilon}$ model and a flow around a 2D square cylinder using the modified ${\kappa}-{\varepsilon}$ model (the Kato-Launder model) are simulated. These results show that the logarithmic form effectively keeps adequate balance of turbulence variables and makes the analysis stable during transient or unsteady processes.

• Ocean Systems Engineering
• /
• v.4 no.3
• /
• pp.185-200
• /
• 2014

The objective of the present simulations is to evaluate the applicability of the standard $k-{\varepsilon}$ turbulence model in engineering practice in the subcritical to supercritical flow regimes. Two-dimensional numerical simulations of flow around a circular cylinder at $Re=1{\times}10^5$, $5{\times}10^5$ and $1{\times}10^6$, had been performed using Unsteady Reynolds-Averaged Navier Stokes (URANS) equations with the standard $k-{\varepsilon}$ turbulence model. Solution verification had been studied by evaluating grid and time step size convergence. For each Reynolds number, several meshes with different grid and time step size resolutions were chosen to calculate the hydrodynamic quantities such as the time-averaged drag coefficient, root-mean square value of lift coefficient, Strouhal number, the coefficient of pressure on the downstream point of the cylinder, the separation angle. By comparing the values of these quantities of adjacent grid or time step size resolutions, convergence study has been performed. Solution validation is obtained by comparing the converged results with published numerical and experimental data. The deviations of the values of present simulated quantities from those corresponding experimental data become smaller as Reynolds numbers increases from $1{\times}10^5$ to $1{\times}10^6$. This may show that the standard $k-{\varepsilon}$ model with enhanced wall treatment appears to be applicable for higher Reynolds number turbulence flow.

• Journal of Ship and Ocean Technology
• /
• v.5 no.2
• /
• pp.27-48
• /
• 2001

Turbulent flow calculations are performed for the two modern practical VLCCs with the sable forebody and the slightly different afterbody, i.e. KVLCC and KVLCC2. Three $\textsc{k}-\varepsilon$ turbulence models are tested to investigate the differences caused by the turbulence models. The calculated results around the two VLCC hull forms using O-O grid topology and profile-fitted surface meshes are compared to the measured data from towing tank experiment. The realizable $\textsc{k}-\varepsilon$model provided realistic wake distribution with hook-like shape, while the standard and RNG-based $\textsc{k}-\varepsilon$models failed. It is very encouraging to see that the CFD with relatively simple turbulence closure can tell the difference quantitatively as well as qualitatively for the two hull forms with stern frameline modification.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.741-750
• /
• 1995

The turbulent buoyancy-driven flow in 2-dimensional enclosed cavities heated from the vertical side is numerically calculated for both cases of a Rayleigh number of 5*10$^{10}$ for air and 2.5*10$^{10}$ for water. Three different turbulence models are considered : standard k-.epsilon. model of Ozoe and low-Reynolds-number model of Lam and Bremhorst, and another low-Reynolds-number model of Davidson. The results indicate that the use of low-Reynolds number models is recommended for the indoor airflow computation, and the results from Davidson model are reasonably close to the reported experimental data. A sensitivity study shows that the amounts of wall-heat transfer and the velocity profiles with the Lam and Bremhorst model largely depend on the choice of the wall function for .epsilon..