• Water Engineering Research
• /
• v.3 no.4
• /
• pp.247-258
• /
• 2002

This paper presents a non-isotropic turbulence modeling of flows over bedforms. The Reynolds stress model is used for the turbulence closure. In the model, Launder, Reece, and Rodi's model and Hanjalic and Launder's model are employed f3r the pressure strain correlation term and the diffusion term, respectively. The mean flow and turbulence structures are simulated and compared with profiles measured in the experiments. The numerical solutions from two-equation turbulence models are also provided for comparisons. The Reynolds stress model yields the separation length of eddy similar to the other numerical results. Using the developed model, the resistance coefficients are also estimated for the flows at different Froude numbers. Karim's (1999) relationship is used to determine the bedform geometry. It is found that the values of the form drag and the skin friction are very similar to those obtained by the other turbulence models. meaning higher values of the form drag and lower values of the skin friction compared with the empirical formulas.

• Journal of the Korean Society of Combustion
• /
• v.14 no.2
• /
• pp.32-41
• /
• 2009

The transported probability density function model combined with the consistent finite volume (FV) method has been applied to simulate the turbulent bluff-body reacting flows. To realistically account for the non-isotropic turbulence effects on the turbulent bluff-body reacting flows, the present PDF transport approach is based on the joint velocity- turbulent frequency-composition PDF formulation. The evolution of the fluctuating velocity of a particle is modeled by a simplified Langevin equation and the particle turbulence frequency is represented by the modified Jayesh - Pope model. Effects of molecular diffusion are represented by the interaction by exchange with the mean (IEM) mixing model. To validate this hybrid FV/PDF transport model, the numerical results are compared with experimental data for the turbulent bluff-body reacting flows.

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

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

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.3
• /
• pp.833-840
• /
• 2014

This study presents numerical simulation of turbulent flows in a rectangular open-channel that has a width-to-depth ratio of 2 using the source code provided by OpenFOAM. Large eddy simulations are carried out by solving the filtered continuity and momentum equations numerically. For the non-isotropic residual stress term, Smagorinsky's (1963) model is used. The flow in the open-channel whose width-to-depth ratio is 2, from experiment of Tominaga et al. (1989), is simulated numerically. Simulation results are compared with measured data by Tominga et al. (1989) and Nezu and Rodi (1985) and with LES data by Shi et al. (1999). Comparisons revealed that the model simulates the mean flow and turbulence statistics well. Specifically, the model reproduced the inner secondary currents located at the corner of sidewall and free surface successfully. In addition, the vortical component of turbulence intensity shows bulged contours towards the bottom edge.

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

In the present paper, turbulent open-channel flows with alternate vegetated zones are numerically simulated using threedimensional model. The Reynolds-averaged Navier-Stokes Equations are solved with the ${\kappa}-{\varepsilon}$ model. The CFD code developed by Olsen(2004) is used for the present study. For model validation, the partly vegetated channel flows are simulated, and the computed depth-averaged mean velocity and Reynolds stress are compared with measured data in the literature. Comparisons reveal that the present model successfully predicts the mean flow and turbulent structures in vegetated open-channel. However, it is found that the ${\kappa}-{\varepsilon}$ model cannot accurately predict the momentum transfer at the interface between the vegetated zone and the non-vegetated zone. It is because the ${\kappa}-{\varepsilon}$ model is the isotropic turbulence model. Next, the open channel flows with alternate vegetated zones are simulated. The computed mean velocities are compared well with the previously reported measured data. Good agreement between the simulated results and the experimental data was found. Also, the turbulent flows are computed for different densities of vegetation. It is found that the vegetation curves the flow and the meandering flow pattern becomes more obvious with increasing vegetation density. When the vegetation density is 9.97%, the recirculation flows occur at the locations opposite to the vegetation zones. The impacts of vegetation on the flow velocity and the water surface elevation are also investigated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2386-2396
• /
• 1996

A numerical simulation has been performed for isothermal and reacting flows in an exisymmetric, bluff-body research combustor. The present formulation is based on the density-weighted averaged Navier-Stokes equations together with a k-epsilon. turbulence model and a modified eddy-breakup combustion model. The PISO algorithm is employed for solution of thel Navier-Stokes system. Comparison between measurements and predictions are made for a centerline axial velocities, location of stagnation points, strength of recirculation zone, and temperature profile. Even though the numerical simulation gives acceptable agreement with experimental data in many respects, the present model is defictient in predicting the recoveryt rate of a central near-wake region, the non-isotropic turbulence effects, and variation of turbulent Schmidt number. Several possible explanations for these discrepancies have been discussed.