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.

A Boussinesq Equation-Based Model for Nearshore Wave Breaking

Based on the wave breaking model by Li and Wang (1999), this work is to apply Dally's analytical solution to the wave-height decay instead of the empirical and semi-empirical hypotheses of wave-height distribution within the wave breaking zone. This enhances the applicability of the model. Computational results of shoaling, location of wave breaking, wave-height decay after wave breaking, set-down and set-up for incident regular waves are shown to have good agreement with experimental and field data.

Numerical Simulation of Abnormal Wave Height Change Induced by Excavated Waterway with Boussinesq Equation

A nonlinear numerical model has been set up by use of Boussinesq Equation with finite difference method, and has been applied to the simulation of the abnormal change of wave height induced by excavated waterway. Numerical results demonstrate that the abnormal change of wave height is due to the adding of the reflected wave height induced by excavated waterway to the incident wave height. Because the angle between the incident wave and the axis of the waterway is smaller than the critical angle, the reflected wave produced by the waterway may propagate to the breakwater and may be added with the incident wave, then the abnormal change of wave height before the breakwater may be caused. So the wave reflection caused by the change of water depth cannot be neglected.

Numerical Models of Higher-Order Boussinesq Equations and Comparisons with Laboratory Measurement

Nonlinear water wave propagation passing a submerged shelf is studied experimentally and numerically. The applicability of two different wave propagation models has been investigated. One is higher-order Boussinesq equations derived by Zou (1999) and the other is the classic Boussinesq equations, Physical experiments are conducted, three different front slopes (1:10, 1:5 and 1:2) of the shelf are set up in the experiment and their effects on wave propagation are investigated. Comparisons of numerical results with test data are made, the model of higher-order Boussinesq equations agrees much better with the measurements than the model of the classical Boussinesq equations, The results show that the higher-order Boussinesq equations can also be applied to the steeper slope case although the mild slope assumption is employed in the derivation of the higher order terms of higher order Boussinesq equations.

Direct Numerical Simulation of Interaction Between Wave and PorousBreakwater Based on N-S Equation

In this paper, a numerical model is established. A modified N-S equation is used as a control equation for the wave field and porous flow area. The control equations are discreted and solved by the finite difference method. The free surface is tracked by the VOF method. The pressure field and velocity field of the whole flow area are solved by the reiterative iteration method. Finally, compared with the physical model test results of wave flume, the numerical model established in the present study is validated.

Numerical Calculation for Nonlinear Waves in Water of Arbitrarily Varying Depth with Boussinesq Equations

Based on the high order nonlinear and dispersive wave equation with a dissipative term, a numerical model for nonlinear waves is developed, It is suitable to calculate wave propagation in water areas with an arbitrarily varying bottom slope and a relative depth h/L(0)less than or equal to1. By the application of the completely implicit stagger grid and central difference algorithm, discrete governing equations are obtained. Although the central difference algorithm of second-order accuracy both in time and space domains is used to yield the difference equations, the order of truncation error in the difference equation is the same as that of the third-order derivatives of the Boussinesq equation. In this paper, the correction to the first-order derivative is made, and the accuracy of the difference equation is improved. The verifications of accuracy show that the results of the numerical model are in good agreement with those of analytical Solutions and physical models.

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.

Wave Numerical Model for Shallow Water

The history of forecasting wind waves by wave energy conservation equation Is briefly described. Several currently used wave numerical models for shallow water based on different wave theories are discussed. Wave energy conservation models for the simulation of shallow water waves are introduced, with emphasis placed on the SWAN model, which takes use of the most advanced wave research achievements and has been applied to several theoretical and field conditions. The characteristics and applicability of the model, the finite difference numerical scheme of the action balance equation and its source terms computing methods are described in detail. The model has been verified with the propagation refraction numerical experiments for waves propagating in following and opposing currents; finally, the model is applied to the Haian Gulf area to simulate the wave height and wave period field there, and the results are compared with observed data.

Applications of sediment sudden deposition model based on the third-generation numerical model for shallow water wave

The existing numerical models for nearshore waves are briefly introduced, and the third-generation numerical model for shallow water wave, which makes use of the most advanced productions of wave research and has been adapted well to be used in the environment of seacoast, lake and estuary area, is particularly discussed. The applied model realizes the significant wave height distribution at different wind directions. To integrate the model into the coastal area sediment, sudden deposition mechanism, the distribution of average silt content and the change of sediment sudden deposition thickness over time in the nearshore area are simulated. The academic productions can give some theoretical guidance to the applications of sediment sudden deposition mechanism for stormy waves in the coastal area. And the advancing directions of sediment sudden deposition model are prospected.

Numerical Model of Internal Solitary Wave Evolution on ImpermeableVariable Seabed in A Stratified Two-Layer Fluid System

In the present study a numerical model developed by Lynett and Liu （2002） is modified to include density difference in a stratified two-layer fluid in a three-dimensional internal wave domain. The internal solitary wave （ISW） in the model is assumed to be weakly nonlinear and weakly dispersive, and the viscosity effects at all boundaries are ignored. The governing equations based on the Navier-Stokes and Euler equations are solved for internal solitary wave propagation over variable seabed topography. Theoretical formulations are established, from which analytical solutions are obtained, in addition to numerical results. Wave profiles from previous experimental studies are compared with the numerical results from the present analytical solutions. Numerical models developed on the basis of the present analytical solutions are better than those developed by Lynett and Liu （2002）. The results of numerical modeling agree well with the experimental data.

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.

A Numerical Model for Nonlinear Wave Propagation on Non-uniform Current

