• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1150-1157
• /
• 2003

Nonlinear relationship between Reynolds stresses and the rate of strain of nonlinear k-$\varepsilon$models is evaluated theoretically by using the boundary layer assumptions against the turbulence-driven secondary flows in noncircular ducts and then their prediction performance is validated numerically through the application to the fully developed turbulent flow in a square duct. Typical predicted quantities such as mean axial and secondary velocities, turbulent kinetic energy and Reynolds stresses are compared with available experimental data. The nonlinear k-$\varepsilon$ model adopted in a commercial code is found to be unable to predict accurately duct flows with the prediction level of secondary flows one order less than that of the experiment.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.1-7
• /
• 2001

Wind flow perturbations, recirculations and turbulence generated by buildings often dominate air pollutant distributions around buildings. This paper describes dispersion of contaminants in the vicinity of a building by solving the concentration equation based on previously simulated wind flow field. Turbulence closure is achieved by using the standard k-ε two-equation model. The paper shows application of the CIP method for solving a species concentration equation of contaminant gas around a rectangular building for two different sources under conditions of neutral atmospheric stratification. Results have been compared to the experimental data and the previous numerical results by hybrid scheme. The computational results of concentration profiles by the CIP method agree well with experimental data.

• Journal of Korean Society for Atmospheric Environment
• /
• v.20 no.1
• /
• pp.47-58
• /
• 2004

The $textsc{k}$-$\varepsilon$ algebraic stress model (KEASM) was applied to atmospheric dispersion simulation using the Lagrangian particle dispersion model and was compared with the most popular turbulence closure model in the field of atmospheric simulation, the Mellor-Yamada (MY) model. KEASM has been rarely applied to atmospheric simulation, but it includes the pressure redistribution effect of buoyancy due to heat and momentum fluxes. On the other hand, such effect is excluded from MY model. In the simulation study, the difference in the two turbulence models was reflected to both the turbulent velocity and the Lagrangian time scale. There was little difference in the vertical diffusion coefficient $\sigma$$_{z}$. However, the horizontal diffusion coefficient or calculated by KEASM was larger than that by MY model, coincided with the Pasquill-Gifford (PG) chart. The applicability of KEASM to atmospheric simulations was demonstrated by the simulations.s.

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

• Wind and Structures
• /
• v.30 no.1
• /
• pp.69-83
• /
• 2020

Modelling incompressible, neutrally stratified, barotropic, horizontally homogeneous and steady-state atmospheric boundary layer (ABL) is an important aspect in computational wind engineering (CWE) applications. The ABL flow can be viewed as a balance of the horizontal pressure gradient force, the Coriolis force and the turbulent stress divergence. While much research has focused on the increase of the wind velocity with height, the Ekman layer effects, entailing veering - the change of the wind velocity direction with height, are far less concerned in wind engineering. In this paper, a modified k-ε model is introduced for the ABL simulation considering wind veering. The self-sustainable method is discussed in detail including the precursor simulation, main simulation and near-ground physical quantities adjustment. Comparisons are presented among the simulation results, field measurement values and the wind profiles used in the conventional wind tunnel test. The studies show that the modified k-ε model simulation results are consistent with field measurement values. The self-sustainable method is effective to maintain the ABL physical quantities in an empty domain. The wind profiles used in the conventional wind tunnel test have deficiencies in the prediction of upper-level winds. The studies in this paper support future practical super high-rise buildings design in CWE.

• Journal of Korea Water Resources Association
• /
• v.44 no.3
• /
• pp.189-198
• /
• 2011

The aim of this study is to compare k-$\varepsilon$ and k-$\omega$ closures under the condition of oscillatory layer flow at high Reynolds number. A one dimensional vertical model incorporated with flow momentum equations and turbulence models (k-$\varepsilon$ and k-$\omega$) is applied to the laboratory measurements in the turbulent oscillatory boundary layer. The numerical simulation reveals that both turbulence models calculate similar velocity profiles and turbulent kinetic energy (TKE). In addition, both deliver high accuracy under the condition of negligible spanwise pressure gradient. Therefore, it is recommended in this study to use k-$\varepsilon$ closure, of which numerical coefficients have been calibrated from many studies, for the cases of straight channel, estuary, and coastal environment where the spanwise pressure gradient is not significant.

• Wind and Structures
• /
• v.4 no.6
• /
• pp.469-480
• /
• 2001

The Leipzig Wind Profile is generally known as a typical neutral planetary boundary layer flow. But it became clear from the present research that it was not completely neutral but weakly stable. We examined whether we could simulate the Leipzig Wind Profile by using a ($k-{\varepsilon}$) turbulence model including the equation of potential temperature. By solving analytically the Second Moment Closure Model under the assumption of local equilibrium and under the condition of a stratified flow, we expressed the turbulent diffusion coefficients (both momentum and thermal) as functions of flux Richardson number. Our ($k-{\varepsilon}$) turbulence model which included the equation of potential temperature and the turbulent diffusion coefficients varying with flux Richardson number reproduced the Leipzig Wind Profile.

• Journal of Hydrogen and New Energy
• /
• v.34 no.2
• /
• pp.196-204
• /
• 2023

This study performed the numerical analysis of the internal flow phenomena of 1 kWe-class solid oxide fuel cell (SOFC) stacks with internal manifold type and planar cells using commercial computational fluid dynamics (CFD) software, Star-CCM+. In particular, the locations where the turbulent phenomena occur inside the SOFC stack were investigated. In addition, the laminar flow model and the standard k-ε turbulent model were used to calculate the SOFC stack, separately. And, the calculation results of both laminar and turbulent models were compared. The calculation results showed that turbulent phenomena occurred mainly in the cathode flow. Especially, the turbulent phenomena were found in the cathode inlet/outlet region, and local turbulence occurred in the end plate near the inlet pipe.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.2
• /
• pp.188-198
• /
• 2012

An improvement to the k-${\varepsilon}$ turbulence model is presented and is shown to lead to better agreement with data regarding supersonic base flows. The improvement was achieved by imposing a grid-independent realizability constraint in the Launder-Sharma k-${\varepsilon}$ model. The effects of compressibility were also examined. The numerical results show that the modified Launder-Sharma model leads to some improvement in the prediction of the velocity and turbulent kinetic energy profiles. Compressibility corrections also lead to better agreement in both the turbulent kinetic energy and the Reynolds stress profiles with the experimental data.