Simulation of Fully Nonlinear 3-D Numerical Wave Tank

A fully nonlinear numerical wave tank (NWT) has been simulated by use of a three-dimensional higher order boundary element method (HOBEM) in the time domain. Within the frame of potential flow and the adoption of simply Rankine source, the resulting boundary integral equation is repeatedly solved at each time step and the fully nonlinear free surface boundary conditions are integrated with time to update its position and boundary values. A smooth technique is also adopted in order to eliminate the possible saw-tooth numerical instabilities. The incident wave at the uptank is given as theoretical wave in this paper. The outgoing waves are absorbed inside a damping zone by spatially varying artificial damping on the free surface at the wave tank end. The numerical results show that the NWT developed by these approaches has a high accuracy and good numerical stability.

Development of A Fully Nonlinear Numerical Wave Tank

A fully nonlinear numerical wave tank based on the solution of the σ-transformed Navier-Stokes equation is developed in this study. The numerical wave is generated from the inflow boundary, where the surface elevation and/or velocity are specified by use of the analytical solution or the laboratory data. The Sommerfeld/Orlanski radiation condition in conjunction with an artificial damping zone is applied to reduce wave reflection from the outflow boundary. The whole numerical solution procedures are split into three steps, i.e., advection, diffusion and propagation, and a new method, the Lagrange-Euler Method, instead of the MAC or VOF method, is introduced to solve the free surface elevation at the new time step. Several typical wave cases, including solitary waves, regular waves and irregular waves, are simulated in the wave tank. The robustness and accuracy of the NWT are verified by the good agreement between the numerical results and the linear or nonlinear analytical solutions. This research will be further developed by study of wave-wave, wave-current, wave-structure or wave-jet interaction in the future.

Generation and Properties of Freak Waves in A Numerical Wave Tank

Freak waves are generated based on the mechanism of wave focusing in a 2D numerical wave tank. To set up the nonlinear numerical wave tank, the Boundary Element Method is used to solve potential flow equations incorporated with fully nonlinear free surface boundary conditions. The nonlinear properties of freak waves, such as high frequency components and wave profile asymmetry, are discussed. The kinematic data, which can be useful for the evaluation of the wave forces exerted on structures to avoid underestimation of linear predictions, are obtained, and discussed, from the simulated results of freak waves.

Determining the Onset and Strength of Unforced Wave Breaking in A Numerical Wave Tank

A numerical wave tank is used to investigate the onset and strength of unforced wave breaking, and the waves have three types of initial spectra： constant amplitude spectrum, constant steepness spectrum and Pierson-Moscowitz spectrum. Numerical tests are performed to validate the model results. Then, the onset of wave breaking is discussed with geometric, kinematic, and dynamic breaking criteria. The strength of wave breaking, which is always characterized by the fractional energy loss and breaking strength coefficient, is studied for different spectra. The results show how the energy growth rate is better than the initial wave steepness on estimating the fractional energy losses as well as breaking strength coefficient.

Generation of linear and nonlinear waves in numerical wave tank using clustering technique-volume of fluid method

A two-dimensional （2D） numerical model is developed for the wave sim- ulation and propagation in a wave flume. The fluid flow is assumed to be viscous and incompressible, and the Navier-Stokes and continuity equations are used as the governing equations. The standard k-e model is used to model the turbulent flow. The Navier- Stokes equations are discretized using the staggered grid finite difference method and solved by the simplified marker and cell （SMAC） method. Waves are generated and propagated using a piston type wave maker. An open boundary condition is used at the end of the numerical flume. Some standard tests, such as the lid-driven cavity, the constant unidirectional velocity field, the shearing flow, and the dam-break on the dry bed, are performed to valid the model. To demonstrate the capability and accuracy of the present method, the results of generated waves are compared with available wave theories. Finally, the clustering technique （CT） is used for the mesh generation, and the best condition is suggested.

A Numerical Wave Tank for Nonlinear Waves with Passive Absorption

