• 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 Korea Water Resources Association
• /
• v.48 no.12
• /
• pp.1023-1035
• /
• 2015

This study applies 3-D N-S solver based on PBM (Porous Body Model), LED-WASS-3D ver 2.0 to directly analyze non-linear interaction of seawater-freshwater-coastal aquifer in order to simulate the seawater infiltration into coastal aquifer. This numerical simulation is the first trial in Korea, as well as unusual and new numerical analysis abroad. Firstly, to validate the applied numerical model, the validity and effectiveness was verified for the numerical model by comparing and considering it with the result of laboratory experiment for seawater-freshwater interface in coastal aquifer. And then it simulated the seawater infiltration into coastal aquifer considering the changed levels of seawater and groundwater in order to analyze the distribution characteristics of flow field and seawater-freshwater interface of coastal aquifer as the level difference between seawater and groundwater and rate of seawater level (${\Delta}h/h$) increased. In addition, the characteristics of seawater infiltration were analyzed from the vertical salinity in the coastal aquifer by ${\Delta}h/h$, which cannot be obtained from existing non-diffusion numerical models. Finally, it analyzed the effect of ${\Delta}h/h$ on the seawater infiltration distance in coastal aquifer, which was indexed.

• Fire Science and Engineering
• /
• v.33 no.5
• /
• pp.48-54
• /
• 2019

A Large Eddy Simulation (LES), adopting the Eddy Dissipation Concept (EDC) 1-step model, was successfully performed for backdraft phenomena. The activation energy of the finite chemistry reaction in the EDC 1-step model was adjusted to simulate the backdraft. The prediction of the EDC 1-step model was similar to that of the Mixing-Controlled Fast Chemistry (MCFC) model, except when the backdraft occurred. The EDC 1-step model could be used to simulate the experimental peak pressure, but not the first peak pressure of the backdraft.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.6
• /
• pp.484-497
• /
• 2016

In order to understand the characteristics of the topography change in front of an impermeable breakwater, a coupled model for a two-way analysis of the existing LES-WASS-2D and newly developed morphodynamic model was suggested. A comparison to existing experimental results revealed that the results computed using the 2-D hydro-morphodynamic model were in good agreement with the experimental results for the wave form, pore water pressure in the seabed, and topographical change in front of a submerged breakwater. It was shown that the two-way model suggested in this study is applicable to a morphological change in the seabed around a submerged breakwater. Then, using the numerical results, the topographical changes in front of an impermeable submerged breakwater were examined in relation to partial standing waves. Moreover, the characteristics of the local scour depths in front of them are also discussed in relation to incident wave conditions, sediment qualities, and submerged breakwater shapes.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2737-2742
• /
• 2007

Strong wind flow around a building complex was numerically studied by LES. The original motivation of this work stemmed from the efforts to develop a risk assessment technique for windstorm hazards. Lagrangian-averaged scale-invariant dynamic subgrid-scale model was used for turbulence modeling, and a log-law-based wall model was employed on all the solid surfaces including the ground and the surface of buildings to replace the no-slip condition. The shape of buildings was implemented on the Cartesian grid system by an immersed boundary method. Key flow quantities for the risk assessment such as mean and RMS values of pressure on the surface of the selected buildings are presented. In addition, characteristics of the velocity field at some selected locations vital to safety of human beings is also reported.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2016.02a
• /
• pp.142.1-142.1
• /
• 2016

We present bulk diffusion rates of hydroxide ions in amorphous solid water (ASW) at 135 ~ 160 K. Previous researches showed that the diffusion mechanism of hydroxide is different from one of hydronium ions, and this implies that they have different diffusion rates. In ultra-high vacuum (UHV) chamber, low-energy scattering (LES) was used to measure ion population and temperature-programmed desorption (TPD) was conducted for measuring ASW thicknesses. To determine the diffusion rates, a simple model for $H_2O/NaOH/H_2O$ sandwich films was developed using Fick's second law. The measured surface population of hydroxide ions as a function of time was well fitted to the model, and the rates were well agreed to an Arrhenius equation.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.392-397
• /
• 2001

A numerical study is performed to investigate the interaction between subsonic axial turbine blade boundary layer and periodically oncoming rotor induced wakes. An implicit scheme for solving the compressible Navier-Stokes equation is developed, which adopts a 4th-order compact difference for spatial discretiztion, a 2nd order Crank-Nicolson scheme for temporal discretization and the dynamic eddy viscosity model as the subgrid scale model. The efficiency and the accuracy of the proposed method are verified by applying to some benchmark problems such as laminar cylinder flow, laminar airfoil cascade flow and a transitional flat plate boundary layer flow. Computational results show good agreements with previous experimental and numerical results. Finally, flow through a stator cascade is simulated at $Re = 7.5{\times}10^5$ without free-stream turbulence intensity. The velocity fields and skin friction coefficients in the transitional region show similar trends with previous boundary layer natural transition.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.603-607
• /
• 2007

Recently, the frequency of unexpecting heavy rains has been increased due to abnormal climate and extreme rainfall. There was a limit to analyze one dimension or two dimension stream flow of domestic rivers that was applied simple momentum equation and fixed energy conservation. Therefore, hydrodynamics flow analysis in rivers has been needed three dimensional numerical analysis for correct stream flow interpolation. In this study, CFD model on FLOW-3D was applied to stream flow analysis, which solves three dimension RANS(Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behavior and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as $k-{\backepsilon}$, RNG $k-{\backepsilon}$ and LES. Those numerical analysis results have been illustrated to bends and junctions by the turbulence energy effects, velocity of flow distributions, water level pressure distributions and eddy flows.