• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.6
• /
• pp.263-269
• /
• 2018

Changes in water depth caused by underwater earthquakes and landslides cause sea surface undulations, which in turn propagate to the coast and result in significant damage as wave heights normally increase due to the wave shoaling process. Various types of numerical models have been developed to simulate the generation and propagation of tsunami waves. Most of tsunami models determine the initial surface of the water based on the assumption that the movement of the seabed is immediately and identically transmitted to the sea surface. However, this approach does not take into account the characteristics of underwater earthquakes that occur with time history and spatial variation. Thus, such an incomplete description on the initial generation of tsunami waves is totally reflected in the error during the simulation. In this study, the analytical solution proposed by Hammack (1973) was applied in the tsunami model in order to simulate the generation of initial water surface elevation by the change of water depth with time history and its propagation. The developed solution is expected to identify the relationship among various type of seabed motions, initial surface undulations, and wave speeds of elevated water surfaces.

• Ocean and Polar Research
• /
• v.29 no.1
• /
• pp.9-31
• /
• 2007

Most seismic sea waves in the East Sea originate from earthquakes occurring near the Japanese west coast. While the waves propagate in the East Sea, they are deformed by refraction, diffraction and scattering. Though the Boussinesq equation is most applicable for such wave phenomena, it was not used in numerical modelling of seismic sea waves in the East Sea. To examine characteristics of seismic sea waves in the East Sea, numerical models based on the Boussinesq equation are established and used to simulate recent tsunamis. By considering Ursell parameter and Kajiura parameter, it is proved that Boussinesq equation is a proper equation for seismic sea waves in the East Sea. Two models based on the Boussinesq equation and linear wave equation are executed with the same initial conditions and grid size ($1min{\times}1min$), and the results are compared in various respects. The Boussinesq equation model produced better results than the linear model in respect to wave propagation and concentration of wave energy. It is also certified that the Boussinesq equation model can be used for operational purpose if it is optimized. Another Boussinesq equation model whose grid size is $40sec{\times}30sec$ is set up to simulate the 1983 and 1993 tsunamis. As the result of simulation, new propagation charts of 2 seismic sea waves focused on the Korean east coast are proposed. Even though the 1983 and 1993 tsunamis started at different areas, the propagation paths near the Korean east coast are similar and they can be distinguished into 4 paths. Among these, total energy and propagating time of the waves passing over North Korea Plateau(NKP) and South Korea Plateau(SKP) determine wave height at the Korean east coast. In case of the 1993 tsunami, the wave passing over NKP has more energy than the wave over SKP. In case of the 1983 tsunami, the huge energy of the wave passing over SKP brought about great maximum wave heights at Mukho and Imwon. The Boussinesq equation model established in this study is more useful for simulation of seismic sea waves near the Korean east coast than it is the Japanese coast. To improve understanding of seismic sea waves in shallow water, a coastal area model based on the Boussinesq equation is also required.

• Journal of Coastal Disaster Prevention
• /
• v.1 no.1
• /
• pp.1-9
• /
• 2014

In order to examine the characteristics of tsunami run-up heights on impermeable/permeable slope, a numerical wave tank by upgrading LES-WASS-3D was used in this study. Then, the model were compared with existing hydraulic model test for its verification. The numerical results well reproduced experimental results of solitary wave deformation, propagation and run-up height under various conditions. Also, the numerical simulation with a slope boundary condition has been carried out to understand solitary wave run-up on impermeable/permeable slope. It is shown that the run-up heights on permeable slope is 52.64-63.2% smaller than those on the impermeable slope because of wave energy dissipation inside the porous media. In addition, it is revealed that the numerical results with slope boundary condition agreed well with experimental results in comparison with the results by using stair type boundary condition.

• Structural Engineering and Mechanics
• /
• v.42 no.6
• /
• pp.831-847
• /
• 2012

Impact from water-borne debris during tsunami and flood events pose a potential threat to structures. Debris impact forces specified by current codes and standards are based on rigid body dynamics, leading to forces that are dependent on total debris mass. However, shipping containers and other debris are unlikely to be rigid compared to the walls, columns and other structures that they impact. The application of a simple one-dimensional model to obtain impact force magnitude and duration, based on acoustic wave propagation in a flexible projectile, is explored. The focus herein is on in-air impact. Based on small-scale experiments, the applicability of the model to predict actual impact forces is investigated. The tests show that the force and duration are reasonably well represented by the simple model, but they also show how actual impact differs from the ideal model. A more detailed three-dimensional finite element model is also developed to understand more clearly the physical phenomena involved in the experimental tests. The tests and the FE results reveal important characteristics of actual impact, knowledge of which can be used to guide larger scale experiments and detailed modeling. The one-dimensional model is extended to consider water-driven debris as well. When fluid is used to propel the 1-D model, an estimate of the 'added mass' effect is possible. In this extended model the debris impact force depends on the wave propagation in the two media, and the conditions under which the fluid increases the impact force are discussed.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.6
• /
• pp.72-80
• /
• 2007

The tsunami that resulted from the Central East sea Earthquake, which registered 7.7 on the Richter scale, that occurred over the entire water region in Akita on May. 26, 1983 and the tsunami that was triggered by the Southwest off Hokkaido Earthquake (7.8 on the Richter scale) that occurred in Southwest off Hokkaido on July 12, 1993 are representative cases that led to considerable damage in life and property, not only in Japan but also in Korea. In this study, multi-grid method was used in order to reproduce sufficiently the shoaling effect that occurs as water depth becomes shallow in the shallow water region and moving boundary condition was introduced to consider the runup in the coastal region. For the tsunamis that exerted considerable effect on the East Sea coast of Korea that were caused by the Central East Sea Earthquake in 1983 and the Southwest off Hokkaido Earthquake in 1993, characteristics like water level rise and propagation in the East Sea coast will be examined using numerical simulations. At the same time, these values will be compared with observed values. In addition, maximum water level rise and change in the water level with respect to time that were caused by the tsunamis were examined at each location along the East sea coast. Usefulness of numerical analysis was verified by comparing with observed values.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.6 no.3
• /
• pp.73-86
• /
• 2002

Prognostic characteristics of hypothetical tsunamis in the East Sea are discussed based on numerical model simulations of linear long wave theory and wave ray for 28 source zones of tsunamigenic earthquake selected by the seismic gap theory. As a result, the propagation patterns of tsunamis due to hypothetical earthquake are presented and analyses also lead to selection of the geographical zones with low risk of tsunamis.

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