• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.6 no.1
• /
• pp.1-15
• /
• 2002

Advection-dominated transport problems possess difficulties in the design of numerical methods for solving them. Because of the hyperbolic nature of advective transport, many characteristic numerical methods have been developed such as the classical characteristic method, the Eulerian-Lagrangian method, the transport diffusion method, the modified method of characteristics, the operator splitting method, the Eulerian-Lagrangian localized adjoint method, the characteristic mixed method, and the Eulerian-Lagrangian mixed discontinuous method. In this paper relationships among these characteristic methods are examined. In particular, we show that these sometimes diverse methods can be given a unified formulation. This paper focuses on characteristic finite element methods. Similar examination can be presented for characteristic finite difference methods.

• Structural Engineering and Mechanics
• /
• v.47 no.5
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.1
• /
• pp.33-55
• /
• 1993

선수충격파의 문제를 푸는데 있어서 Boundary Integral Method(BIM)의 여러가지 수치 해석방법이 검토되었으며, 특히 여러가지 Time stepping scheme, Green function, far-field 조건등에 따른 수치해석안정성과 정확성의 상관관계가 연구되었다. von Neumann 안정성해석과 matrix 안정성해석 등을 이용한 선형 안정성해석을 기초로하여, 수치해석방법의 안정성 여부를 체계적으로 조사할 수 있는 parameter(Free Surface Stability number)를 설정하고, 이 parameter의 변화에 따른 비선형 운동해석을 연구하였다. 그 결과 비선형성이 심하지 않은 기진파의 경우에서는 비선형 운동해석의 수치해석 안정성의 선형 수치해석 안정성과 큰 차이가 없음을 알 수 있게 된다.

• Ocean Systems Engineering
• /
• v.2 no.2
• /
• pp.147-158
• /
• 2012

A time-domain simulation of a land-based Oscillating Water Column (OWC) with various irregular waves as a form of PM spectrum is performed by using a two-dimensional fully nonlinear numerical wave tank (NWT) based on the potential theory, mixed Eulerian-Lagrangian (MEL) approach, and boundary element method. The nonlinear free-surface condition inside the OWC chamber was specially devised to describe both the pneumatic effect of the time-varying pressure and the viscous energy loss due to water column motions. The quadratic models for pneumatic pressure and viscous loss are applied to the air and free surface inside the chamber, and their numerical results are compared with those with equivalent linear ones. Various wave spectra are applied to the OWC system to predict the efficiency of wave-energy take-off for various wave conditions. The cases of regular and irregular waves are also compared.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.5
• /
• pp.1-8
• /
• 2007

A long wave induced by a Gaussian-shape submarine landslide is simulated by a 2D fully nonlinear numerical wave tank (NWT). The NWT is based on the boundary element method and the mixed Eulerian/Lagrangian approach. Using the NWT, physical characteristics of land-slide tsunami, including wave generation, propagation, particle kinematics, hydrodynamic pressure, run-up and depression, are simulated for the early stage of long wave generation and propagation. Various sliding mass heights are applied to the developed model for a systematic sensitivity analysis. In particular, the fully nonlinear NWT results are compared with linear results (exact body-boundary conditions with linear free-surface conditions) to identify the nonlinear effects in the respective cases.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.4
• /
• pp.513-528
• /
• 2013

A floating Oscillating Water Column (OWC) wave energy converter, a Backward Bent Duct Buoy (BBDB), was simulated using a state-of-the-art, two-dimensional, fully-nonlinear Numerical Wave Tank (NWT) technique. The hydrodynamic performance of the floating OWC device was evaluated in the time domain. The acceleration potential method, with a full-updated kernel matrix calculation associated with a mode decomposition scheme, was implemented to obtain accurate estimates of the hydrodynamic force and displacement of a freely floating BBDB. The developed NWT was based on the potential theory and the boundary element method with constant panels on the boundaries. The mixed Eulerian-Lagrangian (MEL) approach was employed to capture the nonlinear free surfaces inside the chamber that interacted with a pneumatic pressure, induced by the time-varying airflow velocity at the air duct. A special viscous damping was applied to the chamber free surface to represent the viscous energy loss due to the BBDB's shape and motions. The viscous damping coefficient was properly selected using a comparison of the experimental data. The calculated surface elevation, inside and outside the chamber, with a tuned viscous damping correlated reasonably well with the experimental data for various incident wave conditions. The conservation of the total wave energy in the computational domain was confirmed over the entire range of wave frequencies.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.6
• /
• pp.12-16
• /
• 2009

In this study, the wave profiles and kinematics of highly nonlinear waves at various water depths were calculated using a 2D fully nonlinear Numerical Wave Tank (NWT). The NWT was developed based on the Boundary Element Method (BEM) with the potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme by 4th-order Runge-Kutta time integration. The spatial variation of intermediate-depth waves along the direction of wave propagation was caused by the unintended generation of 2nd-order free waves, which were originally investigated both theoretically and experimentally by Goda (1998). These free waves were induced by the mismatch between the linear motion of wave maker and nonlinear displacement of water particles adjacent to the maker. When the 2nd-order wave maker motion was applied, the spatial modulation of the waves caused by the free waves was not observed. The respective magnitudes of the nonlinear wave components for various water depths were compared. It was found that the high-order wave components greatly increase as the water depth decreases. The wave kinematics at various locations were calculated and compared with the linear and the Stokes 2nd-order theories.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.355-372
• /
• 2011

Based on the fully nonlinear velocity potential theory, the liquid sloshing in a three dimensional tank under random excitation is studied. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing scheme, B-spline curve, is applied to both the longitudinal and transverse directions of the tank to eliminate the possible saw-tooth instabilities. When the tank is undergoing one dimensional regular motion of small amplitude, the calculated results are found to be in very good agreement with linear analytical solution. In the simulation, the normal standing waves, travelling waves and bores are observed. The extensive calculation has been made for the tank undergoing specified random oscillation. The nonlinear effect of random sloshing wave is studied and the effect of peak frequency used for the generation of random oscillation is investigated. It is found that, even as the peak value of spectrum for oscillation becomes smaller, the maximum wave elevation on the side wall becomes bigger when the peak frequency is closer to the natural frequency.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.6
• /
• pp.34-40
• /
• 2010

Bragg reflection of nonlinear waves is simulated by a 2D fully nonlinear numerical wave tank (NWT). The developed NWT was based on the Boundary Element Method (BEM) with potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme with Runge-Kutta 4th-order time integration. A spatial variation of wave elevations and their Fourier amplitudes of each component are compared to investigate the effect of sea bottom ripples and their relative heights. The incident waves over an undulated sea bottom are partially reflected and changed to partial standing waves due to Bragg reflection. The present results are verified with linear calculations and experimental data. It is found that the 1st-order wave component is mainly affected by Bragg reflection and its spatial modulation is significant in front of the bottom ripples.