On the basis of the new type Boussinesq equations （Madsen et al., 2002）, a set of equations explicitly including the effects of currents on waves are derived. A numerical implementation of the present equations in one dimension is described. The numerical model is tested for wave propagation in a wave flume of uniform depth with current present. The present numerical results are compared with those of other researchers. It is validated that the present numerical model can reasonably reflect the nonlinear influences of currents on waves. Moreover, the effects of inputting different incident boundary conditions on the calculated results are studied.

Analysis Model of Dynamic Response on the Circular Caisson Breakwater Under Random Wave

The random wave load is applied to dynamic response analysis of circular caisson breakwater. The motion process of circular caisson breakwater is classified as rotation motion mode and rotation-and-sliding motion mode. The dynamic model system composed of damper-antislider to control the lateral sliding is introduced, and corresponding dynamic equations of two motion modes are established. The formulas to calculate added mass and new conversion relation of the unit rota- tional stiffness coefficient are put forward according to the characteristic of the circular caisson breakwater. An engineering case is calculated by a program compiled in Fortran language using proposed dynamic model and method. The validity of the model is calibrated.

A Comparative Study of Numerical Models for Wave Propagation and Setup on Steep Coral Reefs

Complex factors including steep slopes, intense wave breaking, large bottom friction and remarkable wave setup should be considered while studying wave propagation over coral reefs, and how to simulate wave propagation and setup on coral reefs efficiently has become a primary focus. Several wave models can be used on coral reefs as have been published, but further testing and comparison of the reliability and applicability of these models are needed. A comparative study of four numerical wave models (i.e., FUNWAVE-TVD, Coulwave, NHWAVE and ZZL18) is carried out in this paper. These models’ governing equations and numerical methods are compared and analyzed firstly to obtain their differences and connections;then the simulation effects of the four wave models are tested in four representative laboratory experiments. The results show that all four models can reasonably predict the spectrum transformation. Coulwave, NHWAVE and ZZL18 can predict the wave height variation more accurately;Coulwave and FUNWAVE-TVD tend to underestimate wave setup on the reef top induced by spilling breaker, while NHWAVE and ZZL18 can predict wave setup relatively accurately for all types of breakers;NHWAVE and ZZL18 can predict wave reflection by steep reef slope more accurately. This study can provide evidence for choosing suitable models for practical engineering or establishing new models.

Solitary Wave Propagation Influenced by Submerged Breakwater

The form of Boussinesq equation derived by Nwogu （1993） using velocity at an arbitrary distance and surface elevation as variables is used to simulate wave surface elevation changes. In the numerical experiment, water depth was divided into five layers with six layer interfaces to simulate velocity at each layer interface. Besides, a physical experiment was carried out to validate numerical model and study solitary wave propagation.“Water column collapsing”method （WCCM） was used to generate solitary wave. A series of wave gauges around an impervious breakwater were set-up in the flume to measure the solitary wave shoaling, run-up, and breaking processes. The results show that the measured data and simulated data are in good agreement. Moreover, simulated and measured surface elevations were analyzed by the wavelet transform method. It shows that different wave frequencies stratified in the wavelet amplitude spectrum. Finally, horizontal and vertical velocities of each layer interface were analyzed in the process of solitary wave propagation through submerged breakwater.

Applying Finite Difference Method to Simulate the Performance of a Perforated Breakwater Under Regular Waves

Using a discretized finite difference method, a numerical model was developed to study the interaction of regular waves with a perforated breakwater. Considering a non-viscous, non-rotational fluid, the governing equations of Laplacian velocity potential were developed, and specific conditions for every single boundary were defined. The final developed model was evaluated based on an existing experimental result. The evaluated model was used to simulate the condition for various wave periods from 0.6 to 2 s. The reflection coefficient and transmission coefficient of waves were examined with different breakwater porosities, wave steepnesses, and angular frequencies. The results show that the developed model can suitably present the effect of the structural and hydraulic parameters on the reflection and transmission coefficients. It was also found that with the increase in wave steepness, the reflection coefficient increased logarithmically, while the transmission coefficient decreased logarithmically.

Numerical Simulation on Hydraulic Performances of Quarter Circular Breakwater

Quarter circular breakwater （QCB） is a new-type breakwater developed from senti-circular breakwater （SCB）. The superstructure of QCB is composed of a quarter circular front wall, a horizontal base slab and a vertical rear wall. The width of QCB＇ s base slab is about half that of SCB, which makes QCB suitable to be used on relatively finn soil foundation. The numerical wave flume based on the Reynolds averaged Navier-Stokes equations for impressible viscosity fluid is adopted in this paper to simulate the hydraulic performances of QCB. Since the geometry of both breakwaters is similar and SCB has been studied in depth, the hydraulic performances of QCB are given in comparison with those of SCB.

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.

3-D acoustic wave equation forward modeling with topography

In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties.

Verification of A Numerical Harbour Wave Model

A numerical model for wave propagation in a harbour is verified by use of physical models. The extended time-dependent mild slope equation is employed as the governing equation, and the model is solved by use of ADI method containing the relaxation factor. Firstly, the reflection coefficient of waves in front of rubble-mound breakwaters under oblique incident waves is determined through physical model tests, and it is regarded as the basis for simulating partial reflection boundaries of the numerical model. Then model tests on refraction, diffraction and reflection of waves in a harbour are performed to measure wave height distribution. Comparative results between physical and numerical model tests show that the present numerical model can satisfactorily simulate the propagation of regular and irregular waves in a harbour with complex topography and boundary conditions.