The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this ...The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.展开更多
A new algorithm called homotopy iteration method based on the homotopy function is studied and improved. By the improved homotopy iteration method, Polynomial systems with high Order and deficient can be solved fast a...A new algorithm called homotopy iteration method based on the homotopy function is studied and improved. By the improved homotopy iteration method, Polynomial systems with high Order and deficient can be solved fast and efficiently comparing to the original homotopy iteration method. Numerical examples for the ninepoint path synthesis of four-bar linkages show the advantages and efficiency of the improved homotopy iteration method.展开更多
In this paper we develop two multilevel iteration methods for solving linear systems resulting from the Galerkin method and Tikhonov regularization for linear ill-posed problems. The two algorithms and their convergen...In this paper we develop two multilevel iteration methods for solving linear systems resulting from the Galerkin method and Tikhonov regularization for linear ill-posed problems. The two algorithms and their convergence analyses are presented in an abstract framework.展开更多
The relativistic Duffin-Kemmer-Petiau equation in the presence of Hulthen potential in (1 +2) dimensions for spin-one particles is studied. Hence, the asymptotic iteration method is used for obtaining energy eigenv...The relativistic Duffin-Kemmer-Petiau equation in the presence of Hulthen potential in (1 +2) dimensions for spin-one particles is studied. Hence, the asymptotic iteration method is used for obtaining energy eigenvalues and eigenfunctions.展开更多
Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the f...Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.展开更多
The asymptotic iteration method (AIM) is used to obtain the quasi-exact solutions of the Schr6dinger equation with a deformed well potential. For arbitrary potential parameters, a numerical aspect of AIM is also app...The asymptotic iteration method (AIM) is used to obtain the quasi-exact solutions of the Schr6dinger equation with a deformed well potential. For arbitrary potential parameters, a numerical aspect of AIM is also applied to obtain highly accurate energy eigenvalues. Additionally, the perturbation expansion, based on the AIM approach, is utilized to obtain simple analytic expressions for the energy eigenvalues.展开更多
This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification o...This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.展开更多
The variational iteration method is successfully extended to the case of solving fractional differential equations, and the Lagrange multiplier of the method is identified in a more accurate way. Some diffusion models...The variational iteration method is successfully extended to the case of solving fractional differential equations, and the Lagrange multiplier of the method is identified in a more accurate way. Some diffusion models with fractional derivatives are investigated analytically, and the results show the efficiency of the new Lagrange multiplier for fractional differential equations of arbitrary order.展开更多
Air exploratory discussion of an ancient Chinese algorithm, the Ying Buzu Shu, in about 2nd century BC, known as the rule of double false position in the West is given. In addition to pointing out that the rule of dou...Air exploratory discussion of an ancient Chinese algorithm, the Ying Buzu Shu, in about 2nd century BC, known as the rule of double false position in the West is given. In addition to pointing out that the rule of double false position is actually a translation version of the ancient Chinese algorithm, a comparison with well-known Newton iteration method is also made. If derivative is introduced, the ancient Chinese algorithm reduces to the Newton method. A modification of the ancient Chinese algorithm is also proposed, and some of applications to nonlinear oscillators are illustrated.展开更多
This paper studies a delayed air-sea coupled oscillator describing the physical mechanism of El Nino Southern Oscillation. The approximate expansions of the delayed differential equation's solution are obtained succe...This paper studies a delayed air-sea coupled oscillator describing the physical mechanism of El Nino Southern Oscillation. The approximate expansions of the delayed differential equation's solution are obtained successfully by the modified variational iteration method. The numerical results illustrate the effectiveness and correctness of the method by comparing with the exact solution of the reduced model.展开更多
The purpose is to accurately predict the performance of foil bearing and achieve accurate results in the design of foil bearing structure.A new type of foil bearing with surface microstructure is used as experimental ...The purpose is to accurately predict the performance of foil bearing and achieve accurate results in the design of foil bearing structure.A new type of foil bearing with surface microstructure is used as experimental material.First,the lubrication mechanism of elastic foil gas bearing is analyzed.Then,the numerical solution process of the static bearing capacity and friction torque is analyzed,including the discretization of the governing equation of rarefied gas pressure based on the non-dimensional modified Reynolds equation and the over relaxation iteration method,the grid planning within the calculation range,the static solution of boundary parameters and static solution of the numerical process.Finally,the solution program is analyzed.The experimental data in National Aeronautics and Space Administration(NASA)public literature are compared with the simulation results of this exploration,so as to judge the accuracy of the calculation process.The results show that under the same static load,the difference between the minimum film thickness calculated and the test results is not obvious;when the rotor speed of the bearing is 60000 r/min,the influence of the boundary slip effect increases with the increase of the micro groove depth on the flat foil surface;when the eccentricity or the micro groove depth of the bearing increases,the bearing capacity will be strengthened.When the eccentricity is 6µm and 14µm,the viscous friction torque of the new foil bearing increases significantly with the increase of the depth of the foil micro groove,but when the eccentricity is 22µm,the viscous friction torque does not change with the change of the depth of the foil micro groove.It shows that the bearing capacity and performance of foil bearing are improved.展开更多
In this paper, a new kind of iteration technique for solving nonlinear ordinary differential equations is described and used to give approximate periodic solutions for some well-known nonlinear problems. The most inte...In this paper, a new kind of iteration technique for solving nonlinear ordinary differential equations is described and used to give approximate periodic solutions for some well-known nonlinear problems. The most interesting features of the proposed methods are its extreme simplicity and concise forms of iteration formula for a wide range of nonlinear problems.展开更多
In this study, we applied the variational iteration method to solve the Boussinesq time equation. Bossiness’s article from 1872 introduced the equations that are now known as the Boussinesq equations. Numerical metho...In this study, we applied the variational iteration method to solve the Boussinesq time equation. Bossiness’s article from 1872 introduced the equations that are now known as the Boussinesq equations. Numerical methods are commonly utilized to solve nonlinear equation systems. Several research papers have documented the values of the variational iteration method and its applications for various categories of differential equations. A comparison of the exact and numerical solutions was obtained using the variational iteration method. The variational iteration method shows that the proposed method is very effective and convenient. The results are shown for different specific cases of the problem. The variational iteration method is useful in numerical simulations and approximate analytical solutions, and it is used to resolve nonlinear differential equations in various situations using Maple. For example, the linear Boussinesq equation was resolved using the variational iteration method. By comparing the numerical results, we found that the variable repetition method produced accurate results and was close to the exact solution, allowing it to be widely applied to the Boussinesq equation. This proves the effectiveness of the method and the capability to quickly and effectively obtain the numerical number solution related to the exact solution using the Maple 18 program. Additionally, the outcomes are extremely precise.展开更多
A model for both stochastic jumps and volatility for equity returns in the area of option pricing is the stochastic volatility process with jumps (SVPJ). A major advantage of this model lies in the area of mean revers...A model for both stochastic jumps and volatility for equity returns in the area of option pricing is the stochastic volatility process with jumps (SVPJ). A major advantage of this model lies in the area of mean reversion and volatility clustering between returns and volatility with uphill movements in price asserts. Thus, in this article, we propose to solve the SVPJ model numerically through a discretized variational iteration method (DVIM) to obtain sample paths for the state variable and variance process at various timesteps and replications in order to estimate the expected jump times at various iterates resulting from executing the DVIM as n increases. These jumps help in estimating the degree of randomness in the financial market. It was observed that the average computed expected jump times for the state variable and variance process is moderated by the parameters (variance process through mean reversion), Θ (long-run mean of the variance process), σ (volatility variance process) and λ (constant intensity of the Poisson process) at each iterate. For instance, when = 0.0, Θ = 0.0, σ = 0.0 and λ = 1.0, the state variable cluttered maximally compared to the variance process with less volatility cluttering with an average computed expected jump times of 52.40607869 as n increases in the DVIM scheme. Similarly, when = 3.99, Θ = 0.014, σ = 0.27 and λ = 0.11, the stochastic jumps for the state variable are less cluttered compared to the variance process with maximum volatility cluttering as n increases in the DVIM scheme. In terms of option pricing, the value 52.40607869 suggest a better bargain compared to the value 20.40344029 due to the fact that it yields less volatility rate. MAPLE 18 software was used for all computations in this research.展开更多
A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO mode...A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.展开更多
In the designing of the flexible fender pile, M-method is usually adopted. But M-method is usable only in the case that the value of the force (moment M and shear force Q) exerted on the pile above the soil is known. ...In the designing of the flexible fender pile, M-method is usually adopted. But M-method is usable only in the case that the value of the force (moment M and shear force Q) exerted on the pile above the soil is known. In practical cases, only the percussive energy of the ship can be found, while the exact value of the percussive force cannot be computed directly. In this paper, a method is introduced which can be used to calculate the inner force of flexible fender piles by a computer in these cases. First, according to the known docking dynamic energy of a ship, the absorbed energy can be found from the deformation of the pile and the rubber fender respectively, at the same time the value of the force acting on the fender piles and the reaction of fenders can be calculated. Then the inner force and displacement of the fender piles can also be calculated further by M-method. The whole computing process is realized by a computer program and the results obtained are reasonable. In this paper, the computation, the iteration method and the computer program are all expounded in detail. This method can be applied to the design of fender piles and type selection of rubber fenders in practical projects.展开更多
In this paper, the modified iteration method is further generalized to the study of axisymmetrical postbuckling of thin circular plates and hereby a new approximate analytic solution of the problem is obtained. Furthe...In this paper, the modified iteration method is further generalized to the study of axisymmetrical postbuckling of thin circular plates and hereby a new approximate analytic solution of the problem is obtained. Further utilizations of this method to postbuckling analyses of plates of more complicated structure are expected.展开更多
This paper presents Modified Chebyshev-Picard Iteration(MCPI)methods for long-term integration of the coupled orbit and attitude dynamics.Although most orbit predictions for operational satellites have assumed that th...This paper presents Modified Chebyshev-Picard Iteration(MCPI)methods for long-term integration of the coupled orbit and attitude dynamics.Although most orbit predictions for operational satellites have assumed that the attitude dynamics is decoupled from the orbit dynamics,the fully coupled dynamics is required for the solutions of uncontrolled space debris and space objects with high area-to-mass ratio,for which cross sectional area is constantly changing leading to significant change on the solar radiation pressure and atmospheric drag.MCPI is a set of methods for solution of initial value problems and boundary value problems.The methods refine an orthogonal function approximation of long-time-interval segments of state trajectories iteratively by fusing Chebyshev polynomials with the classical Picard iteration and have been applied to multiple challenging aerospace problems.Through the studies on integrating a torque-free rigid body rotation and a long-term integration of the coupled orbit-attitude dynamics through the effect of solar radiation pressure,MCPI methods are shown to achieve several times speedup over the Runge-Kutta 7(8)methods with several orders of magnitudes of better accuracy.MCPI methods are further optimized by integrating the decoupled dynamics at the beginning of the iteration and coupling the full dynamics when the attitude solutions and orbit solutions are converging during the iteration.The approach of decoupling and then coupling during iterations provides a unique and promising perspective on the way to warm start the solution process for the longterm integration of the coupled orbit-attitude dynamics.Furthermore,an attractive feature of MCPI in maintaining the unity constraint for the integration of quaternions within machine accuracy is illustrated to be very appealing.展开更多
For over half a century,numerical integration methods based on finite difference,such as the Runge-Kutta method and the Euler method,have been popular and widely used for solving orbit dynamic problems.In general,a sm...For over half a century,numerical integration methods based on finite difference,such as the Runge-Kutta method and the Euler method,have been popular and widely used for solving orbit dynamic problems.In general,a small integration step size is always required to suppress the increase of the accumulated computation error,which leads to a relatively slow computation speed.Recently,a collocation iteration method,approximating the solutions of orbit dynamic problems iteratively,has been developed.This method achieves high computation accuracy with extremely large step size.Although efficient,the collocation iteration method suffers from two limitations:(A)the computational error limit of the approximate solution is not clear;(B)extensive trials and errors are always required in tuning parameters.To overcome these problems,the influence mechanism of how the dynamic problems and parameters affect the error limit of the collocation iteration method is explored.On this basis,a parameter adjustment method known as the“polishing method”is proposed to improve the computation speed.The method proposed is demonstrated in three typical orbit dynamic problems in aerospace engineering:a low Earth orbit propagation problem,a Molniya orbit propagation problem,and a geostationary orbit propagation problem.Numerical simulations show that the proposed polishing method is faster and more accurate than the finite-difference-based method and the most advanced collocation iteration method.展开更多
In this paper, a modified Hermitian and skew-Hermitian splitting (MHSS) iteration method for solving the complex linear matrix equation AXB = C has been presented. As the theoretical analysis shows, the MHSS iterati...In this paper, a modified Hermitian and skew-Hermitian splitting (MHSS) iteration method for solving the complex linear matrix equation AXB = C has been presented. As the theoretical analysis shows, the MHSS iteration method will converge under cer- tain conditions. Each iteration in this method requires the solution of four linear matrix equations with real symmetric positive definite coefficient matrices, although the original coefficient matrices are complex and non-Hermitian. In addition, the optimal parameter of the new iteration method is proposed. Numerical results show that MHSS iteration method is efficient and robust.展开更多
文摘The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.
文摘A new algorithm called homotopy iteration method based on the homotopy function is studied and improved. By the improved homotopy iteration method, Polynomial systems with high Order and deficient can be solved fast and efficiently comparing to the original homotopy iteration method. Numerical examples for the ninepoint path synthesis of four-bar linkages show the advantages and efficiency of the improved homotopy iteration method.
基金The NSF(0611005)of Jiangxi Province and the SF(2007293)of Jiangxi Provincial Education Department.
文摘In this paper we develop two multilevel iteration methods for solving linear systems resulting from the Galerkin method and Tikhonov regularization for linear ill-posed problems. The two algorithms and their convergence analyses are presented in an abstract framework.
文摘The relativistic Duffin-Kemmer-Petiau equation in the presence of Hulthen potential in (1 +2) dimensions for spin-one particles is studied. Hence, the asymptotic iteration method is used for obtaining energy eigenvalues and eigenfunctions.
文摘Using the asymptotic iteration method, we obtain the S-wave solution for a short-range three-parameter central potential with 1/r singularity and with a non-orbital barrier. To the best of our knowledge, this is the first attempt at calculating the energy spectrum for this potential, which was introduced by H. Bahlouli and A. D. Alhaidari and for which they obtained the “potential parameter spectrum”. Our results are also independently verified using a direct method of diagonalizing the Hamiltonian matrix in the J-matrix basis.
文摘The asymptotic iteration method (AIM) is used to obtain the quasi-exact solutions of the Schr6dinger equation with a deformed well potential. For arbitrary potential parameters, a numerical aspect of AIM is also applied to obtain highly accurate energy eigenvalues. Additionally, the perturbation expansion, based on the AIM approach, is utilized to obtain simple analytic expressions for the energy eigenvalues.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10771019 and 10826107)
文摘This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.
基金Project supported by the Key Program of the National Natural Science Foundation of China (Grant No. 51134018).
文摘The variational iteration method is successfully extended to the case of solving fractional differential equations, and the Lagrange multiplier of the method is identified in a more accurate way. Some diffusion models with fractional derivatives are investigated analytically, and the results show the efficiency of the new Lagrange multiplier for fractional differential equations of arbitrary order.
文摘Air exploratory discussion of an ancient Chinese algorithm, the Ying Buzu Shu, in about 2nd century BC, known as the rule of double false position in the West is given. In addition to pointing out that the rule of double false position is actually a translation version of the ancient Chinese algorithm, a comparison with well-known Newton iteration method is also made. If derivative is introduced, the ancient Chinese algorithm reduces to the Newton method. A modification of the ancient Chinese algorithm is also proposed, and some of applications to nonlinear oscillators are illustrated.
基金supported by the National Natural Science Foundation of China (Grant Nos.41105063 and 61070041)
文摘This paper studies a delayed air-sea coupled oscillator describing the physical mechanism of El Nino Southern Oscillation. The approximate expansions of the delayed differential equation's solution are obtained successfully by the modified variational iteration method. The numerical results illustrate the effectiveness and correctness of the method by comparing with the exact solution of the reduced model.
文摘The purpose is to accurately predict the performance of foil bearing and achieve accurate results in the design of foil bearing structure.A new type of foil bearing with surface microstructure is used as experimental material.First,the lubrication mechanism of elastic foil gas bearing is analyzed.Then,the numerical solution process of the static bearing capacity and friction torque is analyzed,including the discretization of the governing equation of rarefied gas pressure based on the non-dimensional modified Reynolds equation and the over relaxation iteration method,the grid planning within the calculation range,the static solution of boundary parameters and static solution of the numerical process.Finally,the solution program is analyzed.The experimental data in National Aeronautics and Space Administration(NASA)public literature are compared with the simulation results of this exploration,so as to judge the accuracy of the calculation process.The results show that under the same static load,the difference between the minimum film thickness calculated and the test results is not obvious;when the rotor speed of the bearing is 60000 r/min,the influence of the boundary slip effect increases with the increase of the micro groove depth on the flat foil surface;when the eccentricity or the micro groove depth of the bearing increases,the bearing capacity will be strengthened.When the eccentricity is 6µm and 14µm,the viscous friction torque of the new foil bearing increases significantly with the increase of the depth of the foil micro groove,but when the eccentricity is 22µm,the viscous friction torque does not change with the change of the depth of the foil micro groove.It shows that the bearing capacity and performance of foil bearing are improved.
文摘In this paper, a new kind of iteration technique for solving nonlinear ordinary differential equations is described and used to give approximate periodic solutions for some well-known nonlinear problems. The most interesting features of the proposed methods are its extreme simplicity and concise forms of iteration formula for a wide range of nonlinear problems.
文摘In this study, we applied the variational iteration method to solve the Boussinesq time equation. Bossiness’s article from 1872 introduced the equations that are now known as the Boussinesq equations. Numerical methods are commonly utilized to solve nonlinear equation systems. Several research papers have documented the values of the variational iteration method and its applications for various categories of differential equations. A comparison of the exact and numerical solutions was obtained using the variational iteration method. The variational iteration method shows that the proposed method is very effective and convenient. The results are shown for different specific cases of the problem. The variational iteration method is useful in numerical simulations and approximate analytical solutions, and it is used to resolve nonlinear differential equations in various situations using Maple. For example, the linear Boussinesq equation was resolved using the variational iteration method. By comparing the numerical results, we found that the variable repetition method produced accurate results and was close to the exact solution, allowing it to be widely applied to the Boussinesq equation. This proves the effectiveness of the method and the capability to quickly and effectively obtain the numerical number solution related to the exact solution using the Maple 18 program. Additionally, the outcomes are extremely precise.
文摘A model for both stochastic jumps and volatility for equity returns in the area of option pricing is the stochastic volatility process with jumps (SVPJ). A major advantage of this model lies in the area of mean reversion and volatility clustering between returns and volatility with uphill movements in price asserts. Thus, in this article, we propose to solve the SVPJ model numerically through a discretized variational iteration method (DVIM) to obtain sample paths for the state variable and variance process at various timesteps and replications in order to estimate the expected jump times at various iterates resulting from executing the DVIM as n increases. These jumps help in estimating the degree of randomness in the financial market. It was observed that the average computed expected jump times for the state variable and variance process is moderated by the parameters (variance process through mean reversion), Θ (long-run mean of the variance process), σ (volatility variance process) and λ (constant intensity of the Poisson process) at each iterate. For instance, when = 0.0, Θ = 0.0, σ = 0.0 and λ = 1.0, the state variable cluttered maximally compared to the variance process with less volatility cluttering with an average computed expected jump times of 52.40607869 as n increases in the DVIM scheme. Similarly, when = 3.99, Θ = 0.014, σ = 0.27 and λ = 0.11, the stochastic jumps for the state variable are less cluttered compared to the variance process with maximum volatility cluttering as n increases in the DVIM scheme. In terms of option pricing, the value 52.40607869 suggest a better bargain compared to the value 20.40344029 due to the fact that it yields less volatility rate. MAPLE 18 software was used for all computations in this research.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 90111011 and 10471039), the National Key Basic Research Special Foundation of China (Grant Nos 2003CB415101-03 and 2004CB418304), the Key Basic Research Foundation of the Chinese Academy of Sciences (Grant No KZCX3-SW-221) and in part by E-Institutes of Shanghai Municipal Education Commission (Grant No N.E03004).
文摘A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.
文摘In the designing of the flexible fender pile, M-method is usually adopted. But M-method is usable only in the case that the value of the force (moment M and shear force Q) exerted on the pile above the soil is known. In practical cases, only the percussive energy of the ship can be found, while the exact value of the percussive force cannot be computed directly. In this paper, a method is introduced which can be used to calculate the inner force of flexible fender piles by a computer in these cases. First, according to the known docking dynamic energy of a ship, the absorbed energy can be found from the deformation of the pile and the rubber fender respectively, at the same time the value of the force acting on the fender piles and the reaction of fenders can be calculated. Then the inner force and displacement of the fender piles can also be calculated further by M-method. The whole computing process is realized by a computer program and the results obtained are reasonable. In this paper, the computation, the iteration method and the computer program are all expounded in detail. This method can be applied to the design of fender piles and type selection of rubber fenders in practical projects.
文摘In this paper, the modified iteration method is further generalized to the study of axisymmetrical postbuckling of thin circular plates and hereby a new approximate analytic solution of the problem is obtained. Further utilizations of this method to postbuckling analyses of plates of more complicated structure are expected.
文摘This paper presents Modified Chebyshev-Picard Iteration(MCPI)methods for long-term integration of the coupled orbit and attitude dynamics.Although most orbit predictions for operational satellites have assumed that the attitude dynamics is decoupled from the orbit dynamics,the fully coupled dynamics is required for the solutions of uncontrolled space debris and space objects with high area-to-mass ratio,for which cross sectional area is constantly changing leading to significant change on the solar radiation pressure and atmospheric drag.MCPI is a set of methods for solution of initial value problems and boundary value problems.The methods refine an orthogonal function approximation of long-time-interval segments of state trajectories iteratively by fusing Chebyshev polynomials with the classical Picard iteration and have been applied to multiple challenging aerospace problems.Through the studies on integrating a torque-free rigid body rotation and a long-term integration of the coupled orbit-attitude dynamics through the effect of solar radiation pressure,MCPI methods are shown to achieve several times speedup over the Runge-Kutta 7(8)methods with several orders of magnitudes of better accuracy.MCPI methods are further optimized by integrating the decoupled dynamics at the beginning of the iteration and coupling the full dynamics when the attitude solutions and orbit solutions are converging during the iteration.The approach of decoupling and then coupling during iterations provides a unique and promising perspective on the way to warm start the solution process for the longterm integration of the coupled orbit-attitude dynamics.Furthermore,an attractive feature of MCPI in maintaining the unity constraint for the integration of quaternions within machine accuracy is illustrated to be very appealing.
基金This study was co-supported by the National Key Research and Development Program of China(No.2021YFA0717100)the National Natural Science Foundation of China(Nos.12072270,U2013206).
文摘For over half a century,numerical integration methods based on finite difference,such as the Runge-Kutta method and the Euler method,have been popular and widely used for solving orbit dynamic problems.In general,a small integration step size is always required to suppress the increase of the accumulated computation error,which leads to a relatively slow computation speed.Recently,a collocation iteration method,approximating the solutions of orbit dynamic problems iteratively,has been developed.This method achieves high computation accuracy with extremely large step size.Although efficient,the collocation iteration method suffers from two limitations:(A)the computational error limit of the approximate solution is not clear;(B)extensive trials and errors are always required in tuning parameters.To overcome these problems,the influence mechanism of how the dynamic problems and parameters affect the error limit of the collocation iteration method is explored.On this basis,a parameter adjustment method known as the“polishing method”is proposed to improve the computation speed.The method proposed is demonstrated in three typical orbit dynamic problems in aerospace engineering:a low Earth orbit propagation problem,a Molniya orbit propagation problem,and a geostationary orbit propagation problem.Numerical simulations show that the proposed polishing method is faster and more accurate than the finite-difference-based method and the most advanced collocation iteration method.
文摘In this paper, a modified Hermitian and skew-Hermitian splitting (MHSS) iteration method for solving the complex linear matrix equation AXB = C has been presented. As the theoretical analysis shows, the MHSS iteration method will converge under cer- tain conditions. Each iteration in this method requires the solution of four linear matrix equations with real symmetric positive definite coefficient matrices, although the original coefficient matrices are complex and non-Hermitian. In addition, the optimal parameter of the new iteration method is proposed. Numerical results show that MHSS iteration method is efficient and robust.