• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.289-301
• /
• 2009

In this work, wave run-up heights and resultant wave forces on a vertical revetment due to tsunami (solitary wave) are investigated numerically using a numerical wave tank model called CADMAS-SURF (CDIT, 2001. Research and Development of Numerical Wave Channel (CADMAS-SURF). CDIT library, No. 12, Japan.), which is based on a 2-D Navier-Stokes solver, coupled to a volume of fluid (VOF) method. The third order approximate solution (Fenton, 1972. A ninth-order solution for the solitary wave. J. of Fluid Mech., Vol. 53, No.2, pp.257-271) is used to generate solitary waves and implemented in original CADMAS-SURF code. Numerical results of the wave profiles and forces are in good agreements with available experimental data. Using the numerical results, the regression curves determined from the least-square analysis are proposed, which can be used to determine the maximum wave run-up height and force on a vertical revetment due to tsunami. In addition, the capability of CADMAS-SURF is demonstrated for tsunami wave forces acting on an onshore structure using various configuration computations including the variations of the crown heights of the vertical wall and the position of the onshore structure. Based on the numerical results such as water level, velocity field and wave force, the direct effects of tsunami on an onshore structure are discussed.

• Atmosphere
• /
• v.23 no.3
• /
• pp.293-305
• /
• 2013

Nonhydrostatic effects on convectively forced mesoscale flows in two dimensions are numerically investigated using a nondimensional model. An elevated heating that represents convective heating due to deep cumulus convection is specified in a uniform basic flow with constant stability, and numerical experiments are performed with different values of the nonlinearity factor and nonhydrostaticity factor. The simulation result in a linear system is first compared to the analytic solution. The simulated vertical velocity field is very similar to the analytic one, confirming the high accuracy of nondimensional model's solutions. When the nonhydrostaticity factor is small, alternating regions of upward and downward motion above the heating top appear. On the other hand, when the nonhydrostaticity factor is relatively large, alternating updraft and downdraft cells appear downwind of the main updraft region. These features according to the nonhydrostaticity factor appear in both linear and nonlinear flow systems. The location of the maximum vertical velocity in the main updraft region differs depending on the degrees of nonlinearity and nonhydrostaticity. Using the Taylor-Goldstein equation in a linear, steady-state, invscid system, it is analyzed that evanescent waves exist for a given nonhydrostaticity factor. The critical wavelength of an evanescent wave is given by ${\lambda}_c=2{\pi}{\beta}$, where ${\beta}$ is the nonhydrostaticity factor. Waves whose wavelengths are smaller than the critical wavelength become evanescent. The alternating updraft and downdraft cells are formed by the superposition of evanescent waves and horizontally propagating parts of propagating waves. Simulation results show that the horizontal length of the updraft and downdraft cells is the half of the critical wavelength (${\pi}{\beta}$) in a linear flow system and larger than ${\pi}{\beta}$ in a weakly nonlinear flow system.

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

In order to evaluate the seismic capacity of massive vertical type breakwaters which have intensively been deployed along the coast of South Korea over the last two decades, we carry out the preliminary numerical simulation against the PoHang, GyeongJu, Hachinohe 1, Hachinohe 2, Ofunato, and artificial seismic waves based on the measured time series of ground acceleration. Numerical result shows that significant sliding can be resulted in once non-negligible portion of seismic energy is shifted toward the longer period during its propagation process toward the ground surface in a form of shear wave. It is well known that during these propagation process, shear waves due to the seismic activity would be amplified, and non-negligible portion of seismic energy be shifted toward the longer period. Among these, the shift of seismic energy toward the longer period is induced by the viscosity and internal friction intrinsic in the soil. On the other hand, the amplification of shear waves can be attributed to the fact that the shear modulus is getting smaller toward the ground surface following the descending effective stress toward the ground surface. And the weakened intensity of soil as the number of attacking shear waves are accumulated can also contribute these phenomenon (Das, 1993). In this rationale, we constitute the numerical model using the model by Hardin and Drnevich (1972) for the weakened shear modulus as shear waves go on, and shear wave equation, in the numerical integration of which $Newmark-{\beta}$ method and Modified Newton-Raphson method are evoked to take nonlinear stress-strain relationship into account. It is shown that the numerical model proposed in this study could duplicate the well known features of seismic shear waves such as that a great deal of probability mass is shifted toward the larger amplitude and longer period when shear waves propagate toward the ground surface.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.3
• /
• pp.139-148
• /
• 2000

