Numerical Simulation and Parallel Computing of Acoustic Wave Equation in Isotropic-Heterogeneous Media

In this paper,we consider the numerical implementation of the 2D wave equation in isotropic-heterogeneous media.The stability analysis of the scheme using the von Neumann stability method has been studied.We conducted a study on modeling the propagation of acoustic waves in a heterogeneous medium and performed numerical simulations in various heterogeneous media at different time steps.Developed parallel code using Compute Unified Device Architecture(CUDA)technology and tested on domains of various sizes.Performance analysis showed that our parallel approach showed significant speedup compared to sequential code on the Central Processing Unit(CPU).The proposed parallel visualization simulator can be an important tool for numerous wave control systems in engineering practice.

Simulating the evolution of focused waves by a two-layer Boussinesq-type model

Accurate simulation of the evolution of freak waves by the wave phase focusing method requires accurate linear and nonlinear properties,especially in deep-water conditions.In this paper,we analyze the ability to simulate deep-water focused waves of a two-layer Boussinesq-type(BT)model,which has been shown to have excellent linear and nonlinear performance.To further improve the numerical accuracy and stability,the internal wavegenerated method is introduced into the two-layer Boussinesq-type model.Firstly,the sensitivity of the numerical results to the grid resolution is analyzed to verify the convergence of the model;secondly,the focused wave propagating in two opposite directions is simulated to prove the symmetry of the numerical results and the feasibility of the internal wave-generated method;thirdly,the limiting focused wave condition is simulated to compare and analyze the wave surface and the horizontal velocity of the profile at the focusing position,which is in good agreement with the measured values.Meanwhile the simulation of focused waves in very deep waters agrees well with the measured values,which further demonstrates the capability of the two-layer BT model in simulating focused waves in deep waters.

A Hybrid Dung Beetle Optimization Algorithm with Simulated Annealing for the Numerical Modeling of Asymmetric Wave Equations

In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.

Numerical Modelling of Wash Waves Generated by Ships Moving over An Uneven Bottom

Unsteady wash waves generated by a ship with constant speed moving across an uneven bottom topography are investigated by numerical simulations based on a Mixed Euler–Lagrange(MEL) method. The transition is accomplished by the ship traveling from the depth h1 into the depth h2 via a step bottom. A small tsunami would be created after this transition. However, the unsteady wave-making resistance induced by this new phenomenon has not been well documented by literature. Therefore, the main purpose of the present study is to quantify the effects of an uneven bottom on the unsteady wash waves and wave-making resistance acting on the ship. An upwind differential scheme is commonly used in the Euler method to deal with the convection terms under free-surface condition to prevent waves in the upstream. Evidently, it cannot be applied to the present problem due to upstream waves generated by the ship would be dampened by the upwind scheme. The central differential scheme provides more accurate results,but it is not unconditionally stable. An MEL method is therefore employed to investigate the upstream wave generated by the ship moving over the uneven bottom. Simulation results show that the hydrodynamic interaction between the ship and the uneven bottom could initiate an upstream tsunami, as well as unsteady wave-making resistance on ships.The unsteady wave-making resistance oscillates periodically, and the amplitude and period of the oscillations are highly dependent on speed and water depth.

A Study of Traveling Wave Structures and Numerical Investigation of Two-Dimensional Riemann Problems with Their Stability and Accuracy

The Riemann wave system has a fundamental role in describing waves in various nonlinear natural phenomena,for instance,tsunamis in the oceans.This paper focuses on executing the generalized exponential rational function approach and some numerical methods to obtain a distinct range of traveling wave structures and numerical results of the two-dimensional Riemann problems.The stability of obtained traveling wave solutions is analyzed by satisfying the constraint conditions of the Hamiltonian system.Numerical simulations are investigated via the finite difference method to verify the accuracy of the obtained results.To extract the approximation solutions to the underlying problem,some ODE solvers in FORTRAN software are applied,and outcomes are shown graphically.The stability and accuracy of the numerical schemes using Fourier’s stabilitymethod and error analysis,respectively,to increase the reassurance are investigated.A comparison between the analytical and numerical results is obtained and graphically provided.The proposed methods are effective and practical to be applied for solving more partial differential equations(PDEs).

Numerical Investigations on the Transient Behavior of Sand Waves in Beibu Gulf Under Normal and Extreme Sea Conditions

