A NEW DIRECT DISCONTINUOUS GALERKIN METHOD WITH SYMMETRIC STRUCTURE FOR NONLINEAR DIFFUSION EQUATIONS

In this paper we continue the study of discontinuous Galerkin finite element methods for nonlinear diffusion equations following the direct discontinuous Galerkin （DDG） meth- ods for diffusion problems [17] and the direct discontinuous Galerkin （DDG） methods for diffusion with interface corrections [18]. We introduce a numerical flux for the test func- tion, and obtain a new direct discontinuous Galerkin method with symmetric structure. Second order derivative jump terms are included in the numerical flux formula and explicit guidelines for choosing the numerical flux are given. The constructed scheme has a sym- metric property and an optimal L2 （L2） error estimate is obtained. Numerical examples are carried out to demonstrate the optimal （k ＋ 1）th order of accuracy for the method with pk polynomial approximations for both linear and nonlinear problems, under one-dimensional and two-dimensional settings.     

ON THE OPTIMAL CONVERGENCE RATE OF A ROBIN-ROBIN DOMAIN DECOMPOSITION METHOD

In this work, we solve a long-standing open problem： Is it true that the convergence rate of the Lions＇ Robin-Robin nonoverlapping domain decomposition （DD） method can be constant, independent of the mesh size h？ We closed this old problem with a positive answer. Our theory is also verified by numerical tests.     

A METHOD FOR SOLVING THE INVERSE SCATTERING PROBLEM FOR SHAPE AND IMPEDANCE

The inverse problem considered in this paper is to determine the shape and the impedance of an obstacle from a knowledge of the time-harmonic incident field and the phase and amplitude of the far field pattern of the scattered wave in two-dimension. Single-layer potential is used to approach the scattered waves. An approximation method is presented and the convergence of the proposed method is established. Numerical examples are given to show that this method is both accurate and easy to use.     

VARIABLE MESH FINITE DIFFERENCE METHOD FOR SELF-ADJOINT SINGULARLY PERTURBED TWO-POINT BOUNDARY VALUE PROBLEMS

A numerical method based on finite difference method with variable mesh is given for self-adjoint singularly perturbed two-point boundary value problems. To obtain parameter- uniform convergence, a variable mesh is constructed, which is dense in the boundary layer region and coarse in the outer region. The uniform convergence analysis of the method is discussed. The original problem is reduced to its normal form and the reduced problem is solved by finite difference method taking variable mesh. To support the efficiency of the method, several numerical examples have been considered.     

THE FINITE DIFFERENCE METHOD FOR DISSIPATIVE KLEIN-GORDON-SCHRODINGER EQUATIONS IN THREE SPACE DIMENSIONS

A fully discrete finite difference scheme for dissipative Klein-Gordon-SchrSdinger equations in three space dimensions is analyzed. On the basis of a series of the time-uniform priori estimates of the difference solutions and discrete version of Sobolev embedding the- orems, the stability of the difference scheme and the error bounds of optimal order for the difference solutions are obtained in H2 × H2 ×H1 over a finite time interval. Moreover, the existence of a maximal attractor is proved for a discrete dynamical system associated with the fully discrete finite difference scheme.     

ENERGY ESTIMATES FOR DELAY DIFFUSION-REACTION EQUATIONS

In this paper we consider nonlinear delay diffusion-reaction equations with initial and Dirichlet boundary conditions. The behaviour and the stability of the solution of such initial boundary value problems （IBVPs） are studied using the energy method. Simple numerical methods are considered for the computation of numerical approximations to the solution of the nonlinear IBVPs. Using the discrete energy method we study the stability and convergence of the numerical approximations. Numerical experiments are carried out to illustrate our theoretical results.     

A COMPACT UPWIND SECOND ORDER SCHEME FOR THE EIKONAL EQUATION

We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on： 1. the numerical observation that classical first order monotone upwind schemes for the Eikonal equation yield numerical upwind gradient which is also first order accurate up to singularities; 2. a remark that partial information on the second derivatives of the solution is known and given in the structure of the Eikonal equation and can be used to reduce the size of the stencil. We implement the second order scheme as a correction to the well known sweeping method but it should be applicable to any first order monotone upwind scheme. Care is needed to choose the appropriate stencils to avoid instabilities.     

ADAPTIVE QUADRILATERAL AND HEXAHEDRAL FINITE ELEMENT METHODS WITH HANGING NODES AND CONVERGENCE ANALYSIS

