A finite volume element predictor-corrector method for a class of nonlinear parabolic system of equations is presented and analyzed. Suboptimal L^2 error estimate for the finite volume element predictor-corrector meth...A finite volume element predictor-corrector method for a class of nonlinear parabolic system of equations is presented and analyzed. Suboptimal L^2 error estimate for the finite volume element predictor-corrector method is derived. A numerical experiment shows that the numerical results are consistent with theoretical analysis.展开更多
The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent ...The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.展开更多
The paper introduces a new class of numerical schemes for the approximate solutions of stochastic pantograph equations. As an effective technique to implement implicit stochastic methods, strong predictor-corrector me...The paper introduces a new class of numerical schemes for the approximate solutions of stochastic pantograph equations. As an effective technique to implement implicit stochastic methods, strong predictor-corrector methods （PCMs） are designed to handle scenario simulation of solutions of stochastic pantograph equations. It is proved that the PCMs are strong convergent with order 1/2.Linear M^-stabiiity of stochastic pantograph equationsand the PCMs are researched in the paper. Sufficient conditions of MS-unstability of stochastic pantograph equations and MS-stability of the PCMs are obtained, respectively. Numerical experiments demonstrate these theoretical results.展开更多
A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is ...A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.展开更多
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 ...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.展开更多
A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the ...A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the iteration parameters of the method are tuned in such a way that fast convergence to the value of corrector is achieved.展开更多
The simplified Newton method, at the expense of fast convergence, reduces the work required by Newton method by reusing the initial Jacobian matrix. The composite Newton method attempts to balance the trade-off betwee...The simplified Newton method, at the expense of fast convergence, reduces the work required by Newton method by reusing the initial Jacobian matrix. The composite Newton method attempts to balance the trade-off between expense and fast convergence by composing one Newton step with one simplified Newton step. Recently, Mehrotra suggested a predictor-corrector variant of primal-dual interior point method for linear programming. It is currently the interiorpoint method of the choice for linear programming. In this work we propose a predictor-corrector interior-point algorithm for convex quadratic programming. It is proved that the algorithm is equivalent to a level-1 perturbed composite Newton method. Computations in the algorithm do not require that the initial primal and dual points be feasible. Numerical experiments are made.展开更多
Quasi-periodic responses can appear in a wide variety of nonlinear dynamical systems. To the best of our knowledge, it has been a tough job for years to solve quasi-periodic solutions, even by numerical algorithms. He...Quasi-periodic responses can appear in a wide variety of nonlinear dynamical systems. To the best of our knowledge, it has been a tough job for years to solve quasi-periodic solutions, even by numerical algorithms. Here in this paper, we will present effective and accurate algorithms for quasi-periodic solutions by improving Wilson-θ and Newmark-β methods, respectively. In both the two methods, routinely, the considered equations are rearranged in the form of incremental equilibrium equations with the coefficient matrixes being updated in each time step. In this study, the two methods are improved via a predictor-corrector algorithm without updating the coefficient matrixes, in which the predicted solution at one time point can be corrected to the true one at the next. Numerical examples show that, both the improved Wilson-θ and Newmark-β methods can provide much more accurate quasi-periodic solutions with a smaller amount of computational resources. With a simple way to adjust the convergence of the iterations, the improved methods can even solve some quasi-periodic systems effectively, for which the original methods cease to be valid.展开更多
This paper proposes a continuous block method for the solution of second order ordinary differential equation. Collocation and interpolation of the power series approximate solution are adopted to derive a continuous ...This paper proposes a continuous block method for the solution of second order ordinary differential equation. Collocation and interpolation of the power series approximate solution are adopted to derive a continuous implicit linear multistep method. Continuous block method is used to derive the independent solution which is evaluated at selected grid points to generate the discrete block method. The order, consistency, zero stability and stability region are investigated. The new method was found to compare favourably with the existing methods in term of accuracy.展开更多
This work considers a special case of Lotka-Volterra equations, which means that in the system of two ordinary differential equations, we take the four parameters equal to one. The reason is that we want just to illus...This work considers a special case of Lotka-Volterra equations, which means that in the system of two ordinary differential equations, we take the four parameters equal to one. The reason is that we want just to illustrate the procedure to reduce that system to only one ordinary differential equation, such that we know its analytical solution. This idea will be applied to study the relations between a system of three ordinary differential equations, and a couple of partial differential equations. Lotka-Volterra equations are solved numerically by a fourth-order predictor-corrector method, which is initialized by an improved Euler method with a rather small time step because it is only a second-order algorithm. Then, it is proposed a model with three species, defined by ordinary differential equations.展开更多
基金The Major State Research Program (G1999030803) of China and the NNSF (G10271066, 19972023) of China.
文摘A finite volume element predictor-corrector method for a class of nonlinear parabolic system of equations is presented and analyzed. Suboptimal L^2 error estimate for the finite volume element predictor-corrector method is derived. A numerical experiment shows that the numerical results are consistent with theoretical analysis.
文摘The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.
文摘The paper introduces a new class of numerical schemes for the approximate solutions of stochastic pantograph equations. As an effective technique to implement implicit stochastic methods, strong predictor-corrector methods （PCMs） are designed to handle scenario simulation of solutions of stochastic pantograph equations. It is proved that the PCMs are strong convergent with order 1/2.Linear M^-stabiiity of stochastic pantograph equationsand the PCMs are researched in the paper. Sufficient conditions of MS-unstability of stochastic pantograph equations and MS-stability of the PCMs are obtained, respectively. Numerical experiments demonstrate these theoretical results.
基金supported by the Yunnan Provincial Applied Basic Research Program of China(No. KKSY201207019)
文摘A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.
文摘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.
文摘A parallel diagonally iterated Runge Kutta (PDIRK) method is constructed to solve stiff initial value problems for delay differential equations. The order and stability of this PDIRK method has been analyzed, and the iteration parameters of the method are tuned in such a way that fast convergence to the value of corrector is achieved.
文摘The simplified Newton method, at the expense of fast convergence, reduces the work required by Newton method by reusing the initial Jacobian matrix. The composite Newton method attempts to balance the trade-off between expense and fast convergence by composing one Newton step with one simplified Newton step. Recently, Mehrotra suggested a predictor-corrector variant of primal-dual interior point method for linear programming. It is currently the interiorpoint method of the choice for linear programming. In this work we propose a predictor-corrector interior-point algorithm for convex quadratic programming. It is proved that the algorithm is equivalent to a level-1 perturbed composite Newton method. Computations in the algorithm do not require that the initial primal and dual points be feasible. Numerical experiments are made.
文摘Quasi-periodic responses can appear in a wide variety of nonlinear dynamical systems. To the best of our knowledge, it has been a tough job for years to solve quasi-periodic solutions, even by numerical algorithms. Here in this paper, we will present effective and accurate algorithms for quasi-periodic solutions by improving Wilson-θ and Newmark-β methods, respectively. In both the two methods, routinely, the considered equations are rearranged in the form of incremental equilibrium equations with the coefficient matrixes being updated in each time step. In this study, the two methods are improved via a predictor-corrector algorithm without updating the coefficient matrixes, in which the predicted solution at one time point can be corrected to the true one at the next. Numerical examples show that, both the improved Wilson-θ and Newmark-β methods can provide much more accurate quasi-periodic solutions with a smaller amount of computational resources. With a simple way to adjust the convergence of the iterations, the improved methods can even solve some quasi-periodic systems effectively, for which the original methods cease to be valid.
文摘This paper proposes a continuous block method for the solution of second order ordinary differential equation. Collocation and interpolation of the power series approximate solution are adopted to derive a continuous implicit linear multistep method. Continuous block method is used to derive the independent solution which is evaluated at selected grid points to generate the discrete block method. The order, consistency, zero stability and stability region are investigated. The new method was found to compare favourably with the existing methods in term of accuracy.
文摘This work considers a special case of Lotka-Volterra equations, which means that in the system of two ordinary differential equations, we take the four parameters equal to one. The reason is that we want just to illustrate the procedure to reduce that system to only one ordinary differential equation, such that we know its analytical solution. This idea will be applied to study the relations between a system of three ordinary differential equations, and a couple of partial differential equations. Lotka-Volterra equations are solved numerically by a fourth-order predictor-corrector method, which is initialized by an improved Euler method with a rather small time step because it is only a second-order algorithm. Then, it is proposed a model with three species, defined by ordinary differential equations.
