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.

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.

Air Gap Prediction for Floating Bodies Using a 3D Numerical Wave Tank Approach

Computations for air gap response of a semisubmersible platform based on a 3D numerical wave tank approach are presented.The developed method is in time domain and can consider nonlinearities associated with incident wave and hydrostatic forces exactly in determining the body response, but the interaction hydrodynamics of radiation and diffraction are based on simplified linearization assumptions. The incident wave can be defined by any suitable wave theory and here defined by a fully nonlinear numerical wave model. After verifying the present computations results in its degenerated linearized version against the usual linear 3D Green function–based frequency-domain results for air gap predictions, systematic comparative studies are undertaken between linear and the approximate nonlinear solutions. It is found that nonlinear computations can yield considerably conservative predictions as compared to fully linear calculations, amounting to a difference of up to 30%–40% in the minimum air gap in steep ambient incident waves at high and moderate frequencies.

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.

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.

DYNAMIC MODEL FOR OIL SLICK DISPERSION INTO A WATER COLUMN- A WINDDRIVEN WAVETANK EXPERIMENT

A dynamic experiment for oil dispersion into a water column was performed with a 21 m long, 0.5 m wide, and 1 m high wind-driven wave tank. At wind velocity between 6-12 m/s and with the oil slide kept constant (about 1 um), the rate of the oil content increase in the water column could be approximated from the difference between the dispersion rate (R) of the oil slick and the coagulation rate (R’) of the dispersed oil slick. Assuming the coagulation rate is directly proportional to the concentration of the water dispersed oil slick (i. e. R’ =KC),, the integral form of the dynamic model can be expressed as C=R*[1-exp(-K*t)]/K and parameters R and K can be regressed with a computer. The relative deviation of model results from the experimental data was mainly less than 10%. The oil slick dispersion rate (R) had exponential relationship with the wind velocity (V), and can be fitted with a formula R=A*(U+1)B.The fitted constant of the coagulation rate, K(0.8-3.0* 10-3 min-1) did not have significant

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.

Evaluation of Hydrodynamic Coefficients of Ships Oscillating in A Numerical Wave Tank

A technique for the evaluation of the hydrodynamic coefficients of ships is outlined for ship oscillating in a numerical wave tank, which is established on Computational Fluid Dynamics （CFD） theories. The numerical simulation of ship sections and bodies forced oscillating in the tank are carried out. The added mass and damping coefficients are obtained by the decomposition of the computational results, which agree well with the corresponding ones of potential theories.

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.

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.

Explosion damage effects of aviation kerosene storage tank under strong ignition

In order to study the blast damage effects of aviation kerosene storage tanks,the out-field explosion experiments of 8 m3fixed-roof tanks were carried out.The fragments,shock wave and fireball thermal radiation of the tank in the presence of bottom oil,half oil and full oil,as well as empty tank,were investigated under internal explosion by various TNT charge contents(1.8 kg,3.5 kg and 6.2 kg).The results showed that the tank roof was the only fragment produced,and the damage forms could be divided into three types.The increase of TNT charge content and oil volume enlarged the deformation of the tank,while the hole ratio presented a trend of increase first and then decrease.The H_r,maxand V_(max)values positively increased as increasing the TNT charge content and oil volume(from empty to half oil),but decreased in full oil.The Pmaxvalues had a progressive increase with the increment of TNT charge content,but not the case with the increase in oil volumes.The development of fireball was divided into three stages:tank roof‘towed'flame,jet flow flame tumbling and rising,and jet flow flame extinguishing.The Dmaxand Hf,maxvalues both increased as increasing TNT charge content and oil volumes.The oscillation phenomenon of fireball temperature was observed in the cooling process.The average temperature of fireball surface was positively correlated with TNT charge content,and negatively correlated with oil volumes.

Poststack reverse-time migration using a non-reflecting recursive algorithm on surface relief

Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be one of highly accurate methods. In this paper, we propose a new non-reflecting recursive algorithm for reverse-time migration by introducing the wave impedance function into the acoustic wave equation and the algorithm for the surface relief case is derived from the coordinate transformation principle. Using the exploding reflector principle and the zero-time imaging condition of poststack reverse- time migration, poststack numerical simulation and reverse-time migration with complex conditions can be realized. The results of synthetic and real data calculations show that the method effectively suppresses unwanted internal reflections and also deals with the seismic imaging problems resulting from surface relief. So, we prove that this method has strong adaptability and practicality.

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.

Efficient Generation of Freak Waves in Laboratory

In the present study, Kriebel＇ s method is improved to generate freak waves in laboratory. The improved method superposes a random wave train with two transient wave trains to simulate freak wave events in a wave tank. The freak waves are more nonlinear than what generated with Kriebel＇ s method of the same energy. It can ＇also generate freak waves to satisfy all the qualifications of the adopted definition with less energy than Kriebel＇ s and can hardly influence the significant wave height.

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.

Efficient computation method for two-dimensional nonlinear waves

The theory and simulation of fully-nonlinear waves in a truncated two-dimensional wave tank in time domain are presented. A piston-type wave-maker is used to generate gravity waves into the tank field in finite water depth. A damping zone is added in front of the wave-maker which makes it become one kind of absorbing wave-maker and ensures the prescribed Neumann condition. The efficiency of numerical tank is further enhanced by installation of a sponge layer beach (SLB) in front of downtank to absorb longer weak waves that leak through the entire wave train front. Assume potential flow, the space- periodic irrotational surface waves can be represented by mixed Euler-lagrange particles Solving the integral equation at each time step for new normal velocities, the instantaneous free surface is integrated following time history by use of fourth-order Runge- Kutta method. The double node technique is used to deal with geometric discontinuity at the wave- body intersections. Several precise smoothing methods have been introduced to treat surface point with high curvature. No saw-tooth like instability is observed during the total simulation. The advantage of proposed wave tank has been verified by comparing with linear theoretical solution and other nonlinear results, excellent agreement in the whole range of frequencies of interest has been obtained.

Validation of a CIP-based tank for numerical simulation of free surface flows

A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the present paper.The tangent of hyperbola for interface capturing（THINC） scheme is applied for capturing complex free surfaces.The new model is capable of simulating a flow with violently varied free surface.A series of computations are conducted to assess the developed algorithm and its versatility.These tests include the collapse of water column with and without an obstacle,sloshing in a fixed tank,the generation of regular waves in a tank,the generation of extreme waves in a tank.Excellent agreements are obtained when numerical results are compared with available analytical,experimental,and other numerical results.

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.