To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical c...To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.展开更多
A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficien...A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficients(k)are constructed based on FDM.The rock bridge area was divided through k-means algorithm and the optimal number of clusters was determined by sum of squared errors(SSE)and elbow method.The influence of maximum principal stress and stress change rate as clustering indexes on the clustering results of rock bridges was compared by using Euclidean distance.The results show that using stress change rate as clustering index is more effective.When the joint coalescence coefficient is less than 0.6,there is no significant stress concentration in the middle area of adjacent joints,that is,no generation of rock bridge.In addition,the range of rock bridge is affected by the coalescence coefficient(k),the relative position of joints and the parameters of weak interlayer.展开更多
In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving th...In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.展开更多
This study investigates the complex heat transfer dynamics inmultilayer bifacial photovoltaic(bPV)solar modules under spectrally resolved solar irradiation.A novel numericalmodel is developed to incorporate internal h...This study investigates the complex heat transfer dynamics inmultilayer bifacial photovoltaic(bPV)solar modules under spectrally resolved solar irradiation.A novel numericalmodel is developed to incorporate internal heat generation resulting from optical absorption,grounded in the physical equations governing light-matter interactions within the module’smultilayer structure.The model accounts for reflection and transmission at each interface between adjacent layers,as well as absorption within individual layers,using the wavelength-dependent dielectric properties of constituent materials.These properties are used to calculate the spectral reflectance,transmittance,and absorption coefficients,enabling precise quantification of internal heat sources from irradiance incidents on both the front and rear surfaces of the module.The study further examines the influence of irradiance reflection on thermal behavior,evaluates the thermal impact of various supporting materials placed beneath the module,and analyzes the role of albedo in modifying heat distribution.By incorporating spectrally resolved heat generation across each layer often simplified or omitted in conventional models,the proposed approach enhances physical accuracy.The transient heat equation is solved using a one-dimensional finite difference(FD)method to produce detailed temperature profiles under multiple operating scenarios,including Standard Test Conditions(STC),Bifacial Standard Test Conditions(BSTC),Normal Operating Cell Temperature(NOCT),and Bifacial NOCT(BNOCT).The results offer valuable insights into the interplay between optical and thermal phenomena in bifacial systems,informing the design and optimization of more efficient photovoltaic technologies.展开更多
Laplace–Fourier(L-F)domain finite-difference(FD)forward modeling is an important foundation for L-F domain full-waveform inversion(FWI).An optimal modeling method can improve the efficiency and accuracy of FWI.A fl e...Laplace–Fourier(L-F)domain finite-difference(FD)forward modeling is an important foundation for L-F domain full-waveform inversion(FWI).An optimal modeling method can improve the efficiency and accuracy of FWI.A fl exible FD stencil,which requires pairing and centrosymmetricity of the involved gridpoints,is used on the basis of the 2D L-F domain acoustic wave equation.The L-F domain numerical dispersion analysis is then performed by minimizing the phase error of the normalized numerical phase and attenuation propagation velocities to obtain the optimization coefficients.An optimal FD forward modeling method is finally developed for the L-F domain acoustic wave equation and applied to the traditional standard 9-point scheme and 7-and 9-point schemes,where the latter two schemes are used in discontinuous-grid FD modeling.Numerical experiments show that the optimal L-F domain FD modeling method not only has high accuracy but can also be applied to equal and unequal directional sampling intervals and discontinuous-grid FD modeling to reduce computational cost.展开更多
Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and effic...Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and efficiency of the conventional numerical modeling, I develop a new seismic modeling method by combining the FD scheme with the numerical dispersion suppression neural network(NDSNN). This method involves the following steps. First, a training data set composed of a small number of wavefield snapshots is generated. The wavefield snapshots with the low-accuracy wavefield data and the high-accuracy wavefield data are paired, and the low-accuracy wavefield snapshots involve the obvious numerical dispersion including both the temporal and spatial dispersion. Second, the NDSNN is trained until the network converges to simultaneously suppress the temporal and spatial dispersion.Third, the entire set of low-accuracy wavefield data is computed quickly using FD modeling with the large time step and the coarse grid. Fourth, the NDSNN is applied to the entire set of low-accuracy wavefield data to suppress the numerical dispersion including the temporal and spatial dispersion.Numerical modeling examples verify the effectiveness of my proposed method in improving the computational accuracy and efficiency.展开更多
In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolso...In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.展开更多
Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using t...Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.展开更多
Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of ...Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we have even been able to apply the scheme without any modification and obtain good results.Many useful PDEs may have stiff relaxation source terms for which the finite difference formulation of WENO is shown to provide some genuine advantages.展开更多
For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple boun...For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.展开更多
In this paper,we give improved error estimates for linearized and nonlinear CrankNicolson type finite difference schemes of Ginzburg-Landau equation in two dimensions.For linearized Crank-Nicolson scheme,we use mathem...In this paper,we give improved error estimates for linearized and nonlinear CrankNicolson type finite difference schemes of Ginzburg-Landau equation in two dimensions.For linearized Crank-Nicolson scheme,we use mathematical induction to get unconditional error estimates in discrete L^(2)and H^(1)norm.However,it is not applicable for the nonlinear scheme.Thus,based on a‘cut-off’function and energy analysis method,we get unconditional L^(2)and H^(1)error estimates for the nonlinear scheme,as well as boundedness of numerical solutions.In addition,if the assumption for exact solutions is improved compared to before,unconditional and optimal pointwise error estimates can be obtained by energy analysis method and several Sobolev inequalities.Finally,some numerical examples are given to verify our theoretical analysis.展开更多
Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurat...Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurate and not affected by strong subsurface structure complexity and horizontal velocity variations. The difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method. It can be used directly for forward modeling on models with complex top surfaces and migration without statics preprocessing. We apply a finite difference method based on triangular grids for post-stack reverse-time migration for the first time. Tests on model data verify that the combination of the two methods can achieve near-perfect results in application.展开更多
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.展开更多
To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to...To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.展开更多
The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduct...The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduction in computer memory requirements and computational time. The computational domain is greatly reduced to enable performance in personal computer. At the same time because edges of a boundary and summits are treated well, the computational results is more accurate and more collector.展开更多
In this paper, the authors establish some theorems that can ascertain the zero solutions of systemsx(n+1)=f(n,x n)(1)are uniformly stable,asymptotically stable or uniformly asymptotically stable. In the obtained theo...In this paper, the authors establish some theorems that can ascertain the zero solutions of systemsx(n+1)=f(n,x n)(1)are uniformly stable,asymptotically stable or uniformly asymptotically stable. In the obtained theorems, ΔV is not required to be always negative, where ΔV(n,x n)≡V(n+1,x(n+1)) -V(n,x(n))=V(n+1,f(n,x n))-V(n,x(n)), especially, in Theorem 1, ΔV may be even positive, which greatly improve the known results and are more convenient to use.展开更多
Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is a...Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...展开更多
Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and ...Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.展开更多
The phase volume fraction has an important role in the match of the strength and plasticity of dual phase steel.The different bainite contents(18–53 vol.%)in polygonal ferrite and bainite(PF+B)dual phase steel were o...The phase volume fraction has an important role in the match of the strength and plasticity of dual phase steel.The different bainite contents(18–53 vol.%)in polygonal ferrite and bainite(PF+B)dual phase steel were obtained by controlling the relaxation finish temperature during the rolling process.The effect of bainite volume fraction on the tensile deformability was systematically investigated via experiments and crystal plasticity finite element model(CPFEM)simulation.The experimental results showed that the steel showed optimal strain hardenability and strength–plasticity matching when the bainite reached 35%.The 3D-CPFEM models with the same grain size and texture characters were established to clarify the influence of stress/strain distribution on PF+B dual phase steel with different bainite contents.The simulation results indicated that an appropriate increase in the bainite content(18%–35%)did not affect the interphase strain difference,but increased the stress distribution in both phases,as a result of enhancing the coordinated deformability of two phases and improving the strength–plasticity matching.When the bainite content increased to 53%,the stress/strain difference between the two phases was greatly increased,and plastic damage between the two phases was caused by the reduction of the coordinated deformability.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
基金Supported by Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)+3 种基金National Natural Science Foundation of China(12301556)Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)Basic Research Plan of Shanxi Province(202203021211129)。
文摘To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.
基金supported by the National Natural Science Foundation of China(No.42277175)Guangxi Emergency Management Department 2024 Innovation and Technology Research Project,China(No.2024GXYJ006)+2 种基金Hunan Provincial Department of Natural Resources Geological Exploration Project,China(No.2023ZRBSHZ056)The First National Natural Disaster Comprehensive Risk Survey in Hunan Province,China(No.2022-70)Guizhou Provincial Major Scientific and Technological Program,China(No.2023-425).
文摘A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficients(k)are constructed based on FDM.The rock bridge area was divided through k-means algorithm and the optimal number of clusters was determined by sum of squared errors(SSE)and elbow method.The influence of maximum principal stress and stress change rate as clustering indexes on the clustering results of rock bridges was compared by using Euclidean distance.The results show that using stress change rate as clustering index is more effective.When the joint coalescence coefficient is less than 0.6,there is no significant stress concentration in the middle area of adjacent joints,that is,no generation of rock bridge.In addition,the range of rock bridge is affected by the coalescence coefficient(k),the relative position of joints and the parameters of weak interlayer.
基金supported by the National Natural Science Foundation of China(12201228,12171047)the Fundamental Research Funds for the Central Universities(3034011102)supported by National Key R&D Program of China(2020YFA0713701).
文摘In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.
文摘This study investigates the complex heat transfer dynamics inmultilayer bifacial photovoltaic(bPV)solar modules under spectrally resolved solar irradiation.A novel numericalmodel is developed to incorporate internal heat generation resulting from optical absorption,grounded in the physical equations governing light-matter interactions within the module’smultilayer structure.The model accounts for reflection and transmission at each interface between adjacent layers,as well as absorption within individual layers,using the wavelength-dependent dielectric properties of constituent materials.These properties are used to calculate the spectral reflectance,transmittance,and absorption coefficients,enabling precise quantification of internal heat sources from irradiance incidents on both the front and rear surfaces of the module.The study further examines the influence of irradiance reflection on thermal behavior,evaluates the thermal impact of various supporting materials placed beneath the module,and analyzes the role of albedo in modifying heat distribution.By incorporating spectrally resolved heat generation across each layer often simplified or omitted in conventional models,the proposed approach enhances physical accuracy.The transient heat equation is solved using a one-dimensional finite difference(FD)method to produce detailed temperature profiles under multiple operating scenarios,including Standard Test Conditions(STC),Bifacial Standard Test Conditions(BSTC),Normal Operating Cell Temperature(NOCT),and Bifacial NOCT(BNOCT).The results offer valuable insights into the interplay between optical and thermal phenomena in bifacial systems,informing the design and optimization of more efficient photovoltaic technologies.
基金National Natural Science Foundation of China(no.41604037)Natural Science Foundation of Hubei Province(no.2022CFB125)+2 种基金Open Fund of Key Laboratory of Exploration Technologies for Oil and Gas Resources(Yangtze University)Ministry of Education(no.K2021-09)College Students'Innovation and Entrepreneurship Training Program(no.2019053)。
文摘Laplace–Fourier(L-F)domain finite-difference(FD)forward modeling is an important foundation for L-F domain full-waveform inversion(FWI).An optimal modeling method can improve the efficiency and accuracy of FWI.A fl exible FD stencil,which requires pairing and centrosymmetricity of the involved gridpoints,is used on the basis of the 2D L-F domain acoustic wave equation.The L-F domain numerical dispersion analysis is then performed by minimizing the phase error of the normalized numerical phase and attenuation propagation velocities to obtain the optimization coefficients.An optimal FD forward modeling method is finally developed for the L-F domain acoustic wave equation and applied to the traditional standard 9-point scheme and 7-and 9-point schemes,where the latter two schemes are used in discontinuous-grid FD modeling.Numerical experiments show that the optimal L-F domain FD modeling method not only has high accuracy but can also be applied to equal and unequal directional sampling intervals and discontinuous-grid FD modeling to reduce computational cost.
基金supported by the National Natural Science Foundation of China (grant numbers: 41874160 and 92055213)。
文摘Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and efficiency of the conventional numerical modeling, I develop a new seismic modeling method by combining the FD scheme with the numerical dispersion suppression neural network(NDSNN). This method involves the following steps. First, a training data set composed of a small number of wavefield snapshots is generated. The wavefield snapshots with the low-accuracy wavefield data and the high-accuracy wavefield data are paired, and the low-accuracy wavefield snapshots involve the obvious numerical dispersion including both the temporal and spatial dispersion. Second, the NDSNN is trained until the network converges to simultaneously suppress the temporal and spatial dispersion.Third, the entire set of low-accuracy wavefield data is computed quickly using FD modeling with the large time step and the coarse grid. Fourth, the NDSNN is applied to the entire set of low-accuracy wavefield data to suppress the numerical dispersion including the temporal and spatial dispersion.Numerical modeling examples verify the effectiveness of my proposed method in improving the computational accuracy and efficiency.
基金supported by the Key Laboratory of Road Construction Technology and Equipment(Chang’an University,No.300102253502)the Natural Science Foundation of Shandong Province of China(GrantNo.ZR2022YQ06)the Development Plan of Youth Innovation Team in Colleges and Universities of Shandong Province(Grant No.2022KJ140).
文摘In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.
基金This work was financially supported by the Key Science and Technology Project of Longmen Laboratory(No.LMYLKT-001)Innovation and Entrepreneurship Training Program for College Students of Henan Province(No.202310464050)。
文摘Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.
基金support via NSF grants NSF-19-04774,NSF-AST-2009776,NASA-2020-1241NASA grant 80NSSC22K0628.DSB+3 种基金HK acknowledge support from a Vajra award,VJR/2018/00129a travel grant from Notre Dame Internationalsupport via AFOSR grant FA9550-20-1-0055NSF grant DMS-2010107.
文摘Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we have even been able to apply the scheme without any modification and obtain good results.Many useful PDEs may have stiff relaxation source terms for which the finite difference formulation of WENO is shown to provide some genuine advantages.
文摘For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.
基金Supported by the National Natural Science Foundation of China(Grant No.11571181)the Research Start-Up Foundation of Nantong University(Grant No.135423602051).
文摘In this paper,we give improved error estimates for linearized and nonlinear CrankNicolson type finite difference schemes of Ginzburg-Landau equation in two dimensions.For linearized Crank-Nicolson scheme,we use mathematical induction to get unconditional error estimates in discrete L^(2)and H^(1)norm.However,it is not applicable for the nonlinear scheme.Thus,based on a‘cut-off’function and energy analysis method,we get unconditional L^(2)and H^(1)error estimates for the nonlinear scheme,as well as boundedness of numerical solutions.In addition,if the assumption for exact solutions is improved compared to before,unconditional and optimal pointwise error estimates can be obtained by energy analysis method and several Sobolev inequalities.Finally,some numerical examples are given to verify our theoretical analysis.
基金sponsored by National Natural Science Foundation(40474041)National Symposium of 863(2006AA06Z206)+1 种基金National Symposium of 973(2007CB209605)CNPC Geophysical Key Laboratory of the China University of Petroleum (East China) Research Department
文摘Compared with other migration methods, reverse-time migration is based on a precise wave equation, not an approximation, and performs extrapolation in the depth domain rather than the time domain. It is highly accurate and not affected by strong subsurface structure complexity and horizontal velocity variations. The difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method. It can be used directly for forward modeling on models with complex top surfaces and migration without statics preprocessing. We apply a finite difference method based on triangular grids for post-stack reverse-time migration for the first time. Tests on model data verify that the combination of the two methods can achieve near-perfect results in application.
基金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.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foun- dation of China National Petroleum Corporation.
文摘To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.
文摘The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduction in computer memory requirements and computational time. The computational domain is greatly reduced to enable performance in personal computer. At the same time because edges of a boundary and summits are treated well, the computational results is more accurate and more collector.
文摘In this paper, the authors establish some theorems that can ascertain the zero solutions of systemsx(n+1)=f(n,x n)(1)are uniformly stable,asymptotically stable or uniformly asymptotically stable. In the obtained theorems, ΔV is not required to be always negative, where ΔV(n,x n)≡V(n+1,x(n+1)) -V(n,x(n))=V(n+1,f(n,x n))-V(n,x(n)), especially, in Theorem 1, ΔV may be even positive, which greatly improve the known results and are more convenient to use.
文摘Contact bounce of relay, which is the main cause of electric abrasion and material erosion, is inevitable. By using the mode expansion form, the dynamic behavior of two different reed systems for aerospace relays is analyzed. The dynamic model uses Euler-Bernoulli beam theory for cantilever beam, in which the driving force (or driving moment) of the electromagnetic system is taken into account, and the contact force between moving contact and stationary contact is simulated by the Kelvin-Voigt vis-coelastic...
基金supported by the National Natural Science Foundation of China(No.42372287 and No.U24A20178)the Fundamental Research Funds for the Central Universities CHD(No.2024SHEEAR002)+3 种基金the Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shaanxi Province(No.2020024)the China Postdoctoral Science Foundation(GZC20232955,2024M753472,and 2024MD763937)the Science-Technology Foundation for Young Scientists of Gansu Province,China(No.24JRRA097)the Study of biodiversity survey and limiting factor analysis of Yinkentala(2023ZL01).
文摘Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.
基金supported by the Project of Liaoning Marine Economic Development(Development of high strength pipeline steel for submarine oil and gas transmission)State Key Laboratory of Metal Material for Marine Equipment and Application Funding(No.SKLMEA-K202205).
文摘The phase volume fraction has an important role in the match of the strength and plasticity of dual phase steel.The different bainite contents(18–53 vol.%)in polygonal ferrite and bainite(PF+B)dual phase steel were obtained by controlling the relaxation finish temperature during the rolling process.The effect of bainite volume fraction on the tensile deformability was systematically investigated via experiments and crystal plasticity finite element model(CPFEM)simulation.The experimental results showed that the steel showed optimal strain hardenability and strength–plasticity matching when the bainite reached 35%.The 3D-CPFEM models with the same grain size and texture characters were established to clarify the influence of stress/strain distribution on PF+B dual phase steel with different bainite contents.The simulation results indicated that an appropriate increase in the bainite content(18%–35%)did not affect the interphase strain difference,but increased the stress distribution in both phases,as a result of enhancing the coordinated deformability of two phases and improving the strength–plasticity matching.When the bainite content increased to 53%,the stress/strain difference between the two phases was greatly increased,and plastic damage between the two phases was caused by the reduction of the coordinated deformability.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
