• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.65-72
• /
• 1993

A numerical method for nonlinear free-surface-wave problem is developed in this paper. The final goal of this study is to simulate the towing tank experiment of a ship model and to partially replace the experiment by the numerical model. The exact problem in the scope of potential flow theory is formulated by a variational principle based on the classical Hamilton's principle. A localized finite element method is used in the present numerical computations which made use of the following two notable steps. The first step is an efficient treatment of the numerical radiation condition by using the intermediate nonlinear-to-linear transition buffer subdomain between the fully nonlinear and linear subdomains. The second is the use of a modal analysis in the final stage of the solution procedures, which enables us to reduce the computation time drastically. With these improvements the present method can treat a much larger computational domain than that was possible previously. A pressure patch on the free surface was chosen as an example. From the present computed results we could investigate the effect of nonlinearity on the down-stream wave pattern more clearly than others, because much larger computational domain was treated. We found, specifically, the widening of the Kelvin angle and the increase of the wave numbers and the magnitude of wave profiles.

• Journal of applied mathematics & informatics
• /
• v.13 no.1_2
• /
• pp.53-65
• /
• 2003

Here we present a study of small-amplitude, shallow water waves on sloping beds. The beds considered in this analysis are linear and quadratic in nature. First we start with stating the relevant governing equations and boundary conditions for the theory of water waves. Once the complete prescription of the water-wave problem is available based on some assumptions (like inviscid, irrotational flow), we normalize it by introducing a suitable set of non-dimensional variables and then we scale the variables with respect to the amplitude parameter. This helps us to characterize the various types of approximation. In the process, a summary of equations that represent different approximations of the water-wave problem is stated. All the relevant equations are presented in rectangular Cartesian coordinates. Then we derive the equations and boundary conditions for small-amplitude and shallow water waves. Two specific types of bed are considered for our calculations. One is a bed with constant slope and the other bed has a quadratic form of surface. These are solved by using separation of variables method.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.237-246
• /
• 1999

A finite element model is studied to simulate unsteady free surface flow based on dynamic wave equation and collocation method. The collocation method is used in conjunction with Hermite polynomials, and resulting matrix equations are solved by skyline method. The model is verified by applying to hydraulic jump, nonlinear disturbance propagation and dam-break flow in a horizontal frictionless channel. The computed results are compared with those by Bubnov-Galerkin and Petrov-Galerkin methods. It is also applied to the North Han River to simulate the floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge hydrographs. The suggested model has proven to be one of the promising scheme for simulating the gradually and rapidly varied unsteady flow in open channels.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.11
• /
• pp.1169-1176
• /
• 2011

Liquid film flows are classified into waveless laminar, wavy laminar, and turbulent flows depending on the Reynolds number or the flow stability. Since the wavy motions of the film flows are so intricate and nonlinear, studies on them have largely been experimental. Most numerical approaches have been limited to the waveless flow regime. The various free surface-tracking schemes adopted for this problem were used to more accurately estimate the average film thickness, rather than to capture the unsteady wavy motion. In this study, the wavy motions in laminar wavy liquid film flows with Reynolds numbers of 200-1000 were simulated with various numerical schemes based on the volume of fluid (VOF) method for interface tracking. The results from each numerical scheme were compared with the experimental results in terms of the average film thickness, the wave velocity, and the wave amplitude.

• Journal of Ocean Engineering and Technology
• /
• v.19 no.3
• /
• pp.18-24
• /
• 2005

Wave overtopping is one of the most important processes for the design of seawalls. The term "wave overtopping" is used here to refer to the processes where waves hit a sloping structure run up the slope and, if the crest level of the slope is lower than the highest run up level, overtop the structure. Wave overtopping is dependent on the processes associated with breaking wave. A numerical model based on Navier-Stokes equation and the Marker-density function for predicting wave overtopping of coastal structures is developed in this paper. In order to evaluate the present model, two simulations are tested. One is overflow without waves at vertical seawall, and the other is wave overtopping at sloping seawalls.