In this paper, the interaction of oblique incident waves with a bottom-mounted and fluid-filled flexible membrane structure is investigated in the frame of linear hydro-elastic theory. The static shape of a membrane structure containing the fluid of a specific density is initially unknown and must be calculated before the hydrodynamic analysis. To solve hydrodynamic problem, the fluid domain is divided into the inner and outer region. The inner solution based on discrete membrane dynamic model and simple-source distribution over the entire fluid boundaries is matched to the outer solution ba~ed on an eigenfunction expansion method. The numerical results were compared to a series of Ohyama's experimental results. The measured reflection and tran¬smission coefficients reasonably follow the trend of predicted values. Using the computer program developed, the performance of a bottom-mounted and fluid-filled flexible membrane strocture is tested with various system parameters (membrane shape, internal pressure, density ratio) and wave characteristics (wave frequencies, incident wave angle). It is found that a bottom-mounted and fluid-filled flexible membrane structure can be an effel;tive wave barrier if properly designed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.6
• /
• pp.409-422
• /
• 2019

In this study, we carried out the 3D numerical simulation to investigate the hydraulic characteristics of Synchronous Edge wave known as the driving mechanism of beach cusp using the Tool Box called the ihFoam that has its roots on the OpenFoam. As a wave driver, RANS (Reynolds Averaged Navier-Stokes equation) and mass conservation equation are used. In doing so, we materialized short-crested waves known as the prerequisite for the formation of Synchronous Edge waves by generating two obliquely colliding Cnoidal waves. Numerical results show that as can be expected, flow velocity along the cross section where waves are focused are simulated to be much faster than the one along the cross section where waves are diverged. It is also shown that along the cross section where waves are focused, up-rush is moving much faster than its associated back-wash, but a duration period of up-rush is shortened, which complies the typical characteristics of nonlinear waves. On the other hand, due to the water-merging effect triggered by the redirected flow toward wave-diverging area at the pinacle of run-up, along the cross section where waves are diverged, offshore-ward velocity is larger than shore-ward velocity at the vicinity of shore-line, while at the very middle of shoaling process, the asymmetry of flow velocity leaned toward the shore is noticeably weakened. Considering that these flow characteristics can be found without exception in Synchronous Edge waves, the numerical simulation can be regarded to be successfully implemented. In doing so, new insight about how the boundary layer streaming occur are also developed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.6
• /
• pp.332-349
• /
• 2016

OLAFOAM is the powerful CFD code and is an expanded version of $OpenFOAM^{(R)}$, for wave mechanics simulation. The $OpenFOAM^{(R)}$ does provide many solvers to correspond to each object of the numerical calculation in a variety of fields. OLAFOAM's governing equation bases on VARANS (Volume-Averaged Reynolds-Averaged Navier-Stokes) equation, and the finite volume method is applied to numerical techniques. The program is coded in C++ and run on the Linux operating system. First of all, in this study, OLAFOAM was validated for 1) wave transformation inside porous structure under bore and regular wave conditions, 2) wave transformation by submerged breakwater under regular wave condition, and 3) regular wave transformation and resultant vertical velocity distribution under current by comparison with existing laboratory measurements. Hereafter, this study, which is almost no examination carried out until now, analyzed closely variation characteristics of water surface level, wave height, frequency spectrum, breaking waves, averaged velocity and turbulent kinetic energy around porous submerged breakwater in the wave and current coexisting field for the case of permeable or impermeable rear beach. It was revealed that the wave height fluctuation according to current direction(following or opposing) was closely related to the turbulent kinetic energy, and others.

• Journal of the Korean GEO-environmental Society
• /
• v.13 no.6
• /
• pp.59-66
• /
• 2012

