Fast and Accurate Predictor-Corrector Methods Using Feedback-Accelerated Picard Iteration for Strongly Nonlinear Problems

Although predictor-corrector methods have been extensively applied,they might not meet the requirements of practical applications and engineering tasks,particularly when high accuracy and efficiency are necessary.A novel class of correctors based on feedback-accelerated Picard iteration(FAPI)is proposed to further enhance computational performance.With optimal feedback terms that do not require inversion of matrices,significantly faster convergence speed and higher numerical accuracy are achieved by these correctors compared with their counterparts;however,the computational complexities are comparably low.These advantages enable nonlinear engineering problems to be solved quickly and accurately,even with rough initial guesses from elementary predictors.The proposed method offers flexibility,enabling the use of the generated correctors for either bulk processing of collocation nodes in a domain or successive corrections of a single node in a finite difference approach.In our method,the functional formulas of FAPI are discretized into numerical forms using the collocation approach.These collocated iteration formulas can directly solve nonlinear problems,but they may require significant computational resources because of the manipulation of high-dimensionalmatrices.To address this,the collocated iteration formulas are further converted into finite difference forms,enabling the design of lightweight predictor-corrector algorithms for real-time computation.The generality of the proposed method is illustrated by deriving new correctors for three commonly employed finite-difference approaches:the modified Euler approach,the Adams-Bashforth-Moulton approach,and the implicit Runge-Kutta approach.Subsequently,the updated approaches are tested in solving strongly nonlinear problems,including the Matthieu equation,the Duffing equation,and the low-earth-orbit tracking problem.The numerical findings confirm the computational accuracy and efficiency of the derived predictor-corrector algorithms.

Optimization of Random Feature Method in the High-Precision Regime

Machine learning has been widely used for solving partial differential equations(PDEs)in recent years,among which the random feature method(RFM)exhibits spectral accuracy and can compete with traditional solvers in terms of both accuracy and efficiency.Potentially,the optimization problem in the RFM is more difficult to solve than those that arise in traditional methods.Unlike the broader machine-learning research,which frequently targets tasks within the low-precision regime,our study focuses on the high-precision regime crucial for solving PDEs.In this work,we study this problem from the following aspects:(i)we analyze the coeffcient matrix that arises in the RFM by studying the distribution of singular values;(ii)we investigate whether the continuous training causes the overfitting issue;(ii)we test direct and iterative methods as well as randomized methods for solving the optimization problem.Based on these results,we find that direct methods are superior to other methods if memory is not an issue,while iterative methods typically have low accuracy and can be improved by preconditioning to some extent.

Another SSOR Iteration Method

Kellogg gave a version of the Peaceman-Radford method. In this paper, we introduce a SSOR iteration method which uses Kellogg’s method. The new algorithm has some advantages over the traditional SSOR algorithm. A Cyclic Reduction algorithm is introduced via a decoupling in Kellogg’s method.

Variational Iteration Method for Solving Time Fractional Burgers Equation Using Maple

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.

Variational iteration method for solving the mechanism of the Equatorial Eastern Pacific El Nino-Southern Oscillation

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.

Variational iteration solving method for El Nio phenomenon atmospheric physics of nonlinear model

A class of E1 Niйo atmospheric physics oscillation model is considered. The E1 Niйo atmospheric physics oscillation is an abnormal phenomenon involved in the tropical Pacific ocean-atmosphere interactions. The conceptual oscillator model should consider the variations of both the eastern and western Pacific anomaly patterns. An E1 Niйo atmospheric physics model is proposed using a method for the variational iteration theory. Using the variational iteration method, the approximate expansions of the solution of corresponding problem are constructed. That is, firstly, introducing a set of functional and accounting their variationals, the Lagrange multiplicators are counted, and then the variational iteration is defined, finally, the approximate solution is obtained. From approximate expansions of the solution, the zonal sea surface temperature anomaly in the equatorial eastern Pacific and the thermocline depth anomaly of the sea-air oscillation for E1 Niйo atmospheric physics model can be analyzed. E1 Niйo is a very complicated natural phenomenon. Hence basic models need to be reduced for the sea-air oscillator and are solved. The variational iteration is a simple and valid approximate method.

ON THE MONOTONE CONVERGENCE OF THE PROJECTED ITERATION METHODS FOR LINEAR COMPLEMENTARITY PROBLEMS

Under suitable conditions,the monotone convergence about the projected iteration method for solving linear complementarity problem is proved and the influence of the involved parameter matrix on the convergence rate of this method is investigated.

Doubly Periodic Wave Solutions of Jaulent-Miodek Equations Using Variational Iteration Method Combined with Jacobian-function Method

One of the advantages of the variational iteration method is the free choice of initial guess. In this paper we use the basic idea of the Jacobian-function method to construct a generalized trial function with some unknown parameters. The Jaulent-Miodek equations are used to illustrate effectiveness and convenience of this method, some new explicit exact travelling wave solutions have been obtained, which include bell-type soliton solution, kink-type soliton solutions, solitary wave solutions, and doubly periodic wave solutions.

