The Asymptotic Behavior for Numerical Solution of a Volterra Equation

Long-time asymptotic stability and convergence properties for the numerical solution of a Volterra equation of parabolic type are studied.The methods are based on the first-second order backward difference methods.The memory term is approximated by the comvolution quadrature and the interpolant quadrature.Discretization of the spatial partial differential operators by the finite element method is also considered.     

Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R~2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.     

RESEARCH ANNOUNCEMENTS High order Zumerical Solution of Integral Transport Equation in slab Geometry

There are some common numerical methods for solving neutron transport equation, which including the well-known discrete ordinates method, PN approximation and integral transport methods[1]. There exists certain singularities in the solution of transport equation near the boundary and interface[2]. It gives rise to the difficulty in the construction of high order accurate numerical methods. The numerical solution obtained by now can not attain the second order convergent accuracy[3,4].     

ON NEWTON-HSS METHODS FOR SYSTEMS OF NONLINEAR EQUATIONS WITH POSITIVE-DEFINITE JACOBIAN MATRICES

The Hermitian and skew-Hermitian splitting （HSS） method is an unconditionally convergent iteration method for solving large sparse non-Hermitian positive definite system of linear equations. By making use of the HSS iteration as the inner solver for the Newton method, we establish a class of Newton-HSS methods for solving large sparse systems of nonlinear equations with positive definite Jacobian matrices at the solution points. For this class of inexact Newton methods, two types of local convergence theorems are proved under proper conditions, and numerical results are given to examine their feasibility and effectiveness. In addition, the advantages of the Newton-HSS methods over the Newton-USOR, the Newton-GMRES and the Newton-GCG methods are shown through solving systems of nonlinear equations arising from the finite difference discretization of a two-dimensional convection-diffusion equation perturbed by a nonlinear term. The numerical implemen- tations also show that as preconditioners for the Newton-GMRES and the Newton-GCG methods the HSS iteration outperforms the USOR iteration in both computing time and iteration step.     

矩阵方程$X A^*X^{-q}A=Q~(q\geq 1)$的Hermite正定解

In this paper, Hermitian positive definite solutions of the nonlinear matrix equation X ＋ A^＊X^-qA = Q （q≥1） are studied. Some new necessary and sufficient conditions for the existence of solutions are obtained. Two iterative methods are presented to compute the smallest and the quasi largest positive definite solutions, and the convergence analysis is also given. The theoretical results are illustrated by numerical examples.     

关于非线性矩阵方程$X^s-A^*X^{-t}A=Q$的Hermitian正定解及其扰动估计

In this paper, the Hermitian positive definite solutions of the nonlinear matrix equation X^s - A^＊X^-tA = Q are studied, where Q is a Hermitian positive definite matrix, s and t are positive integers. The existence of a Hermitian positive definite solution is proved. A sufficient condition for the equation to have a unique Hermitian positive definite solution is given. Some estimates of the Hermitian positive definite solutions are obtained. Moreover, two perturbation bounds for the Hermitian positive definite solutions are derived and the results are illustrated by some numerical examples.     

Convergence analysis for the Secant method based on new recurrence relations

A new convergence theorem for the Secant method in Banach spaces based on new recurrence relations is established for approximating a solution of a nonlinear operator equation. It is assumed that the divided difference of order one of the nonlinear operator is Lipschitz continuous. The convergence conditions differ from some existing ones and are easily satisfied. The results of the paper are justified by numerical examples that cannot be handled by earlier works.     

Convergence analysis for the Secant method based on new recurrence relations

A new convergence theorem for the Secant method in Banach spaces based on new recurrence relations is established for approximating a solution of a nonlinear operator equation.It is assumed that the divided difference of order one of the nonlinear operator is Lipschitz continuous.The convergence conditions differ from some existing ones and are easily satisfied.The results of the paper are justified by numerical examples that cannot be handled by earlier works.     

EXPONENTIAL STABILITY FOR NONLINEAR HYBRID STOCHASTIC PANTOGRAPH EQUATIONS AND NUMERICAL APPROXIMATION

The paper develops exponential stability of the analytic solution and convergence in probability of the numerical method for highly nonlinear hybrid stochastic pantograph equation. The classical linear growth condition is replaced by polynomial growth conditions, under which there exists a unique global solution and the solution is almost surely exponentially stable. On the basis of a series of lemmas, the paper establishes a new criterion on convergence in probability of the Euler-Maruyama approximate solution. The criterion is very general so that many highly nonlinear stochastic pantograph equations can obey these conditions. A highly nonlinear example is provided to illustrate the main theory.     

Wavelet Collocation Methods for Viscosity Solutions to Swing Options in Natural Gas Storage

This paper presents the wavelet collocation methods for the numerical ap- proximation of swing options for natural gas storage in a mean reverting market. The model is characterized by the Hamilton-Jacobi-Bellman （HJB） equations which only have the viscosity solution due to the irregularity of the swing option. The differential operator is formulated exactly and efficiently in the second generation interpolating wavelet setting. The convergence and stability of the numerical scheme are studied in the framework of viscosity solution theory. Numerical experiments demonstrate the accuracy and computational efficiency of the methods.     

On the Nonlinear Matrix Equation X+ A *f1（X）A +B*f2（X）B=Q

In this paper, nonlinear matrix equations of the form X + A*f1 (X)A + B*f2 (X)B = Q are discussed. Some necessary and sufficient conditions for the existence of solutions for this equation are derived. It is shown that under some conditions this equation has a unique solution, and an iterative method is proposed to obtain this unique solution. Finally, a numerical example is given to identify the efficiency of the results obtained.     

Numerical Simulation for Convection-Diffusion Problem with Periodic Micro-Structure

