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.

Modeling and Implementation of a Data Security and Protection Medium Using the Generated Key Based on Electromagnetic Wave Propagation Theories

Today, the advent of quantum computers and algorithms is calling into question the semantic security of symmetrical and asymmetrical cryptosystems. The security of objects connected to the network, which must provide a security service and protect the privacy of users by providing protection against attacks such as identity theft, denial of service, eavesdropping and unauthorised access to personal and sensitive data. It is therefore necessary to find a robust method of using the key that is effective in protecting and preventing data tampering. In this paper, we design and implement a security and data protection method using a key generated on the basis of electromagnetic wave propagation theories. Modelling and implementation of a data security and protection method using a key generated on the basis of electromagnetic wave propagation theories.

Three-dimensional acoustic wave equation modeling based on the optimal finite-difference scheme

Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.

Contact Analysis and Modeling of Standing Wave Linear Ultrasonic Motor

A contact model for describing the contact mechanics between the stator and slider of the standing wave linear ultrasonic motor was presented. The proposed model starts from the assumption that the vibration characteristics of the stator is not affected by the contact process. A modified friction models was used to analyze the contact problems. Firstly, the dynamic normal contact force, interface friction force, and steady-state characteristics were analyzed. Secondly, the influences of the contact layer material, the dynamic characteristics of the stator, and the pre-load on motor performance were simulated. Finally, to validate the contact model, a linear ultrasonic motor based on in-plane modes was used as an example. The corresponding results show that a set of simulation of motor performances based on the proposed contact mechanism is in good agreement with experimental results. This model is helpful to understanding the operation principle of the standing wave linear motor and thus contributes to the design of these types of motor.

Revisiting study on Boussinesq modeling of wave transformation over various reef profiles

To better understand the complex process of wave transformation and associated hydrodynamics over various fringing reef profiles, numerical experiments were conducted with a one-dimensional （1D） Boussinesq wave model. The model is based on higher-order Boussinesq equations and a higher-accuracy finite difference method. The dominant energy dissipation in the surf zone, wave breaking, and bottom friction were considered by use of the eddy viscosity concept and quadratic bottom friction law, respectively. Numerical simulation was conducted for a wide range of wave conditions and reef profiles. Good overall agreement between the computed results and the measurements shows that this model is capable of describing wave processes in the fringing reef environment. Numerical experiments were also conducted to track the source of underestimation of setup for highly nonlinear waves. Linear properties （including dispersion and shoaling） are found to contribute little to the underestimation; the low accuracy in nonlinearity and the ad hoc method for treating wave breaking may be the reason for the problem.

Modeling the interaction of an internal solitary wave with a sill

A nonhydrostatic numerical model was developed and numerical experiments performed on the interaction of an internal solitary wave （ISW） with a sill, for a two-layer fluid with a diffusive interface. Based on the blocking parameter （Br）, the flow was classified into three cases： （1） when bottom topography has little influence on the propagation and spatial structure of the ISW （Br〈0.5）, （2） where the ISW is distorted significantly by the blocking effect of the topography （though no wave breaking occurs, （0.5〈Br〈0.7）, and （3） where the ISW is broken as it encounters and passes over the bottom topography （0.7〈Br）. The numerical results obtained here are consistent with those obtained in laboratory experiments. The breaking process of the incident ISW when Br=0.7 was completely reproduced. Dissipation rate was linearly related to the blocking parameter when B,〈0.7, and the maximum dissipation rate could reach about 34% as Br raised to about 1.0. After that, instead of breaking, more reflection happened. Similarly, breaking induced mixing was also most effective during Br around 1.0, and can be up to 0.16.

The Numerical Modeling of 3-D Elastic Wave Equation Using a High-Order,Staggered-Grid, Finite Difference Scheme

This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.

Self-consistent three-dimensional modeling and simulation of large-scale rectangular surface-wave plasma source

A self-consistent and three-dimensional （3D） model of argon discharge in a large-scale rectangular surface-wave plasma （SWP） source is presented in this paper, which is based on the finite-difference time-domain （FDTD） approximation to Maxwell＇s equations self-consistently coupled with a fluid model for plasma evolution. The discharge characteristics at an input microwave power of 1200 W and a filling gas pressure of 50 Pa in the SWP source are analyzed. The simulation shows the time evolution of deposited power density at different stages, and the 3D distributions of electron density and temperature in the chamber at steady state. In addition, the results show that there is a peak of plasma density approximately at a vertical distance of 3 cm from the quartz window.

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.

Modeling of random wave transformation with strong wave-induced coastal currents

The propagation and transformation of multi-directional and uni-directional random waves over a coast with complicated bathymetric and geometric features are studied experimentally and numerically. Laboratory investigation indicates that wave energy convergence and divergence cause strong coastal currents to develop and inversely modify the wave fields. A coastal spectral wave model, based on the wave action balance equation with diffraction effect （WABED）, is used to simulate the transformation of random waves over the complicated bathymetry. The diffraction effect in the wave model is derived from a parabolic approximation of wave theory, and the mean energy dissipation rate per unit horizontal area due to wave breaking is parameterized by the bore-based formulation with a breaker index of 0.73. The numerically simulated wave field without considering coastal currents is different from that of experiments, whereas model results considering currents clearly reproduce the intensification of wave height in front of concave shorelines.

