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...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.展开更多
A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an o...A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an optimization problem with several parameters determined by applying a generic algorithm. The optimized schemes are analyzed carefully from the aspects of the eigenvalue distribution, the ε-pseudospectra, the short time behavior, and the Fourier analysis. Numerical experiments for the Euler equations are used to show the effectiveness of the final recommended scheme.展开更多
This paper proposed several new types of finite-difference methods for the shallow water equation in absolute coordinate system and put forward an effective two-step predictor-corrector method, a compact and iterative...This paper proposed several new types of finite-difference methods for the shallow water equation in absolute coordinate system and put forward an effective two-step predictor-corrector method, a compact and iterative algorithm for five diagonal matrix. Then the iterative method was used for a multi-grid procedure for shallow water equation. A t last, an initial-boundary value problem was considered, and the numerical results show that the linear sinusoidal wave would successively evolve into conoidal wave.展开更多
Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to ...Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.展开更多
The best finite-difference scheme for the Helmholtz equation is suggested. A method of solving obtained finite-difference scheme is developed. The efficiency and accuracy of method were tested on several examples.
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 met...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.展开更多
In this paper, we present splitting schemes for the two-level Bloch model. After proposing two ways to split the Bloch equation, we show that it is possible in each case togenerate exact numerical solutions of the obt...In this paper, we present splitting schemes for the two-level Bloch model. After proposing two ways to split the Bloch equation, we show that it is possible in each case togenerate exact numerical solutions of the obtained sub-equations. These exact solutionsinvolve matrix exponentials which can be expensive to compute. Here, for 2×2 matriceswe develop equivalent formulations which reduce the computational cost. These splittingschemes are nonstandard ones and conserve all the physical properties (Hermicity, positiveness and trace) of Bloch equations. In addition, they are explicit, making effectivetheir implementation when coupled with the Maxwell’s equations.展开更多
The approach to optimization of finite-difference(FD)schemes for the linear advection equation(LAE)is proposed.The FD schemes dependent on the scalar dimensionless parameter are considered.The parameter is included in...The approach to optimization of finite-difference(FD)schemes for the linear advection equation(LAE)is proposed.The FD schemes dependent on the scalar dimensionless parameter are considered.The parameter is included in the expression,which approximates the term with spatial derivatives.The approach is based on the considering of the dispersive and dissipative characteristics of the schemes as the functions of the parameter.For the proper choice of the parameter,these functions are minimized.The approach is applied to the optimization of two-step schemes with an asymmetric approximation of time derivative and with various approximations of the spatial term.The cases of schemes from first to fourth approximation orders are considered.The optimal values of the parameter are obtained.Schemes with the optimal values are applied to the solution of test problems with smooth and discontinuous initial conditions.Also,schemes are used in the FD-based lattice Boltzmann method(LBM)for modeling of the compressible gas flow.The obtained numerical results demonstrate the convergence of the schemes and decaying of the numerical dispersion.展开更多
The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order sy...The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.展开更多
An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
An explicit multi-conservation finite-difference scheme for solving the spherical shallow-water-wave equation set of barotropic atmosphere has been proposed. The numerical scheme is based on a special semi-discrete fo...An explicit multi-conservation finite-difference scheme for solving the spherical shallow-water-wave equation set of barotropic atmosphere has been proposed. The numerical scheme is based on a special semi-discrete form of the equations that conserves four basic physical integrals including the total energy, total mass, total potential vorticity and total enstrophy. Numerical tests show that the new scheme performs closely like but is much more time-saving than the implicit multi-conservation scheme.展开更多
The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearitie...The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearities.To discretize the continuous problem in space,we utilize a spectral Galerkin method.Subsequently,we introduce a nonlinear-tamed exponential integrator scheme,resulting in a fully discrete scheme.Within the framework of semigroup theory,this study provides precise estimations of the Sobolev regularity,L^(∞) regularity in space,and Hölder continuity in time for the mild solution,as well as for its semi-discrete and full-discrete approximations.Building upon these results,we establish moment boundedness for the numerical solution and obtain strong convergence rates in both spatial and temporal dimensions.A numerical example is presented to validate the theoretical findings.展开更多
As blockchain technology rapidly evolves,smart contracts have seen widespread adoption in financial transactions and beyond.However,the growing prevalence of malicious Ponzi scheme contracts presents serious security ...As blockchain technology rapidly evolves,smart contracts have seen widespread adoption in financial transactions and beyond.However,the growing prevalence of malicious Ponzi scheme contracts presents serious security threats to blockchain ecosystems.Although numerous detection techniques have been proposed,existing methods suffer from significant limitations,such as class imbalance and insufficient modeling of transaction-related semantic features.To address these challenges,this paper proposes an oversampling-based detection framework for Ponzi smart contracts.We enhance the Adaptive Synthetic Sampling(ADASYN)algorithm by incorporating sample proximity to decision boundaries and ensuring realistic sample distributions.This enhancement facilitates the generation of high-quality minority class samples and effectively mitigates class imbalance.In addition,we design a Contract Transaction Graph(CTG)construction algorithm to preserve key transactional semantics through feature extraction from contract code.A graph neural network(GNN)is then applied for classification.This study employs a publicly available dataset from the XBlock platform,consisting of 318 verified Ponzi contracts and 6498 benign contracts.Sourced from real Ethereum deployments,the dataset reflects diverse application scenarios and captures the varied characteristics of Ponzi schemes.Experimental results demonstrate that our approach achieves an accuracy of 96%,a recall of 92%,and an F1-score of 94%in detecting Ponzi contracts,outperforming state-of-the-art methods.展开更多
Clouds play an important role in global atmospheric energy and water vapor budgets, and the low cloud simulations suffer from large biases in many atmospheric general circulation models. In this study, cloud microphys...Clouds play an important role in global atmospheric energy and water vapor budgets, and the low cloud simulations suffer from large biases in many atmospheric general circulation models. In this study, cloud microphysical processes such as raindrop evaporation and cloud water accretion in a double-moment six-class cloud microphysics scheme were revised to enhance the simulation of low clouds using the Global-Regional Integrated Forecast System(GRIST)model. The validation of the revised scheme using a single-column version of the GRIST demonstrated a reasonable reduction in liquid water biases. The revised parameterization simulated medium-and low-level cloud fractions that were in better agreement with the observations than the original scheme. Long-term global simulations indicate the mitigation of the originally overestimated low-level cloud fraction and cloud-water mixing ratio in mid-to high-latitude regions,primarily owing to enhanced accretion processes and weakened raindrop evaporation. The reduced low clouds with the revised scheme showed better consistency with satellite observations, particularly at mid-and high-latitudes. Further improvements can be observed in the simulated cloud shortwave radiative forcing and vertical distribution of total cloud cover. Annual precipitation in mid-latitude regions has also improved, particularly over the oceans, with significantly increased large-scale and decreased convective precipitation.展开更多
In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is de...In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.展开更多
This study proposes a class of augmented subspace schemes for the weak Galerkin(WG)finite element method used to solve eigenvalue problems.The augmented subspace is built with the conforming linear finite element spac...This study proposes a class of augmented subspace schemes for the weak Galerkin(WG)finite element method used to solve eigenvalue problems.The augmented subspace is built with the conforming linear finite element space defined on the coarse mesh and the eigen-function approximations in the WG finite element space defined on the fine mesh.Based on this augmented subspace,solving the eigenvalue problem in the fine WG finite element space can be reduced to the solution of the linear boundary value problem in the same WG finite element space and a low dimensional eigenvalue problem in the augmented sub-space.The proposed augmented subspace techniques have the second order convergence rate with respect to the coarse mesh size,as demonstrated by the accompanying error esti-mates.Finally,a few numerical examples are provided to validate the proposed numerical techniques.展开更多
In this paper,we propose and analyze two second-order accurate finite difference schemes for the one-dimensional heat equation with concentrated capacity on a computa-tional domain=[a,b].We first transform the target ...In this paper,we propose and analyze two second-order accurate finite difference schemes for the one-dimensional heat equation with concentrated capacity on a computa-tional domain=[a,b].We first transform the target equation into the standard heat equation on the domain excluding the singular point equipped with an inner interface matching(IIM)condition on the singular point x=ξ∈(a,b),then adopt Taylor’s ex-pansion to approximate the IIM condition at the singular point and apply second-order finite difference method to approximate the standard heat equation at the nonsingular points.This discrete procedure allows us to choose different grid sizes to partition the two sub-domains[a,ξ]and[ξ,b],which ensures that x=ξ is a grid point,and hence the pro-posed schemes can be generalized to the heat equation with more than one concentrated capacities.We prove that the two proposed schemes are uniquely solvable.And through in-depth analysis of the local truncation errors,we rigorously prove that the two schemes are second-order accurate both in temporal and spatial directions in the maximum norm without any constraint on the grid ratio.Numerical experiments are carried out to verify our theoretical conclusions.展开更多
This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called...This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.展开更多
Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of t...Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.展开更多
基金supported by the National Natural Science Foundation of China(No.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘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.
基金Project supported by the National Natural Science Foundation of China(Nos.11601517,11502296,61772542,and 61561146395)the Basic Research Foundation of National University of Defense Technology(No.ZDYYJ-CYJ20140101)
文摘A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an optimization problem with several parameters determined by applying a generic algorithm. The optimized schemes are analyzed carefully from the aspects of the eigenvalue distribution, the ε-pseudospectra, the short time behavior, and the Fourier analysis. Numerical experiments for the Euler equations are used to show the effectiveness of the final recommended scheme.
文摘This paper proposed several new types of finite-difference methods for the shallow water equation in absolute coordinate system and put forward an effective two-step predictor-corrector method, a compact and iterative algorithm for five diagonal matrix. Then the iterative method was used for a multi-grid procedure for shallow water equation. A t last, an initial-boundary value problem was considered, and the numerical results show that the linear sinusoidal wave would successively evolve into conoidal wave.
文摘Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.
文摘The best finite-difference scheme for the Helmholtz equation is suggested. A method of solving obtained finite-difference scheme is developed. The efficiency and accuracy of method were tested on several examples.
基金This research is supported by the National Natural Science Foundation of China(NSFC)under contract no.42274147.
文摘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.
文摘In this paper, we present splitting schemes for the two-level Bloch model. After proposing two ways to split the Bloch equation, we show that it is possible in each case togenerate exact numerical solutions of the obtained sub-equations. These exact solutionsinvolve matrix exponentials which can be expensive to compute. Here, for 2×2 matriceswe develop equivalent formulations which reduce the computational cost. These splittingschemes are nonstandard ones and conserve all the physical properties (Hermicity, positiveness and trace) of Bloch equations. In addition, they are explicit, making effectivetheir implementation when coupled with the Maxwell’s equations.
文摘The approach to optimization of finite-difference(FD)schemes for the linear advection equation(LAE)is proposed.The FD schemes dependent on the scalar dimensionless parameter are considered.The parameter is included in the expression,which approximates the term with spatial derivatives.The approach is based on the considering of the dispersive and dissipative characteristics of the schemes as the functions of the parameter.For the proper choice of the parameter,these functions are minimized.The approach is applied to the optimization of two-step schemes with an asymmetric approximation of time derivative and with various approximations of the spatial term.The cases of schemes from first to fourth approximation orders are considered.The optimal values of the parameter are obtained.Schemes with the optimal values are applied to the solution of test problems with smooth and discontinuous initial conditions.Also,schemes are used in the FD-based lattice Boltzmann method(LBM)for modeling of the compressible gas flow.The obtained numerical results demonstrate the convergence of the schemes and decaying of the numerical dispersion.
基金supported by the National Natural Science Foundation of China(Grant Nos.60931002 and 61101064)the Universities Natural Science Foundation of Anhui Province,China(Grant Nos.KJ2011A002 and 1108085J01)
文摘The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
基金the National Key Development and Planning Project for the Basic Research (973) (Grant No.2005CB321703)the Science Funds for Creative Research Groups (Grant No.40221503)
文摘An explicit multi-conservation finite-difference scheme for solving the spherical shallow-water-wave equation set of barotropic atmosphere has been proposed. The numerical scheme is based on a special semi-discrete form of the equations that conserves four basic physical integrals including the total energy, total mass, total potential vorticity and total enstrophy. Numerical tests show that the new scheme performs closely like but is much more time-saving than the implicit multi-conservation scheme.
基金partially supported by the National Natural Science Foundation of China(Grant No.12071073)financial support by the Jiangsu Provincial Scientific Research Center of Applied Mathematics(Grant No.BK20233002).
文摘The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise,which possesses both Burgers-type and cubic nonlinearities.To discretize the continuous problem in space,we utilize a spectral Galerkin method.Subsequently,we introduce a nonlinear-tamed exponential integrator scheme,resulting in a fully discrete scheme.Within the framework of semigroup theory,this study provides precise estimations of the Sobolev regularity,L^(∞) regularity in space,and Hölder continuity in time for the mild solution,as well as for its semi-discrete and full-discrete approximations.Building upon these results,we establish moment boundedness for the numerical solution and obtain strong convergence rates in both spatial and temporal dimensions.A numerical example is presented to validate the theoretical findings.
基金supported by the Key Project of Joint Fund of the National Natural Science Foundation of China“Research on Key Technologies and Demonstration Applications for Trusted and Secure Data Circulation and Trading”(U24A20241)the National Natural Science Foundation of China“Research on Trusted Theories and Key Technologies of Data Security Trading Based on Blockchain”(62202118)+4 种基金the Major Scientific and Technological Special Project of Guizhou Province([2024]014)Scientific and Technological Research Projects from the Guizhou Education Department(Qian jiao ji[2023]003)the Hundred-Level Innovative Talent Project of the Guizhou Provincial Science and Technology Department(Qiankehe Platform Talent-GCC[2023]018)the Major Project of Guizhou Province“Research and Application of Key Technologies for Trusted Large Models Oriented to Public Big Data”(Qiankehe Major Project[2024]003)the Guizhou Province Computational Power Network Security Protection Science and Technology Innovation Talent Team(Qiankehe Talent CXTD[2025]029).
文摘As blockchain technology rapidly evolves,smart contracts have seen widespread adoption in financial transactions and beyond.However,the growing prevalence of malicious Ponzi scheme contracts presents serious security threats to blockchain ecosystems.Although numerous detection techniques have been proposed,existing methods suffer from significant limitations,such as class imbalance and insufficient modeling of transaction-related semantic features.To address these challenges,this paper proposes an oversampling-based detection framework for Ponzi smart contracts.We enhance the Adaptive Synthetic Sampling(ADASYN)algorithm by incorporating sample proximity to decision boundaries and ensuring realistic sample distributions.This enhancement facilitates the generation of high-quality minority class samples and effectively mitigates class imbalance.In addition,we design a Contract Transaction Graph(CTG)construction algorithm to preserve key transactional semantics through feature extraction from contract code.A graph neural network(GNN)is then applied for classification.This study employs a publicly available dataset from the XBlock platform,consisting of 318 verified Ponzi contracts and 6498 benign contracts.Sourced from real Ethereum deployments,the dataset reflects diverse application scenarios and captures the varied characteristics of Ponzi schemes.Experimental results demonstrate that our approach achieves an accuracy of 96%,a recall of 92%,and an F1-score of 94%in detecting Ponzi contracts,outperforming state-of-the-art methods.
基金National Natural Science Foundation of China(42375153,42105153,42205157)Development of Science and Technology at Chinese Academy of Meteorological Sciences(2023KJ038)。
文摘Clouds play an important role in global atmospheric energy and water vapor budgets, and the low cloud simulations suffer from large biases in many atmospheric general circulation models. In this study, cloud microphysical processes such as raindrop evaporation and cloud water accretion in a double-moment six-class cloud microphysics scheme were revised to enhance the simulation of low clouds using the Global-Regional Integrated Forecast System(GRIST)model. The validation of the revised scheme using a single-column version of the GRIST demonstrated a reasonable reduction in liquid water biases. The revised parameterization simulated medium-and low-level cloud fractions that were in better agreement with the observations than the original scheme. Long-term global simulations indicate the mitigation of the originally overestimated low-level cloud fraction and cloud-water mixing ratio in mid-to high-latitude regions,primarily owing to enhanced accretion processes and weakened raindrop evaporation. The reduced low clouds with the revised scheme showed better consistency with satellite observations, particularly at mid-and high-latitudes. Further improvements can be observed in the simulated cloud shortwave radiative forcing and vertical distribution of total cloud cover. Annual precipitation in mid-latitude regions has also improved, particularly over the oceans, with significantly increased large-scale and decreased convective precipitation.
基金The Natural Science Foundation of Xinjiang Uygur Autonomous Region of China“RBF-Hermite difference scheme for the time-fractional kdv-Burgers equation”(2024D01C43)。
文摘In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.
基金partly supported by the Beijing Natural Science Foundation(Grant No.Z200003)by the National Natural Science Foundation of China(Grant Nos.12331015,12301475,12301465)+1 种基金by the National Center for Mathematics and Interdisciplinary Science,Chinese Academy of Sciencesby the Research Foundation for the Beijing University of Technology New Faculty(Grant No.006000514122516).
文摘This study proposes a class of augmented subspace schemes for the weak Galerkin(WG)finite element method used to solve eigenvalue problems.The augmented subspace is built with the conforming linear finite element space defined on the coarse mesh and the eigen-function approximations in the WG finite element space defined on the fine mesh.Based on this augmented subspace,solving the eigenvalue problem in the fine WG finite element space can be reduced to the solution of the linear boundary value problem in the same WG finite element space and a low dimensional eigenvalue problem in the augmented sub-space.The proposed augmented subspace techniques have the second order convergence rate with respect to the coarse mesh size,as demonstrated by the accompanying error esti-mates.Finally,a few numerical examples are provided to validate the proposed numerical techniques.
基金supported by the National Natural Science Foundation of China(Grant No.11571181)by the Natural Science Foundation of Jiangsu Province(Grant No.BK20171454).
文摘In this paper,we propose and analyze two second-order accurate finite difference schemes for the one-dimensional heat equation with concentrated capacity on a computa-tional domain=[a,b].We first transform the target equation into the standard heat equation on the domain excluding the singular point equipped with an inner interface matching(IIM)condition on the singular point x=ξ∈(a,b),then adopt Taylor’s ex-pansion to approximate the IIM condition at the singular point and apply second-order finite difference method to approximate the standard heat equation at the nonsingular points.This discrete procedure allows us to choose different grid sizes to partition the two sub-domains[a,ξ]and[ξ,b],which ensures that x=ξ is a grid point,and hence the pro-posed schemes can be generalized to the heat equation with more than one concentrated capacities.We prove that the two proposed schemes are uniquely solvable.And through in-depth analysis of the local truncation errors,we rigorously prove that the two schemes are second-order accurate both in temporal and spatial directions in the maximum norm without any constraint on the grid ratio.Numerical experiments are carried out to verify our theoretical conclusions.
文摘This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.
基金This research was supported by the National Nature Science Foundation of China （No. 41074100） and the Program for NewCentury Excellent Talents in the University of the Ministry of Education of China （No. NCET- 10-0812）.
文摘Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.
