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.     

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.     

SUPERCONVERGENCE OF GRADIENT RECOVERY SCHEMES ON GRADED MESHES FOR CORNER SINGULARITIES

For the linear finite element solution to the Poisson equation, we show that supercon- vergence exists for a type of graded meshes for corner singularities in polygonal domains. In particular, we prove that the L^2-projection from the piecewise constant field △↓UN to the continuous and piecewise linear finite element space gives a better approximation of △↓U in the Hi-norm. In contrast to the existing superconvergence results, we do not assume high regularity of the exact solution.     

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.     

A POSTERIORI ERROR ESTIMATES OF A NON-CONFORMING FINITE ELEMENT METHOD FOR PROBLEMS WITH ARTIFICIAL BOUNDARY CONDITIONS

An a posteriori error estimator is obtained for a nonconforming finite element approximation of a linear elliptic problem, which is derived from a corresponding unbounded domain problem by applying a nonlocal approximate artificial boundary condition. Our method can be easily extended to obtain a class of a posteriori error estimators for various conforming and nonconforming finite element approximations of problems with different artificial boundary conditions. The reliability and efficiency of our a posteriori error estimator are rigorously proved and are verified by numerical examples.     

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.     

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.     

AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES

This paper presents an efficient moving problem with an optimal control constrained mesh method to solve a nonlinear singular condition. The physical problem is governed by a new model of turbulent flow in circular tubes proposed by Luo et al. using Prandtl＇s mixing-length theory. Our algorithm is formed by an outer iterative algorithm for handling the optimal control condition and an inner adaptive mesh redistribution algorithm for solving the singular governing equations. We discretize the nonlinear problem by using a upwinding approach, and the resulting nonlinear equations are solved by using the Newton- Raphson method. The mesh is generated and the grid points are moved by using the arc-length equidistribution principle. The numerical results demonstrate that proposed algorithm is effective in capturing the boundary layers associated with the turbulent model.     

ERROR REDUCTION IN ADAPTIVE FINITE ELEMENT APPROXIMATIONS OF ELLIPTIC OBSTACLE PROBLEMS

We consider an adaptive finite element method （AFEM） for obstacle problems associated with linear second order elliptic boundary value problems and prove a reduction in the energy norm of the discretization error which leads to R-linear convergence. This result is shown to hold up to a consistency error due to the extension of the discrete multipliers （point functionals） to H^-1 and a possible mismatch between the continuous and discrete coincidence and noncoincidence sets. The AFEM is based on a residual-type error estimator consisting of element and edge residuals. The a posteriori error analysis reveals that the significant difference to the unconstrained case lies in the fact that these residuals only have to be taken into account within the discrete noncoincidence set. The proof of the error reduction property uses the reliability and the discrete local efficiency of the estimator as well as a perturbed Galerkin orthogonality. Numerical results are given illustrating the performance of the AFEM.     

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.     

A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF THE CAHN-HILLIARD EQUATION AND THE HELE-SHAW FLOW

This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut ＋ △（ε△Au-ε^-1f（u）） = 0. It is shown that the a posteriori error bounds depends on ε^-1 only in some low polynomial order, instead of exponential order. Using these a posteriori error estimates, we construct at2 adaptive algorithm for computing the solution of the Cahn- Hilliard equation and its sharp interface limit, the Hele-Shaw flow. Numerical experiments are presented to show the robustness and effectiveness of the new error estimators and the proposed adaptive algorithm.     

A MIXED FINITE ELEMENT METHOD ON A STAGGERED MESH FOR NAVIER-STOKES EQUATIONS

In this paper, we introduce a mixed finite element method on a staggered mesh for the numerical solution of the steady state Navier-Stokes equations in which the two components of the velocity and the pressure are defined on three different meshes. This method is a conforming quadrilateral Q1 × Q1 - P0 element approximation for the Navier-Stokes equations. First-order error estimates are obtained for both the velocity and the pressure. Numerical examples are presented to illustrate the effectiveness of the proposed method.     

NUMERICAL BOUNDARY CONDITIONS FOR THE FAST SWEEPING HIGH ORDER WENO METHODS FOR SOLVING THE EIKONAL EQUATION

