• Water for future
• /
• v.29 no.2
• /
• pp.199-207
• /
• 1996

To evaluate the usage of $\kappa-\imath$ turbulence closure for the analysis of thermal discharge behavior, a two-dimensional depth-integrated numerical model is developed. The developed model is applied to a steady flow in an open channel with simle geometry and the numerical results agree well with existing experimental data. The adequate simulation of recirculation, reattachment, and excess temperature rise at downstream of the outlet in the channel attributes to the correct calculation of turbulent eddy viscosity and diffusivity by $\kappa-\imath$ turbulence model. For an accurate prediction of thermal discharge behavior, the introduction of buoyancy production term, the modification of source/sink, and the correct input of turbulence constants of the $\kappa-\imath$ turbulence model are required.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.4
• /
• pp.214-224
• /
• 2021

Recently, as the technology of the supercavitating underwater vehicle is improved, the necessity of research for maneuvering characteristics of the supercavitating underwater vehicle has emerged. In this study, as a preliminary step to analyzing the maneuverability of a supercavitating underwater vehicle, the characteristics of cavity shapes and hydrodynamic forces generated in a supercavitating underwater vehicle with an angle of attack were evaluated numerically. First, the geometry was designed by modifying the shape of the existing supercavitating underwater vehicle. The continuity and the Navier-stokes equations are numerically solved, and turbulent eddy viscosity is solved by the k-ω SST model. The results present the characteristics of cavity shape and the hydrodynamic forces of the designed geometry with an angle of attack.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.3
• /
• pp.357-366
• /
• 1986

The flow characteristics of two-dimensional turbulent offset jet which is discharged parallel to a solid wall has been studied experimentally and numerically. In the experiment, 3-hole pitot tube and 2 channel constant temperature hot-wire anemometer are used to measure local mean velocity, turbulence intensity and Reynolds stress while scannivalve is used to measure the wall pressure distribution. It is confirmed experimentally that local mean velocity is closely related to wall pressure distribution. It is also verified that for large Reynolds numbers and fixed step height there exists a similarity in the distribution of wall pressure coefficient. The maximum values of turbulence intensity occur in the top and bottom mixing layers and the magnitude of Reynolds stress becomes large in the lower mixing layer than in the top mixing layer due to the effect of streamline curvature and entrainment. In the numerical analysis, standard k-.epsilon. model based on eddy viscosity model and Leschziner and Rodi model based on algebraic stress model are adopted. The numerical analyses predict shorter reattachment lengths than the experiment, and this difference is judged to be due mainly to the problem of turbulence model constants and numerical algorithm. This also causes the inconsistency between the two results for other turbulence quantities in the recirculation region and impingement region, which constitutes a subject of a continued future study.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.363-371
• /
• 2008

The aim of this study is to compare the numerical results obtained by 2-D and 3-D models which are used to examine the wave field around a permeable submerged breakwater. At first, the numerical model, which is able to consider the flow through a porous medium with inertial, laminar and turbulent resistance terms and determine the eddy viscosity with LES turbulent model, is used and validated by comparing with existing experimental data. And then, the numerical test on the wave field around a permeable submerged breakwater is performed. It is revealed from the numerical results that, at the onshore side of the submerged breakwater, the wave height by 2-D analysis is higher than that by 3-D analysis. Also, the time-averaged mean flow around a submerged breakwater is discussed in detail.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.3B
• /
• pp.325-335
• /
• 2010

The aim of this study is to survey the effects of laying drainage layer in sandy beach on beach stabilization. At first, the numerical model developed by Hur and Lee (2007), which is able to consider the flow through a porous medium with inertia, laminar and turbulent resistance terms, i.e. simulate directly WAve Structure Seabed/Sandy beach interaction and can determine the eddy viscosity with LES turbulent model in 3-D wave field (LES-WASS-3D), is validated by comparing with existing experimental data. And then, numerical simulation is carried out to examine the characteristics of wave-sandy beach interaction for a beach with/without drainage layer. From the numerical results, it is shown that mean ground-water level around a foreshore decreases and offshore-ward flow over a seabed reduces in case of a beach with drainage layer. Moreover, the effects of cross profile of drainage layer and incident wave condition on mean ground-water level around a foreshore are also discussed as well the distribution of wave setup around the foreshore.

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

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.20-28
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic mode1. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach arc used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vortex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the k-w SST turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.5B
• /
• pp.483-489
• /
• 2011

In this study, the coupled-numerical model has been newly developed to investigate numerically scouring and deposition around a coastal structure like a submerged breakwater using a numerical wave model and a lagrangian particle model for sand transport. As a numerical wave model, LES-WASS-2D (Hur and Choi, 2008) is adopted. The model is able to consider the flow through a porous midium with inertial, laminar and turbulent resistance term and determine the eddy viscosity with LES turbulence model. Distinct element method (Cundall and Strack, 1979), which is able to apply to many dynamical analysis of particulate media, as a lagrangian particle model for sand transport is newly coupled to the numerical wave model. The numerical simulation has been carried out to examine the scour problem in front of an impermeable submerged breakwater using the newly coupled-numerical model. The numerical results has been compared qualitatively with an existing experimental data and then its applicability has been discussed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.2B
• /
• pp.249-259
• /
• 2008

The present paper studies the effect of the slope gradient of a fully permeable submerged breakwater using a newly developed numerical model that is able to consider the flow through a porous midium with inertial, laminar and turbulent resistance terms, i.e. simulate directly WAve-Structure (submerged breakwater)-Sand seabed interaction and can determine the eddy viscosity with LES turbulence model in 2-Dimensional wave field (LES-WASS-2D). The developed model was validated through the comparison with an existing experimental data, and further used for various numerical experiments in oder to investigate the complicated hydrodynamics on the varying slope gradient of permeable submerged breakwater. We found an acceptable phenomenon, as we expect intuitively, that reflection and transmission coefficients decrease simultaneously as slope gradient decrease. In addition, the breaking point, the circulation flow and mean vorticity around a submerged breakwater are throughly discussed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.575-589
• /
• 2008

In the present work, we investigate the hydrodynamic behavior of a turbulent bore, such as tsunami bore and tidal bore, generated by the removal of a gate with water impounded on one side. The bore generation system is similar to that used in a general dam-break problem. In order to the numerical simulation of the formation and propagation of a bore, we consider the incompressible flows of two immiscible fluids, liquid and gas, governed by the Navier-Stokes equations. The interface tracking between two fluids is achieved by the volume-of-fluid (VOF) technique and the M-type cubic interpolated propagation (MCIP) scheme is used to solve the Navier-Stokes equations. The MCIP method is a low diffusive and stable scheme and is generally extended the original one-dimensional CIP to higher dimensions, using a fractional step technique. Further, large eddy simulation (LES) closure scheme, a cost-effective approach to turbulence simulation, is used to predict the evolution of quantities associated with turbulence. In order to verify the applicability of the developed numerical model to the bore simulation, laboratory experiments are performed in a wave tank. Comparisons are made between the numerical results by the present model and the experimental data and good agreement is achieved.