Maximum initial primary wave model for low-Froude-number reservoir landslides based on wave theory

The impulse waves induced by large-reservoir landslides can be characterized by a low Froude number.However,systematic research on predictive models specifically targeting the initial primary wave is lacking.Taking the Shuipingzi 1#landslide that occurred in the Baihetan Reservoir area of the Jinsha River in China as an engineering example,this study established a large-scale physical model(with dimensions of 30 m×29 m×3.5 m at a scale of 1:150)and conducted scaled experiments on 3D landslide-induced impulse waves.During the process in which a sliding mass displaced and compressed a body of water to generate waves,the maximum initial wave amplitude was found to be positively correlated with the sliding velocity and the volume of the landslide.With the increase in the water depth,the wave amplitude initially increased and then decreased.The duration of pressure exertion by the sliding mass at its maximum velocity directly correlated with an elevated wave amplitude.Based on the theories of low-amplitude waves and energy conservation,while considering the energy conversion efficiency,a predictive model for the initial wave amplitude was derived.This model could fit and validate the functions of wavelength and wave velocity.The accuracy of the initial wave amplitude was verified using physical experiment data,with a prediction accuracy for the maximum initial wave amplitude reaching 90%.The conversion efficiency(η)directly determined the accuracy of the estimation formula.Under clear conditions for landslide-induced impulse wave generation,estimating the value ofηthrough analogy cases was feasible.This study has derived the landslide-induced impulse waves amplitude prediction formula from the standpoints of wave theory and energy conservation,with greater consideration given to the intrinsic characteristics in the formation process of landslide-induced impulse waves,thereby enhancing the applicability and extensibility of the formula.This can facilitate the development of empirical estimation methods for landslide-induced impulse waves toward universality.

Contrasts of bimodal tropical instability waves(TIWs)-induced wind stress perturbations in the Pacific Ocean among observations,ocean models,and coupled climate models

The coupling between wind stress perturbations and sea surface temperature(SST)perturbations induced by tropical instability waves(TIWs)in the Pacific Ocean has been revealed previously and proven crucial to both the atmosphere and ocean.However,an overlooked fact by previous studies is that the loosely defined“TIWs”actually consist of two modes,including the Yanai wave-based TIW on the equator(hereafter eTIW)and the Rossby wave-based TIW off the equator(hereafter vTIW).Hence,the individual feedbacks of the wind stress to the bimodal TIWs remain unexplored.In this study,individual coupling relationships are established for both eTIW and v TIW,including the relationship between the TIW-induced SST perturbations and two components of wind stress perturbations,and the relationship between the TIW-induced wind stress perturbation divergence(curl)and the downwind(crosswind)TIW-induced SST gradients.Results show that,due to different distributions of eTIW and vTIW,the coupling strength induced by the eTIW is stronger on the equator,and that by the vTIW is stronger off the equator.The results of any of eTIW and vTIW are higher than those of the loosely defined TIWs.We further investigated how well the coupling relationships remained in several widely recognized oceanic general circulation models and fully coupled climate models.However,the coupling relationships cannot be well represented in most numerical models.Finally,we confirmed that higher resolution usually corresponds to more accurate simulation.Therefore,the coupling models established in this study are complementary to previous research and can be used to refine the oceanic and coupled climate models.

Comparisons of Wave Force Model Effects on the Structural Responses and Fatigue Loads of a Semi-Submersible Floating Wind Turbine

The selection of wave force models will significantly impact the structural responses of floating wind turbines.In this study,comparisons of wave force model effects on the structural responses and fatigue loads of a semi-submersible floating wind turbine(SFWT)were conducted.Simulations were performed by employing the Morison equation(ME)with linear or second-order wave kinematics and potential flow theory(PFT)with first-or second-order wave forces.A comparison of regular waves,irregular waves,and coupled wind/waves analyses with the experimental data showed that many of the simulation results and experimental data are relatively consistent.However,notable discrepancies are found in the response amplitude operators for platform heave,tower base bending moment,and tension in mooring lines.PFT models give more satisfactory results of heave but more significant discrepan-cies in tower base bending moment than the ME models.In irregular wave analyses,low-frequency resonances were captured by PFT models with second-order difference-frequency terms,and high-frequency resonances were captured by the ME models or PFT models with second-order sum-frequency terms.These force models capture the response frequencies but do not reasonably predict the response amplitudes.The coupled wind/waves analyses showed more satisfactory results than the wave-only analyses.However,an important detail to note is that this satisfactory result is based on the overprediction of wind-induced responses.

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.

RANS Modeling of Solitary Wave Propagation over A Submerged Rectangular Breakwater

