Dual-mixed finite element approximation of Stokes and nonlinear Stokes problems using trace-free velocity gradients

In this work a finite element method for a dual-mixed approximation of Stokes and nonlinear Stokes problems is studied. The dual-mixed structure, which yields a twofold saddle point problem, arises in a formulation of this problem through the introduction of unknown variables with relevant physical meaning. The method approximates the velocity, its gradient, and the total stress tensor, but avoids the explicit computation of the pressure, which can be recovered through a simple postprocessing technique. This method improves an existing approach for these problems and uses Raviart-Thomas elements and discontinuous piecewise polynomials for approximating the unknowns. Existence, uniqueness, and error results for the method are given, and numerical experiments that exhibit the reduced computational cost of this approach are presented.     

Continuous-discontinuous finite element method for convection-diffusion problems with characteristic layers

We study convergence properties of a numerical method for convection-diffusion problems with characteristic layers on a layer-adapted mesh. The method couples standard Galerkin with an h-version of the nonsymmetric discontinuous Galerkin finite element method with bilinear elements. In an associated norm, we derive the error estimate as well as the supercloseness result that are uniform in the perturbation parameter. Applying a post-processing operator for the discontinuous Galerkin method, we construct a new numerical solution with enhanced convergence properties.     

An interpolation matched interface and boundary method for elliptic interface problems

An interpolation matched interface and boundary (IMIB) method with second-order accuracy is developed for elliptic interface problems on Cartesian grids, based on original MIB method proposed by Zhou et al. [Y. Zhou, G. Wei, On the fictious-domain and interpolation formulations of the matched interface and boundary method, J. Comput. Phys. 219 (2006) 228-246]. Explicit and symmetric finite difference formulas at irregular grid points are derived by virtue of the level set function. The difference scheme using IMIB method is shown to satisfy the discrete maximum principle for a certain class of problems. Rigorous error analyses are given for the IMIB method applied to one-dimensional (1D) problems with piecewise constant coefficients and two-dimensional (2D) problems with singular sources. Comparison functions are constructed to obtain a sharp error bound for 1D approximate solutions. Furthermore, we compare the ghost fluid method (GFM), immersed interface method (IIM), MIB and IMIB methods for 1D problems. Finally, numerical examples are provided to show the efficiency and robustness of the proposed method.     

A new fitted operator finite difference method to solve systems of evolutionary reaction-diffusion equations

Abstract

In recent years, fitted operator finite difference methods (FOFDMs) have been developed for numerous types of singularly perturbed ordinary differential equations. The construction of most of these methods differed though the final outcome remained similar. The most crucial aspect was how the difference operator was designed to approximate the differential operator in question. Very often the approaches for constructing these operators had limited scope in the sense that it was difficult to extend them to solve even simple one-dimensional singularly perturbed partial differential equations. However, in some of our most recent work, we have successfully designed a class of FOFDMs and extended them to solve singularly perturbed time-dependent partial differential equations. In this paper, we design and analyze a robust FOFDM to solve a system of coupled singularly perturbed parabolic reaction-diffusion equations. We use the backward Euler method for the semi-discretization in time. An FOFDM is then developed to solve the resulting set of boundary value problems. The proposed method is analyzed for convergence. Our method is uniformly convergent with order one and two, respectively, in time and space, with respect to the perturbation parameters. Some numerical experiments supporting the theoretical investigations are also presented.     

Finite element approximation of multi-scale elliptic problems using patches of elements

In this paper we present a method for the numerical solution of elliptic problems with multi-scale data using multiple levels of not necessarily nested grids. The method consists in calculating successive corrections to the solution in patches whose discretizations are not necessarily conforming. This paper provides proofs of the results published earlier (see C. R. Acad. Sci. Paris, Ser. I 337 (2003) 679–684), gives a generalization of the latter to more than two domains and contains extensive numerical illustrations. New results including the spectral analysis of the iteration operator and a numerical method to evaluate the constant of the strengthened Cauchy-Buniakowski-Schwarz inequality are presented. Supported by CTI Project 6437.1 IWS-IW.     

