Seismic inversion with generalized Radon transform based on local second-order approximation of scattered field in acoustic media

Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms （Born/ single scattering approximation） are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation （heterogeneous） medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps： firstly, to handle the forward scattering by taking into account the second- order Born approximation, which is related to generalized Radon transform （GRT） about quadratic scattering poten- tial; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation （3c/c0） of background media up to 10 %, andits inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the per- turbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a trans- mission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process.

Evolutionary Safe Padé Approximation Scheme for Dynamical Study of Nonlinear Cervical Human Papilloma Virus Infection Model

This study proposes a structure-preserving evolutionary framework to find a semi-analytical approximate solution for a nonlinear cervical cancer epidemic(CCE)model.The underlying CCE model lacks a closed-form exact solution.Numerical solutions obtained through traditional finite difference schemes do not ensure the preservation of the model’s necessary properties,such as positivity,boundedness,and feasibility.Therefore,the development of structure-preserving semi-analytical approaches is always necessary.This research introduces an intelligently supervised computational paradigm to solve the underlying CCE model’s physical properties by formulating an equivalent unconstrained optimization problem.Singularity-free safe Padérational functions approximate the mathematical shape of state variables,while the model’s physical requirements are treated as problem constraints.The primary model of the governing differential equations is imposed to minimize the error between approximate solutions.An evolutionary algorithm,the Genetic Algorithm with Multi-Parent Crossover(GA-MPC),executes the optimization task.The resulting method is the Evolutionary Safe PadéApproximation(ESPA)scheme.The proof of unconditional convergence of the ESPA scheme on the CCE model is supported by numerical simulations.The performance of the ESPA scheme on the CCE model is thoroughly investigated by considering various orders of non-singular Padéapproximants.

Determining Hubbard U of VO_(2) by the quasi-harmonic approximation

Vanadium dioxide VO_(2) is a strongly correlated material that undergoes a metal-to-insulator transition around 340 K.In order to describe the electron correlation effects in VO_(2), the DFT+U method is commonly employed in calculations.However, the choice of the Hubbard U parameter has been a subject of debate and its value has been reported over a wide range. In this paper, taking focus on the phase transition behavior of VO_(2), the Hubbard U parameter for vanadium oxide is determined by using the quasi-harmonic approximation(QHA). First-principles calculations demonstrate that the phase transition temperature can be modulated by varying the U values. The phase transition temperature can be well reproduced by the calculations using the Perdew–Burke–Ernzerhof functional combined with the U parameter of 1.5eV. Additionally,the calculated band structure, insulating or metallic properties, and phonon dispersion with this U value are in line with experimental observations. By employing the QHA to determine the Hubbard U parameter, this study provides valuable insights into the phase transition behavior of VO_(2). The findings highlight the importance of electron correlation effects in accurately describing the properties of this material. The agreement between the calculated results and experimental observations further validates the chosen U value and supports the use of the DFT+U method in studying VO_(2).

A multi-scale second-order autoregressive recursive filter approach for the sea ice concentration analysis

To effectively extract multi-scale information from observation data and improve computational efficiency,a multi-scale second-order autoregressive recursive filter(MSRF)method is designed.The second-order autoregressive filter used in this study has been attempted to replace the traditional first-order recursive filter used in spatial multi-scale recursive filter(SMRF)method.The experimental results indicate that the MSRF scheme successfully extracts various scale information resolved by observations.Moreover,compared with the SMRF scheme,the MSRF scheme improves computational accuracy and efficiency to some extent.The MSRF scheme can not only propagate to a longer distance without the attenuation of innovation,but also reduce the mean absolute deviation between the reconstructed sea ice concentration results and observations reduced by about 3.2%compared to the SMRF scheme.On the other hand,compared with traditional first-order recursive filters using in the SMRF scheme that multiple filters are executed,the MSRF scheme only needs to perform two filter processes in one iteration,greatly improving filtering efficiency.In the two-dimensional experiment of sea ice concentration,the calculation time of the MSRF scheme is only 1/7 of that of SMRF scheme.This means that the MSRF scheme can achieve better performance with less computational cost,which is of great significance for further application in real-time ocean or sea ice data assimilation systems in the future.

Saddlepoint Approximation Method in Reliability Analysis:A Review

