• Journal of Ocean Engineering and Technology
• /
• v.33 no.3
• /
• pp.280-288
• /
• 2019

A submerged breakwater is one of the coastal structures used to reduce wave energy and coastal erosion. However, a submerged breakwater has a negative aspect in that a strong rip current occurring around an open inlet due to a difference in mean water levels at the front and rear sides of the structure leads to scouring. Such scouring has a bad effect on its stability. In order to eliminate this kind of demerit, this study investigated an artificial reef of the overflow type with openings. We also developed a program where the flows around the artificial reef of the overflow type could be analyzed numerically. An unstructured grid system was used to cover the various geometries, and the level set method was applied to treat the movement of the free surface. To verify these numerical schemes, hydraulic physical tests were performed on the submerged breakwater and double breaking type artificial reef. Then, the wave height and velocity distribution around the reef were examined using the experimental results. Comparisons between the results of hydraulic and numerical tests showed reasonable agreement.

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

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.20 no.1
• /
• pp.37-44
• /
• 2017

Complexity of multiphase flows due to existence of more than two interface including free-surface in one system, cannot be simulated easily. Since more than two fluids affect to flows and disturb interface, non-linearities such as instabilities can be appeared. Among several instabilities on multiphase flows, one of representative is Rayleigh-taylor instability. In order to examine in importance of density disparity, several cases with numerous Atwood number are set. Moreover, investigation of influence on initial disturbance were also considered. Moving particle simulation (MPS) method, which was employed in this paper, was not widely used for multiphase problem. In this study, by adding new particle interaction models such as self-buoyance correction, surface tension, and boundary condition at interface models, MPS were developed having more strength of physics and robust. By applying newly developed multiphase MPS, considered cases are performed and compared each other. Additionally, though existence of disagreement of magnitude of rising velocity between theoretical values from linear potential theory and that of numerical simulation, agreement of tendency can be proved of similarity of result. the discordance of magnitude can be explained due to non-linear effects on numerical simulation which was not considered in theoretical result.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.725-733
• /
• 1990

An analysis of laminar mixed convection flow adjacent to the inclined flat surface which is subjected to a uniform temperature in a uniform free stream is performed. Nonsimilar boundary layed equation are derived by using the mixed convection parameters such that smooth transition from the purely forced convection limit to the purely free convection limit is possible. The governing equations are solved by a finite difference method using the coupled box scheme of sixth order. Numerical results are presented for prandtl numbers of 0.7 and 7 with the angle of inclination ranging from 0 to 90 degree from the vertical. The velocity distributions for the buoyancy assisting flow exhibit a significant overshoot above the free stream value in the region of intense mixed convection and the velocity field is found to be more sensitive to the buoyancy effect than the temperature field. The separation point near the wall was obtained for the buoyancy opposing flow. The local Nusselt number increases for buoyancy assisting flow and decreases for opposing flow with increasing value of the local Grashoff number in the mixed convection parameter. For large Prandtl number, the Nusselt number and the friction factor decrease significantly near the separation point. Present numerical predictions are in good agreement with recent experimental results by Ramachandran.

• Journal of Korea Water Resources Association
• /
• v.43 no.10
• /
• pp.911-920
• /
• 2010

The closed water bodies, such as reservoirs and lakes, could be contaminated by an inflow of pollutants in the upstream as well as a stratification caused by seasonal natural phenomena. The vertical circulation particularly plays an important role in reduction of environmental pollutants. The factors of the vertical circulation are the temperature, wind, thermal diffusivity and sunlight. The wind is probably the most significant factor among them. Thus, it is necessary to describe the validation and application of a three-dimensional numerical model of wind-induced circulation in a thermally stratified flow. In this paper, a three-dimensional numerical model for the thermally stratified flows is presented. The model is conducted in three steps to calculate the velocity components from the momentum equations in x- and y- axis directions, the elevations from the free surface equation and the temperature from the scalar transport equation. Numerical predictions are compared with available analytical solutions for the sloshing free surface movement in a rectangular basin. The numerical results generally show a reasonable agreement with analytical solutions. And the model is applied to the circulation for the wind induced flow in a thermally stratification. Consequently, the developed model is validated by two verifications and phenomena of the internal flow.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.12
• /
• pp.1017-1026
• /
• 2016

The present paper describes the two off-surface visualization techniques and their application examples to vortex flows. One of the two visualization techniques is the classical smoke-wire technique, and the other is the visualization technique using the micro water-droplets generated by the home-style ultrasonic humidifier. The smoke-wire technique has the limit of air flow speed (about 5 m/sec for 0.07 mm-diameter wire) and the pollution problem, but it produces very fine and clear streak line sheet. It is applied to visualize the wing-tip vortices of a 3-dimensional wing. The micro water-droplet technique has the larger limit of air flow speed (above 10 m/sec) and is free from pollution and toxic problems compared to the smoke-wire technique. It is successfully applied to visualize the complex vortex system of a double-delta wing with an apex strake.

