• Journal of the Korean Statistical Society
• /
• v.28 no.2
• /
• pp.151-165
• /
• 1999

This paper considers the trimmed least squares estimator of the autoregression parameter in the unstable AR(1) model: X\ulcorner=ØX\ulcorner+$\varepsilon$\ulcorner, where $\varepsilon$\ulcorner are iid random variables with mean 0 and variance $\sigma$$^2$> 0, and Ø is the real number with │Ø│=1. The trimmed least squares estimator for Ø is defined in analogy of that of Welsh(1987). The limiting distribution of the trimmed least squares estimator is derived under certain regularity conditions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.231-237
• /
• 2009

This study presents a numerical model to simulate density currents developing two dimensionally. The ${\kappa}-{\varepsilon}$ model is used for the turbulence closure. Elliptic flow equations are solved by the finite volume method. In order to investigate the applicability of the numerical model, discontinuous density currents are simulated numerically. The vortices due to the instability at the interface are simulated, showing a good agreement with the experimental visualizations in the literature. It is also investigated that the transition from slumping phase to inertial phase occurs when a bore generated at the end wall overtakes the front. However, the propagation of the density current is retarded compared with the experimental results. Two-dimensional modeling seems to have an effect on underestimating the front velocity of the density current.

• Wind and Structures
• /
• v.4 no.3
• /
• pp.177-196
• /
• 2001

At present the most popular turbulence models used for engineering solutions to flow problems are the $k-{\varepsilon}$ and Reynolds stress models. The shortcoming of these models based on the isotropic eddy viscosity concept and Reynolds averaging in flow fields of the type found in the field of Wind Engineering are well documented. In view of these shortcomings this paper presents the implementation of a non-linear model and its evaluation for flow around a building. Tests were undertaken using the classical bluff body shape, a surface mounted cube, with orientations both normal and skewed at $45^{\circ}$ to the incident wind. Full-scale investigations have been undertaken at the Silsoe Research Institute with a 6 m surface mounted cube and a fetch of roughness height equal to 0.01 m. All tests were originally undertaken for a number of turbulence models including the standard, RNG and MMK $k-{\varepsilon}$ models and the differential stress model. The sensitivity of the CFD results to a number of solver parameters was tested. The accuracy of the turbulence model used was deduced by comparison to the full-scale predicted roof and wake recirculation zone lengths. Mean values of the predicted pressure coefficients were used to further validate the turbulence models. Preliminary comparisons have also been made with available published experimental and large eddy simulation data. Initial investigations suggested that a suitable turbulence model should be able to model the anisotropy of turbulent flow such as the Reynolds stress model whilst maintaining the ease of use and computational stability of the two equations models. Therefore development work concentrated on non-linear quadratic and cubic expansions of the Boussinesq eddy viscosity assumption. Comparisons of these with models based on an isotropic assumption are presented along with comparisons with measured data.

• 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 Ocean Engineering and Technology
• /
• v.22 no.3
• /
• pp.24-33
• /
• 2008

In this study, flow and creation and dissipation of vorticity around rectangular floating breakwater is investigated both experimentally and numerically. The PIV system(Particle image velocimetry) is employed to obtain the velocity field in the vorticity of rectangular structure. The numerical model, combined with ${\kappa}-{\varepsilon}$ turbulence model and the VOF method based on RANS equation, is used to analyze the turbulence structure. In the results of this study, the vorticity is found around conner of rectangular structure at all time domain, and creation and dissipation of vorticity are closely related to wave period. Separation points of phase of vortex due to flow separation for longer period waves are faster then for shorter period waves.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.2
• /
• pp.157-166
• /
• 1999

Numerical calculations of turbulent natural convection in an enclosure of the 20 kYA oil-immersed transformer model are presented. The transformer is modelled as two concentric cylinders with different heights and diameters. The thermal boundary layers are well represented in the temperature distributions along the wall of the transformer model. The flow stratification between the hot and cold walls can not be seen in the transformer model. The turbulence eddy viscosity has its maximum at the center of the core and its maximum values at the top of the core are larger than those at the bottom of the core.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.13-14
• /
• 2003

Two different ${\kappa}-{\varepsilon}-\overline{{\nu}{\nu}}-f$ turbulence models together with the two-layer model are evaluated for natural convection in a rectangular cavity. The numerical problem and accuracy of the turbulence models are discussed. The original $\overline{{\nu}{\nu}}-f$ model suffers from the numerical stiffness problem when used with the segregate solution procedure like the SIMPLE algorithm, and a remedy for this problem is proposed. It is shown that original $\overline{{\nu}{\nu}}-f$ model best predicts the mean velocity, Reynolds stresses and the turbulent heat flux while the modified $\overline{{\nu}{\nu}}-f$ model (N=6) overpredicts the turbulent quantities.

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.317-328
• /
• 2002

A numerical investigation on the turbulent flows over a three-dimensional steep hill is presented. The numerical model developed for the present work is based on the finite volume method and the SIMPLE algorithm with a non-staggered grid system. Standard $k-{\varepsilon}$ model and Shih's non-linear model are tested for the validation of the prediction accuracy in the 3D separated flow. Comparisons of the mean velocity and turbulence profiles between the numerical predictions and the measurements show good agreement. The Shih's non-linear model is found to predict mean flow and turbulence better than the Standard $k-{\varepsilon}$. Flow patterns have also been examined to explain the difference in the cavity zone between 2D and 3D hills.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.5
• /
• pp.48-57
• /
• 1994

The specific power output and thermal effeciency of any two-stroke engine are dependent on its scavenging behavior. Among the many factors which influence on the scavenging process, the cylinder head shape is one of the important factor. Hence in this study three different type models of cylinder head shape which are the cylindrical, the spherical and the arbitrary shape are studied to show the effects of the turbulent scavenging process in the cylinder with one inlet port, two side ports and one exhaust port. A modified version of KIVA-II which strip out of or add planes of cells across the mesh above the piston for flow simulation of two-stroke engine is used. The $k-{\varepsilon}$ turbulent model is used. The results show that the flow in a two-stroke engine cylinder of the spherical head shape among the three different type model is a desirable for efficient scavenging.