The escalating need for reliability analysis(RA)and reliability-based design optimization(RBDO)within engineering challenges has prompted the advancement of saddlepoint approximationmethods(SAM)tailored for such problems.This article offers a detailed overview of the general SAM and summarizes the method characteristics first.Subsequently,recent enhancements in the SAM theoretical framework are assessed.Notably,the mean value first-order saddlepoint approximation(MVFOSA)bears resemblance to the conceptual framework of the mean value second-order saddlepoint approximation(MVSOSA);the latter serves as an auxiliary approach to the former.Their distinction is rooted in the varying expansion orders of the performance function as implemented through the Taylor method.Both the saddlepoint approximation and third-moment(SATM)and saddlepoint approximation and fourth-moment(SAFM)strategies model the cumulant generating function(CGF)by leveraging the initial random moments of the function.Although their optimal application domains diverge,each method consistently ensures superior relative precision,enhanced efficiency,and sustained stability.Every method elucidated is exemplified through pertinent RA or RBDO scenarios.By juxtaposing them against alternative strategies,the efficacy of these methods becomes evident.The outcomes proffered are subsequently employed as a foundation for contemplating prospective theoretical and practical research endeavors concerning SAMs.The main purpose and value of this article is to review the SAM and reliability-related issues,which can provide some reference and inspiration for future research scholars in this field.

Gamma Approximation Based Multi-Antenna Covert Communication Detection

Covert communication technology makes wireless communication more secure,but it also provides more opportunities for illegal users to transmit harmful information.In order to detect the illegal covert communication of the lawbreakers in real time for subsequent processing,this paper proposes a Gamma approximation-based detection method for multi-antenna covert communication systems.Specifically,the Gamma approximation property is used to calculate the miss detection rate and false alarm rate of the monitor firstly.Then the optimization problem to minimize the sum of the missed detection rate and the false alarm rate is proposed.The optimal detection threshold and the minimum error detection probability are solved according to the properties of the Lambert W function.Finally,simulation results are given to demonstrate the effectiveness of the proposed method.

Relaxed Stability Criteria for Time-Delay Systems:A Novel Quadratic Function Convex Approximation Approach

This paper develops a quadratic function convex approximation approach to deal with the negative definite problem of the quadratic function induced by stability analysis of linear systems with time-varying delays.By introducing two adjustable parameters and two free variables,a novel convex function greater than or equal to the quadratic function is constructed,regardless of the sign of the coefficient in the quadratic term.The developed lemma can also be degenerated into the existing quadratic function negative-determination(QFND)lemma and relaxed QFND lemma respectively,by setting two adjustable parameters and two free variables as some particular values.Moreover,for a linear system with time-varying delays,a relaxed stability criterion is established via our developed lemma,together with the quivalent reciprocal combination technique and the Bessel-Legendre inequality.As a result,the conservatism can be reduced via the proposed approach in the context of constructing Lyapunov-Krasovskii functionals for the stability analysis of linear time-varying delay systems.Finally,the superiority of our results is illustrated through three numerical examples.

An Overview of Sequential Approximation in Topology Optimization of Continuum Structure

This paper offers an extensive overview of the utilization of sequential approximate optimization approaches in the context of numerically simulated large-scale continuum structures.These structures,commonly encountered in engineering applications,often involve complex objective and constraint functions that cannot be readily expressed as explicit functions of the design variables.As a result,sequential approximation techniques have emerged as the preferred strategy for addressing a wide array of topology optimization challenges.Over the past several decades,topology optimization methods have been advanced remarkably and successfully applied to solve engineering problems incorporating diverse physical backgrounds.In comparison to the large-scale equation solution,sensitivity analysis,graphics post-processing,etc.,the progress of the sequential approximation functions and their corresponding optimizersmake sluggish progress.Researchers,particularly novices,pay special attention to their difficulties with a particular problem.Thus,this paper provides an overview of sequential approximation functions,related literature on topology optimization methods,and their applications.Starting from optimality criteria and sequential linear programming,the other sequential approximate optimizations are introduced by employing Taylor expansion and intervening variables.In addition,recent advancements have led to the emergence of approaches such as Augmented Lagrange,sequential approximate integer,and non-gradient approximation are also introduced.By highlighting real-world applications and case studies,the paper not only demonstrates the practical relevance of these methods but also underscores the need for continued exploration in this area.Furthermore,to provide a comprehensive overview,this paper offers several novel developments that aim to illuminate potential directions for future research.

Nonlinear Filtering With Sample-Based Approximation Under Constrained Communication:Progress, Insights and Trends

