• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.1
• /
• pp.44-51
• /
• 1996

The accurate calculation of run-up heights of long waves along the coastline is important in the view of engineering. In this paper the run-up heights of long waves are estimated by using the cnoidal wave theory which also covers both sinusoidal and solitary waves. However, the generation and the calculation of run-up heights of cnoidal waves are difficult both in laboratory and numerical experiments. In this study, the maximum run-up heights of cnoidal waves on steep slopes are computed by using the boundary integral equation model. It has been shown that the run-up heights of cnoidal waves are less than those of solitary waves, while they are larger than those of sinusoidal waves having the same wavelengths and heights. The variation of run-up heights of cnoidal waves is not a monotonic function of the wavelength. However, the run-up heights of cnoidal waves asymptotically approach that of a solitary wave as the wavelength approaches infinity. The calculated run-up heights agreed reasonably with experimental data.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.4
• /
• pp.390-403
• /
• 2017

Smoothed Particle Hydrodynamics (SPH) method has a good adaptability for the simulation of free surface flow problems. There are two forms of SPH. One is weak compressible SPH and the other one is incompressible SPH (ISPH). Compared with the former one, ISPH method performs better in many cases. ISPH based on Rankine source solution can perform better than traditional ISPH, as it can use larger stepping length by avoiding the second order derivative in pressure Poisson equation. However, ISPH_R method needs to solve the sparse linear matrix for pressure Poisson equation, which is one of the most expensive parts during one time stepping calculation. Iterative methods are normally used for solving Poisson equation with large particle numbers. However, there are many iterative methods available and the question for using which one is still open. In this paper, three iterative methods, CGS, Bi-CGstab and GMRES are compared, which are suitable and typical for large unsymmetrical sparse matrix solutions. According to the numerical tests on different cases, still water test, dam breaking, violent tank sloshing, solitary wave slamming, the GMRES method is more efficient than CGS and Bi-CGstab for ISPH method.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.10a
• /
• pp.450-455
• /
• 1995

A three-dimensional numerical tool is developed to calculate the potential distribution, electric field, and conductance for any types of conductance probes immersed in the wavy liquid film with various shapes of its free surface. The tool is validated against various analytical solutions. It is applied to find out the characteristics of the wire-wire probe, the flush-wire probe and the flush-flush probe in terms of resolution, linearity, and sensitivity. The wire-wire probe shows high resolution and excellent linearity for various film thickness, but comparably low sensitivity for low film thickness fixed. The flush-wire probe shows good linearity and high sensitivity for varying film thickness, but resolution degrading with an increase in film thickness. In order to check the applicability of the three types of probes in the real situation, the Korteweg-de Vries(KdV) two-dimensional solitary wave is simulated. The wire-wire probe is strongly affected by the installation direction of the two wires; when the wires are installed perpendicularly to the flow direction, the wire-wire probe shows large distortion of the solitary wave. In order to measure the transverse profile of waves, the wire-wire probes and the flush-wire probes are required to be separately installed 2mm and 2mm, respectively.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.12 no.3
• /
• pp.133-141
• /
• 2007

A finite element model is developed for the extended Boussinesq equations that is capable of simulating the dynamics of long and short waves. Galerkin weighted residual method and the introduction of auxiliary variables for 3rd spatial derivative terms in the governing equations are used for the model development. The Adams-Bashforth-Moulton Predictor Corrector scheme is used as a time integration scheme for the extended Boussinesq finite element model so that the truncation error would not produce any non-physical dispersion or dissipation. This developed model is applied to the problems of solitary wave propagation. Predicted results is compared to available analytical solutions and laboratory measurements. A good agreement is observed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.1
• /
• pp.79-92
• /
• 2011

In this study, wave pressure of short-period gravity waves and tsunami acting on the upright section of the horizontal-slit type caisson placed on the impermeable or permeable seabed, which is a well-known permeable breakwater with a good wave controlling ability, are investigated via numerical simulations. Further, the permeable seabed was modeled as the porous media with porosity of 0.4. Using the numerical results, the effects of the seabed conditions on the wave pressure on the front wall and inside wall of the chamber have been studied. In the numerical simulations, short-period gravity waves and tsunami(solitary wave or bore) with the same amplitude to the gravity wave are considered. A numerical wave tank is used, which is able to consider a gas-liquid two-phase flow in the same calculation zone. Numerical results show that the wave pressure of the tsunami was 3~5 times higher than the short-period gravity waves acting on the front wall and it was 2~4 times higher than the short-period gravity waves acting on the inner wall.