In this article, the convection dominated convection-diffusion problems with the periodic micro-structure are discussed. A two-scale finite element scheme based on the homogenization technique for this kind of problems is provided. The error estimates between the exact solution and the approximation solution, of the homogenized equation or the two-scale finite element scheme are analyzed. It is shown that the scheme provided in this article is convergent for any fixed diffusion coefficient 5, and it may be convergent independent of δ under some conditions. The numerical results demonstrating the theoretical results are presented in this article.     

OSCILLATION CRITERIA FOR DELAY DIFFERENTIAL EQUATIONS WITH IMPULSES

The present paper is devoted to the investigation of the oscillation of second order nonlinear delay differential equation (DDE) with impulses, sufficient conditions of oscillation for the solutions of DDE with impulses are obtained by using differential inequality. Some known results are improved and generalized.     

A CLASS OF SINGULARLY PERTURBED INITIAL BOUNDARY PROBLEM FOR REACTION DIFFUSION EQUATION

The singularly perturbed initial boundary value problem for a class of reaction diffusion equation is considered. Under appropriate conditions, the existence-uniqueness and the asymptotic behavior of the solution are showed by using the fixed-point theorem.     

ERROR ANALYSIS FOR A FAST NUMERICAL METHOD TO A BOUNDARY INTEGRAL EQUATION OF THE FIRST KIND

For two-dimensional boundary integral equations of the first kind with logarithmic kernels, the use of the conventional boundary element methods gives linear systems with dense matrix. In a recent work [J. Comput. Math., 22 （2004）, pp. 287-298], it is demonstrated that the dense matrix can be replaced by a sparse one if appropriate graded meshes are used in the quadrature rules. The numerical experiments also indicate that the proposed numerical methods require less computational time than the conventional ones while the formal rate of convergence can be preserved. The purpose of this work is to establish a stability and convergence theory for this fast numerical method. The stability analysis depends on a decomposition of the coefficient matrix for the collocation equation. The formal orders of convergence observed in the numerical experiments are proved rigorously.     

GL method for solving equations in math-physics and engineering

In this paper, we propose a GL method for solving the ordinary and the partial differential equation in mathematical physics and chemics and engineering. These equations govern the acustic, heat, electromagnetic, elastic, plastic, flow, and quantum etc. macro and micro wave field in time domain and frequency domain. The space domain of the differential equation is infinite domain which includes a finite inhomogeneous domain. The inhomogeneous domain is divided into finite sub domains. We present the solution of the differential equation as an explicit recursive sum of the integrals in the inhomogeneous sub domains. Actualy, we propose an explicit representation of the inhomogeneous parameter nonlinear inversion. The analytical solution of the equation in the infinite homogeneous domain is called as an initial global field. The global field is updated by local scattering field successively subdomaln by subdomain. Once all subdomains are scattered and the updating process is finished in all the sub domains, the solution of the equation is obtained. We call our method as Global and Local field method, in short , GL method. It is different from FEM method, the GL method directly assemble inverse matrix and gets solution. There is no big matrix equation needs to solve in the GL method. There is no needed artificial boundary and no absorption boundary condition for infinite domain in the GL method. We proved several theorems on relationships between the field solution and Green＇s function that is the theoretical base of our GL method. The numerical discretization of the GL method is presented. We proved that the numerical solution of the GL method convergence to the exact solution when the size of the sub domain is going to zero. The error estimation of the GL method for solving wave equation is presented. The simulations show that the GL method is accurate, fast, and stable for solving elliptic, parabolic, and hyperbolic equations. The GL method has advantages and wide applications in the 3D electromagnetic （EM）     

OSCILLATION AND GLOBAL ATTRACTIVITY OF IMPULSIVE PERIODIC DELAY RESPIRATORY DYNAMICS MODEL

This paper studies the nonlinear delay impulsive respiratory dynamics model. The model describes the sudden changes of the concentration of CO2 in the blood of the mammal. It is proved that the model has a unique positive periodic solution. Some sufficient conditions for oscillation of all positive solutions about the positive periodic solution are established and also some sufficient conditions for the global attractivity of the periodic solution are obtained.     

TWO-GRID CHARACTERISTIC FINITE VOLUME METHODS FOR NONLINEAR PARABOLIC PROBLEMS＊

In this work, two-grid characteristic finite volume schemes for the nonlinear parabolic problem are considered. In our algorithms, the diffusion term is discretized by the finite volume method, while the temporal differentiation and advection terms are treated by the characteristic scheme. Under some conditions about the coefficients and exact solution, optimal error estimates for the numerical solution are obtained. Furthermore, the two- grid characteristic finite volume methods involve solving a nonlinear equation on coarse mesh with mesh size H, a large linear problem for the Oseen two-grid characteristic finite volume method on a fine mesh with mesh size h = O（H2） or a large linear problem for the Newton two-grid characteristic finite volume method on a fine mesh with mesh size h = 0（I log hll/2H3）. These methods we studied provide the same convergence rate as that of the characteristic finite volume method, which involves solving one large nonlinear problem on a fine mesh with mesh size h. Some numerical results are presented to demonstrate the efficiency of the proposed methods.     

D-CONVERGENCE OF RUNGE-KUTTA METHODS FOR STIFF DELAY DIFFERENTIAL EQUATIONS

1. IntroductionWhen considering the applicability of numerical methods for the solution of the delay differential equation (DDE) y'(t) = f(t, y(t), y(t - T)), it is necessary to analyze the error behaviourof the methods. In fact, many papers have investigated the local and global error behaviour ofDDE solvers (cL[1,2,14]). These error analyses are based on the assumption that the fUnctionf(t,y,z) satisfies Lipschitz conditions in both the last two variables. They are suitable fornonstiff …