ADAPTIVE FINITE ELEMENT APPROXIMATION FOR A CLASS OF PARAMETER ESTIMATION PROBLEMS

In this paper, we study adaptive finite element discretisation schemes for a class of parameter estimation problem. We propose to efficient algorithms for the estimation problem use adaptive multi-meshes in developing We derive equivalent a posteriori error estimators for both the state and the control approximation, which particularly suit an adaptive multi-mesh finite element scheme. The error estimators are then implemented and tested with promising numerical 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.     

THE ELLIPSOID ARTIFICIAL BOUNDARY METHOD FOR THREE-DIMENSIONAL UNBOUNDED DOMAINS

The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived explicitly in terms of an infinite series. Then the well-posedness of the coupled variational problem is obtained. It is found that error estimates derived depend on the mesh size, truncation term and the location of the artificial boundary. Three numerical examples are presented to demonstrate the effectiveness and accuracy of the proposed method.     

STABILITY FOR IMPOSING ABSORBING BOUNDARY CONDITIONS IN THE FINITE ELEMENT SIMULATION OF ACOUSTIC WAVE PROPAGATION＊

It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation coupled with the first and the second order absorbing boundary conditions. The computational scheme is also developed. The approach allows the absorbing boundary conditions to be naturally imposed, which makes it easier for us to construct high order schemes for the absorbing boundary conditions. A thirdorder Lagrange finite element method with mass lumping is applied to obtain the spatial discretization of the wave equation. The resulting scheme is stable and is very efficient since no matrix inversion is needed at each time step. Moreover, we have shown both abstract and explicit conditional stability results for the fully-discrete schemes. The results are helpful for designing computational parameters in computations. Numerical computations are illustrated to show the efficiency and accuracy of our method. In particular, essentially no boundary reflection is seen at the artificial boundaries.     

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 RATES FOR DIFFERENCE SCHEMES FOR POLYHEDRAL NONLINEAR PARABOLIC EQUATIONS

We build finite difference schemes for a class of fully nonlinear parabolic equations. The schemes are polyhedral and grid aligned. While this is a restrictive class of schemes, a wide class of equations are well approximated by equations from this class. For regular （C2,α） solutions of uniformly parabolic equations, we also establish of convergence rate of O（α）. A case study along with supporting numerical results is included.     

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.     

A Shishkin mesh for a singularly perturbed Riccati equation

In this paper a singularly perturbed Riccati initial value problem is examined. Parameter explicit bounds on the solution and its derivatives are given. A numerical method composed of an implicit difference operator and a piecewise-uniform Shishkin mesh is constructed. A theoretical parameter independently bound on the errors in the numerical approximations is established. Numerical results are presented which are in agreement with the theoretical error bound.     

A PARAMETER-UNIFORM TAILORED FINITE POINT METHOD FOR SINGULARLY PERTURBED LINEAR ODE SYSTEMS＊

In scientific applications from plasma to chemical kinetics, a wide range of temporal scales can present in a system of differential equations. A major difficulty is encountered due to the stiffness of the system and it is required to develop fast numerical schemes that are able to access previously unattainable parameter regimes. In this work, we consider an initial-final value problem for a multi-scale singularly perturbed system of linear ordi- nary differential equations with discontinuous coefficients. We construct a tailored finite point method, which yields approximate solutions that converge in the maximum norm, uniformly with respect to the singular perturbation parameters, to the exact solution. A parameter-uniform error estimate in the maximum norm is also proved. The results of numerical experiments, that support the theoretical results, are reported.     

THE SECOND-ORDER OPTIMALITY CONDITIONS FOR VARIABLE PROGRAMMING

We study in this paper the continuity of the objective function for variable program- ming. In particular, we study the second-order optimality conditions for unconstrained and constrained variable programming. Some new second-order sufficient and necessary conditions are obtained.     

SPLITTING SCHEMES FOR A NAVIER-STOKES-CAHN-HILLIARD MODEL FOR TWO FLUIDS WITH DIFFERENT DENSITIES

In this work, we focus on designing efficient numerical schemes to approximate a ther- modynamically consistent Navier-Stokes/Cahn-Hilliard problem given in [3] modeling the mixture of two incompressible fluids with different densities. The model is based on a diffuse-interface phase-field approach that is able to describe topological transitions like droplet coalescense or droplet break-up in a natural way. We present a splitting scheme, decoupling computations of the Navier-Stokes part from the Cahn-Hilliard one, which is unconditionally energy-stable up to the choice of the potential approximation. Some nu- merical experiments are carried out to validate the correctness and the accuracy of the scheme, and to study the sensitivity of the scheme with respect to different physical pa- rameters.     

VARIABLE STEP-SIZE IMPLICIT-EXPLICIT LINEAR MULTISTEP METHODS FOR TIME-DEPENDENT PARTIAL DIFFERENTIAL EQUATIONS

Implicit-explicit （IMEX） linear multistep methods are popular techniques for solving partial differential equations （PDEs） with terms of different types. While fixed timestep versions of such schemes have been developed and studied, implicit-explicit schemes also naturally arise in general situations where the temporal smoothness of the solution changes. In this paper we consider easily implementable variable step-size implicit-explicit （VSIMEX） linear multistep methods for time-dependent PDEs. Families of order-p, pstep VSIMEX schemes are constructed and analyzed, where p ranges from 1 to 4. The corresponding schemes are simple to implement and have the property that they reduce to the classical IMEX schemes whenever constant time step-sizes are imposed. The methods are validated on the Burgers＇ equation. These results demonstrate that by varying the time step-size, VSIMEX methods can outperform their fixed time step counterparts while still maintaining good numerical behavior.     

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 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.     

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.     

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.     

Solutions to Boundary Value Problem of Nonlinear Impulsive Differential Equation of Fractional Order

This paper is devoted to study the existence and uniqueness of solutions to a boundary value problem of nonlinear fractional differential equation with impulsive effects.The arguments are based upon Schauder and Banach fixed-point theorems. We improve and generalize the results presented in[B.Ahmad,S.Sivasundaram, Existence results for nonlinear impulsive hybrid boundary value problems involving fractional differential equations,Nonlinear Analysis：Hybrid Systems,3（2009）,251- 258].     

FINITE ELEMENT APPROXIMATIONS FOR SCHROEDINGER EQUATIONS WITH APPLICATIONS TO ELECTRONIC STRUCTURE COMPUTATIONS

In this paper, both the standard finite element discretization and a two-scale finite element discretization for SchrSdinger equations are studied. The numerical analysis is based on the regularity that is also obtained in this paper for the Schroedinger equations. Very satisfying applications to electronic structure computations are provided, too.     

Blow-up vs.Global Finiteness for an Evolution p-Laplace System with Nonlinear Boundary Conditions

In this paper, the authors consider the positive solutions of the system of the evolution p-Laplacian equations
with nonlinear boundary conditions
and the initial data （u0, v0）, where Ω is a bounded domain in Rn with smooth boundary δΩ, p 〉 2, h（·,·） and s（·,· ） are positive C1 functions, nondecreasing in each variable. The authors find conditions on the functions f, g, h, s that prove the global existence or finite time blow-up of positive solutions for every （u0, v0）.     

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.