• Journal of Korean Port Research
• /
• v.13 no.2
• /
• pp.399-408
• /
• 1999

Copeland’s(1985) hyperbolic mild-slope equation including diffraction refraction and reflection in the wave field is used as a governing equation in this study. The result of Maruyama & Kajima(1985) is used to calculate wave direction and that of Watanabe & Maruyama(1986) is used as a energy dissipation formula. Numerical solutions are obtained by the Leap-Frog scheme and compared with Watanabe & Maruyama’s (1984) hydraulic experimental results and numerical simulation results for the detached breakwater. This wave model is applied to a detached breakwater and compared with Watanabe and Maruyama’s (1984) hydraulic model results to check the characteristics of reflected wave field around a detached breakwater. The distribution of wave height and we phase in front of a detached breakwater is more accurate than the Watanabe and Maruyama’s numerical results. The results from our wave model show good agreements with the others and also show nonlinear effects around the detached breakwater. This model is applied to the Gamcheon harbor of pusan. the field observations were carried out at Pusan harbor wave station in 1986-1995 and the results were accepted as a design wave condition in this study. The wave height and wave period was measured by Dong-A university at one station in the Gamcheon harbor in 1996-1997 and used as a calibration criterion. The measured data were used as input data for the numerical simulation and also compared with simulated results. The numerical simulation shows a fairly good results which considering the effect of topographic characteristics and effect of narrow entrance due to two separated breakwaters in Gamcheon harbor. The wave distribution characteristics inside Gamcheon harbor is quite different with the offshore wave direction and wave period.

• Geomechanics and Engineering
• /
• v.21 no.3
• /
• pp.227-236
• /
• 2020

Non-destructive exploration using elastic waves has been widely used to characterize rock mass properties. Wave propagation in jointed rock masses is significantly governed by the characteristics and orientation of discontinuities. The relationship between spatial heterogeneity (i.e., joint spacing) and wavelength for elastic waves propagating through jointed rock masses have been investigated previously. Discontinuous rock masses can be considered as an equivalent continuum material when the wavelength of the propagating elastic wave exceeds the spatial heterogeneity. However, it is unclear how stress-dependent long-wavelength elastic waves propagate through a repetitive rock-joint system with multiple joints. A preliminary numerical simulation was performed in in this study to investigate long-wavelength elastic wave propagation in regularly jointed rock masses using the three-dimensional distinct element code program. First, experimental studies using the quasi-static resonant column (QSRC) testing device are performed on regularly jointed disc column specimens for three different materials (acetal, aluminum, and gneiss). The P- and S-wave velocities of the specimens are obtained under various normal stress levels. The normal and shear joint stiffness are calculated from the experimental results using an equivalent continuum model and used as input parameters for numerical analysis. The spatial and temporal sizes are carefully selected to guarantee a stable numerical simulation. Based on the calibrated jointed rock model, the numerical and experimental results are compared.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2010.10a
• /
• pp.27-27
• /
• 2010

There have been several studies about empirical relation between seismic source parameters(e.g., focal mechanisms, depths, magnitudes, etc.) and T-wave observation. In order to delineate the relation, numerical and theoretical approaches to figure out T-wave excitation mechanism are required. In an attempt to investigate source radiation and wave scattering effects in the oceanic crust on T-wave envelopes, we perform three-dimensional numerical modeling to synthesize T-wave envelopes. We first calculate seismic P- and SV-wave energy on the seafloor using the Direct Simulation Monte Carlo based on the Radiative Transfer Theory, which enables us to take into account both realistic seismic source parameters and wave scattering in heterogeneous media, and then estimate excited T-wave energy by normal mode computation. The numerical simulation has been carried out considering the following different conditions: source types (strike and normal faults), source depths (shallow and deep), and wave propagation through homogeneous and heterogeneous Earth media. From the results of numerical modeling, we confirmed that T-wave envelopes vary according to spatial seismic energy distributions on the seafloor for the various input parameters. Furthermore, the synthesized T-wave envelopes show directional patterns due to anisotropic source radiation, and the slope change of T-wave envelopes caused by focal depth. Seismic wave scattering in the oceanic crust is likely to control the shape of envelopes.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.39-44
• /
• 2008