Error estimates of the lumped mass finite element method for semilinear elliptic problems

We are concerned with the semilinear elliptic problems. We first investigate the L²-error estimate for the lumped mass finite element method. We then use the cascadic multigrid method to solve the corresponding discrete problem. On the basis of the finite element error estimates, we prove the optimality of the proposed multigrid method. We also report some numerical results to support the theory.     

Interior penalty discontinuous approximations of convection–diffusion problems with parabolic layers

Summary A nonsymmetric discontinuous Galerkin finite element method with interior penalties is considered for two–dimensional convection–diffusion problems with regular and parabolic layers. On an anisotropic Shishkin–type mesh with bilinear elements we prove error estimates (uniformly in the perturbation parameter) in an integral norm associated with this method. On different types of interelement edges we derive the values of discontinuity–penalization parameters. Numerical experiments complement the theoretical results.     

Discontinuous Least-Squares finite element method for the Div-Curl problem

In this paper, we consider the div-curl problem posed on nonconvex polyhedral domains. We propose a least-squares method based on discontinuous elements with normal and tangential continuity across interior faces, as well as boundary conditions, weakly enforced through a properly designed least-squares functional. Discontinuous elements make it possible to take advantage of regularity of given data (divergence and curl of the solution) and obtain convergence also on nonconvex domains. In general, this is not possible in the least-squares method with standard continuous elements. We show that our method is stable, derive a priori error estimates, and present numerical examples illustrating the method.     

Analysis of the cell boundary element methods for convection dominated convection-diffusion equations

We introduce two kinds of the cell boundary element (CBE) methods for convection dominated convection-diffusion equations: one is the CBE method with the exact bubble function and the other with inexact bubble functions. The main focus of this paper is on inexact bubble CBE methods. For inexact bubble CBE methods we introduce a family of numerical methods depending on two parameters, one for control of interior layers and the other for outflow boundary layers. Stability and convergence analysis are provided and numerical tests for inexact bubble CBEs with various choices of parameters are presented.     

Superconvergent biquadratic finite volume element method for two-dimensional Poisson’s equations

In this paper, a kind of biquadratic finite volume element method is presented for two-dimensional Poisson’s equations by restricting the optimal stress points of biquadratic interpolation as the vertices of control volumes. The method can be effectively implemented by alternating direction technique. It is proved that the method has optimal energy norm error estimates. The superconvergence of numerical gradients at optimal stress points is discussed and it is proved that the method has also superconvergence displacement at nodal points by a modified dual argument technique. Finally, a numerical example verifies the theoretical results and illustrates the effectiveness of the method.     

Error estimators for advection–reaction–diffusion equations based on the solution of local problems

This paper deals with a posteriori error estimates for advection–reaction–diffusion equations. In particular, error estimators based on the solution of local problems are derived for a stabilized finite element method. These estimators are proved to be equivalent to the error, with equivalence constants eventually depending on the physical parameters. Numerical experiments illustrating the performance of this approach are reported.     

An adaptive stabilized finite element method for the generalized Stokes problem

In this work we present an adaptive strategy (based on an a posteriori error estimator) for a stabilized finite element method for the Stokes problem, with and without a reaction term. The hierarchical type estimator is based on the solution of local problems posed on appropriate finite dimensional spaces of bubble-like functions. An equivalence result between the norm of the finite element error and the estimator is given, where the dependence of the constants on the physics of the problem is explicited. Several numerical results confirming both the theoretical results and the good performance of the estimator are given.     

Characteristics difference method for incompressible miscible flow in porous media

Two-phase ,incompressible miscible flow in porous media is governed by a system ofnonlinear partial differential equations. The pressure equation ,which is e11iptic in appearance ,isdiseretizod by a standard five-points difference method, The concentration equation is treated byan impliclt finite difference method that appbes a form of the method of characterlstics to thetransport terms. A class of biquadlatle interpolation is introduced for the method of chracteristics.Convergence rate is proved to be O(△t h^2)。     