A numerical wave tank with passive absorption for irregular waves is considered in this paper. Waves with spectral shapes corresponding to that of the Mitsuyasu- Bretschneider type are used as the initial condition at one end of the flume, An absorbing boundary is imposed at the other end of the wave flume to minimize reflection. By use of a Lagrangian description for the Surface elevation, and finite difference for approximation of the time derivative, the problem is then solved by the boundary element method, The effects of the absorbing boundary are investigated by varying the values of the absorption coefficient mu, and studying the time histories of the Surface elevations 'recorded' on pre-selected locations.

Two Dimensional Fully Nonlinear Numerical Wave Tank Based on the BEM

The development of a two dimensional numerical wave tank （NWT） with a rocker or piston type wavemaker based on the high order boundary element method （BEM） and mixed Eulerian-Lagrangian （MEL） is examined. The cauchy principle value （CPV） integral is calculated by a special Gauss type quadrature and a change of variable. In addition the explicit truncated Taylor expansion formula is employed in the time-stepping process. A modified double nodes method is assumed to tackle the comer problem, as well as the damping zone technique is used to absorb the propagation of the free surface wave at the end of the tank. A variety of waves are generated by the NWT, for example; a monochromatic wave, solitary wave and irregular wave. The results confirm the NWT model is efficient and stable.

Numerical Wave Tank Based on A Conserved Wave-Absorbing Method

Recently the numerical wave tank has become a widely-used tool to study waves as well as wave-structure interactions, and the wave-absorbing method is very important as its effect on the quality of waves generated. The relaxation method and the derived momentum source method are often utilized, however, the damping weight is constant during calculation and repeated trials are required to obtain an acceptable wave-absorbing effect. To address the abovementioned issues, a conserved wave-absorbing method is developed. The damping weight is determined by solving the mass conservation equation of the absorbing region at every time step. Based on this method, a two-dimensional numerical wave tank is established by using the VB language to simulate various waves by which the validation of this method is evaluated.

Development of two-dimensional numerical wave tank based on lattice Boltzmann method

A 2-D numerical wave tank(NWT)is developed using the lattice Boltzmann method(LBM)and a multi-relaxation-time(MRT)collision model coupled with an algebraic volume of fluid(VOF)scheme for free surface tracking.An external force based on the momentum source function is used to generate the waves,and a zone of porous media is used to absorb the waves.Numerical simulations of the progressive and standing waves show that the NWT can generate stable wave trains in agreement with the analytical solutions and eliminate the re-reflection waves.The NWT is used to simulate two problems encountered in practice,namely:the wave transformation over a submerged breakwater and the wave runup on a sea dike.The numerical predictions are in good agreement with the measured data.

Numerical Investigation on the Generation and Propagation of Irregular Waves in A Two-Dimensional Wave Tank

The modeling of generation and subsequent propagation of irregular waves in a numerical wave flume is performed by mean of the boundary element method. Random waves are generated by a piston-type wave generator at one end of the flume with the Mitsuyasu-Bretschneider spectrum used as the target spectrum for the generation. An artificial absorbing beach is placed at the other end of the flume to minimize wave reflection. Surface fluctuations are described by use of the Lagrangian description, and finite difference is adopted for the approximation of time derivative. To monitor the developments of the waves, a number of pseudo wave gauges are installed along the tank. Through comparison of the spectra from those gauges with the target spectrum, satisfactory results can be obtained from the present numerical scheme.

A third-order KdV solution for internal solitary waves and its application in the numerical wave tank

A third-order KdV solution to the internal solitary wave is derived by a new method based on the weakly nonlinear assumptions in a rigid-lid two-layer system.The solution corrects an error by Mirie and Su(1984).A two-dimensional numerical wave tank has been established with the help of the open source CFD library OpenFOAM and the third-party software waves2Foam.Various analytical solutions,including the first-order to third-order KdV solutions,the eKdV solution and the MCC solution,have been used to initialise the flow fields in the CFD simulations of internal solitary waves.Two groups including 11 numerical cases have been carried out.In the same group,the initial wave amplitudes are the same but the implemented analytical solutions are different.The simulated wave profiles at different moments have been presented.The relative errors in terms of the wave amplitude between the last time step and the initial input have been analysed quantitatively.It is found that the third-order KdV solution results in the most stable internal solitary wave in the numerical wave tank for both small-amplitude and finite-amplitude cases.The finding is significant for the further simulations involving internal solitary waves.

