Analytical Method for the Wave Diffraction of Asymmetrically Arranged Breakwaters

The layout forms of several breakwater structures can be generalized as asymmetrical arrangements in actual engineering.However,the problem of wave diffraction around asymmetrically arranged breakwaters has not been adequately investigated.In this study,we propose an analytical method of wave diffraction for regular waves passing through asymmetrically arranged breakwaters,and we use the Nyström method to obtain the analytical solution numerically.We compared the results of this method with those of previous analytical solutions and with numerical results to demonstrate the validity of our approach.We also provided diffraction coefficient diagrams of breakwaters with different layout forms.Moreover,we described the analytical expression for the problem of diffraction through long-wave incident breakwaters and presented an analysis of the relationship between the diffraction coefficients and the widths of breakwater gates.The analytical method presented in this study contributes to the limited literature on the theory of wave diffraction through asymmetrically arranged breakwaters.

Time-Domain Higher-Order Boundary Element Method for Simulating High Forward-Speed Ship Motions in Waves

The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.

An empirical method for joint inversion of wave and wind parameters based on SAR and wave spectrometer data

Synthetic aperture radar(SAR)and wave spectrometers,crucial in microwave remote sensing,play an essential role in monitoring sea surface wind and wave conditions.However,they face inherent limitations in observing sea surface phenomena.SAR systems,for instance,are hindered by an azimuth cut-off phenomenon in sea surface wind field observation.Wave spectrometers,while unaffected by the azimuth cutoff phenomenon,struggle with low azimuth resolution,impacting the capture of detailed wave and wind field data.This study utilizes SAR and surface wave investigation and monitoring(SWIM)data to initially extract key feature parameters,which are then prioritized using the extreme gradient boosting(XGBoost)algorithm.The research further addresses feature collinearity through a combined analysis of feature importance and correlation,leading to the development of an inversion model for wave and wind parameters based on XGBoost.A comparative analysis of this model with ERA5 reanalysis and buoy data for of significant wave height,mean wave period,wind direction,and wind speed reveals root mean square errors of 0.212 m,0.525 s,27.446°,and 1.092 m/s,compared to 0.314 m,0.888 s,27.698°,and 1.315 m/s from buoy data,respectively.These results demonstrate the model’s effective retrieval of wave and wind parameters.Finally,the model,incorporating altimeter and scatterometer data,is evaluated against SAR/SWIM single and dual payload inversion methods across different wind speeds.This comparison highlights the model’s superior inversion accuracy over other methods.

A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

Structural Strength Study of A Floating Wind Turbine Under Freak Waves Through the CFD-FEA Method

In recent years,the exploitation of offshore wind resources has been attached with greater importance.As a result,semi-submersible floating wind turbines(FWTs)have gradually become a popular research topic,with the structural strength being a research hotspot as it can ensure the safe operation of FWTs.The severe sea conditions of freak waves result in enormous wave heights,extremely fast wave speeds,and concentrated energy.Thus,it is difficult to accurately simulate these effects on the loads of floating wind turbines using the potential flow theory and other theories.In this paper,the structural strength of a floating wind turbine under the action of freak waves is analyzed based on the CFD-FEA coupled method.The effects of the mooring system and the wind load are considered in the time domain,and the CFD method is applied to analyze the wave load of the floating wind turbine under the extreme sea state of freak waves.The strength and motion of the floating wind turbine float structure are analyzed by combining the CFD method and the FEA method,and the analytical results of the mutual transfer of these two methods are taken as the initial quantities for further analysis.The accuracy of the analytical model of the CFD-FEA method is verified by the results of the tank test analysis,and the structural strength analysis under freak wave conditions is carried out for a new type of floating wind turbine.The results of this research provide useful guidance and references for the design and engineering applications of offshore floating wind turbines.

An investigation on the wind profiles and gravity wave dynamics in MLT region based on the meteor radars from the Meridian Project

