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

• Journal of Ocean Engineering and Technology
• /
• v.32 no.4
• /
• pp.244-252
• /
• 2018

This paper provides quantification of the effects of the turbulence model and grid refinement on the analysis of tip vortex flows by using the RANS(Reynolds averaged Navier-Stokes) method. Numerical simulations of the tip vortex flows of the NACA $66_2$-415 elliptic hydrofoil were conducted, and two turbulence models for RANS closure were tested, i.e., the Realizable $k-{\varepsilon}$ model and the Reynolds stress transport model. Numerical results were compared with available experimental data, and it was shown that the data for the Reynolds stress transport model that were computed on the finest grid system had better agreement in reproducing the development and propagation of the tip vortex. The Realizable $k-{\varepsilon}$ model overestimated the turbulence level in the vortex core and showed a diffusive behavior of the tip vortex. The tip vortex cavitation on the hydrofoil and its trajectory also showed good agreement between the current numerical results that were obtained using the Reynolds stress transport model and the results observed in the experiment.

• Journal of Energy Engineering
• /
• v.21 no.3
• /
• pp.228-236
• /
• 2012

Large eddy simulation (LES) is increasingly used as a tool for studying the dynamics of turbulence in combustion chamber flows due to the promise of wider generality and more accurate results compared to Reynolds averaged Navier-Stokes(RANS) models. This study presents the appropriate subgrid-scale(SGS) model in LES for predicting the turbulent flow field in the internal combustion engine. The study of the effects of model and numerical parameters such as discretization scheme, initial condition, time step and SGS model was performed. The results of LES using the SGS model were found to be in the good agreement with experimental data.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.47-56
• /
• 2015

In this paper, we study the effect of various turbulence models by comparing the aerodynamic characteristics and the flow patterns computed for aircraft models. An in-house CFD solver, MSAPv, that solves the three dimensional RANS equations with the turbulence model equations is used. The turbulence models used in this study are the Spalart-Allmaras model, Menter's $k-{\omega}$ SST model, Coakley's $q-{\omega}$ model, and Huang and Coakley's $k-{\varepsilon}$ model. DLR-F6 WB and WBNP configurations are selected for the study. We concentrate on the separated flow pattern variations with the turbulence models at the wing-body junction and the wing-pylon junction as well as drag polar curves.

• Journal of computational fluids engineering
• /
• v.14 no.3
• /
• pp.115-122
• /
• 2009

Numerical methods are applied to simulate the smoke behavior in a ventilated tunnel using large eddy simulation (LES) which is incorporated in FDS (Fire Dynamics Simulator) with proper combustion and radiation model. In this study, present numerical results are compared with data obtained from experiments on pool fires in a ventilated tunnel. The model tunnel is $182m(L){\times}5.4m(W){\times}2.4m(H)$. Two fire scenarios with different ventilation rates are considered with two different fire strengths. The present results are analyzed with those from LES without combustion and radiation model and from RANS ($\kappa-\epsilon$) model as well. Temperature distributions caused by fire in tunnel are compared with each other. It is found that thermal stratification and smoke back-layer can be predicted by FDS and the temperature predictions by FDS show better results than LES without combustion and radiation model. The FDS solver, however, failed to predict correct flow pattern when the high ventilation rate is considered in tunnel because of the defects in the tunnel-inlet turbulence and the near-wall turbulence.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.3
• /
• pp.183-188
• /
• 2017

This paper provides numerical results for the estimation of the efficiency of KRISO energy saving devices in the design stage. A finite volume method is used to solve Reynolds averaged Navier-Stokes (RANS) equations, where the SST k-$\omega$ model is selected for turbulence closure. The propeller rotating motion is determined using a rigid body motion (RBM) scheme, which is called a sliding mesh technique. The numerical analysis focuses on predicting the power reduction by the designed KRISO devices (K-DUCT) under a self-propulsion condition. The present numerical results show good agreement with the available experimental data. Finally, it is concluded that CFD can be a useful method, along with model tests, for assessing the performance of energy saving devices for propulsion efficiency improvement.

• Wind and Structures
• /
• v.28 no.4
• /
• pp.239-253
• /
• 2019