The purpose of this study is to estimate the behavior and the mechanism of debris flow on the slope, which has specially various gradient plane. The numerical simulation was performed by using the Finite Differential Element method (FDM) based on the equation for the mass conservation and momentum conservation. The mechanism of flow type for debris flow is divided into three flow types which are stony debris flow, immature debris flow, and turbulent water flow, respectively. First, flow discharge, water flow depth, sediment volume concentration was investigated by variable input of flow discharge at the straight slope angle and two step inclined plane. As the input of flow discharge was decrease, flow discharge and water flow depth was increased, after the first coming debris flow only reached at the downstream. As the input of flow discharge was increased, the curve of flow discharge and flow depth was highly fluctuated. As the results of RMS ratio, the flow discharge and flow depth was lower two step slope angle than the straight slope angle. Second, the behavior of debris flow was investigated by the four cases of gradient degree at the downstream of slope angle. The band width of flow discharge and flow depth for $14^{\circ}$ between $16^{\circ}$ was higher than other gradient degree, and fluctuation curve was continuously high after 10 seconds.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.3
• /
• pp.180-198
• /
• 2019

In recent years, countries like Europe and Japan have been involved in many researches on the Low-Crested Structure (LCS) which is the method to protect beach erosion and it is regarded as an alternative to the submerged breakwaters, and compiled its results and released the design manual. In the past, studies on LCS have focused on two-dimensional wave transmission and calculating required weight of armor units, and these were mainly examined and discussed based on experiments. In this study, three-dimensional numerical analysis is performed on permeable LCS. The open-source CFD code olaFlow based on the Navier-Stokes momentum equations is applied to the numerical analysis, which is a strongly nonlinear analysis method that enables breaking and turbulence analysis. As a result, the distribution characteristics of the LCS such as water level, water flow, and turbulent kinetic energy were examined and discussed, then they were carefully compared and examined in the case of submerged breakwaters. The study results indicate that there is a difference between the flow patterns of longshore current near the shoreline, the spatial distribution of longshore and on-offshore directions of mean turbulent kinetic energy in case of submerged breakwaters and LCS. It is predicted that the difference in these results leads to the difference in sand movement.

• Geophysics and Geophysical Exploration
• /
• v.24 no.3
• /
• pp.89-97
• /
• 2021

Fault detection in seismic data is well suited to the application of machine learning algorithms. Accordingly, various machine learning techniques are being developed. In recent studies, machine learning models, which utilize synthetic data, are the particular focus when training with deep learning. The use of synthetic training data has many advantages; Securing massive data for training becomes easy and generating exact fault labels is possible with the help of synthetic training data. To interpret real data with the model trained by synthetic data, the synthetic data used for training should be geologically realistic. In this study, we introduce a method to generate realistic synthetic seismic data. Initially, reflectivity models are generated to include realistic fault structures, and then, a one-way wave equation is applied to efficiently generate seismic stack sections. Next, a migration algorithm is used to remove diffraction artifacts and random noise is added to mimic actual field data. A convolutional neural network model based on the U-Net structure is used to verify the generated synthetic data set. From the results of the experiment, we confirm that realistic synthetic data effectively creates a deep learning model that can be applied to field data.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2010

Seismic reflection surveying is one of the most widely used and effective techniques for coal seam structure delineation and risk mitigation for underground longwall mining. However, the ability of the method can be compromised by the presence of volcanic cover. This problem arises within parts of the Bowen and Sydney Basins of Australia and seismic surveying can be unsuccessful. As a consequence, such areas are less attractive for coal mining. Techniques to improve the success of seismic surveying over basalt flows are needed. In this paper, we use elastic wave-equation-based forward modelling techniques to investigate the effects and characteristics of seismic wave propagation under different settings involving changes in basalt properties, its thickness, lateral extent, relative position to the shot position and various forms of inhomogeneity. The modelling results suggests that: 1) basalts with high impedance contrasts and multiple flows generate strong multiples and weak reflectors; 2) thin basalts have less effect than thick basalts; 3) partial basalt cover has less effect than full basalt cover; 4) low frequency seismic waves (especially at large offsets) have better penetration through the basalt than high frequency waves; and 5) the deeper the coal seams are below basalts of limited extent, the less influence the basalts will have on the wave propagation. In addition to providing insights into the issues that arise when seismic surveying under basalts, these observations suggest that careful management of seismic noise and the acquisition of long-offset seismic data with low-frequency geophones have the potential to improve the seismic results.