In this paper a series of numerical simulations are performed to investigate the vortex shedding mechanism for a solitary wave propagating over a submerged breakwater by use of Reynolds averaged Navier-Stokes （RAINS） model combined with a k-ε model. Flows of different Reynolds numbers up to Re = 1.4 × 10^5 corresponding to varying incident wave heights are considered in which the characteristic fluid velocity is represented by the maximum horizontal velocity above the submerged breakwater. For the verification of the accuracy of the numerical model, the incident waves and the velocity field in the vicinity of the breakwater are compared with experimental data. The result shows that the model is capable of describing vortex shedding for a solitary wave propagating over a rectangular submerged breakwater. Key features of vortex generation, evolution and dissipation are investigated. It is found that the vortex shedding and their evolution due to separated boundary layer over the breakwater are strongly related to the Reynolds number. A considerable number of vortices and complicated vortex pattern are observed as the Reynolds number increases.

Numerical Modeling of Water Wave Interaction with A Soft Mud Bed

A vertical 2-D numerical model is presented for simulating the interaction between water waves and a soft mud bed. Taking into account nonlinear rheology, a semi-empirical rheological model is applied to this water-mud model, reflecting the combined visco-elasto-plastic properties of soft mud under such oscillatory external forces as water waves. In order to increase the resolution of the flow in the neighborhood of both sides of the inter-surface, a logarithmic grid in the vertical direction is employed for numerical treatment. Model verifications are given through comparisons between the calculated and the measured mud mass transport velocities as well as wave height changes.

A stable staggered-grid finite-difference scheme for acoustic modeling beyond conventional stability limit

Staggered-grid finite-difference(SGFD)schemes have been widely used in acoustic wave modeling for geophysical problems.Many improved methods are proposed to enhance the accuracy of numerical modeling.However,these methods are inevitably limited by the maximum Courant-Friedrichs-Lewy(CFL)numbers,making them unstable when modeling with large time sampling intervals or small grid spacings.To solve this problem,we extend a stable SGFD scheme by controlling SGFD dispersion relations and maximizing the maximum CFL numbers.First,to improve modeling stability,we minimize the error between the FD dispersion relation and the exact relation in the given wave-number region,and make the FD dispersion approach a given function outside the given wave-number area,thus breaking the conventional limits of the maximum CFL number.Second,to obtain high modeling accuracy,we use the SGFD scheme based on the Remez algorithm to compute the FD coefficients.In addition,the hybrid absorbing boundary condition is adopted to suppress boundary reflections and we find a suitable weighting coefficient for the proposed scheme.Theoretical derivation and numerical modeling demonstrate that the proposed scheme can maintain high accuracy in the modeling process and the value of the maximum CFL number of the proposed scheme can exceed that of the conventional SGFD scheme when adopting a small maximum effective wavenumber,indicating that the proposed scheme improves stability during the modeling.

Simulation of wave scattering over a floating platform in the ocean with a coupled CFD-IBM model

A numerical study of linear wave scattering over a floating platform has been simulated by an efficient numericalmodel in this letter.The non-hydrostatic model is used to simulate the free surface and the uneven bottom.For thesolid body modelling,the immersed boundary method(IBM)is implemented by introducing a virtual boundaryforce into the momentum equations to emulate the boundary conditions.This implementation enhances theability of the model to simulate interactions between waves and floating structures.A numerical case involvingwave interactions with a floating platform is studied to validate the numerical model.By simulating the wavepropagation,the numerical model captures the variation of the wave scattering very well,which verifies theperformance of the numerical model and the robust strategy of the IBM.

Predicting hepatocellular carcinoma: A new non-invasive model based on shear wave elastography

BACKGROUND Integrating conventional ultrasound features with 2D shear wave elastography(2D-SWE)can potentially enhance preoperative hepatocellular carcinoma(HCC)predictions.AIM To develop a 2D-SWE-based predictive model for preoperative identification of HCC.METHODS A retrospective analysis of 884 patients who underwent liver resection and pathology evaluation from February 2021 to August 2023 was conducted at the Oriental Hepatobiliary Surgery Hospital.The patients were divided into the modeling group(n=720)and the control group(n=164).The study included conventional ultrasound,2D-SWE,and preoperative laboratory tests.Multiple logistic regression was used to identify independent predictive factors for RESULTS In the modeling group analysis,maximal elasticity(Emax)of tumors and their peripheries,platelet count,cirrhosis,and blood flow were independent risk indicators for malignancies.These factors yielded an area under the curve of 0.77(95%confidence interval:0.73-0.81)with 84%sensitivity and 61%specificity.The model demonstrated good calibration in both the construction and validation cohorts,as shown by the calibration graph and Hosmer-Lemeshow test(P=0.683 and P=0.658,respectively).Additionally,the mean elasticity(Emean)of the tumor periphery was identified as a risk factor for microvascular invasion(MVI)in malignant liver tumors(P=0.003).Patients receiving antiviral treatment differed significantly in platelet count(P=0.002),Emax of tumors(P=0.033),Emean of tumors(P=0.042),Emax at tumor periphery(P<0.001),and Emean at tumor periphery(P=0.003).CONCLUSION 2D-SWE’s hardness value serves as a valuable marker for enhancing the preoperative diagnosis of malignant liver lesions,correlating significantly with MVI and antiviral treatment efficacy.