• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.1-10
• /
• 1997

Viscous flow around actual ship is calculated by an use of RANS equations. The propriety of this computing method, usefulness to hull form design and the scale effect which is the effect of viscous flow depending on the scale of ship model are investigated. Reynolds stress is modelled by using k-${\varepsilon}$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the arbitrary 3-dimensional shape of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM. In the calculation of pressure, SIMPLE method is adopted and the solution of the discretized equation is obtained by the line-by-line method with the use of TDMA The calculations of two ships, 4410 TEU container carrier and 50,000 DWT class bulk carrier, are performed at model and actual ship scale. The results are compared and discussed with the model test results which are viscous resistance, nominal wake distribution at propeller plane and limiting streamline on the hull surface. They describe the effect of stem form and the scale effect very well. In particular, the calculated nominal wake distribution and limiting streamline are agreed qualitatively with the experiments and the viscous resistance values are estimated within ${\pm}5%$ difference from the resistance tests.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.1
• /
• pp.81-96
• /
• 2007

Nonlinear shoaling characteristics over surf zone are numerically investigated based on spatially averaged NavierStokes equation. We also test the validity of gradient model for turbulent stresses due to wave breaking using the data acquainted during SUPERTANK LABORATORY DATA COLLECTION PROJECT(Krauss et al., 1992). It turns out that the characteristics length scale of breaking induced current is not negligible, which firmly stands against ever popular gradient model, ${\kappa}-{\varepsilon}$ model, but favors Large Eddy Simulation with finer grid. Based on these observations, we model the residual stress of spatially averaged NavierStokes equation after Lagrangian Dynamic Smagorinsky(Meneveau et al., 1996). We numerically integrate newly proposed wave equations using SPH with Gaussian kernel function. Severely deformed water surface profile, free falling water particle, queuing splash after landing of water particle on the free surface and wave finger due to structured vortex on rear side of wave crest(Narayanaswamy and Dalrymple, 2002) are successfully duplicated in the numerical simulation of wave propagation over uniform slope beach, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.8 no.3
• /
• pp.189-199
• /
• 2010

In a high-level waste (HLW) repository, heat is generated by the radioactive decay of the waste. This can affect the safety of the repository because the surrounding environment can be changed by the heat transfer through the rock. Thus, it is important to determine the heat transfer coefficient of the atmosphere in the underground repository. In this study, the heat transfer coefficient was estimated by measuring the indoor environmental factors in the Korea Atomic Energy Research Institute Underground Research Tunnel (KURT) under forced convection. For the experiment, a heater of 5 kw capacity, 2 meters long, was inserted through the tunnel wall in the heating section of KURT in order to heat up the inside of the rock to $90^{\circ}C$, and fresh air was provided by an air supply fan connected to the outside of the tunnel. The results showed that the average air velocity in the heating section after the provision of the air from outside of the tunnel was 0.81 m/s with the Reynolds number of 310,000~340,000. The seasonal heat transfer coefficient in the heating section under forced convection was $7.68\;W/m^2{\cdot}K$ in the summer and $7.24\;W/m^2{\cdot}K$ in the winter.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.2
• /
• pp.1533-1540
• /
• 2015

In order to demonstrate the flow properties of the river bed and the design of hydraulic structures, the estimation of friction velocity is essentially required. However, existing friction velocity equations such as Log method and Power law have trouble to estimate the friction velocity because a boundary condition and various hydraulic properties are changed constantly in near the wall. In the present study, therefore, a new friction velocity equation that can minimize the parameters and reduce an error was suggested. To verify accuracy and reliability for the proposed equation, Clauser method, $\sqrt{gRI}$ method, reynolds stress method by Dr. Song were compared with the proposed method by estimated entropy parameter M for each channel. Consequently, the results show that uniform flow condition as well as non-uniform flow condition with highly accuracy nearly matched in case of accelerating non-uniform condition of $R^2=0.9621$, Decelerating Non Uniform condition of $R^2=0.9274$, Uniform condition of $R^2=0.8865$.

• Journal of the Korean Institute of Gas
• /
• v.2 no.1
• /
• pp.28-39
• /
• 1998