In this paper we study the convergence of adaptive finite element methods for the gen- eral non-attine equivalent quadrilateral and hexahedral elements on 1-irregular meshes with hanging nodes. Based on several basic ingredients, such as quasi-orthogonality, estimator reduction and D6fler marking strategy, convergence of the adaptive finite element methods for the general second-order elliptic partial equations is proved. Our analysis is effective for all conforming Qm elements which covers both the two- and three-dimensional cases in a unified fashion.     

CONVERGENCE OF GRADIENT METHOD WITH MOMENTUM FOR BACK-PROPAGATION NEURAL NETWORKS

In this work, a gradient method with momentum for BP neural networks is considered. The momentum coefficient is chosen in an adaptive manner to accelerate and stabilize the learning procedure of the network weights. Corresponding convergence results are proved.     

ERROR ESTIMATES FOR THE RECURSIVE LINEARIZATION OF INVERSE MEDIUM PROBLEMS

This paper is devoted to the mathematical analysis of a general recursive linearization algorithm for solving inverse medium problems with multi-frequency measurements. Under some reasonable assumptions, it is shown that the algorithm is convergent with error estimates. The work is motivated by our effort to analyze recent significant numerical results for solving inverse medium problems. Based on the uncertainty principle, the recursive linearization allows the nonlinear inverse problems to be reduced to a set of linear problems and be solved recursively in a proper order according to the measurements. As an application, the convergence of the recursive linearization algorithm [Chen, Inverse Problems 13（1997）, pp.253-282] is established for solving the acoustic inverse scattering problem.     

ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS

The main purpose of this work is to provide a novel numerical approach for the Volterra integral equations based on a spectral approach. A Legendre-collocation method is proposed to solve the Volterra integral equations of the second kind. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical results confirm the theoretical prediction of the exponential rate of convergence. The result in this work seems to be the first successful spectral approach （with theoretical justification） for the Volterra type equations.     

MIXED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR ONE-DIMENSIONAL FULLY NONLINEAR SECOND ORDER ELLIPTIC AND PARABOLIC EQUATIONS

This paper is concerned with developing accurate and efficient numerical methods for one-dimensional fully nonlinear second order elliptic and parabolic partial differential equations （PDEs）. In the paper we present a general framework for constructing high order interior penalty discontinuous Galerkin （IP-DG） methods for approximating viscosity solutions of these fully nonlinear PDEs. In order to capture discontinuities of the second order derivative uxx of the solution u, three independent functions p1,p2 and p3 are introduced to represent numerical derivatives using various one-sided limits. The proposed DG frame- work, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a nonlinear problem into a mostly linear system of equations where the nonlinearity has been modified to include multiple values of the second order derivative uxz. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes. In addition to the nonstandard mixed formulation setting, another main idea is to replace the fully nonlinear differential operator by a numerical operator which is consistent with the differential operator and satisfies certain monotonicity （called g-monotonicity） properties. To ensure such a g-monotonicity, the crux of the construction is to introduce the numerical moment, which plays a critical role in the proposed DG frame- work. The g-monotonicity gives the DG methods the ability to select the mathematically ＂correct＂ solution （i.e., the viscosity solution） among all possible solutions. Moreover, the g-monotonicity allows for the possible development of more efficient nonlinear solvers as the special nonlinearity of the algebraic systems can be explored to decouple the equations. This paper also presents and analyzes numerical results for several numerical test problems which are used to guage the accuracy and efficiency of the proposed DG methods.     

UNIFORMLY-STABLE FINITE ELEMENT METHODS FOR DARCY-STOKES-BRINKMAN MODELS

In this paper, we consider 2D and 3D Darcy-Stokes interface problems. These equations are related to Brinkman model that treats both Darcy＇s law and Stokes equations in a single form of PDE but with strongly discontinuous viscosity coefficient and zerothorder term coefficient. We present three different methods to construct uniformly stable finite element approximations. The first two methods are based on the original weak formulations of Darcy-Stokes-Brinkman equations. In the first method we consider the existing Stokes elements. We show that a stable Stokes element is also uniformly stable with respect to the coefficients and the jumps of Darcy-Stokes-Brinkman equations if and only if the discretely divergence-free velocity implies almost everywhere divergence-free one. In the second method we construct uniformly stable elements by modifying some well-known H（div）-conforming elements. We give some new 2D and 3D elements in a unified way. In the last method we modify the original weak formulation of Darcy-Stokes- Brinkman equations with a stabilization term. We show that all traditional stable Stokes elements are uniformly stable with respect to the coefficients and their jumps under this new formulation.     