Research on Numerical Wave Tank Based on the Improved Moving-Particle Semi-Implicit Method with Large Eddy Simulation

Moving-particle semi-implicit(MPS) method is a new mesh-free numerical method based on Lagrangian particle. In this paper, MPS method is applied to the study on numerical wave tank. For the purpose of simulating numerical wave, we combine the MPS method with large eddy simulation(LES) which can simulate the turbulence in the flow. The intense pressure fluctuation is a significant shortcoming in MPS method. So, we improve the original MPS method by using a new pressure Poisson equation to ease the pressure fluctuation. Divergencefree condition representing fluid incompressible is used to calculate pressure smoothly. Then, area-time average technique is used to deal with the calculation. With these improvements, the modified MPS-LES method is applied to the simulation of numerical wave. As a contrast, we also use the original MPS-LES method to simulate the wave in a numerical wave tank. The result shows that the new method is better than the original MPS-LES method.

Numerical Study of Air Chamber for Oscillating Water Column Wave Energy Convertor

Oscillating Water Column （OWC） wave energy converting system is one of the most widely used facilities all over the world. The air chamber is utilized to convert the wave energy into the pneumatic energy. The numerical wave tank based on the two-phase VOF model is established in the present study toinvestigate the operating performance of OWC air chamber. The RANS equations, standard k-ε turbulence model and dynamic mesh technology are employed in the numerical model. The effects of incident wave conditions and shape parameters on the wave energy converting efficiency are studied and the capability of the present numerical wave tank on the corresponding engineering application is validated.

Numerical investigation of water-entry of flatted-bottom seafloor mining tool in ocean waves

A numerical wave load model based on two-phase(water-air) Reynolds-Averaged Navier Stokes(RANS) type equations is used to evaluate hydrodynamic forces exerted on flatted-bottom seafloor mining tool during its entering ocean waves of deploying process.The discretization of the RANS equations is achieved by a finite volume approach(FV).The volume of fluid method(VOF) is employed to track the complicated free surface.A numerical wave tank is built to generate the ocean waves which are suitable for deploying seafloor mining tool.A typical deploying condition is employed to reflect the process of flatted-bottom body impacting with waves,and the pressure distribution of bottom is also presented.Four different lowering velocities are applied to obtain the time histories of maximum pressure of bottom,and it can be concluded that the pressure coefficient decreases with water velocity increasing,which is similar with ordinary water entry case.The numerical results clearly demonstrate the characteristics of flatted-bottom body entering ocean waves.

Nonlinear numerical simulation on extreme-wave kine-matics

A fully nonlinear numerical model based on a time-domain higher-order boundary element method （HOBEM） is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here.

A simple method of depressing numerical dissipation effects during wave simulation within the Euler model

Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numerical wave tanks suffer from issues known as the numerical dissipation and numerical dispersion. The former causes wave energy to be slowly dissipated and the latter shifts wave frequencies during wave propagation. This paper proposes a simple method of depressing numerical dissipation effects on the basis of solving Euler equations using the finite difference method(FDM). The wave propagation solutions are solved analytically taking into account the influence of the damping terms. The main idea of the method is to append a source term to the momentum equation, whose strength is determined by how strong the numerical damping effect is. The method is verified by successfully depressing numerical effects during the simulation of regular linear waves, Stokes waves and irregular waves. By applying the method, wave energy is able to be close to its initial value after long distance of travel.

Numerical Simulation of Waves Generated by Seafloor Movements