The nonlinear filtering problem has enduringly been an active research topic in both academia and industry due to its ever-growing theoretical importance and practical significance.The main objective of nonlinear filtering is to infer the states of a nonlinear dynamical system of interest based on the available noisy measurements. In recent years, the advance of network communication technology has not only popularized the networked systems with apparent advantages in terms of installation,cost and maintenance, but also brought about a series of challenges to the design of nonlinear filtering algorithms, among which the communication constraint has been recognized as a dominating concern. In this context, a great number of investigations have been launched towards the networked nonlinear filtering problem with communication constraints, and many samplebased nonlinear filters have been developed to deal with the highly nonlinear and/or non-Gaussian scenarios. The aim of this paper is to provide a timely survey about the recent advances on the sample-based networked nonlinear filtering problem from the perspective of communication constraints. More specifically, we first review three important families of sample-based filtering methods known as the unscented Kalman filter, particle filter,and maximum correntropy filter. Then, the latest developments are surveyed with stress on the topics regarding incomplete/imperfect information, limited resources and cyber security.Finally, several challenges and open problems are highlighted to shed some lights on the possible trends of future research in this realm.

The Convergence Rate of Fréchet Distribution under the Second-Order Regular Variation Condition

In this article we consider the asymptotic behavior of extreme distribution with the extreme value index γ>0 . The rates of uniform convergence for Fréchet distribution are constructed under the second-order regular variation condition.

RBMDO Using Gaussian Mixture Model-Based Second-Order Mean-Value Saddlepoint Approximation

Actual engineering systems will be inevitably affected by uncertain factors.Thus,the Reliability-Based Multidisciplinary Design Optimization(RBMDO)has become a hotspot for recent research and application in complex engineering system design.The Second-Order/First-Order Mean-Value Saddlepoint Approximate(SOMVSA/-FOMVSA)are two popular reliability analysis strategies that are widely used in RBMDO.However,the SOMVSA method can only be used efficiently when the distribution of input variables is Gaussian distribution,which significantly limits its application.In this study,the Gaussian Mixture Model-based Second-Order Mean-Value Saddlepoint Approximation(GMM-SOMVSA)is introduced to tackle above problem.It is integrated with the Collaborative Optimization(CO)method to solve RBMDO problems.Furthermore,the formula and procedure of RBMDO using GMM-SOMVSA-Based CO(GMM-SOMVSA-CO)are proposed.Finally,an engineering example is given to show the application of the GMM-SOMVSA-CO method.

Study of(e,2e) process on potassium at 6 e V-60 eV above threshold in a second-order Born approximation

The standard distorted wave Born approximation （DWBA） method has been extended to second-order Born amplitude in order to describe the multiple interactions between the projectile and the atomic target. Second-order DWBA calculations have been preformed to investigate the triple differential cross sections （TDCS） of coplanar doubly symmetric （e, 2e） collisions for the alkali target potassium at excess energies of 6 eV-60 eV. Compared with the previous first-order DWBA calculations, the present theoretical model improves the degree of agreement with experiments, especially for the backward scattering angle region of TDCS. This indicates that the present second-order Born term is capable of giving a reasonable correction to the DWBA model in studying coplanar symmetric （e, 2e） problems in low and intermediate energy ranges.

Study of scattering from time-varying Gerstners sea surface using second-order small slope approximation

Backscattered fields from one-dimensional time-varying Gerstners sea surface are calculated utilising the secondorder small slope approximation. It is well known that spectral properties of the backscattered echoes relate to the velocity of the small elementary scatterers on sea surface profiles. Therefore, modeling Doppler spectra from the ocean requires an accurate description of the sea surface motion. The profile of nonlinear Gerstners sea surface shows verticalskewness of sea waves, it is sharper at the crest and flatter at the trough than linear waves, and its maximum slope position is closer to the crest than to the trough. Furthermore, the horizontal component of the small elementary scatterers orbit velocity on the sea surface, which yields noticeable influence on Doppler spectra, can be obtained conveniently by Gerstners sea surface model. In this study the characteristics of Doppler spectra of backscattered fields from time-varying Gerstners sea surface are investigated and the dependences of the Doppler frequency and the Doppler bandwidth on the parameters, such as the wind speed, the radar frequency, the incident angle, etc. are discussed. It is shown that the Doppler bandwidth of microwave scattered fields from Gerstners sea surface is considerably broadened. For the case of high frequency backscattered fields, the values of the higher-order spectrum peaks are larger than those obtained by linear sea surface.

Dynamics of Rabi model under second-order Born Oppenheimer approximation

We apply the second-order Born-Oppenheimer （BO） approximation to investigate the dynamics of the Rabi model, which describes the interaction between a two-level system and a single bosonic mode beyond the rotating wave approxi- mation. By comparing with the numerical results, we find that our approach works well when the frequency of the two-level system is much smaller than that of the bosonic mode.

Doppler spectral analysis for time-evolving sea surfaces using second-order small slope approximation