Block method for problems on L-shaped domains

An extremely accurate, exponentially convergent solution is presented for both symmetric and non-symmetric Laplacian problems on L-shaped domains by using one-block version of the block method (BM). A simple and highly accurate formula for computing the stress intensity factor is given. Comparisons with various results in the literature are included.     

A posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods

We consider some (anisotropic and piecewise constant) diffusion problems in domains of R², approximated by a discontinuous Galerkin method with polynomials of any fixed degree. We propose an a posteriori error estimator based on gradient recovery by averaging. It is shown that this estimator gives rise to an upper bound where the constant is one up to some additional terms that guarantee reliability. The lower bound is also established. Moreover these additional terms are negligible when the recovered gradient is superconvergent. The reliability and efficiency of the proposed estimator is confirmed by some numerical tests.     

A GLOBAL-LOCAL FINITE ELEMENT METHOD IN SPACE-TIME FOR A HYPERBOLIC PROBLEM

ABSTRACT

A finite element method for solving the wave equation with couples boundary conditions is presented. In this approach finite elements are applied globally with respect to space and simultaneously but locally with respect to time. This gives rise to a single-step method in time. The method is a practical and economic one and the numerical results obtained compare favourably with the available analytic solution.     

Superconvergence of mixed covolume method for elliptic problems on triangular grids

In this paper, we consider the superconvergence of a mixed covolume method on the quasi-uniform triangular grids for the variable coefficient-matrix Poisson equations. The superconvergence estimates between the solution of the mixed covolume method and that of the mixed finite element method have been obtained. With these superconvergence estimates, we establish the superconvergence estimates and the L^∞$L^{\infty }$-error estimates for the mixed covolume method for the elliptic problems. Based on the superconvergence of the mixed covolume method, under the condition that the triangulation is uniform, we construct a post-processing method for the approximate velocity which improves the order of approximation of the approximate velocity.     

Iterative schemes for mixed finite element methods with applications to elasticity and compressible flow problems

Summary Iterative schemes for mixed finite element methods are proposed and analyzed in two abstract formulations. The first one has applications to elliptic equations and incompressible fluid flow problems, while the second has applications to linear elasticity and compressible Stokes problems. These schemes are constructed through iteratively penalizing the mixed finite element scheme, of which iterated penalty method and augmented Lagrangian method are special cases. Convergence theorems are demonstrated in abstract formulations in Hilbert spaces, and applications to individual physical problems are considered as examples. Theoretical analysis and computational experiments both show that the proposed schemes have very fast convergence; a few iterations are normally enough to reduce the iterative error to a prescribed precision. Numerical examples with continuous and discontinuous coefficients are presented.     

DOMAIN DECOMPOSITION METHOD FOR PARABOLIC EQUATION ON GENERAL DOMAIN

The nonoxerlapping domain deoomposition method for parabolic partial differential equation on general domain is considered. A kind of domain decomposition that uses the finite element procedure ks given. The problem.over the domains can be implemented on parallel computer. Convergence analysis is also presented.     

Fourier analysis of schwarz alternating methods for piecewise hermite bicubic orthogonal spline collocation

The rates of convergence of two Schwarz alternating methods are analyzed for the iterative solution of a discrete problem which arises when orthogonal spline collocation with piecewise Hermite bicubics is applied to the Dirichlet problem for Poisson's equation on a rectangle. In the first method, the rectangle is divided into two overlapping subrectangles, while three overlapping subrectangles are used in the second method. Fourier analysis is used to obtain explicit formulas for the convergence factors by which theH ¹-norm of the errors is reduced in one iteration of the Schwarz methods. It is shown numerically that while these factors depend on the size of overlap, they are independent of the partition stepsize. Results of numerical experiments are presented which confirm the established rates of convergence of the Schwarz methods.This research was supported in part by funds from the National Science Foundation grant CCR-9103451.