Waves generated by vertical seafloor movements are simulated by use of a fully nonlinear two-dimensional numerical wave tank. In the souree region, the seafloor lifts to a designated height by a generation function. The numerical tests show that the linear theory is only valid for estimating the wave behaviors induced by the seafloor movements with a small amplitude, and the fully nonlinear numerical model should be adopted in the simulation of the wave generation by the large amplitude seafloor movements. Without the background surface waves, many numerical tests on the stable maximum elevations η0^max are carried out by beth the linear theory and the fully nonlinear model. The results of two models are compared and analyzed. For the fully nonlinear model, the influences of the amplitudes and the horizontal lengths on η^max are stronger than that of the characteristic duration times. Furthermore, results reveal that there are significant differences be- tween the linear theory and the fully nonlinear model. When the influences of the background surface waves are considered, the corresponding numerical analyses reveal that with the fully nonlinear model the η0^max near-linearly varies with the wave amplitudes of the surface waves, and the η0^max has significant dependences on the wave lengths and the wave phases of the surface waves. In addition, the differences between the linear theory and the fully nonlinear model are still obvious, and these differences are significantly affected by the wave parameters of the background surface waves, such as the wave amplitude, the wave length and the wave phase.

Numerical Investigation of the Effect of Current on Wave Focusing

In the present study, a numerical wave tank is developed to simulate the nonlinear wave-current interactions based on High Order Spectral （HOS） method. The influences of current on wave focusing are investigated by use of numerical model. The current is assumed to be constant in space. Focused waves with different amplitudes and frequency spectra are simulated with and without current. The focused wave characteristics, such as surface elevation, the maximum crest and frequency spectrum, with different current are compared. The results show that the opposing current increases the maximum crest and the energy transform during wave focusing process, and vice versa for the following current.

Comparative Experimental and Numerical Study of Wave Loads on A Monopile Structure Using Different Turbulence Models

This study numerically and experimentally investigates the effects of wave loads on a monopile-type offshore wind turbine placed on a 1:25 slope at different water depths as well as the effect of choosing different turbulence models on the efficiency of the numerical model.The numerical model adopts a two-phase flow by solving Unsteady Reynolds-Averaged Navier−Stokes(URANS)equations using the Volume Of Fluid(VOF)method and three differentk-ωturbulence models.Typical environmental conditions from the East China Sea are studied.The wave run-up and the wave loads applied on the monopile are investigated and compared with relevant experimental data as well as with mathematical predictions based on relevant theories.The numerical model is well validated against the experimental data at model scale.The use of different turbulence models results in different predictions on the wave height but less differences on the wave period.The baseline k-ωturbulence model and Shear-Stress Transport(SST)k-ωturbulence model exhibit better performance on the prediction of hydrodynamic load,at a model-scale water depth of 0.42 m,while the laminar model provides better results for large water depths.The SST turbulence model performs better in predicting wave run-up for water depth 0.42 m,while the laminar model and standard k-ωmodel perform better at water depth 0.52 m and 0.62 m,respectively.

RANS-VOF Solver for Solitary Wave Run-up on A Circular Cylinder

Simulation of solitary wave run-up on a vertical circular cylinder is carried out in a viscous numerical wave tank developed based on the open source codes Open FOAM. An incompressible two-phase flow solver naoe-FOAM-SJTU is used to solve the Reynolds-Averaged Navier–Stokes（RANS） equations with the SST k ？？ turbulence model. The PISO algorithm is utilized for the pressure-velocity coupling. The air-water interface is captured via Volume of Fluid（VOF） technique. The present numerical model is validated by simulating the solitary wave run-up and reflected against a vertical wall, and solitary wave run-up on a vertical circular cylinder. Comparisons between numerical results and available experimental data show satisfactory agreement. Furthermore, simulations are carried out to study the solitary wave run-up on the cylinder with different incident wave height H and different cylinder radius a. The relationships of the wave run-up height with the incident wave height H, cylinder radius a are analyzed. The evolutions of the scattering free surface and vortex shedding are also presented to give a better understanding of the process of nonlinear wave-cylinder interaction.