• Journal of the korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.21-31
• /
• 1980

Two phase boundary layer equations of laminar filmwise condensation are solved by an approximate integral method under the following condition; saturated vapour flows vertically downward over a cooled surface of uniform temperature, the condensate film is so thin that the inertia and convection terms are neglected. The following conclusions are drawn under the above assumptions. 1. free convection In case of the linear temperature profile in a liquid film, numerical results for the average coefficients of heat transfer may be expressed as N $u_{m}$=4/3,(G $r_{l}$ /4.H)$^{1}$4/ and in case of the quadratic profile, numerical results may be expressed as N $u_{m}$=2/1.682,(G $r_{l}$ /H)$^{1}$4/. 2. Forced convection When the temperature profile is assumed to be linear in a liquid film, numerical results fir the average heat transfer coefficients may be expressed as N $u_{m}$=(A, R $e_{l}$ /H)$^{1}$2/. This expression is compared with the experimental results hitherto reported; For theoretical Nusselt number (N $u_{m}$)$_{th}$<2*10$^{4}$, the experimental Nusselt number (N $u_{m}$)$_{exp}$ is on the average larger than theoretical Nusselt number (N $u_{m}$)$_{th}$ by 30%. For (N $u_{m}$)$_{th}$>2*10$^{4}$, experimental Nusselt number (N $u_{m}$)$_{exp}$ is about 1.6 times as large as theoretical Nusselt number (N $u_{m}$)$_{th}$. These large deviation may be caused by the presence of turbulence in the liquid film. In case of the quadratic temperature profile in a liquid film, numerical results for the average coefficients of heat transfer may be expressed as N $u_{m}$'=(2,A,Re/H)$^{1}$2/. This formular shows that theoretical Nusselt number (N $u_{m}$)$_{th}$ is larger than experimental Nusselt number (N $u_{m}$)$_{exp}$ by 60%. It is speculated that when the temperature difference between cooled surface and saturated vapour is small, temperature profile in a liquid film is quadratic.quadratic.. quadratic.quadratic..atic..

• Journal of Ocean Engineering and Technology
• /
• v.31 no.6
• /
• pp.388-396
• /
• 2017

A Lagrangian approach based computational fluid dynamics (CFD) was used to simulate large and/or sharp deformations and fragmentations of interfaces, including free surfaces, through tracing each particle with physical quantities. According to the concept of the particle-based CFD method, it is possible to apply it to both fluid particles and solid particles such as sand, gravel, and rock. However, the presence of more than two different phases in the same domain can make it complicated to calculate the interaction between different phases. In order to solve multiphase problems, particle interaction models for multiphase problems, including surface tension, buoyancy-correction, and interface boundary condition models, were newly adopted into the moving particle semi-implicit (MPS) method. The newly developed MPS method was used to simulate a typical validation problem involving dam breaking. Because the soil and other particles, excluding the water, may have different viscosities, various viscosity coefficients were applied in the simulations for validation. The newly developed and validated MPS method was used to simulate the mobile beds induced by broken dam flows. The effects of the viscosity on soil particles were also investigated.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.7 no.4
• /
• pp.536-547
• /
• 1997

For various angular velocities of crucible and crystal, the characteristics of melt flows, temperatures and concentrations of oxygen are numerically studied in the Czochralski furnace with a uniform axial magnetic field. Buoyancy effect due to the heating of crucible wall and thermocapillary effect due to the temperature gradient at the free surface, can be differentiably suppressed by the centrifugal forces due to the rotations of the crucible and crystal. The most important factor which yields the centrifugal forces is the rotation velocity of the crucible, that influences the fields of velocities, temperatures and concentrations. In the case that the crucible rotation velocity is not high, the rotations of the crystal gives rise to the centrifugal forces effectively.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.20-26
• /
• 2015

This paper introduces a computational method for analysis of the 6-DOF motions of a ship in waves using an overset grid technique which consists of inner and outer domains for representing body motions and numerical wave tank, respectively. High order interpolation scheme is employed to increase numerical accuracy over the interface where physical values, such as velocities and pressure, interact between the inner and outer domains. The numerical schemes and algorithm are addressed in the present paper. An application to motion of KCS container carrier in head waves is presented, and the comparison of responses on heave and pitch motions shows good agreement with those of model tests.