• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

• The Journal of the Acoustical Society of Korea
• /
• v.40 no.4
• /
• pp.320-329
• /
• 2021

In this paper, an empirical model of air bubble size to be applied to an air masker for reduction of underwater radiation noise is presented. The proposed model improves the divergence problem under the low-speed flow condition of the existing model derived using Rayleigh's jet instability model and simple continuity condition by introducing a jet flow velocity of air. The jet flow velocity of air is estimated using the bubble size where the liquid is quiescent. In a medium without flow, the size of the bubble is estimated by an empirical method where bubble formation regime is divided into a laminar-flow range, a transition range, and a turbulent-flow range based on the Reynolds number of the injected air. The proposed bubble size model is confirmed to be in good agreement with the Computational Fluid Dynamics (CFD) analysis result and the experimental results of the existing literature. Using the acoustic inversion method, the air bubble population is estimated from the insertion loss measured during the air injection experiment of the air- masker model in a large cavitation tunnel. The results of the experiments and the bubble size model are compared in the paper.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.7
• /
• pp.922-930
• /
• 2020

Marine accidents involving fishing boats, caused by a loss of stability, have been increasing over the last decade. One of the main reasons for these accidents is a sudden wind attacks. In this regard, the wind loads acting on the ship hull need to be estimated accurately for safety assessments of the motion and maneuverability of the ship. Therefore, this study aims to develop a computational model for the inlet boundary condition and to numerically estimate the wind load acting on a fishing boat. In particular, wind loads acting on a fishing boat at the wind speed profile boundary condition were compared with the numerical results obtained under uniform wind speed. The wind loads were estimated at intervals of 15° over the range of 0° to 180°, and i.e., a total of 13 cases. Furthermore, a numerical mesh model was developed based on the results of the mesh dependency test. The numerical analysis was performed using the RANS-based commercial solver STAR-CCM+ (ver. 13.06) with the k-ω turbulent model in the steady state. The wind loads for surge, sway, and heave motions were reduced by 39.5 %, 41.6 %, and 46.1 % and roll, pitch, and yaw motions were 48.2 %, 50.6 %, and 36.5 %, respectively, as compared with the values under uniform wind speed. It was confirmed that the developed inlet boundary condition describing the wind speed gradient with respect to height features higher accuracy than the boundary condition of uniform wind speed. The insights obtained in this study can be useful for the development of a numerical computation method for ships.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.27 no.1
• /
• pp.145-152
• /
• 2021

The purpose of this study was to determine the optimal distances between pipes to minimize the pressure loss and turbulent intensity. This was accomplished by investigating the distances between sleeve-jointed pipes and the flow changes in pipes based on variations in the Reynolds (Re) number when installing adjusting pieces for the pipes. When the thickness tp of the sleeve-jointed piping was fixed at 5 mm and the pipe lengths Lp were 10, 50, 100, and 200 mm, the correlations with the velocity of the sleeve-jointed part, pressure distribution, length of the reattachment point in the recirculation area, and Re number were analyzed. The flow characteristic of the sleeve-jointed part from a laminar to a turbulent flow region was determined by setting the Re range to 200 ≤ Re ≤ 5,000. This was done by utilizing Ansys Fluent 18.1, which is a commercial program. The enlargement and contraction ratios of the sleeve-jointed part were 1.2 and 0.83, respectively, and the turbulent intensity of the sleeve downstream edge and pressure change both increased as the Re number increased while Lp remained constant. The fact that the flow on the sleeve wall surface was disturbed by tp resulted in losses in velocity energy. Therefore, the edge of the sleeve-jointed part was also effected. When Lp was 10 mm or less, the turbulent intensity of the edge part did not change significantly as the Re number increased. The reattachment point in the recirculation area did not appear at Lp of 10 mm or less and was not affected by the vortex. In the case of 3,000 ≤ Re, the reattachment length of the wall surface of the sleeve-jointed part was nearly constant as Lp increased.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.7
• /
• pp.1244-1251
• /
• 2022

Marine accidents due to a loss of stability, have been gradually increasing over the last decade. Measures must be taken on the roll reduction of a ship. Amongst the measures, building an anti-roll tank in a ship is recognized as the most simple and effective way to reduce the roll motion. Therefore, this study aims to develop a computational model for a U-tube type anti-roll tank and to validate it by experiment. In particular, to validate the developed computational model, the height of the free surface in the tank was measured in the experiment. To develop a computational model, the mesh dependency test was carried out. Further, the effects of a turbulence model, time step size, and the number of iterations on the numerical solution were analyzed. In summary, a U-tube type anti-roll tank simulation had to be performed accurately with conditions of a realizable k-𝜖 turbulence model, 10^-2s time step size, and 15 iterations. In validation, the two cases of measured data from the experiment were compared with the numerical results. In the present study, STAR-CCM+ (ver. 17.02), a RANS-based commercial solver was used.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.643-651
• /
• 2008

In the present paper, turbulent open-channel flows over longitudinal bedforms are numerically simulated. The Reynolds- averaged Navier-Stokes equations in curvilinear coordinates are solved with the non-linear $k-{\varepsilon}$ model by Speziale( 1987). First, the developed model is applied to rectangular open channel flows for purposes of model validation and parameter sensitivity studies. It is found that the parameters $C_D$ and $C_E$ are important to the intensity of secondary currents and the level of turbulent anisotropy, respectively. It is found that the non-linear $k-{\varepsilon}$ model can hardly reproduce the turbulence anisotropy near the free surface. However, the overall pattern of the secondary currents by the present model is seen to coincide with measured data. Then, numerical simulations of turbulent flows over longitudinal bedforms are performed, and the simulated results are compared with the experimental data in the literature. The simulated secondary currents clearly show upflows and downflows over the ridges and troughs, respectively. The numerical results of secondary currents, streamwise mean velocity, and turbulence structures compare favorably with the measured data. However, it is observed that the secondary currents towards the troughs were significantly weak compared with the measured data.