• Journal of Korea Water Resources Association
• /
• v.50 no.9
• /
• pp.647-652
• /
• 2017

Many islands are scattered around the southern area of the Korean Peninsula and they may be very vulnerable to unexpected tsunami attacks. During the East Japan Tsunami Event occurred on March 11, 2011, many islands located at the southern area were affected by tsunamis. In this study, maximum run-up heights of solitary waves on a circular island with asymmetrical crest lengths investigated by using a numerical model based on the shallow-water theory. The obtained results could be used by local authorities to establish a defense plan against unexpected tsunami invasion.

• Journal of applied mathematics & informatics
• /
• v.28 no.1_2
• /
• pp.383-395
• /
• 2010

In the present paper, we construct the traveling wave solutions involving parameters of nonlinear evolution equations in the mathematical physics via the (3+1)- dimensional potential- YTSF equation, the (3+1)- dimensional generalized shallow water equation, the (3+1)- dimensional Kadomtsev- Petviashvili equation, the (3+1)- dimensional modified KdV-Zakharov- Kuznetsev equation and the (3+1)- dimensional Jimbo-Miwa equation by using a simple method which is called the ($\frac{G'}{G}$)- expansion method, where $G\;=\;G(\xi)$ satisfies a second order linear ordinary differential equation. When the parameters are taken special values, the solitary waves are derived from the travelling waves. The travelling wave solutions are expressed by hyperbolic, trigonometric and rational functions.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.174-180
• /
• 2002

A virtual reality technology for multipurpose numerical simulation is developed to reproduce and investigate a variety of ocean environmental problems in a 3D-Numerical Wave Tank. The governing equations for solving incompressible fluid motion are Navier-Stokes equation and continuity equation, and the Marker-Density function technique is adopted to implement the fully-nonlinear free-surface kinematic condition. The marine environmental situations, i.e. waves, currents, wind, etc., are reproduced by use of multi-segmented wavemaker on the basis of the so-called "snake-principle". In this paper, some numerical reproduction techniques for regular and irregular waves, multi-directional waves, Bull's-eye wave, wave-current, and solitary wave are presented, and a model test in motion with large amplitude of roll angle is conducted in the developed 3D-NWT, using a overlaid grid system.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.4
• /
• pp.202-216
• /
• 2015

In the narrow water channel, which has been frequently deployed in the artificial canal in the South Korea due to the lack of available land, solitary wave type ship induced waves can occur. In order to test this hypothetical view, we carried out the numerical simulation. Numerical model consists of Navier-Stokes Equations and VOF, and the verification is implemented using the data by PIANC (1987) and the analytical model derived in this study. It was shown that numerically simulated front wave height are much larger than the one by PIANC (1987), and the fluctuation of free surface near the channel bank persists much longer (around 20s). For the case of stern waves, numerically simulated wave height are somewhat smaller than the data by PIANC (1987). These results seriously deviates from the general characteristics of ship induced waves observed in the wide water channels, and leads us to conclude that ship induced waves is severely affected by the width of water channel. It was also shown that the currents from the channel banks toward a ship, and currents from the ship toward the channel banks are alternatively occurring due to reflection at the channel banks. The velocity of currents reaches its maximum at 0.90 m/s, and these values are sustained through the entire depth. which implies that severe scourings at the channel bottom can be underway.

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

Internal solitary waves between two-stratified fluid layers are examined. Hamilton's principle is used to derive a governing equation. The wave speed-wave height relation is obtained and this agrees well with the experiment. The wave profiles obtained also agree well with the experimental result. There exist minimum and maximum wave speeds depending on the depth ratio and the density ratio of two layers. The KdV equation overpredicts the wave speed for a given wave height and underpredict the wave height for a given wave speed.