In this study, a morphodynamic numerical model is established with the Regional Ocean Modeling System(ROMS)to investigate the transient behavior of sand waves under realistic sea conditions. The simulation of sand wave evolution comprises two steps: 1) a regional-scale model is configured first to simulate the ocean hydrodynamics, i.e., tides and tidal currents, and 2) the transient behavior of sand waves is simulated in a small computational domain under the time-variant currents extracted from the large model. The evolution of sand waves on the continental shelf in the Beibu Gulf is specifically investigated. The numerical results of the two-year evolution of sand waves under normal sea conditions compare well with the field survey data. The transient behavior of sand waves in individual months shows that the sand waves are more stable in April and October than that in other months, which can be selected as the windows for seabed operations. The effects of sediment properties, including settling velocity, critical shear stress and surface erosion rate, on sand wave evolution are also analyzed. Then, the typhoon-induced currents are further superimposed on the tidal currents as the extreme weather conditions. Sand waves with the average wavelength generally have more active behavior than smaller or larger sand waves. The characteristics of the evolution of sand waves in an individual typhoon process are quite different for different hydrodynamic combinations. For the storm conditions, i.e., the real combination and maximum combination cases, the sand waves experience a significant migration together with a damping in height due to the dominant suspended sediment transport. For the mild conditions, i.e., the pure tidal current and minimum combination cases, the sand waves migrate less, but the heights continue growing due to the dominant bedload transport.

Linear and Nonlinear Stokes Waves Theory: Numerical Hydrodynamic and Energy Studies

The increase of wave energy in electricity production is an objective shared by many countries to meet growing demand and global warming. To analyze devices capable of converting the energy of sea waves into electrical energy, it is important to master the various theories of gravity waves and generation. We will in our work consider a numerical waves tank for an amplitude A=0.5, a wavelength λ=0.25 , an average height He=10 and a Froude number fixed at 1 × 105. Numerical wave channel analysis is used to reproduce the natural phenomenon of wave propagation in an experimental model. Wave makers are usually used to generate waves in the channel. In theory, the influence of an incident wave can be considered, as in the case of our study. In this study, the evolution of the hydrodynamic parameters and the energy transported in one wavelength can be determined by calculation. A change of variable will be done in this work to facilitate the writing of the boundary conditions at the free surface and at the bottom. The nonlinear Stokes theory will be studied in this case in order to provide hydrodynamic solutions through the Navier-Stokes equations to finally deduce the energetic results. To do this, the finite difference method will be used for the hydrodynamic results such as the velocity potential and the free surface elevation and the trapezium method of Newton for the energetic results. Thus, we will determine the energetic potential according to the decrease in the slope of the tank. To do this, we will take as values of beta representing the inverse of the slope of the tank, β=100, β=105, β=110 and β=105. .

Numerical Simulation of Nonlinear Three-Dimensional Waves in Water of Arbitrary Varying Topography

The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.

2-D Composite Model for Numerical Simulations of Nonlinear Waves

A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.

A numerical study of generation and propagation of type-a and type-b internal solitary waves in the northern South China Sea

Numerical simulations based on a high-resolution three-dimensional MIT general circulation model(MITgcm)using realistic topography and tidal forcing are conducted to investigate the generation and propagation of the so-called type-a waves(large-amplitude rank-ordered wave packets)and type-b waves(isolated wave packets)in the northern South China Sea.At first,we summarized and analyzed the generation and propagation characteristics of these waves.Then,energy budget at the Luzon Strait is calculated.Energy generation has three local maxima every day,of which the largest one corresponds to the emergence of the type-a wave.Energy flux at the west boundary of the Luzon Strait shows two local maxima each day.The larger one is consistent with the generation of the type-a wave and the smaller one is in correspondence with the generation of the type-b wave.Sensitivity experiments are designed to explore the role of the east and west ridge of the Luzon Strait on the generation and propagation of the type-a and type-b waves.It is found that the east ridge is indispensable on the generation of the type-a wave while the west ridge has little contribution.The west ridge diminishes the type-a waves'amplitude but hardly changes their propagation speed.The type-b waves also come from perturbation signals which originate from the east ridge and are enhanced in amplitude and reduced in propagation speed by the west ridge.

Numerical simulation of sediment transport in coastal waves and wave-induced currents

Prediction of coastal sediment transport is of particularly importance for analyzing coast erosion accurately and solving the corresponding coast protection engineering problems.The present study provided a numerical scheme for sediment transport in coastal waves and wave-induced currents.In the scheme,the sand transport model was implemented with wave refraction-diffraction model and near-shore current model.Coastal water wave was simulated by using the parabolic mild-slope equation in which wave refraction,diffraction and breaking effects are considered.Wave-induced current was simulated by using the nonlinear shallow water equations in which wave provides radiation stresses for driving current.Then,sediment transport in waves and wave-induced currents was simulated by using the two-dimensional suspended sediment transport equations for suspended sediment and the bed-load transport equation for bed load.The numerical scheme was validated by experiment results from the Large-scale Sediment Transport Facility at the US Army Corps of Engineer Research and Development Center in Vicksburg.The numerical results showed that the present scheme is an effective tool for modeling coastal sediment transport in waves and near-shore currents.

A VOF-based numerical model for breaking waves in surf zone