Accurate numericalsimulation of wind field over complex terrain is an important prerequisite for wind resource assessment. In this study, numerical simulation of wind field over complex terrain was further carried out by taking the complex terrain around Siu Ho Wan station in Hong Kong as an example. By artificially expanding the original digital model data, Gambit and ICEM CFD software were used to create high-precision complex terrain model with high-quality meshing. The equilibrium atmospheric boundary layer simulation based on RANS turbulence model was carried out in a flat terrain domain, and the approximate inflow boundary conditions for the wind field simulation over complex terrain were established. Based on this, numerical simulations of wind field over complex terrain under different inflow wind directions were carried out. The numerical results were compared with the wind tunnel test and field measurement data for land and sea fetches. The results show that the numerical results are in good agreement with the wind tunnel data and the field measurement data which can verify the accuracy and reliability of the numerical simulation. The near ground wind field over complex terrain is complex and affected obviously by the terrain, and the wind field characteristics should be fully understood by numerical simulation when carrying out engineering application on it.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.4
• /
• pp.206-216
• /
• 2017

The temporal distribution of the turbulence kinetic energy (TKE) and the vertical component of Reynolds stresses ($-{\bar{u^{\prime}w^{\prime}}}$) was measured during one wave period under high wave energy conditions. The wave data were obtained at Hujeong Beach in the east coast of Korea at January 14~18 of 2017 when an extratropical cyclone was developed in the East Sea. Among the whole thousands of waves measured during the period, hundreds of regular waves that had with similar pattern were selected for the analysis in order to give three representing mean wave patterns using the ensemble average technique. The turbulence properties were then estimated based on the selected wave data. It is interesting to find out that $-{\bar{u^{\prime}w^{\prime}}}$ has one clear peak near the time of flow reversal while TKE has two peaks at the corresponding times of maximum cross-shore velocity magnitudes. The distinguished pattern of Reynolds stress indicates that vertical fluxes of such properties as suspended sediments may be enhanced at the time when the horizontal flow direction is reversed to disturb the flows, supporting the turbulence convection process proposed by Nielsen (1992). The characteristic patterns of turbulence properties are examined using the CADMAS-SURF Reynolds-Averaged Navier-Stokes (RANS) model. Although the model can reasonably simulate the distribution of TKE pattern, it fails to produce the $-{\bar{u^{\prime}w^{\prime}}}$ peak at the time of flow reversal, which indicates that the application of RANS model is limited in the prediction of some turbulence properties such as Reynolds stresses.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.7
• /
• pp.703-714
• /
• 2010

This paper describes the design optimization of a rotating rectangular channel with staggered arrays of pin-fins by Kriging metamodeling technique. Two non-dimensional variables, the ratio of the height to the diameter of the pin-fins and the ratio of the spacing between the pin-fins to the diameter of the pin-fins are chosen as the design variables. The objective function that is a linear combination of heat transfer and friction loss related terms with a weighting factor is selected for the optimization. To construct the Kriging model, objective function values at 20 training points generated by Latin hypercube sampling are evaluated by a three-dimensional Reynolds-averaged Navier-Stokes (RANS) analysis method with the SST turbulence model. The Kriging model predicts the objective function value that agrees well with the value calculated by the RANS analysis at the optimum point. The objective function is reduced by 11% by the optimization of the channel.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.104-110
• /
• 2008

Detached Eddy Simulation (DES) is applied to an axisymmetric base flow at supersonic mainstream. DES is a hybrid approach to modeling turbulence that combines the best features of the Reynolds-averaged Navier-Stokes RANS) and large-eddy simulation (LES) approaches. In the Reynolds-averaged mode, the model is currently based on either the Spalart-Allmaras (S-A) turbulence model. In the large eddy simulation mode, it is based on the Smagorinski subgrid scale model. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology with less computational cost than that of pure LES and monotone integrated large-eddy simulation (MILES) approaches. The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region, small-eddy motions inside the recirculating region. Comparison of the results shows that it is necessary to resolve approaching boundary layers and free shear-layer velocity profiles from the base edge correctly for the accurate prediction of base flows. The consideration of an empirical constant CDES for a compressible flow analysis may suggest that the optimal value of empirical constant CDES may be larger in the flows with strong compressibility than in incompressible flows.