A RELAXATION RESULT OF FUNCTIONALS IN THE SPACE SBD（Ω）

In this article,the authors obtain an integral representation for the relaxation of the functional
F（x,u,Ω）：={∫^f（x,u（x）,εu（x））dx Ω if u∈W^1,1（Ω,R^N）, ＋∞ otherwise, in the space of functions of bounded deformation,with respect to L^1-convergence.Here Eu represents the absolutely continuous part of the symmetrized distributional derivative Eu.f（x,p,ξ）satisfying weak convexity assumption.     

Optimal integrability of some system of integral equations

We obtain the optimal integrability for positive solutions of the Euler-Lagrange system of the weighted Hardy-Littlewood-Sobolev inequality in R^n ：{u（x）=1/｜x｜^α｜∫R^n v（y）^q｜y｜^β｜x-y｜^λdy,v（x）=1/｜x｜^β∫R^n u（y）^p｜y｜^α｜x-y｜^λdy.C. Jin and C. Li [Calc. Var. Partial Differential Equations, 2006, 26： 447-457] developed some very interesting method for regularity lifting and obtained the optimal integrability for p, q 〉 1. Here, based on some new observations, we overcome the difficulty there, and derive the optimal integrability for the case of p, q ≥1 and pq ≠1. This integrability plays a key role in estimating the asymptotic behavior of positive solutions when ｜x｜ →0 and when ｜x｜→∞.     

ON THE ORDER OF SUMMABILITY OF THE FOURIER INVERSION FORMULA

In this article we show that the order of the point value, in the sense of Lojasiewicz, of a tempered distribution and the order of summability of the pointwise Fourier inversion formula are closely related. Assuming that the order of the point values and certain order of growth at infinity are given for a tempered distribution, we estimate the order of summability of the Fourier inversion formula. For Fourier series, and in other cases, it is shown that if the distribution has a distributional point value of order k, then its Fourier series is e.v. Cesaro summable to the distributional point value of order k＋1. Conversely, we also show that if the pointwise Fourier inversion formula is e.v. Cesaro summable of order k, then the distribution is the （k ＋ 1）-th derivative of a locally integrable function, and the distribution has a distributional point value of order k ＋ 2. We also establish connections between orders of summability and local behavior for other Fourier inversion problems.     

HIGHLY OSCILLATORY DIFFUSION-TYPE EQUATIONS

We explore new asymptotic-numeric solvers for partial differential equations with highly oscillatory forcing terms. Such methods represent the solution as an asymptotic series, whose terms can be evaluated by solving non-oscillatory problems and they guarantee high accuracy at a low computational cost. We consider two forms of oscillatory forcing terms, namely when the oscillation is in time or in space： each lends itself to different treatment. Numerical examples highlight the salient features of the new approach.     

NUMERICAL RECONSTRUCTION OF SMALL PERTURBATIONS IN THE ELECTROMAGNETIC COEFFICIENTS OF A DIELECTRIC MATERIAL

The aim of this paper is to solve numerically the inverse problem of reconstructing small amplitude perturbations in the magnetic permeability of a dielectric material from partial or total dynamic boundary measurements. Our numerical algorithm is based on the resolution of the time-dependent Maxwell equations, an exact controllability method and Fourier inversion for localizing the perturbations. Two-dimensional numerical experiments illustrate the performance of the reconstruction method for different configurations even in the case of limited-view data.     

Uncertainty Principles for the Generalized Fourier Transform Associated with Spherical Mean Operator

The aim of this paper is to establish an extension of qualitative and quantitative uncertainty principles for the Fourier transform connected with the spherical mean operator.     

A RELAXED HSS PRECONDITIONER FOR SADDLE POINT PROBLEMS FROM MESHFREE DISCRETIZATION＊

In this paper, a relaxed Hermitian and skew-Hermitian splitting （RHSS） preconditioner is proposed for saddle point problems from the element-free Galerkin （EFG） discretization method. The EFG method is one of the most widely used meshfree methods for solving partial differential equations. The RHSS preconditioner is constructed much closer to the coefficient matrix than the well-known HSS preconditioner, resulting in a RHSS fixed-point iteration. Convergence of the RHSS iteration is analyzed and an optimal parameter, which minimizes the spectral radius of the iteration matrix is described. Using the RHSS pre- conditioner to accelerate the convergence of some Krylov subspace methods （like GMRES） is also studied. Theoretical analyses show that the eigenvalues of the RHSS precondi- tioned matrix are real and located in a positive interval. Eigenvector distribution and an upper bound of the degree of the minimal polynomial of the preconditioned matrix are obtained. A practical parameter is suggested in implementing the RHSS preconditioner. Finally, some numerical experiments are illustrated to show the effectiveness of the new preconditioner.     

ON THE SEPARABLE NONLINEAR LEAST SQUARES PROBLEMS

Separable nonlinear least squares problems are a special class of nonlinear least squares problems, where the objective functions are linear and nonlinear on different parts of variables. Such problems have broad applications in practice. Most existing algorithms for this kind of problems are derived from the variable projection method proposed by Golub and Pereyra, which utilizes the separability under a separate framework. However, the methods based on variable projection strategy would be invalid if there exist some constraints to the variables, as the real problems always do, even if the constraint is simply the ball constraint. We present a new algorithm which is based on a special approximation to the Hessian by noticing the fact that certain terms of the Hessian can be derived from the gradient. Our method maintains all the advantages of variable projection based methods, and moreover it can be combined with trust region methods easily and can be applied to general constrained separable nonlinear problems. Convergence analysis of our method is presented and numerical results are also reported.     

A LOCAL MULTILEVEL PRECONDITIONER FOR THE ADAPTIVE MORTAR FINITE ELEMENT METHOD＊

In this paper, we propose a local multilevel preconditioner for the mortar finite element approximations of the elliptic problems. With some mesh assumptions on the interface, we prove that the condition number of the preconditioned systems is independent of the large jump of the coefficients but depends on the mesh levels around the cross points. Some numericM experiments are presented to confirm our theoreticM results.     

SUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS

We propose and analyze a C^0 spectral element method for a model eigenvalue problem with discontinuous coefficients in the one dimensional setting. A super-geometric rate of convergence is proved for the piecewise constant coefficients case and verified by numerical tests. Furthermore, the asymptotical equivalence between a Gauss-Lobatto collocation method and a spectral Galerkin method is established for a simplified model.     

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.     

ON SMOOTH LU DECOMPOSITIONS WITH APPLICATIONS TO SOLUTIONS OF NONLINEAR EIGENVALUE PROBLEMS

We study the smooth LU decomposition of a given analytic functional A-matrix A（A） and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A（A）, and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.     

COUPLING OF FINITE ELEMENT AND BOUNDARY ELEMENT METHODS FOR THE SCATTERING BY PERIODIC CHIRAL STRUCTURES

Consider a time-harmonic electromagnetic plane wave incident on a biperiodic structure in R^3. The periodic structure separates two homogeneous regions. The medium inside the structure is chiral and nonhomogeneous. In this paper, variational formulations coupling finite element methods in the chiral medium with a method of integral equations on the periodic interfaces are studied. The well-posedness of the continuous and discretized problems is established. Uniform convergence for the coupling variational approximations of the model problem is obtained.     

POISSON PRECONDITIONING FOR SELF-ADJOINT ELLIPTIC PROBLEMS

In this paper, we formulate interface problem and Neumann elliptic boundary value problem into a form of linear operator equations with self-adjoint positive definite op- erators. We prove that in the discrete level the condition number of these operators is independent of the mesh size. Therefore, given a prescribed error tolerance, the classical conjugate gradient algorithm converges within a fixed number of iterations. The main computation task at each iteration is to solve a Dirichlet Poisson boundary value problem in a rectangular domain, which can be furnished with fast Poisson solver. The overall computational complexity is essentially of linear scaling.     

EFFICIENT NUMERICAL ALGORITHMS FOR THREE-DIMENSIONAL FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS

This paper detailedly discusses the locally one-dimensional numerical methods for ef- ficiently solving the three-dimensional fractional partial differential equations, including fractional advection diffusion equation and Riesz fractional diffusion equation. The second order finite difference scheme is used to discretize the space fractional derivative and the Crank-Nicolson procedure to the time derivative. We theoretically prove and numerically verify that the presented numerical methods are unconditionally stable and second order convergent in both space and time directions. In particular, for the Riesz fractional dif- fusion equation, the idea of reducing the splitting error is used to further improve the algorithm, and the unconditional stability and convergency are also strictly proved and numerically verified for the improved scheme.     

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.     

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.     

DISSIPATIVE NUMERICAL METHODS FOR THE HUNTER-SAXTON EQUATION

In this paper, we present further development of the local discontinuous Galerkin （LDG） method designed in [21] and a new dissipative discontinuous Galerkin （DG） method for the HuntermSaxton equation. The numerical fluxes for the LDG and DG methods in this paper are based on the upwinding principle. The resulting schemes provide additional energy dissipation and better control of numerical oscillations near derivative singularities. Stability and convergence of the schemes are proved theoretically, and numerical simulation results are provided to compare with the scheme in [21].     

FINITE ELEMENT AND DISCONTINUOUS GALERKIN METHOD FOR STOCHASTIC HELMHOLTZ EQUATION IN TWO- AND THREE-DIMENSIONS

In this paper, we consider the finite element method and discontinuous Galerkin method for the stochastic Helmholtz equation in R^d （d = 2, 3）. Convergence analysis and error estimates are presented for the numerical solutions. The effects of the noises on the accuracy of the approximations are illustrated. Numerical experiments are carried out to verify our theoretical results.     

ON THE ERROR ESTIMATES OF A NEW OPERATE SPLITTING METHOD FOR THE NAVIER-STOKES EQUATIONS

In this paper, a new operator splitting scheme is introduced for the numerical solution of the incompressible Navier-Stokes equations. Under some mild regularity assumptions on the PDE solution, the stability of the scheme is presented, and error estimates for the velocity and the pressure of the proposed operator splitting scheme are given.