This paper introduces a numerical model for studying the evolution of a periodic wave train, shoaling, and breaking in surf zone. The model can solve the Reynolds averaged Navier-Stokes （RANS） equations for a mean flow, and the k-ε equations for turbulence kinetic energy k and turbulence dissipation rate ε. To track a free surface, the volume of fluid （VOF） function, satisfying the advection equation was introduced. In the numerical treatment, third-order upwind difference scheme was applied to the convection terms of the RANS equations in order to reduce the effect of numerical viscosity. The shoaling and breaking processes of a periodic wave train on gently sloping beaches were modeled. The computed wave heights of a sloping beach and the distribution of breaking wave pressure on a vertical wall were compared with laboratory data.

Numerical study on water waves and wave-induced longshore currents in Obaky coastal water

In this paper, the water waves and wave-induced longshore currents in Obak6y coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic mild-slope equation for coastal water waves and the nonlinear shallow water equation for the wave-induced currents. The wave transformation under the effects of shoaling, refraction, diffraction and breaking is considered, and the wave provides radiation stresses for driving currents in the model. The numerical results for the water wave-induced longshore currents were validated by the measured data to demonstrate the efficiency of the numerical model. Then the water waves and longshore currents induced by the waves from main directions were numerically simulated and analyzed based on the numerical results. The numerical results show that the movement of the longshore currents was different while the wave proDaRated to a coastal zone from different directions.

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 STUDY OF THE INFLUENCE OF WAVES AND TIDE-SURGE INTERACTION ON TIDE-SURGES IN THE BOHAI SEA

The author’s combined numerical model consisting of a third generation shallow water wave model and a 3 D tide surge model with wave dependent surface wind stress were used to study the influence of waves on tide surge motion. For the typical weather case, in this study, the magnitude and mechanism of the influence of waves on tide surges in the Bohai Sea were revealed for the first time. The results showed that although consideration of the wave dependent surface wind stresses raise slightly the traditional surface wind stress, due to the accumulated effects, the computed results are improved on the whole. Storm level maximum modulation can reach 0.4 m. The results computed by the combined model agreed well 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 second order volume of fluid (VOF) scheme for numerical simulation of 2-D breaking waves

Among all environmental forces acting on ocean structures and marine vessels, those resulting from wave impacts are likely to yield the highest loads. Being highly nonlinear, transient and complex, a theoretical analysis of their impact would be impossible without numerical simulations. In this paper, a pressure-split two-stage numerical algorithm is proposed based on Volume Of Fluid （VOF） methodology. The algorithm is characterized by introduction of two pressures at each half and full cycle time step, and thus it is a second-order accurate algorithm in time. A simplified second-order Godunov-type solver is used for the continuity equations. The method is applied to simulation of breaking waves in a 2-D water tank, and a qualitative comparison with experimental photo observations is made. Quite consistent results are observed between simulations and experiments. Commercially available software and Boundary Integral Method （BIM） have also been used to simulate the same problem. The results from present code and BIM are in good agreement with respect to breaking location and timing, while the results obtained from the comrnercial software which is only first-order accurate in time has clearly showed a temporal and spatial lag, verifying the need to use a higher order numerical scheme.

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.

Numerical study on wave loads and motions of two ships advancing in waves by using three-dimensional translating-pulsating source

A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course in waves. Two experiments are carried out respectively to mea- sure the wave loads and the free motions for a pair of side-by- side arranged ship models advancing with an identical speed in head regular waves. For comparison, each model is also tested alone. Predictions obtained by the present solution are found in favorable agreement with the model tests and are more accurate than the traditional method based on the three dimensional pulsating （3DP） source Green function. Numer- ical resonances and peak shift can be found in the 3DP pre- dictions, which result from the wave energy trapped in the gap between two ships and the extremely inhomogeneous wave load distribution on each hull. However, they can be eliminated by 3DTP, in which the speed affects the free sur- face and most of the wave energy can be escaped from the gap. Both the experiment and the present prediction show that hydrodynamic interaction effects on wave loads and free motions are significant. The present solver may serve as a validated tool to predict wave loads and motions of two ves- sels under replenishment at sea, and may help to evaluate the hydrodynamic interaction effects on the ships safety in replenishment operation.

Numerical Simulation of Freak Waves Based on the Four-Order Nonlinear Schrdinger Equation

A numerical wave model based on the modified four-order nonlinear Schoedinger （NKS） equation in deep water is developed to simulate freak waves. A standard split-step, pseudo-spectral method is used to solve NLS equation. The validation of the model is firstly verified, and then the simulation of freak waves is perforated by changing sideband condi- tions. Results show that freak waves entirely consistent with the definition in the evolution of wave trains are obtained. The possible occurrence mechanism of freak waves is discussed and the relevant characteristics are also analyzed.