• Journal of Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.

• Journal of Korea Water Resources Association
• /
• v.47 no.8
• /
• pp.703-715
• /
• 2014

This study aims to investigating the effect of Schmidt number (${\sigma}_c$) on sediment suspension and hydrodynamics calculation. The range of ${\sigma}_c$ is also studied based on the flux Richardson number ($Ri_f$) and gradient Richardson number ($Ri_g$). Numerical experiments are carried out by 1 dimensional vertical model. Both cohesive and non-cohesive sediments are tested under the conditions of pure current and oscillatory flow. The turbulence damping effect due to sediment suspension is examined considering ${\sigma}_c$ as a constant for the damping effect. The results of this study show the consistent effect of ${\sigma}_c$ on sediment suspension regardless of hydrodynamic condition. It is also found that the model overestimates the flow velocity and turbulent kinetic energy when the damping effect is not considered. Under the conditions of $Ri_f$ and $Ri_g$ causing density stratification, it is known that the vertical mixing of sediment is reasonably calculated in the range of ${\sigma}_c$ from 0.3 to 0.5.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.48-56
• /
• 2002

The flow characteristics around KRISO 3600TEU container ship model have been experimentally investigated in a circulating water channel. The instantaneous velocity vectors were measured using 2-frame PIV measurement system. The mean velocity fields and turbulent statistics including turbulent kinetic energy and vorticity were obtained by ensemble-averaging 400 instantaneous velocity fields. The free stream velocity was fixed at 0.6m/s and the corresponding Reynolds number was $9{\times}10^5$. The test sections were divided into two regions, three transverse sections of the wake region(Station -0.5767, -1, -3) and five longitudinal sections of the wake((Z/(B/2)=0, 0.1, 0.2, 0.4, 0.6). In the wake region, large-scale longitudinal vortices of nearly same strength are symmetric with respect to the wake centerline and a relatively weak secondary vortex is formed near the waterline. With going downstream, the strength of longitudinal vortex is decreased and the wake region expands.

• Journal of Korea Water Resources Association
• /
• v.57 no.8
• /
• pp.493-507
• /
• 2024

RANS-based CFD analysis is widely applied in various engineering fields, including practical hydraulic engineering, due to its high computational efficiency. However, problems of non-physical behavior in the analysis of two phase flow, such as free surfaces, have long been raised. The two-equation turbulence models used in general RANS-based analysis were developed for single phase flow and simulate unrealistically high turbulence energy at the interface where there are abrupt changes in fluid density. To solve this issue, one of the methods recently developed is the buoyancy-modified turbulence model, which has been partially validated in coastal engineering, but has not been applied to open channel flows. In this study, the applicability of the buoyancy-modified turbulence model is evaluated using the VOF method in the open-source program OpenFoam. The results of the uniform flow showed that both the buoyancy-modified k-𝜖 model and the buoyancy-modified k-ω SST model effectively simulated the reduction of turbulence energy near the free surface. Specifically, the buoyancy-modified k-ω SST model accurately simulated the vertical velocity distribution. Additionally, the model is applied to dam-break flows to examine cases with significant surface variation and cavity formation. The simulation results show that the buoyancy-modified turbulence models produce varying results depending on the VOF method and shows non-physical behavior different from experimental results. While the buoyancy-modified turbulence model is applicable in cases with stable surface shapes, it still has limitations in general application when there are rapid changes in the free surface. It is concluded that appropriate adjustments to the turbulence model are necessary for flows with rapid surface changes or cavity formation.

• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• The Journal of the Acoustical Society of Korea
• /
• v.34 no.1
• /
• pp.36-45
• /
• 2015

Recently, a lot of studies have been made about the methods used to generate turbulent velocity fields stochastically in order to effectively predict broadband flow noise. Among them, the FRPM (Fast Random Particle Mesh) method which generates turbulence with specific statistical properties using turbulence kinetic energy and dissipation obtained from the steady solution of the RANS (Reynolds Averaged Navier-Stokes) equations has been successfully applied. However, the FRPM method cannot be applied to the flow noise problems involving intrinsic unsteady characteristics such as centrifugal fan. In this paper, to effectively predict the broadband noise generated by centrifugal fan, U-FRPM (unsteady FRPM) method is developed by extending the FRPM method to be combined with the unsteady numerical solutions of the unsteady RANS equations to generate the turbulence considered as broadband noise sources. Firstly, an unsteady flow field is obtained from the unsteady RANS equations through CFD (Computational Fluid Dynamics). Then, noise sources are generated using the U-FRPM method combined with acoustic analogy. Finally, the linear propagation model which is realized through BEM (Boundary Element Method) is combined with the generated sources to predict broadband noise at the listeners' position. The proposed technique is validated to compare its prediction result with the measured data.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.7
• /
• pp.930-938
• /
• 2018

Underwater noise radiated from submarines is directly related to the probability of being detected by the sonar of an enemy vessel. Therefore, minimizing the noise of a submarine is essential for improving survival outcomes. For modern submarines, as the speed and size of a submarine increase and noise reduction technology is developed, interest in flow noise around the hull has been increasing. In this study, a noise analysis technique was developed to predict flow noise generated around a submarine shape considering the free surface effect. When a submarine is operated near a free surface, turbulence-induced noise due to the turbulence of the flow and bubble noise from breaking waves arise. First, to analyze the flow around a submarine, VOF-based incompressible two-phase flow analysis was performed to derive flow field data and the shape of the free surface around the submarine. Turbulence-induced noise was analyzed by applying permeable FW-H, which is an acoustic analogy technique. Bubble noise was derived through a noise model for breaking waves based on the turbulent kinetic energy distribution results obtained from the CFD results. The analysis method developed was verified by comparison with experimental results for a submarine model measured in a Large Cavitation Tunnel (LCT).

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.11 no.4
• /
• pp.192-205
• /
• 2009

The unexpected wind over the Mt. Hwawang on 9 February 2009 was deadly when many spectators were watching a traditional event to burn dried grasses and the fire went out of control due to the wind. We analyzed the fatal wind based on wind flow simulations over a digitized complex terrain of the mountain with a localized heating area using a three dimensional computational fluid dynamics model, CFD_NIMR_SNU (Computational Fluid Dynamics_National Institute of Meteorological Research_Seoul National University). Three levels of fire intensity were simulated: no fire, $300^{\circ}C$ and $600^{\circ}C$ of surface temperature at the site on fire. The surface heat accelerated vertical wind speed by as much as $0.7\;m\;s^{-1}$ (for $300^{\circ}C$) and $1.1\;m\;s^{-1}$ (for $600^{\circ}C$) at the center of the fire. Turbulent kinetic energy was increased by the heat itself and by the increased mechanical force, which in turn was generated by the thermal convection. The heating together with the complex terrain and strong boundary wind induced the unexpected high wind conditions with turbulence at the mountain. The CFD_NIMR_SNU model provided valuable analysis data to understand the consequences of the fatal mountain fire. It is suggested that the place of fire was calm at the time of the fire setting due to the elevated terrain of the windward side. The suppression of wind was easily reversed when there was fire, which caused updraft of hot air by the fire and the strong boundary wind. The strong boundary wind in conjunction with the fire event caused the strong turbulence, resulting in many fire casualties. The model can be utilized in turbulence forecasting over a small area due to surface fire in conjunction with a mesoscale weather model to help fire prevention at the field.