A Unification of the Concepts of the Variational Iteration, Adomian Decomposition and Picard Iteration Methods;and a Local Variational Iteration Method

This paper compares the variational iteration method(VIM),the Adomian decomposition method(ADM)and the Picard iteration method(PIM)for solving a system of first o rder n onlinear o rdinary d ifferential e quations(ODEs).A unification of the concepts underlying these three methods is attempted by considering a very general iterative algorithm for VIM.It is found that all the three methods can be regarded as special cases of using a very general matrix of Lagrange multipliers in the iterative algorithm of VIM.The global variational iteration method is briefly reviewed,and further recast into a Local VIM,which is much more convenient and capable of predicting long term complex dynamic responses of nonlinear systems even if they are chaotic.

Solutions of the Duffin-Kemmer-Petiau equation in the presence of Hulthn potential in(1+2) dimensions for unity spin particles using the asymptotic iteration method

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.

Methods of Variations and Their Applications

In this review article, we begin with reviewing Calculus of variations giving few examples on its use to solve a large number of problems in geometry, physics, and other branches of knowledge. Afterwards, we direct our attention to different methods of variations which evolved during the last century and which include their use in eigenvalue problems and in finite difference methods and those adopted in classical and quantum mechanics. The methods used in evaluating products and quotients of functionals are also discussed along with variational iteration methods. Later on, a good number of applications in different areas are presented and discussed;then a concluding discussion is given.

A New Multidimensional Time Series Forecasting Method Based on the EOF Iteration Scheme

In this paper a new .mnultidimensional time series forecasting scheme based on the empirical orthogonal function (EOF) stepwise iteration process is introduced. The scheme is tested in a series of forecast experiments of Nino3 SST anomalies and Tahiti-Darwin SO index. The results show that the scheme is feasible and ENSO predictable.

Multilevel Iteration Methods for Solving Linear Operator Equations of the First Kind

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.

Multistage Numerical Picard Iteration Methods for Nonlinear Volterra Integral Equations of the Second Kind

Using the Picard iteration method and treating the involved integration by numerical quadrature formulas, we propose a numerical scheme for the second kind nonlinear Volterra integral equations. For enlarging the convergence region of the Picard iteration method, multistage algorithm is devised. We also introduce an algorithm for problems with some singularities at the limits of integration including fractional integral equations. Numerical tests verify the validity of the proposed schemes.

Sign-Changing Solutions for Superlinear Kirchhoff Type Problem via the Nehari Method

In this paper, we consider a class of Kirchhoff type problem with superlinear nonlinearity. A sign-changing solution with exactly two nodal domains will be obtained by combining the Nehari method and an iterative technique.

Integration of the Coupled Orbit-Attitude Dynamics Using Modified Chebyshev-Picard Iteration Methods

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.

On Minimization of Upper Bound for the Convergence Rate of the QHSS Iteration Method

For an upper bound of the spectral radius of the QHSS (quasi Hermitian and skew-Hermitian splitting) iteration matrix which can also bound the contraction factor of the QHSS iteration method,we give its minimum point under the conditions which guarantee that the upper bound is strictly less than one. This provides a good choice of the involved iteration parameters,so that the convergence rate of the QHSS iteration method can be significantly improved.

A Study of Some Nonlinear Partial Differential Equations by Using Adomian Decomposition Method and Variational Iteration Method

In this paper, a numerical solution of nonlinear partial differential equation, Benjamin-Bona-Mahony (BBM) and Cahn-Hilliard equation is presented by using Adomain Decomposition Method (ADM) and Variational Iteration Method (VIM). The results reveal that the two methods are very effective, simple and very close to the exact solution.

On the Application of Adomian Decomposition Method to Special Equations in Physical Sciences

The current manuscript makes use of the prominent iterative procedure, called the Adomian Decomposition Method (ADM), to tackle some important special differential equations. The equations of curiosity in this study are the singular equations that arise in many physical science applications. Thus, through the application of the ADM, a generalized recursive scheme was successfully derived and further utilized to obtain closed-form solutions for the models under consideration. The method is, indeed, fascinating as respective exact analytical solutions are accurately acquired with only a small number of iterations.

Existence and numerical approximation of a solution to frictional contact problem for electro-elastic materials

In this paper,a frictional contact problem between an electro-elastic body and an electrically conductive foundation is studied.The contact is modeled by normal compliance with finite penetration and a version of Coulomb’s law of dry friction in which the coefficient of friction depends on the slip.In addition,the effects of the electrical conductivity of the foundation are taken into account.This model leads to a coupled system of the quasi-variational inequality of the elliptic type for the displacement and the nonlinear variational equation for the electric potential.The existence of a weak solution is proved by using an abstract result for elliptic variational inequalities and a fixed point argument.Then,a finite element approximation of the problem is presented.Under some regularity conditions,an optimal order error estimate of the approximate solution is derived.Finally,a successive iteration technique is used to solve the problem numerically and a convergence result is established.