The Fast Sweeping Method based on equivalent slowness is an improved algorithm for solving seismic wave traveltime,aiming to address the significant errors caused by the large wavefront curvature near the source point...The Fast Sweeping Method based on equivalent slowness is an improved algorithm for solving seismic wave traveltime,aiming to address the significant errors caused by the large wavefront curvature near the source point in traditional methods.The core of this method lies in transferring part of the complex curvature information of the traveltime eld to the distance term of the analytical solution,thereby simplifying the calculation process of the traveltime eld.Specically,by dening the equivalent slowness as the ratio of traveltime to the straight-line distance from a point to the source point,using a nite dierence scheme to discretely solve the equivalent slowness,and then converting it into traveltime values,this approach can signicantly enhance computational accuracy and efciency.In this study,numerical examples are employed to verify the effectiveness and superiority of the Fast Sweeping Method based on equivalent slowness in handling complex seismic wave traveltime calculations.展开更多
Accurate traveltime computation in complex geometries is crucial for seismological applications such as traveltime tomography and migration.However,eikonal solvers relying on finite-differences encounter the source si...Accurate traveltime computation in complex geometries is crucial for seismological applications such as traveltime tomography and migration.However,eikonal solvers relying on finite-differences encounter the source singularity,leading to numerical errors that propagate throughout the computational domain and compromise traveltime accuracy.Solving the factored eikonal equation has proven effective in addressing this singularity,but it has not yet been applied to triangular meshes,which are significant for modeling complex geological structures and irregular topographies.To address this challenge,we propose a new coordinate-transformation-based fast sweeping method(FsMcT)designed for triangular meshes.FSMCT maps each triangular element in the physical domain into a canonical reference triangle through coordinate transformation,and constructs an upwind finite-difference scheme on the reference triangle to solve the eikonal equation.The framework enables a straight forward extension to the multiplicatively/additively factored eikonal equations,overcoming the long-standing limitation that the conventional fast sweeping method(FsM)cannot solve the factored eikonal equation on triangular meshes,thereby effectively mitigating source singularity.Furthermore,FSMcT solves the eikonal equation on the reference triangle,while existing FSM typically require further subdivision of obtuse triangles.Numerical simulations demonstrate that FsMcT accurately solves the factored eikonal equation on triangular meshes,substantially reducing the errors caused by source singularity and achieving significantly higher accuracy than the conventional FSM.展开更多
A gradient-based optimization method for producing a contoured beam by using a single-fed reflector antenna is presented. First, a quick and accurate pattern approximation formula based on physical optics(PO) is adopt...A gradient-based optimization method for producing a contoured beam by using a single-fed reflector antenna is presented. First, a quick and accurate pattern approximation formula based on physical optics(PO) is adopted to calculate the gradients of the directivity with respect to reflector's nodal displacements. Because the approximation formula is a linear function of nodal displacements, the gradient can be easily derived. Then, the method of the steepest descent is adopted, and an optimization iteration procedure is proposed. The iteration procedure includes two loops: an inner loop and an outer loop. In the inner loop, the gradient and pattern are calculated by matrix operation, which is very fast by using the pre-calculated data in the outer loop. In the outer loop, the ideal terms used in the inner loop to calculate the gradient and pattern are updated, and the real pattern is calculated by the PO method. Due to the high approximation accuracy, when the outer loop is performed once, the inner loop can be performed many times, which will save much time because the integration is replaced by matrix operation. In the end, a contoured beam covering the continental United States(CONUS) is designed, and simulation results show the effectiveness of the proposed algorithm.展开更多
A steady thermo-hydraulic model of the helical tube steam generator was first constructed to study the coupled heat transfer process between the primary and secondary sides based on a discrete modeling method,and obta...A steady thermo-hydraulic model of the helical tube steam generator was first constructed to study the coupled heat transfer process between the primary and secondary sides based on a discrete modeling method,and obtain the heat flux density distribution along the steam generator.Then,taking the obtained coupled heat flux density distribution as the thermal boundary condition input,considering the dynamic variation of physical properties on the secondary side,a dynamic model based on the time-domain method suitable for two-phase flow instability among parallel multiple channels of the steam generator was constructed.Finally,taking the lead-bismuth fast reactor as an example,flow instability of the steam generator was analyzed under an inlet lead-bismuth temperature of 320℃~480℃ and an inlet water temperature of 160℃~240℃.It was found that flow instability is less likely to occur under coupled heat conditions,compared with that under uniform or linear distribution.Flow excursion is prone to occur under low inlet temperature of the primary or secondary side.As the inlet lead bismuth temperature increases from 320℃ to 480℃,average heat flux significantly increases by 2.5 times,and the non-uniformity of heat flux distribution increases of 49%.Meanwhile,the density wave oscillation amplitude gradually increases,and system stability weakens.展开更多
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%.展开更多
Fast beam migration(FBM),characterized by its super-high efficiency in velocity model building,consists of three main steps:beam forming,beam propagation,and image forming.The super-high efficiency is achieved by beam...Fast beam migration(FBM),characterized by its super-high efficiency in velocity model building,consists of three main steps:beam forming,beam propagation,and image forming.The super-high efficiency is achieved by beam forming,as it needs only to be performed once for one dataset and is independent of velocity,and the other two steps take relatively little time.However,compared to the beam-propagation and image-forming steps,the beam-forming step is still quite time-consuming owing to the high-dimensional computing problem of estimating the source and receiver slope orientation of a beam.Furthermore,previous methods for estimating the source and receiver slope orientation of a beam struggled to deal with intersecting events,leading to poor imaging results for complex subsurface structures,such as unconformities or faults,where events often intersect.We propose the use of a three-step multimodal optimization method based on the neighborhood crowding differential evolution(NCDE)algorithm to estimate the source and receiver slope orientation of a beam during the beam-forming step,which can quickly and accurately obtain slope orientations when events intersect.We first test the three-step multimodal optimization algorithm on a 3D super-gather and provide the parameter criteria.We then apply the FBM based on the three-step multimodal optimization algorithm to the Marmousi 2 and 3D SEG/EAGE salt models.Both results demonstrate that the proposed method can image intersecting events well and that the imaging quality of complex zones is improved.We also apply the proposed method to a 2D offshore seismic dataset containing abundant intersecting events,which validates the practicality of the proposed method.展开更多
This paper presents a fast algorithm for solving the scattering problem from an open rectangular cavity embedded in the ground plane.The computational region is chosen as the union of two rectangular regions:one is a ...This paper presents a fast algorithm for solving the scattering problem from an open rectangular cavity embedded in the ground plane.The computational region is chosen as the union of two rectangular regions:one is a region above the ground,the other one is a region containing the cavity.The finite difference scheme is constructed in each region.An intermediate layer of the mesh is shared by both regions,which is the key of the algorithm.A cyclic reduction method is employed to solve the difference equation in the region above the ground.Then the numerical solution on the cavity aperture can be obtained.The numerical experiments are provided to verify the feasibility of the proposed algorithm.展开更多
Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse ef...Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse effects on the final location results.To achieve a high-accuracy location in a complex cavern-containing structure,this study develops an MS location method using the fast marching method(FMM)with a second-order difference approach(FMM2).Based on the established velocity model with three-dimensional(3D)discrete grids,the realization of the MS location can be achieved by searching the minimum residual between the theoretical and actual first arrival times.Moreover,based on the calculation results of FMM2,the propagation paths from the MS sources to MS sensors can be obtained using the linear interpolation approach and the Runge–Kutta method.These methods were validated through a series of numerical experiments.In addition,our proposed method was applied to locate the recorded blasting and MS events that occurred during the excavation period of the underground caverns at the Houziyan hydropower station.The location results of the blasting activities show that our method can effectively reduce the location error compared with the results based on the uniform velocity model.Furthermore,the obtained MS location was verified through the occurrence of shotcrete fractures and spalling,and the monitoring results of the in-situ multipoint extensometer.Our proposed method can offer a more accurate rock fracture location and facilitate the delineation of damage zones inside the surrounding rock mass.展开更多
A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be effi...3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position inf...In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.展开更多
In 2D fast multipole method for scattering problems,square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel,and numerical integration error is introduced...In 2D fast multipole method for scattering problems,square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel,and numerical integration error is introduced. Taking advantage of the relationship between Euler-Maclaurin formula and trapezoidal quadrature rule,and the relationship between trapezoidal and square quadrature rule,sharp computable bound with analytical form on the error of numerical integration of Bessel integral identity by square quadrature rule is derived in this paper. Numerical experiments are presented at the end to demonstrate the accuracy of the sharp computable bound on the numerical integration error.展开更多
In actual power systems,most of the high-voltage buses of the transformers are zero injection buses without load or generation.Power injections into these buses are strictly 0,so based on Kirchhoff's current law(K...In actual power systems,most of the high-voltage buses of the transformers are zero injection buses without load or generation.Power injections into these buses are strictly 0,so based on Kirchhoff's current law(KCL),equality constraints should be used to handle these buses in a state estimation model.It is a challenge to ensure that these zero injection constraints can be strictly satisfied without losing computational efficiency.展开更多
Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternati...Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions.A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems.Hence,they are easy to be applied to a general hyperbolic system.To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains,inverse Lax-Wendroff(ILW)procedures were developed as a very effective approach in the literature.In this paper,we combine a fifthorder fixed-point fast sweeping WENO method with an ILW procedure to solve steadystate solution of hyperbolic conservation laws on complex computing regions.Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids.Numerical results show highorder accuracy and good performance of the method.Furthermore,the method is compared with the popular third-order total variation diminishing Runge-Kutta(TVD-RK3)time-marching method for steady-state computations.Numerical examples show that for most of examples,the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.展开更多
For accurate trajectory tracking and obstacle avoidance in finite time of a nonholonomic mobile robot,a trajectory tracking controller based on global fast terminal sliding mode method is proposed,which has the advant...For accurate trajectory tracking and obstacle avoidance in finite time of a nonholonomic mobile robot,a trajectory tracking controller based on global fast terminal sliding mode method is proposed,which has the advantages of chattering-free and adjustable convergence time.First of all,the kinematics model of the robot is established in mobile carrier coordinates.Secondly,the global structure including terminal attractor and exponential convergence of the fast terminal sliding mode trajectory tracking controller is proved by Lyapunov stability theory,ensuring that the trajectory and heading angle tracking error converges to a smaller zero range in finite time.Finally,the artificial potential field obstacle avoidance method is introduced to make the robot not only track the reference trajectory strictly,but also avoid the obstacles.The simulation results show that the proposed method can achieve a stable tracking control in finite time for a given reference trajectory.展开更多
This paper is dedicated to applying the Fourier amplitude sensitivity test(FAST)method to the problem of mixed extension and inflation of a circular cylindrical tube in the presence of residual stresses.The metafuncti...This paper is dedicated to applying the Fourier amplitude sensitivity test(FAST)method to the problem of mixed extension and inflation of a circular cylindrical tube in the presence of residual stresses.The metafunctions and the Ishigami function are considered in the sensitivity analysis(SA).The effects of the input variables on the output variables are investigated,and the most important parameters of the system under the applied pressure and axial force such as the axial stretch and the azimuthal stretch are determined.展开更多
基金supported by the Key Project of the Spark Program of Earthquake Science and Technology in Hebei Province(DZ2025091100001 and DZ2025081200001)University Stability Support Program of Shenzhen-General Project(20231128113233002).
摘要The Fast Sweeping Method based on equivalent slowness is an improved algorithm for solving seismic wave traveltime,aiming to address the significant errors caused by the large wavefront curvature near the source point in traditional methods.The core of this method lies in transferring part of the complex curvature information of the traveltime eld to the distance term of the analytical solution,thereby simplifying the calculation process of the traveltime eld.Specically,by dening the equivalent slowness as the ratio of traveltime to the straight-line distance from a point to the source point,using a nite dierence scheme to discretely solve the equivalent slowness,and then converting it into traveltime values,this approach can signicantly enhance computational accuracy and efciency.In this study,numerical examples are employed to verify the effectiveness and superiority of the Fast Sweeping Method based on equivalent slowness in handling complex seismic wave traveltime calculations.
基金supported by the National Natural Science Foundation of China(42325403,42074162)Xin Fu was additionally supported by the National Natural Science Foundation of China(No.42504116)+1 种基金the Deep Earth Probe and Mineral Resources Exploration-National Science and Technology Major Project(No.2024zD1004201)the Young Expert of Taishan ScholarsProject(No.tsqn202408095).
摘要Accurate traveltime computation in complex geometries is crucial for seismological applications such as traveltime tomography and migration.However,eikonal solvers relying on finite-differences encounter the source singularity,leading to numerical errors that propagate throughout the computational domain and compromise traveltime accuracy.Solving the factored eikonal equation has proven effective in addressing this singularity,but it has not yet been applied to triangular meshes,which are significant for modeling complex geological structures and irregular topographies.To address this challenge,we propose a new coordinate-transformation-based fast sweeping method(FsMcT)designed for triangular meshes.FSMCT maps each triangular element in the physical domain into a canonical reference triangle through coordinate transformation,and constructs an upwind finite-difference scheme on the reference triangle to solve the eikonal equation.The framework enables a straight forward extension to the multiplicatively/additively factored eikonal equations,overcoming the long-standing limitation that the conventional fast sweeping method(FsM)cannot solve the factored eikonal equation on triangular meshes,thereby effectively mitigating source singularity.Furthermore,FSMcT solves the eikonal equation on the reference triangle,while existing FSM typically require further subdivision of obtuse triangles.Numerical simulations demonstrate that FsMcT accurately solves the factored eikonal equation on triangular meshes,substantially reducing the errors caused by source singularity and achieving significantly higher accuracy than the conventional FSM.
基金supported by the National Natural Science Foundation of China(51805399)the Fundamental Research Funds for the Central Universities(JB180403)+2 种基金the Chinese Academy of Sciences(CAS)"Light of West China" Program(2017-XBQNXZ-B-024)the National Basic Research Program of China(973 Program)(2015CB857100)the Operation,Maintenance and Upgrading Fund for Astronomical Telescopes and Facility Instruments,budgeted from the Ministry of Finance of China(MOF)and administrated by the CAS
摘要A gradient-based optimization method for producing a contoured beam by using a single-fed reflector antenna is presented. First, a quick and accurate pattern approximation formula based on physical optics(PO) is adopted to calculate the gradients of the directivity with respect to reflector's nodal displacements. Because the approximation formula is a linear function of nodal displacements, the gradient can be easily derived. Then, the method of the steepest descent is adopted, and an optimization iteration procedure is proposed. The iteration procedure includes two loops: an inner loop and an outer loop. In the inner loop, the gradient and pattern are calculated by matrix operation, which is very fast by using the pre-calculated data in the outer loop. In the outer loop, the ideal terms used in the inner loop to calculate the gradient and pattern are updated, and the real pattern is calculated by the PO method. Due to the high approximation accuracy, when the outer loop is performed once, the inner loop can be performed many times, which will save much time because the integration is replaced by matrix operation. In the end, a contoured beam covering the continental United States(CONUS) is designed, and simulation results show the effectiveness of the proposed algorithm.
基金supported by Natural Science Basic Research Program of Shaanxi Province(Grant Nos.2025JC-YBMS-475,2025JC-YBMS-420)Joint Funds of the National Natural Science Foundation of China(Grant No.U20B2036)National Natural Science Foundation of China(Grant No.52274064).
摘要A steady thermo-hydraulic model of the helical tube steam generator was first constructed to study the coupled heat transfer process between the primary and secondary sides based on a discrete modeling method,and obtain the heat flux density distribution along the steam generator.Then,taking the obtained coupled heat flux density distribution as the thermal boundary condition input,considering the dynamic variation of physical properties on the secondary side,a dynamic model based on the time-domain method suitable for two-phase flow instability among parallel multiple channels of the steam generator was constructed.Finally,taking the lead-bismuth fast reactor as an example,flow instability of the steam generator was analyzed under an inlet lead-bismuth temperature of 320℃~480℃ and an inlet water temperature of 160℃~240℃.It was found that flow instability is less likely to occur under coupled heat conditions,compared with that under uniform or linear distribution.Flow excursion is prone to occur under low inlet temperature of the primary or secondary side.As the inlet lead bismuth temperature increases from 320℃ to 480℃,average heat flux significantly increases by 2.5 times,and the non-uniformity of heat flux distribution increases of 49%.Meanwhile,the density wave oscillation amplitude gradually increases,and system stability weakens.
基金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%.
基金Natural Science Foundation of China(42304128 and 42074150)the National Key Research and Development Program of China(2023YFC2906704-5 and 2023YFC370790)+1 种基金the Mount Tai Industry Leading Talent Project Special Fund Support(tscx202312018)the Jinan Science and Technology Innovation Development Plan(Social Livelihood Special Project)(202131001)。
摘要Fast beam migration(FBM),characterized by its super-high efficiency in velocity model building,consists of three main steps:beam forming,beam propagation,and image forming.The super-high efficiency is achieved by beam forming,as it needs only to be performed once for one dataset and is independent of velocity,and the other two steps take relatively little time.However,compared to the beam-propagation and image-forming steps,the beam-forming step is still quite time-consuming owing to the high-dimensional computing problem of estimating the source and receiver slope orientation of a beam.Furthermore,previous methods for estimating the source and receiver slope orientation of a beam struggled to deal with intersecting events,leading to poor imaging results for complex subsurface structures,such as unconformities or faults,where events often intersect.We propose the use of a three-step multimodal optimization method based on the neighborhood crowding differential evolution(NCDE)algorithm to estimate the source and receiver slope orientation of a beam during the beam-forming step,which can quickly and accurately obtain slope orientations when events intersect.We first test the three-step multimodal optimization algorithm on a 3D super-gather and provide the parameter criteria.We then apply the FBM based on the three-step multimodal optimization algorithm to the Marmousi 2 and 3D SEG/EAGE salt models.Both results demonstrate that the proposed method can image intersecting events well and that the imaging quality of complex zones is improved.We also apply the proposed method to a 2D offshore seismic dataset containing abundant intersecting events,which validates the practicality of the proposed method.
基金Supported by the National Natural Science Foundation of China(12101205)the Natural Science Foundation of Heilongjiang Province of China(PL2024A010)。
摘要This paper presents a fast algorithm for solving the scattering problem from an open rectangular cavity embedded in the ground plane.The computational region is chosen as the union of two rectangular regions:one is a region above the ground,the other one is a region containing the cavity.The finite difference scheme is constructed in each region.An intermediate layer of the mesh is shared by both regions,which is the key of the algorithm.A cyclic reduction method is employed to solve the difference equation in the region above the ground.Then the numerical solution on the cavity aperture can be obtained.The numerical experiments are provided to verify the feasibility of the proposed algorithm.
基金the Key Program of National Natural Science Foundation of China(52039007)for providing financial support.
摘要Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse effects on the final location results.To achieve a high-accuracy location in a complex cavern-containing structure,this study develops an MS location method using the fast marching method(FMM)with a second-order difference approach(FMM2).Based on the established velocity model with three-dimensional(3D)discrete grids,the realization of the MS location can be achieved by searching the minimum residual between the theoretical and actual first arrival times.Moreover,based on the calculation results of FMM2,the propagation paths from the MS sources to MS sensors can be obtained using the linear interpolation approach and the Runge–Kutta method.These methods were validated through a series of numerical experiments.In addition,our proposed method was applied to locate the recorded blasting and MS events that occurred during the excavation period of the underground caverns at the Houziyan hydropower station.The location results of the blasting activities show that our method can effectively reduce the location error compared with the results based on the uniform velocity model.Furthermore,the obtained MS location was verified through the occurrence of shotcrete fractures and spalling,and the monitoring results of the in-situ multipoint extensometer.Our proposed method can offer a more accurate rock fracture location and facilitate the delineation of damage zones inside the surrounding rock mass.
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
摘要A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
基金supported by NSFC(Nos.41274120,41404085,and 41504084)
摘要3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
摘要A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
基金Projects(40804027,41074085) supported by the National Natural Science Foundation of ChinaProject(09JJ3048) supported by the Natural Science Foundation of Hunan Province,ChinaProject(200805331082) supported by the Research Fund for the Doctoral Program of Higher Education,China
摘要In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.
基金the National Natural Science Foundation of China (No. 11074170)the Independent Research Program of State Key Laboratory of Machinery System and Vibration (SKLMSV) (No. MSV-MS-2008-05)the Visiting Scholar Program of SKLMSV (No. MSV-2009-06)
摘要In 2D fast multipole method for scattering problems,square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel,and numerical integration error is introduced. Taking advantage of the relationship between Euler-Maclaurin formula and trapezoidal quadrature rule,and the relationship between trapezoidal and square quadrature rule,sharp computable bound with analytical form on the error of numerical integration of Bessel integral identity by square quadrature rule is derived in this paper. Numerical experiments are presented at the end to demonstrate the accuracy of the sharp computable bound on the numerical integration error.
摘要In actual power systems,most of the high-voltage buses of the transformers are zero injection buses without load or generation.Power injections into these buses are strictly 0,so based on Kirchhoff's current law(KCL),equality constraints should be used to handle these buses in a state estimation model.It is a challenge to ensure that these zero injection constraints can be strictly satisfied without losing computational efficiency.
基金Research was supported by the NSFC Grant 11872210Research was supported by the NSFC Grant 11872210 and Grant No.MCMS-I-0120G01+1 种基金Research supported in part by the AFOSR Grant FA9550-20-1-0055NSF Grant DMS-2010107.
摘要Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions.A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems.Hence,they are easy to be applied to a general hyperbolic system.To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains,inverse Lax-Wendroff(ILW)procedures were developed as a very effective approach in the literature.In this paper,we combine a fifthorder fixed-point fast sweeping WENO method with an ILW procedure to solve steadystate solution of hyperbolic conservation laws on complex computing regions.Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids.Numerical results show highorder accuracy and good performance of the method.Furthermore,the method is compared with the popular third-order total variation diminishing Runge-Kutta(TVD-RK3)time-marching method for steady-state computations.Numerical examples show that for most of examples,the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.
基金National Natural Science Foundation of China(No.61673042)Shanxi Province Science Foundation for Youths(No.201701D221123)。
摘要For accurate trajectory tracking and obstacle avoidance in finite time of a nonholonomic mobile robot,a trajectory tracking controller based on global fast terminal sliding mode method is proposed,which has the advantages of chattering-free and adjustable convergence time.First of all,the kinematics model of the robot is established in mobile carrier coordinates.Secondly,the global structure including terminal attractor and exponential convergence of the fast terminal sliding mode trajectory tracking controller is proved by Lyapunov stability theory,ensuring that the trajectory and heading angle tracking error converges to a smaller zero range in finite time.Finally,the artificial potential field obstacle avoidance method is introduced to make the robot not only track the reference trajectory strictly,but also avoid the obstacles.The simulation results show that the proposed method can achieve a stable tracking control in finite time for a given reference trajectory.
摘要This paper is dedicated to applying the Fourier amplitude sensitivity test(FAST)method to the problem of mixed extension and inflation of a circular cylindrical tube in the presence of residual stresses.The metafunctions and the Ishigami function are considered in the sensitivity analysis(SA).The effects of the input variables on the output variables are investigated,and the most important parameters of the system under the applied pressure and axial force such as the axial stretch and the azimuthal stretch are determined.