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

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.10 s.253
• /
• pp.1235-1241
• /
• 2006

A new hydrocode which is still under development using Lagrangian, Eulerian and arbitrary Lagrangian-Eulerian operators, has been described. The three operators are implemented into a single framework by incorporating the sequential three stages of Lagrangian, remesh and remap stages. Several numerical schemes used for each operator are discussed briefly in this paper. In order to evaluate the characteristics of each operator, the Taylor Impact Test has been simulated using each operator and the results are compared. Currently the code is 1st order accuracy in the material interface tracking algorithm and can not handle multimaterial in the mixed cell. The areas of possible enhancement of the code are also discussed.

• 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 the Korean Society of Propulsion Engineers
• /
• v.14 no.2
• /
• pp.1-9
• /
• 2010

LES(Large eddy simulation) of breakup and droplet atomization of a liquid jet into cross turbulent flow was performed. Two phase flow of gas and liquid phases were modeled by the mixed numerical scheme of both Eulerian and Lagrangian methods for gas and liquid droplet respectively. The breakup process of a liquid column and droplets was observed by implementing the blob-KH wave breakup model. The penetration depth into cross flow was comparable with experimental data for several variants of the liquid-gas momentum flux ratio by varying liquid injection velocity. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.

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

• 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의 변화에 따른 비선형 운동해석을 연구하였다. 그 결과 비선형성이 심하지 않은 기진파의 경우에서는 비선형 운동해석의 수치해석 안정성의 선형 수치해석 안정성과 큰 차이가 없음을 알 수 있게 된다.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.388-397
• /
• 2010

A two-dimensional time-domain, potential-theory-based fully nonlinear numerical wave tank (NWT) was developed by using boundary element method and the mixed Eulerian-Lagrangian (MEL) approach for free-surface node treatment. The NWT was applied to prediction of primary wave energy conversion efficiency of a bottom-mounted oscillating water column (OWC) wave power device. The nonlinear free-surface condition inside the chamber was specially devised to represent the pneumatic pressure due to airflow velocity and viscous energy loss at the chamber entrance due to wave column motion. The newly developed NWT technique was verified through comparison with given experimental results. The maximum energy extraction was estimated with various chamber-air duct volume ratios.

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