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.

Analytical Approach to Differential Equations with Piecewise Continuous Arguments via Modified Piecewise Variational Iteration Method

In the present article, we apply the modified piecewise variational iteration method to obtain the approximate analytical solutions of the differential equations with piecewise continuous arguments. This technique provides a sequence of functions which converges to the exact solution of the problem. Moreover, this method reduces the volume of calculations because it does not need discretization of the variables, linearization or small perturbations. The results seem to show that the method is very reliable and convenient for solving such equations.

The Appearance of Noise Terms in Modified Adomian Decomposition Method for Quadratic Integral Equations

In this paper, we apply the modified Adomian Decomposition Method to get the numerical solutions of Quadratic integral equations. The appearance of noise terms in Decomposition Method was investigated. The method was described along with several examples.

Numerical Solution of Generalized Abel’s Integral Equation by Variational Iteration Method

In this paper, a user friendly algorithm based on the variational iteration method (VIM) is proposed to solve singular integral equations with generalized Abel’s kernel. It is observed that an approximate solutions yn(x) converges to the exact solution irrespective of the initial choice y0 (x). Illustrative numerical examples are given to demonstrate the efficiency and simplicity of the method in solving these types of singular integral equations.

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.

Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

Iterative convergence of boundary-volume integral equation method

The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.

Modified Efficient Families of Two and Three-Step Predictor-Corrector Iterative Methods for Solving Nonlinear Equations

In this paper, we present and analyze modified families of predictor-corrector iterative methods for finding simple zeros of univariate nonlinear equations, permitting near the root. The main advantage of our methods is that they perform better and moreover, have the same efficiency indices as that of existing multipoint iterative methods. Furthermore, the convergence analysis of the new methods is discussed and several examples are given to illustrate their efficiency.

The Alternating Group Explicit Iterative Method for the Regularized Long-Wave Equation

An Alternating Group Explicit (AGE) iterative method with intrinsic parallelism is constructed based on an implicit scheme for the Regularized Long-Wave (RLW) equation. The method can be used for the iteration solution of a general tridiagonal system of equations with diagonal dominance. It is not only easy to implement, but also can directly carry out parallel computation. Convergence results are obtained by analysing the linear system. Numerical experiments show that the theory is accurate and the scheme is valid and reliable.

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.

Solving mKdV-sinh-Gordon equation by a modified variable separated ordinary differential equation method

By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ODE method is developed for solving the mKdV-sinh-Gordon equation. As a result, many explicit and exact solutions including some new formal solutions are successfully picked up for the mKdV-sinh-Gordon equation by this approach.

ON MONOTONE ITERATIVE METHOD FOR INITIAL VALUE PROBLEMS OF NONLINEAR SECOND-ORDER INTEGRO-DIFFERENTIAL EQUATIONS IN BANACH SPACES

Using the monotone iterative method and Monch Fixed point theorem, the existence of solutions and coupled minimal and maximal quasisolutions of initial value problems for mixed monotone second-order integro-differential equations in Banach spaces are studied. Some existence theorems of solutions and coupled minimal and maximal quasisolutions are obtained.

Modified Burgers' equation by the local discontinuous Galerkin method

In this paper,we present the local discontinuous Galerkin method for solving Burgers＇ equation and the modified Burgers＇ equation.We describe the algorithm formulation and practical implementation of the local discontinuous Galerkin method in detail.The method is applied to the solution of the one-dimensional viscous Burgers＇ equation and two forms of the modified Burgers＇ equation.The numerical results indicate that the method is very accurate and efficient.

OPTIMIZATION OF ADAPTIVE DIRECT METHOD FOR APPROXIMATE SOLUTION OF INTEGRAL EQUATIONS OF SEVERAL VARIABLES

This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.

A Method for Solving Fredholm Integral Equations of the First Kind Based on Chebyshev Wavelets

In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.

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 modified Tikhonov regularization method for a Cauchy problem of a time fractional diffusion equation

In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.

A Fourth-Order Modified Method for the Cauchy Problem of the Modified Helmholtz Equation

This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0<x≤1,y∈R.The Cauchy data at x = 0 is given and the solution is then sought for the interval 0<x≤1.This problem is highly ill-posed and the solution(if it exists) does not depend continuously on the given data. In this paper,we propose a fourth-order modified method to solve the Cauchy problem. Convergence estimates are presented under the suitable choices of regularization parameters and the a priori assumption on the bounds of the exact solution.Numerical implementation is considered and the numerical examples show that our proposed method is effective and stable.

Exact Traveling Wave Solutions for the System of Shallow Water Wave Equations and Modified Liouville Equation Using Extended Jacobian Elliptic Function Expansion Method

In this work, an extended Jacobian elliptic function expansion method is proposed for constructing the exact solutions of nonlinear evolution equations. The validity and reliability of the method are tested by its applications to the system of shallow water wave equations and modified Liouville equation which play an important role in mathematical physics.