• Special Issue of the Society of Naval Architects of Korea
• /
• 2005.06a
• /
• pp.22-30
• /
• 2005

Sloshing loads, produced by the violent liquid free-surface motions inside the cargo tank have become an important design parameter in ship building industry since there have been demands for the increased sizes of the cargo containment system of LNG carriers. In this study, sloshing impact pressure acting on the shell of the spherical cargo tank of an LNG carrier as well as dynamic pressure and flow behavior around the pump tower located at the center of the tank have been calculated. Comparative numerical sloshing simulations for a spherical LNG tank using 2-D LR.FLUIDS which is based on the finite difference method and 3-D MSC.DYTRAN which is capable of calculating nonlinear fluid-structure interaction have been carried out. A method of calculating sloshing-induced dynamic loads and the subsequent structural strength analysis for pump tower of a spherical LNG carrier using MSC. DYTRAN and MSC.NASTRAN have been presented.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.6
• /
• pp.369-381
• /
• 2017

It has been presumed that highly nonlinear skewed waves frequently observed in a surf zone could significantly influence the transmission behaviour via a porous wave breaker due to its larger inertia force than its nonlinear counterparts of zero skewness [Cnoidal waves]. In this study, in order to confirm this perception, a numerical simulation has been implemented for 6 waves the skewness of that range from 1.02 to 1.032. A numerical simulation are based on the Tool Box called as the ihFoam that has its roots on the OpenFoam. Skewed waves are guided by the shoal of 1:30 slope, and the flow in the porous media are analyzed by adding the additional damping term into the RANS (Reynolds Averaged Navier-Stokes equation). Numerical results show that the highly nonlinear skewed waves are of higher transmitted ratio than its counterparts due to its stronger inertia force. In this study, in order to see whether or not the damping at the porous structure has an effect on the wave celerity, we also derived the dispersive relationships of Nonlinear Shallow Water Eq. [NSW] with damping at the porous structure being accounted. The newly derived dispersive relationships shows that the phase lag between the damping friction and the free surface elevation due to waves significantly influence the wave celerity.

• Journal of Korea Water Resources Association
• /
• v.36 no.6
• /
• pp.937-947
• /
• 2003

A shoreline, which has no the water depth, moves continuously as waves rise up and recede. Therefore, a special boundary treatment is required to track properly the movements of the shoreline in numerical modeling of the behavior of tsunamis or tides near a coastal zone. In this study, convective terms in nonlinear shallow-water equations are discretized explicitly by using a second-order upwind scheme to describe a moving shoreline more accurately. An oscillatory flow motion in a circular paraboloidal basin has been employed to validate the performance of the developed numerical model. Computed results of instantaneous free surface displacements are compared with those of analytical solutions and existing numerical solutions. The run-up heights in the vicinity of a circular island have also been calculated and obtained numerical results have been shown against available laboratory measurements. A good agreement has been observed.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.399-413
• /
• 2020

The numerical simulation of liquid sloshing in the three-dimensional tanks under horizontal excitation and roll excitation was carried out, and the inhibition effect of different baffles on the sloshing phenomenon was investigated. The numerical calculations were carried out by the nonlinear Boundary Element Method (BEM) with Green's theorem based on the potential flow, which was conducted with the governing equation corresponding to the boundaries of each region. The validity of the method was verified by comparing with experimental values and published literatures. The horizontal baffle, the vertical baffle and the T-shaped baffle in the sloshing tanks were investigated respectively, and the baffles' position, dimension and the liquid depth were provided and discussed in detail. It is drawn that the baffle shape plays a non-negligible role in the tank sloshing. The vertical baffle is a more effective way to reduce the sloshing amplitude when the tank is under a horizontal harmonic excitation while the horizontal baffle is a more effective way when the tank is under a roll excitation. The amplitude of free surface elevation at right tank wall decreases with the increasing of the horizontal baffle length and the vertical baffle height. Although the T-shaped baffle has the best suppression effect on tank sloshing under horizontal excitation, it has limited suppression effect under roll excitation and will complicate the sloshing phenomenon when changing baffle height.