The meteor radar can detect the zenith angle,azimuth,radial velocity,and altitude of meteor trails so that one can invert the wind profiles in the mesosphere and low thermosphere(MLT)region,based on the Interferometric and Doppler techniques.In this paper,the horizontal wind field,gravity wave(GW)disturbance variance,and GW fluxes are analyzed through the meteor radar observation from 2012−2022,at Mohe(53.5°N,122.4°E)and Zuoling(30.5°N,114.6°E)stations of the(Chinese)Meridian Project.The Lomb−Scargle periodogram method has been utilized to analyze the periodic variations for time series with observational data gaps.The results show that the zonal winds at both stations are eastward dominated,while the meridional winds are southward dominated.The variance of GW disturbances in the zonal and meridional directions increases gradually with height,and there is a strong pattern of annual variation.The zonal momentum flux of GW changes little with height,showing weak annual variation.The meridional GW flux varies gradually from northward to southward with height,and the annual periodicity is stronger.For both stations,the maximum values of zonal and meridional wind occur close to the peak heights of GW flux,with opposite directions.This observational evidence is consistent with the filtering theory.The horizontal wind velocity,GW flux,and disturbance variance of the GW at Mohe are overall smaller than those at Zuoling,indicating weaker activities in the MLT at Mohe.The power spectral density(PSD)calculated by the Lomb−Scargle periodogram shows that there are 12-month period and 6-month period in horizontal wind field,GW disturbance variance and GW flux at both stations,and especially there is also a 4-month cycle in the disturbance variance.The PSD of the 12-month and 6-month cycles exhibits maximum values below 88 km and above 94 km.

A weak formulation of heterogenous viscoacoustic wave propagation in infinite domain

The accurate simulation of wave propagation in real media requires properly taking the attenuation into account,which leads to wave dissipation together with its causal companion,wave dispersion.In this study,to obtain a weak formulation of heterogenous viscoacoustic wave propagation in an infinite domain,the viscoacoustic medium is first characterized by its frequency-dependent complex bulk compliance instead of the classically used complex bulk modulus.Then,a mechanical model using serially connected standard linear solids(SSLS)is built to obtain the rational approximation of the complex bulk compliance whose parameters are calculated via an adapted nonlinear optimization method.Utilizing the obtained bulk compliance-based constitutive relation,a novel second-order viscoacoustic wave equation in the frequency domain is derived,of which the weak formulation can be physically explained as the virtual work equation and can thus be discretized using a continuous spectral element method in space.Additionally,a new method is introduced to address the convolution terms involved in the inverse Fourier transform,whose accurate time integration can then be achieved using an explicit time scheme,which avoids the transient growth that exists in the classical method.The resulting full time-space decoupling scheme can handle wave propagation in arbitrary heterogeneous media.Moreover,to treat the wave propagation in an infinite domain,a perfectly matched layer in weak formulation is derived for the truncation of the infinite domain via complex coordinate stretching of the virtual work equation.With only minor modification,the resulting perfectly matched layer can be implemented using the same time scheme as for the wave equation inside the truncated domain.The accuracy,numerical stability,and versatility of the new proposed scheme are demonstrated with numerical examples.

Joint 3D traveltime calculation based on fast marching method and wavefront construction

3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.

Multi-scale seismic full waveform inversion in the frequency-domain with a multi-grid method

Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.

NUMERICAL INVESTIGATION OF TOROIDAL SHOCK WAVES FOCUSING USING DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD

A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

Calculation of Valence Subband Structures for Strained Quantum-Wells by Plane Wave Expansion Method Within 6×6 Luttinger-Kohn Model

The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6×6 Luttinger Kohn model.The effect of the number and period of plane waves used for expansion on the stability of energy eigenvalues is examined.For practical calculation,it should choose the period large sufficiently to ensure the envelope functions vanish at the boundary and the number of plane waves large enough to ensure the energy eigenvalues keep unchanged within a prescribed range.

3D elastic wave equation forward modeling based on the precise integration method

The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.

Seismic wave input method for three-dimensional soil-structure dynamic interaction analysis based on the substructure of artificial boundaries

The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.

Numerical simulation of solitary and random wave propagation through vegetation based on VOF method

A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.

Focused Wave Properties Based on A High Order Spectral Method with A Non-Periodic Boundary

In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.

ASYMPTOTIC METHOD OF TRAVELLING WAVE SOLUTIONS FOR A CLASS OF NONLINEAR REACTION DIFFUSION EQUATION

In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for the corresponding problem is obtained.

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.

Study on Measurement Methods of Multi-Directional Waves—Wave Gauge Array Method

Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orientations, etc. of the three array arrangements, i. e., linear array, T-type array and pentagon array, on the resolution of the directional spreading of waves, are investigated experimentally. This study can be used as a reference in the experimental study and the field measurement of directional waves.