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.     

FINITE ELEMENTS WITH LOCAL PROJECTION STABILIZATION FOR INCOMPRESSIBLE FLOW PROBLEMS

In this paper we review recent developments in the analysis of finite element methods for incompressible flow problems with local projection stabilization （LPS）. These methods preserve the favourable stability and approximation properties of classical residual-based stabilization （RBS） techniques but avoid the strong coupling of velocity and pressure in the stabilization terms. LPS-methods belong to the class of symmetric stabilization techniques and may be characterized as variational multiscale methods. In this work we summarize the most important a priori estimates of this class of stabilization schemes developed in the past 6 years. We consider the Stokes equations, the Oseen linearization and the NavierStokes equations. Furthermore, we apply it to optimal control problems with linear（ized） flow problems, since the symmetry of the stabilization leads to the nice feature that the operations ＂discretize＂ and ＂optimize＂ commute.     

Regularity of the inverse of a homeomorphism with finite inner distortion

Let f:Ω→ f(Ω) R~n be a W~(1,1)-homeomorphism with L~1-integrable inner distortion.We show that finiteness of min{lip_f(x),k_f/(x)},for every x ∈Ω\E,implies that f~(-1) ∈ W~(1,n)and has finite distortion,provided that the exceptional set E has cr-finite H~1-measure.Moreover,/ has finite distortion,differentiable a.e.and the Jacobian J_f 0 a.e.     

ON THE P1 POWELL-SABIN DIVERGENCE-FREE FINITE ELEMENT FOR THE STOKES EQUATIONS

The stability of the P1-P0 mixed-element is established on general Powell-Sabin triangular grids. The piecewise linear finite element solution approximating the velocity is divergence-free pointwise for the Stokes equations. The finite element solution approximating the pressure in the Stokes equations can be obtained as a byproduct if an iterative method is adopted for solving the discrete linear system of equations. Numerical tests are presented confirming the theory on the stability and the optimal order of convergence for the P1 Powell-Sabin divergence-free finite element method.     

A Quasilinear Parabolic System with Nonlocal Boundary Conditions and Localized Sources

This paper investigates the properties of solutions to a quasilinear parabolic system with nonlocal boundary conditions and localized sources. Conditions for the existence of global or blow-up solutions are given. Global blow-up property and blow- up rate estimates are also derived.     

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.     

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.     

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

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.     

A TAILORED FINITE POINT METHOD FOR THE HELMHOLTZ EQUATION WITH HIGH WAVE NUMBERS IN HETEROGENEOUS MEDIUM

In this paper, we propose a tailored-finite-point method for the numerical simulation of the Helmholtz equation with high wave numbers in heterogeneous medium. Our finite point method has been tailored to some particular properties of the problem, which allows us to obtain approximate solutions with the same behaviors as that of the exact solution very naturally. Especially, when the coefficients are piecewise constant, we can get the exact solution with only one point in each subdomain. Our finite-point method has uniformly convergent rate with respect to wave number k in L^2-norm.     

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.     

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.     

A POSTERIORI ESTIMATOR OF NONCONFORMING FINITE ELEMENT METHOD FOR FOURTH ORDER ELLIPTIC PERTURBATION PROBLEMS

In this paper, we consider the nonconforming finite element approximations of fourth order elliptic perturbation problems in two dimensions. We present an a posteriori error estimator under certain conditions, and give an h-version adaptive algorithm based on the error estimation. The local behavior of the estimator is analyzed as well. This estimator works for several nonconforming methods, such as the modified Morley method and the modified Zienkiewicz method, and under some assumptions, it is an optimal one. Numerical examples are reported, with a linear stationary Cahn-HiUiard-type equation as a model problem.     

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.     

UNIFORMLY CONVERGENT NONCONFORMING ELEMENT FOR 3-D FOURTH ORDER ELLIPTIC SINGULAR PERTURBATION PROBLEM

In this paper, using a bubble function, we construct a cuboid element to solve the fourth order elliptic singular perturbation problem in three dimensions. We prove that the nonconforming CO-cuboid element converges in the energy norm uniformly with respect to the perturbation parameter.     

A NUMERICAL STUDY FOR THE PERFORMANCE OF THE WENO SCHEMES BASED ON DIFFERENT NUMERICAL FLUXES FOR THE SHALLOW WATER EQUATIONS

In this paper we investigate the performance of the weighted essential non-oscillatory （WENO） methods based on different numerical fluxes, with the objective of obtaining better performance for the shallow water equations by choosing suitable numerical fluxes. We consider six numerical fluxes, i.e., Lax-Friedrichs, local Lax-Friedrichs, Engquist-Osher, Harten-Lax-van Leer, HLLC and the first-order centered fluxes, with the WENO finite volume method and TVD Runge-Kutta time discretization for the shallow water equations. The detailed numerical study is performed for both one-dimensional and two-dimensional shallow water equations by addressing the property, and resolution of discontinuities. issues of CPU cost, accuracy, non-oscillatory     

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.     

Contact Finite Determinacy of Relative Map Germs

The strong contact finite determinacy of relative map germs is studied by means of classical singularity theory. We first give the definition of a strong relative contact equivalence （or KS,T equivalence） and then prove two theorems which can be used to distinguish the contact finite determinacy of relative map germs, that is, f is finite determined relative to KS,T if and only if there exists a positive integer k, such that M^5（n）ε（S;n）p ∩→ TKs,T（f）.     

Reducing Subspaces of Toeplitz Operators on N_φ-type Quotient Modules on the Torus

In this paper, we prove that the Toeplitz operator with finite Blaschke product symbol Sψ（z） on Nφ has at least m non-trivial minimal reducing subspaces, where m is the dimension of H^2（Гω）⊙φ（ω）H^2（Гω）. Moreover, the restriction of Sψ（z） on any of these minimal reducing subspaces is unitary equivalent to the Bergman shift Mz.