Stress wave propagation plays an important role in many engineering problems for reducing industrial noise and vibrations. In this paper, the dispersion-corrected finite element model is proposed for reducing the dispersion error in simulation of stress wave propagation. At eliminating the numerical dispersion error arising from the numerical simulation of stress wave propagation, numerical dispersion characteristics of the wave equation based finite element model are analyzed and some dispersion control scheme are proposed. The validity of the dispersion correction techniques is demonstrated by comparing the numerical solutions obtained using the present techniques.

• Journal of Mechanical Science and Technology
• /
• v.19 no.12
• /
• pp.2263-2269
• /
• 2005

The ship wave phenomena in the restricted waterway were investigated by a numerical analysis. The Euler and continuity equations were employed for the present study. The boundary fitted and moving grid system was adopted to enhance the computational efficiency. The convective terms in the governing equations and the kinematic free surface boundary condition were solved by the Constrained Interpolated Profile (CIP) algorithm in order to solve accurately wave heights in far field as well as near field. The advantage of the CIP method was verified by the comparison of the computed results by the CIP and the Maker and Cell (MAC) method. The free surface flow simulation around Wigley hull was performed and compared with the experiment for the sake of the validation of the numerical method. The present numerical scheme was applied to the free surface simulation for various canal sections in order to understand the effect of the sectional shape of waterways on the ship waves. The wave heights on the side wall and the shape of the wave patterns with their characteristics of flow are discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.675-678
• /
• 2008

The effect of wave and current interactions on jet-like current flowing against waves was investigated based on numerical simulations. The numerical simulations are conducted by a combined model system of REF/DIF(a wave model) plus SHORECIRC(a current model) and a Boussinesq equation model, FUNWAVE. In the simulations, regular and irregular waves refracted due to the jet-like opposing current were focused along the core region of current, and the jet-like current was earlier spreaded when the waves had larger wave heights. The numerical results show that the rapid change of wave height distribution in transverse direction near current inlet plays a significant role to spread the jet-like current. In other words, the gradients of radiation stress forcing in transverse direction have a more significant effect on the jet-like current than its relatively small gradients forcing in flowing direction, which tend to accelerate the current, do. In conclusion, it is indispensible to take into account the interaction effect of wave transformation and current characteristics when waves meet the opposing jet-like current such as river mouth.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2000.04a
• /
• pp.140-145
• /
• 2000

The spectral method which is composed of an eigenfunction expansion of free modes in the wave number domain is used to produce two dimensional unsteady inviscid wave simulation such as progressive waves in a numerical pneumatic wave tank. A spatial and time dependent free surface elevation and the potential are calculated by integrating ODE derived from fully nonlinear kinematic and dynamic free surface boundary condition at each time step. The nonlinear characteristics in the waves by this method were notable as increasing wave steepness. This method is very useful and powerful in terms of saving computational time caused by rapid convergence exponentially with increasing number of nodes, even preserving accurate numerical results. Moreover, it will given us many possibilities to apply to naval and ocean engineering fields.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.504-512
• /
• 2003

Simulation programs have been developed and used as an attempt to improve the accuracy of Non-Destructive Evaluation(NDE) techniques. The applicability of these programs is very limited, however, because it is difficult to describe the delicacy of the propagation of stress waves. To investigate the applicability of the finite element analysis for stress wave-based NDE techniques numerical simulation for Impact-Echo method and SASW method is performed. The numerical studies are performed to determine the essential parameters such as contact time of impact load, mesh size and time step size. These studies show that the choice of parameter is very important for improving the accuracy and confidence of the numerical procedure and, thereby, the applicability of the numerical analysis for stress wave-based NDE techniques

• Journal of Ocean Engineering and Technology
• /
• v.26 no.6
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.