High order fast sweeping methods have been developed recently in the literature to solve static Hamilton-Jacobi equations efficiently. Comparing with the first order fast sweeping methods, the high order fast sweeping methods are more accurate, but they often require additional numerical boundary treatment for several grid points near the boundary because of the wider numerical stencil. It is particularly important to treat the points near the inflow boundary accurately, as the information would flow into the computational domain and would affect global accuracy. In the literature, the numerical solution at these boundary points are either fixed with the exact solution, which is not always feasible, or computed with a first order discretization, which could reduce the global accuracy. In this paper, we discuss two strategies to handle the inflow boundary conditions. One is based on the numerical solutions of a first order fast sweeping method with several different mesh sizes near the boundary and a Richardson extrapolation, the other is based on a Lax-Wendroff type procedure to repeatedly utilizing the PDE to write the normal spatial derivatives to the inflow boundary in terms of the tangential derivatives, thereby obtaining high order solution values at the grid points near the inflow boundary. We explore these two approaches using the fast sweeping high order WENO scheme in [18] for solving the static Eikonal equation as a representative example. Numerical examples are given to demonstrate the performance of these two approaches.     

A COARSENING ALGORITHM ON ADAPTIVE GRIDS BY NEWEST VERTEX BISECTION AND ITS APPLICATIONS

In this paper, an efficient and easy-to-implement coarsening algorithm is proposed for adaptive grids obtained using the newest vertex bisection method in two dimemsions. The new coarsening algorithm does not require storing the binary refinement tree explicitly. Instead, the structure is implicitly contained in a special ordering of triangular elements. Numerical experiments demonstrate that the proposed coarsening algorithm is efficient when applied for multilevel preconditioners and mesh adaptivity for time-dependent problems.     

AN ANISOTROPIC NONCONFORMING FINITE ELEMENT METHOD FOR APPROXIMATING A CLASS OF NONLINEAR SOBOLEV EQUATIONS

An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.     

REDUCED BASIS METHOD FOR PARAMETRIZED ELLIPTIC ADVECTION-REACTION PROBLEMS

In this work we consider the Reduced Basis method for the solution of parametrized advection-reaction partial differential equations. For the generation of the basis we adopt a stabilized finite element method and we define the Reduced Basis method in the ＂primal- dual＂ formulation for this stabilized problem. We provide a priori Reduced Basis error estimates and we discuss the effects of the finite element approximation on the Reduced Basis error. We propose an adaptive algorithm, based on the a posteriori Reduced Basis error estimate, for the selection of the sample sets upon which the basis are built; the idea leading this algorithm is the minimization of the computational costs associated with the solution of the Reduced Basis problem. Numerical tests demonstrate the efficiency, in terms of computational costs, of the ＂primal-dual＂ Reduced Basis approach with respect to an ＂only primal＂ one. Parametrized advection-reaction partial differential equations, Reduced Basis method, ＂primal-dual＂ reduced basis approach, Stabilized finite element method, a posteriori error estimation.     

CONVERGENCE OF A MIXED FINITE ELEMENT FOR THE STOKES PROBLEM ON ANISOTROPIC MESHES

The main aim of this paper is to study the convergence properties of a low order mixed finite element for the Stokes problem under anisotropic meshes. We discuss the anisotropic convergence and superconvergence independent of the aspect ratio. Without the shape regularity assumption and inverse assumption on the meshes, the optimal error estimates and natural superconvergence at central points are obtained. The global superconvergence for the gradient of the velocity and the pressure is derived with the aid of a suitable postprocessing method. Furthermore, we develop a simple method to obtain the superclose properties which improves the results of the previous works .     

EVALUATION OF SINGULAR AND NEARLY SINGULAR INTEGRALS IN THE BEM WITH EXACT GEOMETRICAL REPRESENTATION＊

The geometries of many problems of practical interest are created from circular or ellip- tic arcs. Arc boundary elements can represent these boundaries exactly, and consequently, errors caused by representing such geometries using polynomial shape functions can be removed. To fully utilize the geometry of circular boundary, the non-singular boundary integral equations （BIEs） and a general nonlinear transformation technique available for arc elements are introduced to remove or damp out the singular or nearly singular proper- ties of the integral kernels. Several benchmark 2D elastostatic problems demonstrate that the present algorithm can effectively handle singular and nearly singular integrals occur- ring in the boundary element method （BEM） for boundary layer effect and thin-walled structural problems. Owing to the employment of exact geometrical representation, only a small number of elements need to be divided along the boundary and high accuracy can be achieved without increasing other more computational efforts.     

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.