• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.3
• /
• pp.201-207
• /
• 2011

In this paper, we focused on computing the higher-harmonic components of the transmitted wave passing over a submerged circular cylinder to show that it is causing a horizontal negative drift force. As numerical models, a circular cylinder held fixed under free surface in deep water is adopted. As the submergence of a circular cylinder decreases and the incident wavelength becomes longer, the higher-harmonic components of the transmitted wave starts to increase. An increase of the higher-harmonic components of the transmitted wave makes the horizontal drift force be negative. It is also found that the higher-harmonic amplitudes averaged over the transmitted wave region become larger with the increase of wave steepness and wavelength as well as the decrease of submergence depth.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• Journal of the Korean Mathematical Society
• /
• v.49 no.3
• /
• pp.585-604
• /
• 2012

In this paper we consider the evolution of the rolling stone with a rotationally symmetric nonconvex compact initial surface ${\Sigma}_0$ under the Gauss curvature flow. Let $X:S^n{\times}[0,\;{\infty}){\rightarrow}\mathbb{R}^{n+1}$ be the embeddings of the sphere in $\mathbb{R}^{n+1}$ such that $\Sigma(t)=X(S^n,t)$ is the surface at time t and ${\Sigma}(0)={\Sigma}_0$. As a consequence the parabolic equation describing the motion of the hypersurface becomes degenerate on the interface separating the nonconvex part from the strictly convex side, since one of the curvature will be zero on the interface. By expressing the strictly convex part of the surface near the interface as a graph of a function $z=f(r,t)$ and the non-convex part of the surface near the interface as a graph of a function $z={\varphi}(r)$, we show that if at time $t=0$, $g=\frac{1}{n}f^{n-1}_{r}$ vanishes linearly at the interface, the $g(r,t)$ will become smooth up to the interface for long time before focusing.

• Journal of Ocean Engineering and Technology
• /
• v.35 no.4
• /
• pp.296-304
• /
• 2021

Research on the development of marine renewable energy is actively in progress. Various studies are being conducted on the development of wave energy converters. In this study, a numerical analysis of wave-energy extraction performance was performed according to the body shape and scale of the sloped oscillating water column (OWC) wave energy converter (WEC), which can be connected with the breakwater. The sloped OWC WEC was modeled in the time domain using a two-dimensional fully nonlinear numerical wave tank. The nonlinear free surface condition in the chamber was derived to represent the pneumatic pressure owing to the wave column motion and viscous energy loss at the chamber entrance. The free surface elevations in the sloped chamber were calculated at various incident wave periods. For verification, the results were compared with the 1:20 scaled model test. The maximum wave energy extraction was estimated with a pneumatic damping coefficient. To calculate the energy extraction of the actual size WEC, OWC models approximately 20 times larger than the scale model were calculated, and the viscous damping coefficient according to each size was predicted and applied. It was verified that the energy, owing to the airflow in the chamber, increased as the incident wave period increased, and the maximum efficiency of energy extraction was approximately 40% of the incident wave energy. Under the given incident wave conditions, the maximum extractable wave power at a chamber length of 5 m and a skirt draft of 2 m was approximately 4.59 kW/m.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.135-140
• /
• 2006

Numerical computations are made to predict wave loads on a vertical cylinder in an extreme wave. To generate the extreme wave, a frequency-focused unidirectional wave is adopted in three-dimensional numerical wave tank. The mathematical formulation is wide in the scope of the potential theory with fully nonlinear free surface conditions. As a numerical method, finite element method based on variational principle is applied. Comparisons between the present numerical results and the previous computation data. show a good agreement.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.174-180
• /
• 2002

A virtual reality technology for multipurpose numerical simulation is developed to reproduce and investigate a variety of ocean environmental problems in a 3D-Numerical Wave Tank. The governing equations for solving incompressible fluid motion are Navier-Stokes equation and continuity equation, and the Marker-Density function technique is adopted to implement the fully-nonlinear free-surface kinematic condition. The marine environmental situations, i.e. waves, currents, wind, etc., are reproduced by use of multi-segmented wavemaker on the basis of the so-called "snake-principle". In this paper, some numerical reproduction techniques for regular and irregular waves, multi-directional waves, Bull's-eye wave, wave-current, and solitary wave are presented, and a model test in motion with large amplitude of roll angle is conducted in the developed 3D-NWT, using a overlaid grid system.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.32-39
• /
• 1998

A numerical analysis for large amplitude motions of submerged circular cylinder is presented. The method is based on potential theory and two-dimensional motions in regular harmonic waves are tented as an initial value problem. The fully nonlinear free surface boundary condition is assumed in an inner domain and this solution is matched along an assumed an assumed common boundary to a linear solution in outer domain. Calculations of the large amplitude motion of a submerged circular cylinder are directly simulated in time domain. It is shown that relative motion between the body and fluid particle gives a significant effect on the lift and drift motions.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.484-491
• /
• 2007

Numerical experiments using a numerical wave tank have been performed to verier the nonlinear interaction between monochromatic waves and a submerged horizontal plate. As a model for numerical wave tank, we used a higher-order Boundary Element Method(BEM) based on fully nonlinear potential flow theory and CADMAS-SURF for solving Navier Stokes equations and exact free surface conditions. Both nonlinear models are able to predict the higher harmonic generation in the shallow water region over a submerged horizontal plate. CADMAS-SURF, which involves the viscous effect, can evaluate the higher harmonic generation by flow separation and vortices at the each ends of plate. The comparison of reflection and transmission coefficients with experimental results(Patarapanich and Cheong, 1989) at different lengths and submergence depths of a horizontal plate are presented with a good agreement. It is found that the transfer of energy from the incident fundamental waves to higher harmonics becomes larger as the submergence depth ratio decreases and the length ratio increases.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1445-1451
• /
• 1994

To calculate the wave pressure acting on coastal structures under the design wave condition, it is often necessary to numerically reproduce the big standing wave profiles close to wave breaking condition. For this, the governing equation and all nonlinear terms occurring in boundary conditions should be effectively considered in the numerical wave profile. In particular, the velocity square term in the free surface boundary condition is very important. A boundary element method is applied here to calculate the standing wave profile with the velocity square term fully treated by Newton iterative method. In order to check the validity of the method, the numerical wave profiles are compared to ones calculated by the perturbation method, the Fourier approximation method and the hydraulic experiment.