• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.43-48
• /
• 1998

The development of unsteady three-dimensional incompressible viscous solver based on unsteady physical curvilinear coordinate system is presented. A 12-point finite analytic scheme based on local uniform grid spacing is extended for nonuniform grid spacing. The formulation of a condition is suggested to avoid the oscillation of the series summations produced by the application of the method of separation of variables. SIMPLER and pressure Poisson equation techniques are used for solving a velocity-pressure coupled problem. The matrix is solved using the Generalized Minimal RESidual (GMRES) method to enhance the convergence rate of unsteady flow solver and the Kinematic boundary condition of a free surface flow. It is demonstrated that the numerical solutions of these equations are less mesh sensitive.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.3
• /
• pp.157-167
• /
• 2015

In order to provide simple and accurate wave theory in design of offshore structure, an analytical approximation is introduced in this paper. The solution is limited to flat bottom having a constant water depth. Water is considered as inviscid, incompressible and irrotational. The solution satisfies the continuity equation, bottom boundary condition and non-linear kinematic free surface boundary condition exactly. Error for dynamic condition is quite small. The solution is suitable in description of breaking waves. The solution is presented with closed form and dispersion relation is also presented with closed form. In the last century, there have been two main approaches to the nonlinear problems. One of these is perturbation method. Stokes wave and Cnoidal wave are based on the method. The other is numerical method. Dean's stream function theory is based on the method. In this paper, power series method was considered. The power series method can be applied to certain nonlinear differential equations (initial value problems). The series coefficients are specified by a nonlinear recurrence inherited from the differential equation. Because the non-linear wave problem is a boundary value problem, the power series method cannot be applied to the problem in general. But finite number of coefficients is necessary to describe the wave profile, truncated power series is enough. Therefore the power series method can be applied to the problem. In this case, the series coefficients are specified by a set of equations instead of recurrence. By using the set of equations, the nonlinear wave problem has been solved in this paper.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.2
• /
• pp.13-29
• /
• 1996

This paper presents calculated results of the free-surface flow around a Series 60($C_B=0.6$) model. Three-dimensional Navier-Sotkes equations are solved and Baldwin-Lomax algebraic turbulence model is adopted to simulate the high Reynolds-number flow. To reduce computational efforts, velocity components near the wall are extrapolated with a the solved by using the Implicit Approximate Factorization method[2]. The successive-over-relaxation method is used for solving pressure-Poisson equation when obtaining the pressure field projecting the divergence-free velocity field. To simulate the free-surface flows more precisely, the numerical scheme solving the equation for the kinematic boundary condition is very important. In this paper, there numerical schemes are employed and the results are compared with the available experimental data.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.2
• /
• pp.101-107
• /
• 2022

Dean (1965) proposed the use of the root mean square error (RMSE) in the dynamic free surface boundary condition (DFSBC) and kinematic free-surface boundary condition (KFSBC) as an error evaluation criterion for wave theories. There are well known wave theories with RMSE more than 1%, such as Airy theory, Stokes theory, Dean's stream function theory, Fenton's theory, and trochodial theory for deep-water waves. However, none of them can be applied for deep-water breaking waves. The purpose of this study is to provide a closed-form solution for deep-water waves with RMSE less than 1% even for breaking waves. This study is based on a previous study (Shin, 2016), and all flow fields were simplified for deep-water waves. For a closed-form solution, all Fourier series coefficients and all related parameters are presented with Newton's polynomials, which were determined by curve fitting data (Shin, 2016). For verification, a wave in Miche's limit was calculated, and, the profiles, velocities, and the accelerations were compared with those of 5^th-order Stokes theory. The results give greater velocities and acceleration than 5^th-order Stokes theory, and the wavelength depends on the wave height. The results satisfy the Laplace equation, bottom boundary condition (BBC), and KFSBC, while Stokes theory satisfies only the Laplace equation and BBC. RMSE in DFSBC less than 7.25×10^-2% was obtained. The series order of the proposed method is three, but the series order of 5^th-order Stokes theory is five. Nevertheless, this study provides less RMSE than 5^th-order Stokes theory. As a result, the method is suitable for offshore structural design.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.34 no.4
• /
• pp.93-99
• /
• 2022

Analytical solutions for regular waves generated by bottom wave makers in a 3-dimensional wave basin were derived in this study. Bottom wave makers which have triangular, rectangular and combination of two shapes were adopted. The 3-dimensional velocity potential was derived based on the linear wave theory with the bottom moving boundary condition, kinematic and dynamic free surface boundary conditions in a wave basin. Then, analytical solutions of 3-dimensional particle velocities and free surface displacement were derived from the velocity potential. The solutions showed physically valid results for regular waves generated by bottom wave makers in a wave basin. The analytical solution for obliquely propagating wave generation from bottom wave maker which works like a snake was also derived. Numerical results of the solution agree well with theoretically predicted results.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.6
• /
• pp.7-15
• /
• 2003

A Digital Wave Tank simulation technique, based on a finite-difference method and a modified marker-and-cell (MAC) algorithm, is applied in order to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach, Ohkushiri Island, and to predict maximum wove run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain, and the boundary values are updated at each time step, by a finite-difference time-marching scheme in the frame of the rectangular coordinate system. The fully nonlinear, kinematic, free-surface condition is satisfied by the modified marker-density function technique. The near shore Tsunami is assumed to be a solitary wave, and is generated from the numerical wave-maker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods, based on the shallow-water wave theory.

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

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.05a
• /
• pp.231-239
• /
• 2003

A Digital Wave Tank simulation technique based on a finite-difference method and a modified marker-and-cell (MAC) algorithm is applied to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach and Ohkushiri island, and to predict maximum wave run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain and the boundary values updated at each time step by a finite-difference time-marching scheme in the frame of rectangular coordinate system. The fully nonlinear kinematic free-surface condition is satisfied by the modified marker-density function technique. The Nearshore Tsunami is assumed to be a solitary wave and generated from the numerical wavemaker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods based on the shallow-water wave theory.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.2
• /
• pp.33-39
• /
• 1982

It is one of the basic problems of naval architecture and ocean engineering how to describe the wave kinematics normally under the assumption of an ideal fluid. At present, there are many wave theories available for design purposes. These can be classified into two groups: One is the analytic theory and the other is the numerical theory. This paper briefly introduces the stream function method of R.G. Dean which belongs to the latter group and shows its numerical evaluations exemplified for two cases: One is applied to observed waves and the other is for design waves. In the former case, the wave profiles are calculated by the stream function method and compared with those of the observed waves and also with the results of R.G. Dean. They show good agreement. In the latter case, the wave kinematics and wave loads on a column of diameter 1m are calculated by the stream function method and these are compared with those resulted from the 5th-order gravity wave theory. As a result of comparison the values by the stream function method are slightly larger than those by the 5th-order gravity wave theory but the difference are negligible. From this it is concluded that the stream function method is very useful. And as characteristics of the numerical theories, the stream function method of R.G. Dean can be easily extended to the higher order terms and can include easily the current velocity and the pressure distribution on the free surface. In addition, when the data of observed wave profile are given, this method can reproduced the observed wave profile as closely as possible so that this method seems to describe the ocean wave more realistically. And from standpoint of a mathematical principle the stream function method exactly satisfies the kinematic free-surface boundary condition.