THE GENERALIZED LOCAL HERMITIAN AND SKEW-HERMITIAN SPLITTING ITERATION METHODS FOR THE NON-HERMITIAN GENERALIZED SADDLE POINT PROBLEMS

For large and sparse saddle point problems, Zhu studied a class of generalized local Hermitian and skew-Hermitian splitting iteration methods for non-Hermitian saddle point problem [M.-Z. Zhu, Appl. Math. Comput. 218 （2012） 8816-8824 ]. In this paper, we further investigate the generalized local Hermitian and skew-Hermitian splitting （GLHSS） iteration methods for solving non-Hermitian generalized saddle point problems. With different choices of the parameter matrices, we derive conditions for guaranteeing the con- vergence of these iterative methods. Numerical experiments are presented to illustrate the effectiveness of our GLHSS iteration methods as well as the preconditioners.     

THE OPTIMAL CONVERGENCE ORDER OF THE DISCONTINUOUS FINITE ELEMENT METHODS FOR FIRST ORDER HYPERBOLIC SYSTEMS

In this paper, a discontinuous finite element method for the positive and symmetric, first-order hyperbolic systems （steady and nonsteady state） is constructed and analyzed by using linear triangle elements, and the O（h^2）-order optimal error estimates are derived under the assumption of strongly regular triangulation and the Ha-regularity for the exact solutions. The convergence analysis is based on some superclose estimates of the interpolation approximation. Finally, we discuss the Maxwell equations in a two-dimensional domain, and numerical experiments are given to validate the theoretical results.     

STABILIZED FEM FOR CONVECTION-DIFFUSION PROBLEMS ON LAYER-ADAPTED MESHES

The application of a standard Galerkin finite element method for convection-diffusion problems leads to oscillations in the discrete solution, therefore stabilization seems to be necessary. We discuss several recent stabilization methods, especially its combination with a Galerkin method on layer-adapted meshes. Supercloseness results obtained allow an improvement of the discrete solution using recovery techniques.     

A TWO-SCALE HIGHER-ORDER FINITE ELEMENT DISCRETIZATION FOR SCHROEDINGER EQUATION

In this paper, a two-scale higher-order finite element discretization scheme is proposed and analyzed for a Schroedinger equation on tensor product domains. With the scheme, the solution of the eigenvalue problem on a fine grid can be reduced to an eigenvalue problem on a much coarser grid together with some eigenvalue problems on partially fine grids. It is shown theoretically and numerically that the proposed two-scale higher-order scheme not only significantly reduces the number of degrees of freedom but also produces very accurate approximations.     

MULTIGRID METHOD FOR FLUID DYNAMIC PROBLEMS

This paper covers the dynamics problems. The review and some aspects of main development stages of using Multigrid method for fluid multigrid technics are presented. Some approaches for solving Navier-Stokes equations and convection- diffusion problems are considered.     

ON CONTRACTION AND SEMI-CONTRACTION FACTORS OF GSOR METHOD FOR AUGMENTED LINEAR SYSTEMS

The generalized successive overrelaxation （GSOR） method was presented and studied by Bai, Parlett and Wang [Numer. Math. 102（2005）, pp.1-38] for solving the augmented system of linear equations, and the optimal iteration parameters and the corresponding optimal convergence factor were exactly obtained. In this paper, we further estimate the contraction and the semi-contraction factors of the GSOR method. The motivation of the study is that the convergence speed of an iteration method is actually decided by the contraction factor but not by the spectral radius in finite-step iteration computations. For the nonsingular augmented linear system, under some restrictions we obtain the contraction domain of the parameters involved, which guarantees that the contraction factor of the GSOR method is less than one. For the singular but consistent augmented linear system, we also obtain the semi-contraction domain of the parameters in a similar fashion. Finally, we use two numerical examples to verify the theoretical results and the effectiveness of the GSOR method.     

A FEM-BEM FORMULATION FOR AN EXTERIOR QUASILINEAR ELLIPTIC PROBLEM IN THE PLANE

In this paper, the finite element method and the boundary element method are combined to solve numerically an exterior quasilinear elliptic problem. Based on an appropriate transformation and the Fourier series expansion, the exact quasilinear artificial boundary conditions and a series of the corresponding approximations for the given problem are presented. Then the original problem is reduced into an equivalent problem defined in a bounded computational domain. We provide error estimate for the Galerkin method. Numerical results are presented to illustrate the theoretical results.