The present paper is devoted to the convergence control and accelerating the traditional Decomposition Methodof Adomian (ADM). By means of perturbing the initial or early terms of the Adomian iterates by adding aparam...The present paper is devoted to the convergence control and accelerating the traditional Decomposition Methodof Adomian (ADM). By means of perturbing the initial or early terms of the Adomian iterates by adding aparameterized term, containing an embedded parameter, new modified ADM is constructed. The optimal value ofthis parameter is later determined via squared residual minimizing the error. The failure of the classical ADM is alsoprevented by a suitable value of the embedded parameter, particularly beneficial for the Duan–Rach modification ofthe ADM incorporating all the boundaries into the formulation. With the presented squared residual error analysis,there is no need to check out the results against the numerical ones, as usually has to be done in the traditional ADMstudies to convince the readers that the results are indeed converged to the realistic solutions. Physical examplesselected from the recent application of ADM demonstrate the validity, accuracy and power of the presented novelapproach in this paper. Hence, the highly nonlinear equations arising from engineering applications can be safelytreated by the outlined method for which the classical ADM may fail or be slow to converge.展开更多
The key objective of this paper is to improve the approximation of a sufficiently smooth nonperiodic function defined on a compact interval by proposing alternative forms of Fourier series expansions. Unlike in classi...The key objective of this paper is to improve the approximation of a sufficiently smooth nonperiodic function defined on a compact interval by proposing alternative forms of Fourier series expansions. Unlike in classical Fourier series, the expansion coefficients herein are explicitly dependent not only on the function itself, but also on its derivatives at the ends of the interval. Each of these series expansions can be made to converge faster at a desired polynomial rate. These results have useful implications to Fourier or harmonic analysis, solutions to differential equations and boundary value problems, data compression, and so on.展开更多
In the author’s recent publications, a parametric system biorthogonal to the corresponding segment of the exponential Fourier system was unusually effective. On its basis, it was discovered that knowledge of a finite...In the author’s recent publications, a parametric system biorthogonal to the corresponding segment of the exponential Fourier system was unusually effective. On its basis, it was discovered that knowledge of a finite number of Fourier coefficients of function f from an infinite-dimensional set of elementary functions allows f to be accurately restored (the phenomenon of over-convergence). Below, parametric biorthogonal systems are constructed for classical trigonometric Fourier series, and the corresponding phenomena of over-convergence are discovered. The decisive role here was played by representing the space L2 as an orthogonal sum of two corresponding subspaces. As a result, fast parallel algorithms for reconstructing a function from its truncated trigonometric Fourier series are proposed. The presented numerical experiments confirm the high efficiency of these convergence accelerations for smooth functions. In conclusion, the main results of the work are summarized, and some prospects for the development and generalization of the proposed approaches are discussed.展开更多
In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various accelerati...In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.展开更多
The current paper considers the problem of recovering a function using a limited number of its Fourier coefficients. Specifically, a method based on Bernoulli-like polynomials suggested and developed by Krylov, Lanczo...The current paper considers the problem of recovering a function using a limited number of its Fourier coefficients. Specifically, a method based on Bernoulli-like polynomials suggested and developed by Krylov, Lanczos, Gottlieb and Eckhoff is examined. Asymptotic behavior of approximate calculation of the so-called "jumps" is studied and asymptotic L2 constants of the rate of convergence of the method are computed.展开更多
Convergence acceleration of the classical trigonometric interpolation by the Eckhoff method is considered, where the exact values of the "jumps" are approximated by solution of a system of linear equations. The accu...Convergence acceleration of the classical trigonometric interpolation by the Eckhoff method is considered, where the exact values of the "jumps" are approximated by solution of a system of linear equations. The accuracy of the "jump" approximation is explored and the corresponding asymptotic error of interpolation is derived. Numerical results validate theoretical estimates.展开更多
In this paper, a new method, so called A-method, is given for the convergence analysis of the MQ-algorithm. And the finer relaxation parameter θA is obtained. The numerical results show that our new method has the ou...In this paper, a new method, so called A-method, is given for the convergence analysis of the MQ-algorithm. And the finer relaxation parameter θA is obtained. The numerical results show that our new method has the outstanding effect of accelerating convergence. Moreover, the relaxation parameter θA is the optimum in a point of view.展开更多
In this paper, we address one of the issues in the frequency assignment problem for cellular mobile networks in which we intend to minimize the interference levels when assigning frequencies from a limited frequency s...In this paper, we address one of the issues in the frequency assignment problem for cellular mobile networks in which we intend to minimize the interference levels when assigning frequencies from a limited frequency spectrum. In order to satisfy the increasing demand in such cellular mobile networks, we use a hybrid approach consisting of a Particle Swarm Optimization(PSO) combined with a Tabu Search(TS) algorithm. This approach takes both advantages of PSO efficiency in global optimization and TS in avoiding the premature convergence that would lead PSO to stagnate in a local minimum. Moreover, we propose a new efficient, simple, and inexpensive model for storing and evaluating solution's assignment. The purpose of this model reduces the solution's storage volume as well as the computations required to evaluate thesesolutions in comparison with the classical model. Our simulation results on the most known benchmarking instances prove the effectiveness of our proposed algorithm in comparison with previous related works in terms of convergence rate, the number of iterations, the solution storage volume and the running time required to converge to the optimal solution.展开更多
The paper considers the Krylov-Lanczos and the Eckhoff approximations for recovering a bivariate function using limited number of its Fourier coefficients. These approximations are based on certain corrections associa...The paper considers the Krylov-Lanczos and the Eckhoff approximations for recovering a bivariate function using limited number of its Fourier coefficients. These approximations are based on certain corrections associated with jumps in the partial derivatives of the approximated function. Approximation of the exact jumps is accomplished by solution of systems of linear equations along the idea of Eckhoff. Asymptotic behaviors of the approximate jumps and the Eckhoff approximation are studied. Exact constants of the asymptotic errors are computed. Numerical experiments validate theoretical investigations.展开更多
In this paper, we show that the coupled modified Kd V equations possess rich mathematical structures and some remarkable properties. The connections between the system and skew orthogonal polynomials,convergence accel...In this paper, we show that the coupled modified Kd V equations possess rich mathematical structures and some remarkable properties. The connections between the system and skew orthogonal polynomials,convergence acceleration algorithms and Laurent property are discussed in detail.展开更多
The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of s...The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of slip surface is given, and the formula-tions of FERA based on incremental tangent stiffness method and modified Aitken accelerating algorithm are developed. The limited step length iteration method (LSLIM) is adopted to calculate the reliability index. The nonlinear FERA code using the SSA technique is developed and the main flow chart is illustrated. Numerical examples are used to demonstrate the efficiency and robustness of this method. It is found that the accelerating convergence algorithm proposed in this study proves to be very efficient for it can reduce the iteration number greatly, and LSLIM is also efficient for it can assure the convergence of the iteration of the reliability index.展开更多
The Hankel transform is widely used to solve various engineering and physics problems,such as the representation of electromagnetic field components in the medium,the representation of dynamic stress intensity factors...The Hankel transform is widely used to solve various engineering and physics problems,such as the representation of electromagnetic field components in the medium,the representation of dynamic stress intensity factors,vibration of axisymmetric infinite membrane and displacement intensity factors which all involve this type of integration.However,traditional numerical integration algorithms cannot be used due to the high oscillation characteristics of the Bessel function,so it is particularly important to propose a high precision and efficient numerical algorithm for calculating the integral of high oscillation.In this paper,the improved Gaver-Stehfest(G-S)inverse Laplace transform method for arbitrary real-order Bessel function integration is presented by using the asymptotic characteristics of the Bessel function and the accumulation of integration,and the optimized G-S coefficients are given.The effectiveness of the algorithm is verified by numerical examples.Compared with the linear transformation accelerated convergence algorithm,it shows that the G-S inverse Laplace transform method is suitable for arbitrary real order Hankel transform,and the time consumption is relatively stable and short,which provides a reliable calculation method for the study of electromagnetic mechanics,wave propagation,and fracture dynamics.展开更多
文摘The present paper is devoted to the convergence control and accelerating the traditional Decomposition Methodof Adomian (ADM). By means of perturbing the initial or early terms of the Adomian iterates by adding aparameterized term, containing an embedded parameter, new modified ADM is constructed. The optimal value ofthis parameter is later determined via squared residual minimizing the error. The failure of the classical ADM is alsoprevented by a suitable value of the embedded parameter, particularly beneficial for the Duan–Rach modification ofthe ADM incorporating all the boundaries into the formulation. With the presented squared residual error analysis,there is no need to check out the results against the numerical ones, as usually has to be done in the traditional ADMstudies to convince the readers that the results are indeed converged to the realistic solutions. Physical examplesselected from the recent application of ADM demonstrate the validity, accuracy and power of the presented novelapproach in this paper. Hence, the highly nonlinear equations arising from engineering applications can be safelytreated by the outlined method for which the classical ADM may fail or be slow to converge.
文摘The key objective of this paper is to improve the approximation of a sufficiently smooth nonperiodic function defined on a compact interval by proposing alternative forms of Fourier series expansions. Unlike in classical Fourier series, the expansion coefficients herein are explicitly dependent not only on the function itself, but also on its derivatives at the ends of the interval. Each of these series expansions can be made to converge faster at a desired polynomial rate. These results have useful implications to Fourier or harmonic analysis, solutions to differential equations and boundary value problems, data compression, and so on.
文摘In the author’s recent publications, a parametric system biorthogonal to the corresponding segment of the exponential Fourier system was unusually effective. On its basis, it was discovered that knowledge of a finite number of Fourier coefficients of function f from an infinite-dimensional set of elementary functions allows f to be accurately restored (the phenomenon of over-convergence). Below, parametric biorthogonal systems are constructed for classical trigonometric Fourier series, and the corresponding phenomena of over-convergence are discovered. The decisive role here was played by representing the space L2 as an orthogonal sum of two corresponding subspaces. As a result, fast parallel algorithms for reconstructing a function from its truncated trigonometric Fourier series are proposed. The presented numerical experiments confirm the high efficiency of these convergence accelerations for smooth functions. In conclusion, the main results of the work are summarized, and some prospects for the development and generalization of the proposed approaches are discussed.
基金supported by the National Natural Science Foundation of China(Nos.11775126,11545013,11605101)the Young Elite Scientists Sponsorship Program by CAST(No.2016QNRC001)+1 种基金Science Challenge Project by MIIT of China(No.TZ2018001)Tsinghua University,Initiative Scientific Research Program。
文摘In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.
文摘The current paper considers the problem of recovering a function using a limited number of its Fourier coefficients. Specifically, a method based on Bernoulli-like polynomials suggested and developed by Krylov, Lanczos, Gottlieb and Eckhoff is examined. Asymptotic behavior of approximate calculation of the so-called "jumps" is studied and asymptotic L2 constants of the rate of convergence of the method are computed.
基金Supported in part by grant PS 1867 from the Armenian National Science and Education Fund (ANSEF) based in New York, USA
文摘Convergence acceleration of the classical trigonometric interpolation by the Eckhoff method is considered, where the exact values of the "jumps" are approximated by solution of a system of linear equations. The accuracy of the "jump" approximation is explored and the corresponding asymptotic error of interpolation is derived. Numerical results validate theoretical estimates.
文摘In this paper, a new method, so called A-method, is given for the convergence analysis of the MQ-algorithm. And the finer relaxation parameter θA is obtained. The numerical results show that our new method has the outstanding effect of accelerating convergence. Moreover, the relaxation parameter θA is the optimum in a point of view.
文摘In this paper, we address one of the issues in the frequency assignment problem for cellular mobile networks in which we intend to minimize the interference levels when assigning frequencies from a limited frequency spectrum. In order to satisfy the increasing demand in such cellular mobile networks, we use a hybrid approach consisting of a Particle Swarm Optimization(PSO) combined with a Tabu Search(TS) algorithm. This approach takes both advantages of PSO efficiency in global optimization and TS in avoiding the premature convergence that would lead PSO to stagnate in a local minimum. Moreover, we propose a new efficient, simple, and inexpensive model for storing and evaluating solution's assignment. The purpose of this model reduces the solution's storage volume as well as the computations required to evaluate thesesolutions in comparison with the classical model. Our simulation results on the most known benchmarking instances prove the effectiveness of our proposed algorithm in comparison with previous related works in terms of convergence rate, the number of iterations, the solution storage volume and the running time required to converge to the optimal solution.
文摘The paper considers the Krylov-Lanczos and the Eckhoff approximations for recovering a bivariate function using limited number of its Fourier coefficients. These approximations are based on certain corrections associated with jumps in the partial derivatives of the approximated function. Approximation of the exact jumps is accomplished by solution of systems of linear equations along the idea of Eckhoff. Asymptotic behaviors of the approximate jumps and the Eckhoff approximation are studied. Exact constants of the asymptotic errors are computed. Numerical experiments validate theoretical investigations.
基金supported by National Natural Science Foundation of China(Grant Nos.11331008,11201469,11571358 and 11601237)the China Postdoctoral Science Foundation Funded Project(Grant Nos.2012M510186 and 2013T60761)the Hong Kong Research Grant Council(Grant No.GRF HKBU202512)
文摘In this paper, we show that the coupled modified Kd V equations possess rich mathematical structures and some remarkable properties. The connections between the system and skew orthogonal polynomials,convergence acceleration algorithms and Laurent property are discussed in detail.
基金supported by the National Natural Science Foundation of China (No. 50748033)the Specific Foundation for PhD of Hefei University of Technology (No. 2007GDBJ044), China
文摘The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of slip surface is given, and the formula-tions of FERA based on incremental tangent stiffness method and modified Aitken accelerating algorithm are developed. The limited step length iteration method (LSLIM) is adopted to calculate the reliability index. The nonlinear FERA code using the SSA technique is developed and the main flow chart is illustrated. Numerical examples are used to demonstrate the efficiency and robustness of this method. It is found that the accelerating convergence algorithm proposed in this study proves to be very efficient for it can reduce the iteration number greatly, and LSLIM is also efficient for it can assure the convergence of the iteration of the reliability index.
基金Supported by the National Natural Science Foundation of China(42064004,12062022,11762017,11762016)
文摘The Hankel transform is widely used to solve various engineering and physics problems,such as the representation of electromagnetic field components in the medium,the representation of dynamic stress intensity factors,vibration of axisymmetric infinite membrane and displacement intensity factors which all involve this type of integration.However,traditional numerical integration algorithms cannot be used due to the high oscillation characteristics of the Bessel function,so it is particularly important to propose a high precision and efficient numerical algorithm for calculating the integral of high oscillation.In this paper,the improved Gaver-Stehfest(G-S)inverse Laplace transform method for arbitrary real-order Bessel function integration is presented by using the asymptotic characteristics of the Bessel function and the accumulation of integration,and the optimized G-S coefficients are given.The effectiveness of the algorithm is verified by numerical examples.Compared with the linear transformation accelerated convergence algorithm,it shows that the G-S inverse Laplace transform method is suitable for arbitrary real order Hankel transform,and the time consumption is relatively stable and short,which provides a reliable calculation method for the study of electromagnetic mechanics,wave propagation,and fracture dynamics.