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.展开更多
3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of ...Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.展开更多
This paper proposes a hybrid method based on the forward-backward method (FBM) and the reciprocity theorem (RT) for evaluating the scattering field from dielectric rough surface with a 2D target above it. Here, th...This paper proposes a hybrid method based on the forward-backward method (FBM) and the reciprocity theorem (RT) for evaluating the scattering field from dielectric rough surface with a 2D target above it. Here, the equivalent electric/magnetic current densities on the rough surface as well as the scattering field from it are numerically calculated by FBM, and the scattered field from the isolated target is obtained utilizing the method of moments (MOM). Meanwhile, the rescattered coupling interactions between the target and the surface are evaluated employing the combination of FBM and RT. Our hybrid method is first validated by available MOM results. Then, the functional dependences of bistatic and monostatic scattering from the target above rough surface upon the target altitude, incident and scattering angles are numerically simulated and discussed. This study presents a numerical description for the scattering mechanism associated with rescattered coupling interactions between a target and an underlying randomly rough surface.展开更多
Fatigue crack propagation characteristics of a diesel engine crankshaft are studied by measuring the fatigue crack growth rate using a frequency sweep method on a resonant fatigue test rig. Based on the phenomenon tha...Fatigue crack propagation characteristics of a diesel engine crankshaft are studied by measuring the fatigue crack growth rate using a frequency sweep method on a resonant fatigue test rig. Based on the phenomenon that the system frequency will change when the crack becomes large, this method can be directly applied to a complex component or structure. Finite element analyses (FEAs) are performed to calibrate the relation between the frequency change and the crack size, and to obtain the natural frequency of the test rig and the stress intensity factor (SIF) of growing cracks. The crack growth rate i.e. da/dN-AK of each crack size is obtained by combining the testing-time monitored data and FEA results. The results show that the crack growth rate of engine crankshaft, which is a component with complex geometry and special surface treatment, is quite different from that of a pure material. There is an apparent turning point in the Paris's crack partition. The cause of the fatigue crack growth is also discussed.展开更多
With the development of distribution automation system, the centralized meter reading system has been adopted more and more extensively, which provides real-time electricity consumption data of end-users, and conseque...With the development of distribution automation system, the centralized meter reading system has been adopted more and more extensively, which provides real-time electricity consumption data of end-users, and consequently lays foundation for operating condition on-line analysis of distribution network. In this paper, a modified back/forward sweep method, which directly uses real-time electricity consumption data acquired from the centralized meter reading system, is proposedto realize voltage analysis based on 24-hour electricity consumption data of a typical transformer district. Furthermore, the calculated line losses are verified through data collected from the energy metering of the distribution transformer, illustrating that the proposed method can be applied in analyzing voltage level and discovering unknown energy losses, which will lay foundation for on-line analysis, calculation and monitoring of power distribution network.展开更多
A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-or...A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-order central difference scheme at the middle interface are used. Maximum norm error estimate for the procedure is derived. Then an iterative method based on domain decomposition is presented for the numerical scheme and the convergence of the given method is established. Then numerical experiments are presented to support the theoretical analysis.展开更多
A general type of forward-backward doubly stochastic differential equations (FBDSDEs) is studied. It extends many important equations that have been well studied, including stochastic Hamiltonian systems. Under some...A general type of forward-backward doubly stochastic differential equations (FBDSDEs) is studied. It extends many important equations that have been well studied, including stochastic Hamiltonian systems. Under some much weaker monotonicity assumptions, the existence and uniqueness of measurable solutions are established with a incthod of continuation. Furthermore, the continuity and differentiability of the solutions to FBDSDEs depending on parameters is discussed.展开更多
Introduction Sensitive linear sweep voltammetry to have been applied to the catalytic kinetic analysis may be helpful for the application of some organic reagents to the analysis of inorganic ones by voltammetry, The ...Introduction Sensitive linear sweep voltammetry to have been applied to the catalytic kinetic analysis may be helpful for the application of some organic reagents to the analysis of inorganic ones by voltammetry, The catalytic reaction-linear sweep voltammetric method for the determination of noble metals, such as Ag, Os, Ru and Ir has been developed, but no report about the method for the determination of rhodium has been published so far.展开更多
To improve the transmission performance of XCTD channel, this paper proposes a method to measure directly and fit the channel transmission characteristics by using frequency sweeping method. Sinusoidal signals with a ...To improve the transmission performance of XCTD channel, this paper proposes a method to measure directly and fit the channel transmission characteristics by using frequency sweeping method. Sinusoidal signals with a frequency range of 100 Hz to 10 k Hz and an interval of 100 Hz are used to measure transmission characteristics of channels with lengths of 300 m, 800 m, 1300 m, and 1800 m. The correctness of the fitted channel characteristics by transmitting square wave, composite waves of different frequencies, and ASK modulation are verified. The results show that when the frequency of the signal is below 1500 Hz, the channel has very little effect on the signal. The signal compensated for amplitude and phase at the receiver is not as good as the uncompensated signal.Alternatively, when the signal frequency is above 1500 Hz, the channel distorts the signal. The quality of signal compensated for amplitude and phase at receiver is better than that of the uncompensated signal. Thus, we can select the appropriate frequency for XCTD system and the appropriate way to process the received signals. Signals below1500 Hz can be directly used at the receiving end. Signals above 1500 Hz are used after amplitude and phase compensation at the receiving end.展开更多
基金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.
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
文摘3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
文摘Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.
基金Project supported by the National Natural Science Foundation of China (Grant No 60571058)the National Defense Foundation of China
文摘This paper proposes a hybrid method based on the forward-backward method (FBM) and the reciprocity theorem (RT) for evaluating the scattering field from dielectric rough surface with a 2D target above it. Here, the equivalent electric/magnetic current densities on the rough surface as well as the scattering field from it are numerically calculated by FBM, and the scattered field from the isolated target is obtained utilizing the method of moments (MOM). Meanwhile, the rescattered coupling interactions between the target and the surface are evaluated employing the combination of FBM and RT. Our hybrid method is first validated by available MOM results. Then, the functional dependences of bistatic and monostatic scattering from the target above rough surface upon the target altitude, incident and scattering angles are numerically simulated and discussed. This study presents a numerical description for the scattering mechanism associated with rescattered coupling interactions between a target and an underlying randomly rough surface.
文摘Fatigue crack propagation characteristics of a diesel engine crankshaft are studied by measuring the fatigue crack growth rate using a frequency sweep method on a resonant fatigue test rig. Based on the phenomenon that the system frequency will change when the crack becomes large, this method can be directly applied to a complex component or structure. Finite element analyses (FEAs) are performed to calibrate the relation between the frequency change and the crack size, and to obtain the natural frequency of the test rig and the stress intensity factor (SIF) of growing cracks. The crack growth rate i.e. da/dN-AK of each crack size is obtained by combining the testing-time monitored data and FEA results. The results show that the crack growth rate of engine crankshaft, which is a component with complex geometry and special surface treatment, is quite different from that of a pure material. There is an apparent turning point in the Paris's crack partition. The cause of the fatigue crack growth is also discussed.
文摘With the development of distribution automation system, the centralized meter reading system has been adopted more and more extensively, which provides real-time electricity consumption data of end-users, and consequently lays foundation for operating condition on-line analysis of distribution network. In this paper, a modified back/forward sweep method, which directly uses real-time electricity consumption data acquired from the centralized meter reading system, is proposedto realize voltage analysis based on 24-hour electricity consumption data of a typical transformer district. Furthermore, the calculated line losses are verified through data collected from the energy metering of the distribution transformer, illustrating that the proposed method can be applied in analyzing voltage level and discovering unknown energy losses, which will lay foundation for on-line analysis, calculation and monitoring of power distribution network.
基金Supported by National Science Foundation of China(Grant 10871179)the National Basic Research Programme of China(Grant 2008CB717806)the Department of Education of Zhejiang Province(GrantY200803559).
文摘A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-order central difference scheme at the middle interface are used. Maximum norm error estimate for the procedure is derived. Then an iterative method based on domain decomposition is presented for the numerical scheme and the convergence of the given method is established. Then numerical experiments are presented to support the theoretical analysis.
基金supported by the National Natural Science Foundation of China (No. 10771122)the NaturalScience Foundation of Shandong Province of China (No. Y2006A08)the National Basic ResearchProgram of China (973 Program) (No. 2007CB814900)
文摘A general type of forward-backward doubly stochastic differential equations (FBDSDEs) is studied. It extends many important equations that have been well studied, including stochastic Hamiltonian systems. Under some much weaker monotonicity assumptions, the existence and uniqueness of measurable solutions are established with a incthod of continuation. Furthermore, the continuity and differentiability of the solutions to FBDSDEs depending on parameters is discussed.
基金Supported by the National Natural Science Foundation of China
文摘Introduction Sensitive linear sweep voltammetry to have been applied to the catalytic kinetic analysis may be helpful for the application of some organic reagents to the analysis of inorganic ones by voltammetry, The catalytic reaction-linear sweep voltammetric method for the determination of noble metals, such as Ag, Os, Ru and Ir has been developed, but no report about the method for the determination of rhodium has been published so far.
基金financially supported by the National Key Research and Development Program of China(Grant No.2016YFC1400400)
文摘To improve the transmission performance of XCTD channel, this paper proposes a method to measure directly and fit the channel transmission characteristics by using frequency sweeping method. Sinusoidal signals with a frequency range of 100 Hz to 10 k Hz and an interval of 100 Hz are used to measure transmission characteristics of channels with lengths of 300 m, 800 m, 1300 m, and 1800 m. The correctness of the fitted channel characteristics by transmitting square wave, composite waves of different frequencies, and ASK modulation are verified. The results show that when the frequency of the signal is below 1500 Hz, the channel has very little effect on the signal. The signal compensated for amplitude and phase at the receiver is not as good as the uncompensated signal.Alternatively, when the signal frequency is above 1500 Hz, the channel distorts the signal. The quality of signal compensated for amplitude and phase at receiver is better than that of the uncompensated signal. Thus, we can select the appropriate frequency for XCTD system and the appropriate way to process the received signals. Signals below1500 Hz can be directly used at the receiving end. Signals above 1500 Hz are used after amplitude and phase compensation at the receiving end.