The second-order small slope approximation （SSA2） method is introduced to study the Doppler characteristics from time-evolving sea surfaces. Simulation results show better agreement between the SSA2 model and the numerical method for both vertical and horizontal polarizations, meaning that SSA2 gives a satisfactory prediction of the spectral difference between two po- larizations; while such discrepancy cannot be captured using the lowest-order SSA （SSA1） model. In particular, the Doppler shifts and spectral widths for different incident angles, wind directions and polarizations are analyzed, demonstrating correct variations with respect to such parameters. Those observations prove that the SSA2 provides an efficient and relatively fast tool for sea surface Doppler spectral analysis.

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems

The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.

Integrating a novel irrigation approximation method with a process-based remote sensing model to estimate multi-years'winter wheat yield over the North China Plain

Accurate estimation of regional winter wheat yields is essential for understanding the food production status and ensuring national food security.However,using the existing remote sensing-based crop yield models to accurately reproduce the inter-annual and spatial variations in winter wheat yields remains challenging due to the limited ability to acquire irrigation information in water-limited regions.Thus,we proposed a new approach to approximating irrigations of winter wheat over the North China Plain(NCP),where irrigation occurs extensively during the winter wheat growing season.This approach used irrigation pattern parameters(IPPs)to define the irrigation frequency and timing.Then,they were incorporated into a newly-developed process-based and remote sensing-driven crop yield model for winter wheat(PRYM–Wheat),to improve the regional estimates of winter wheat over the NCP.The IPPs were determined using statistical yield data of reference years(2010–2015)over the NCP.Our findings showed that PRYM–Wheat with the optimal IPPs could improve the regional estimate of winter wheat yield,with an increase and decrease in the correlation coefficient(R)and root mean square error(RMSE)of 0.15(about 37%)and 0.90 t ha–1(about 41%),respectively.The data in validation years(2001–2009 and 2016–2019)were used to validate PRYM–Wheat.In addition,our findings also showed R(RMSE)of 0.80(0.62 t ha–1)on a site level,0.61(0.91 t ha–1)for Hebei Province on a county level,0.73(0.97 t ha–1)for Henan Province on a county level,and 0.55(0.75 t ha–1)for Shandong Province on a city level.Overall,PRYM–Wheat can offer a stable and robust approach to estimating regional winter wheat yield across multiple years,providing a scientific basis for ensuring regional food security.

Approximate Controllability of Second-Order Neutral Stochastic Differential Equations with Infinite Delay and Poisson Jumps

The modelling of risky asset by stochastic processes with continuous paths, based on Brow- nian motions, suffers from several defects. First, the path continuity assumption does not seem reason- able in view of the possibility of sudden price variations （jumps） resulting of market crashes. A solution is to use stochastic processes with jumps, that will account for sudden variations of the asset prices. On the other hand, such jump models are generally based on the Poisson random measure. Many popular economic and financial models described by stochastic differential equations with Poisson jumps. This paper deals with the approximate controllability of a class of second-order neutral stochastic differential equations with infinite delay and Poisson jumps. By using the cosine family of operators, stochastic analysis techniques, a new set of sufficient conditions are derived for the approximate controllability of the above control system. An example is provided to illustrate the obtained theory.

EXISTENCE OF SOLUTION AND APPROXIMATE CONTROLLABILITY OF A SECOND-ORDER NEUTRAL STOCHASTIC DIFFERENTIAL EQUATION WITH STATE DEPENDENT DELAY

This paper has two sections which deals with a second order stochastic neutral partial differential equation with state dependent delay. In the first section the existence and uniqueness of mild solution is obtained by use of measure of non-compactness. In the second section the conditions for approximate controllability are investigated for the distributed second order neutral stochastic differential system with respect to the approximate controllability of the corresponding linear system in a Hilbert space. Our method is an extension of co-author N. Sukavanam’s novel approach in [22]. Thereby, we remove the need to assume the invertibility of a controllability operator used by authors in [5], which fails to exist in infinite dimensional spaces if the associated semigroup is compact. Our approach also removes the need to check the invertibility of the controllability Gramian operator and associated limit condition used by the authors in [20], which are practically difficult to verify and apply. An example is provided to illustrate the presented theory.

Closed-form steady-state solutions for forced vibration of second-order axially moving systems

Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for the forced vibration when the damping effect and the coupling effect of multiple second-order models are considered.In this paper, Green's function method based on the Laplace transform is used to obtain closed-form solutions for the forced vibration of second-order axially moving systems. By taking the axially moving damping string system and multi-string system connected by springs as examples, the detailed solution methods and the analytical Green's functions of these second-order systems are given. The mode functions and frequency equations are also obtained by the obtained Green's functions. The reliability and convenience of the results are verified by several examples. This paper provides a systematic analytical method for the dynamic analysis of second-order axially moving systems, and the obtained Green's functions are applicable to different second-order systems rather than just string systems. In addition, the work of this paper also has positive significance for the study on the forced vibration of high-order systems.