This study investigated the natural convection and radiation in a rectangular enclosure with ceiling vent experimentally and numerically. A heat source is located on the center of the bottom surface. The analysis was peformed a pure convection and is combination of natural convection and radiation. The shape of the considered two dimensional model is a square whose center of ceiling($30\%$) is opened. The numerical simulations are carried out for the pure natural convection case and the combined heat transfer case by using the SIMPLE algorithm. For the turbulent flow, Reynolds stresses are closed by the standard $k-{\epsilon}$ model and the wall function is used to determine the wall boundary conditions. The experiment was performed on the same geometrical shape as the computations. The radiative heat transfer is analized by the S-N discrete ordinates method. The results of pure natural convection are compared with those of combined heat transfer by the velocity vectors, stream lines, isothermal lines. The results obtained are as follows 1. Comparing the results of pure convection with those of the combined convection-radiation through the shape of stream lines, isothermal lines are similar to each other. 2. The temperature fields obtained by numerical method are compared to those obtained by experimental one, and it is found that they are showed mean relative error $8.5\%$. 3. Visualization bt smoke is similar to computational results.

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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2018.05a
• /
• pp.71-71
• /
• 2018

Sea pollution accidents have been occurred due to the increase of marine ship traffic. Oil spill from the hull hole induced by tanker collision results in the huge sea pollution. Proper and prompt reaction on such oil spill disaster is needed to minimize the damage. Thru-hull emergency wood plug is typically used to manually close small holes, while it is required to develop some mechanical devices for closing large holes in the hull due to huge fluid pressure. Accurate estimation of oil discharge and velocity from such holes are important to develop proper device to control hull hole damage. High resolution CFD modeling investigation on the configurations of hull hole of 7.5 m initial depth and 30 cm diameter, which was observed in the oil spill accident of the Hebei Sprit off the west coast of Korea in 2007, has been carried out to compute the oil spill velocity distribution in terms of flow depth. Friction loss due to the viscous flow and the discharge coefficient of crude oil with specific gravity SG = 0.85 and viscosity of $4-12cP(mPa{\cdot}s)$ at the temperature of $20^{\circ}C-100^{\circ}C$ are presented in terms of Reynolds number based on the results of high-resolution CFD modeling.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.20 no.6
• /
• pp.760-767
• /
• 2014

An experimental study is carried out to investigate the effect of jet stream in the gab of rectangular cylinders with different length in a parallel arrangement by using PIV method in a circulating water channel. The height(h) of the rectangular cylinder and the gap between the cylinder is 10mm, and the width(B) which is 300mm. The length of the model for flow direction was applied to 30mm, 60mm, 90mm & 120mm, The aspect ratio of a model on the basis of height(H=30mm) is 1, 2, 3 and 4. Reynolds number $Re=1.4{\times}10^4$, $Re=2.0{\times}10^4$, $Re=2.9{\times}10^4$ based on the height(H) of model for the distance of tidal distributions as of water depth have been applied during the whole experiments. The measurement area was set to 5H rear of the cylinder. As a result, Vortex size in the wake area were increased as velocity increased. and high aspect ratio increased through-flow velocity component in the near wake. Velocity deficit increased highly after near-wake area and low aspect ratio.

• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1834-1845
• /
• 2021

Liquid lead-bismuth cooled fast reactor is one of the most promising reactor types among the fourth-generation nuclear energy systems. The flow and heat transfer characteristics of lead-bismuth eutectic (LBE) are completely different from ordinary fluids due to its special thermal properties, causing that the traditional Reynolds analogy is no longer recommended and appropriate. More accurate turbulence flow and heat transfer model for the liquid metal lead-bismuth should be developed and applied in CFD simulation. In this paper, a specific CFD solver for simulating the flow and heat transfer of liquid lead-bismuth based on the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model was developed based on the open source platform OpenFOAM. Then the advantage of proposed model was demonstrated and validated against a set of experimental data. Finally, the simulation of LBE turbulent flow and heat transfer in a 7-pin wire-wrapped rod bundle with the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model was carried out. The influence of wire on the flow and heat transfer characteristics and the three-dimensional distribution of key thermal hydraulic parameters such as temperature, cross-flow velocity and Nusselt number were studied and presented. Compared with the traditional SED model with a constant Pr_t = 1.5 or 2.0, the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model is more accurate on predicting the turbulence flow and heat transfer of liquid lead-bismuth. The average relative error of the k - 𝜀 - k_𝜃 - 𝜀_𝜃 model is reduced by 11.1% at most under the simulation conditions in this paper. This work is meaningful for the thermal hydraulic analysis and structure design of fuel assembly in the liquid lead-bismuth cooled fast reactor.

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