Discrete maximum principles for finite element solutions of nonlinear elliptic problems with mixed boundary conditions

Summary. One of the most important problems in numerical simulations is the preservation of qualitative properties of solutions of the mathematical models by computed approximations. For problems of elliptic type, one of the basic properties is the (continuous) maximum principle. In our work, we present several variants of the maximum principles and their discrete counterparts for (scalar) second-order nonlinear elliptic problems with mixed boundary conditions. The problems considered are numerically solved by the continuous piecewise linear finite element approximations built on simplicial meshes. Sufficient conditions providing the validity of the corresponding discrete maximum principles are presented. Geometrically, they mean that the employed meshes have to be of acute or nonobtuse type, depending of the type of the problem. Finally some examples of real-life problems, where the preservation of maximum principles plays an important role, are presented.The first author was supported by the Hungarian Research Fund OTKA under grant no. F034840The second author was supported by the Agora Center under Grant InBCT of TEKES, Finland, and by the Academy Research Fellowship no. 208628 from the Academy of FinlandMathematics Subject Classification (2000): 35B50, 35J65, 65N30, 65N50     

Discrete minimum and maximum principles for finite element approximations of non-monotone elliptic equations

Summary. Uniform lower and upper bounds for positive finite-element approximations to semilinear elliptic equations in several space dimensions subject to mixed Dirichlet-Neumann boundary conditions are derived. The main feature is that the non-linearity may be non-monotone and unbounded. The discrete minimum principle provides a positivity-preserving approximation if the discretization parameter is small enough and if some structure conditions on the non-linearity and the triangulation are assumed. The discrete maximum principle also holds for degenerate diffusion coefficients. The proofs are based on Stampacchias truncation technique and on a variational formulation. Both methods are settled on careful estimates on the truncation operator.Mathematics Subject Classification (2000): 65N30, 65N12     

Error bounds for finite element solutions of mildly nonlinear elliptic boundary value problems

Summary The equivalence in a Hilbert space of variational and weak formulations of linear elliptic boundary value problems is well known. This same equivalence is proved here for mildly nonlinear problems where the right hand side of the differential equation involves the solution function. A finite element approximation to the solution of the weak problem ina finite dimensional subspace of the original Hilbert space is defined. An inequality bounding the error in this approximation over all functions of the space is derived, and in particular this holds for an interpolant to the weak solution. Thus this inequality, together with previously known, interpolation error bounds, produces a bound on the finite element solution to this nonlinear problem. An example of a mildly nonlinear Poisson problem is given.     

A nonconforming mixed finite element for second-order elliptic problems

In a recent work, Hiptmair [Mathematisches Institut, M9404, 1994] has constructed and analyzed a family of nonconforming mixed finite elements for second-order elliptic problems. However, his analysis does not work on the lowest order elements. In this article, we show that it is possible to construct a nonconforming mixed finite element for the lowest order case. We prove the convergence and give estimates of optimal order for this finite element. Our proof is based on the use of the properties of the so-called nonconforming bubble function to control the consistency terms introduced by the nonconforming approximation. We further establish an equivalence between this mixed finite element and the nonconforming piecewise quadratic finite element of Fortin and Soulie [J. Numer. Methods Eng., 19, 505–520, 1983]. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 445–457, 1997     

Sign-changing solutions for some fourth-order nonlinear elliptic problems

In this paper, we consider the existence and multiplicity of sign-changing solutions for some fourth-order nonlinear elliptic problems and some existence and multiple are obtained. The weak solutions are sought by means of sign-changing critical theorems.     

NontriviaL solutions of some nonlinear elliptic problems

This paper is concerned with a class of semilinear elliptic Dirichlet problems approximating degenerate equations. The aim is to prove the existence of at least 4^k?1 nontrivial solutions when the degeneration set consists of k distinct connected components     

Existence of solutions for some fourth-order nonlinear elliptic problems

In this paper, we consider the existence of positive, negative and sign-changing solutions for some fourth order semilinear elliptic boundary value problems. We present new results on invariant sets of the gradient flows of the corresponding variational functionals. The structure of the invariant sets will be built into minimax procedures to construct the sign-changing solutions.     

Discrete maximum principles for FE solutions of nonstationary diffusion‐reaction problems with mixed boundary conditions

In this article, we derive and discuss sufficient conditions for providing validity of the discrete maximum principle for nonstationary diffusion‐reaction problems with mixed boundary conditions, solved by means of simplicial finite elements and the θ time discretization method. The theoretical analysis is supported by numerical experiments. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Multilevel iterative methods for mixed finite element discretizations of elliptic problems

Summary For solving second order elliptic problems discretized on a sequence of nested mixed finite element spaces nearly optimal iterative methods are proposed. The methods are within the general framework of the product (multiplicative) scheme for operators in a Hilbert space, proposed recently by Bramble, Pasciak, Wang, and Xu [5,6,26,27] and make use of certain multilevel decomposition of the corresponding spaces for the flux variable.     

Interior penalty preconditioners for mixed finite element approximations of elliptic problems

It is established that an interior penalty method applied to second-order elliptic problems gives rise to a local operator which is spectrally equivalent to the corresponding nonlocal operator arising from the mixed finite element method. This relation can be utilized in order to construct preconditioners for the discrete mixed system. As an example, a family of additive Schwarz preconditioners for these systems is constructed. Numerical examples which confirm the theoretical results are also presented.

     

Two-grid methods of finite element solutions for semi-linear elliptic interface problems

Numerical Algorithms - In this paper, we present two efficient two-grid algorithms for solving two-dimensional semi-linear elliptic interface problems using finite element method. To linearize the...     

Two-grid finite volume element method for linear and nonlinear elliptic problems

Two-grid finite volume element discretization techniques, based on two linear conforming finite element spaces on one coarse and one fine grid, are presented for the two-dimensional second-order non-selfadjoint and indefinite linear elliptic problems and the two-dimensional second-order nonlinear elliptic problems. With the proposed techniques, solving the non-selfadjoint and indefinite elliptic problem on the fine space is reduced into solving a symmetric and positive definite elliptic problem on the fine space and solving the non-selfadjoint and indefinite elliptic problem on a much smaller space; solving a nonlinear elliptic problem on the fine space is reduced into solving a linear problem on the fine space and solving the nonlinear elliptic problem on a much smaller space. Convergence estimates are derived to justify the efficiency of the proposed two-grid algorithms. A set of numerical examples are presented to confirm the estimates. The work is supported by the National Natural Science Foundation of China (Grant No: 10601045).     

A posteriori error estimates for mixed finite element approximations of elliptic problems

We derive residual based a posteriori error estimates of the flux in L ²-norm for a general class of mixed methods for elliptic problems. The estimate is applicable to standard mixed methods such as the Raviart–Thomas–Nedelec and Brezzi–Douglas–Marini elements, as well as stabilized methods such as the Galerkin-Least squares method. The element residual in the estimate employs an elementwise computable postprocessed approximation of the displacement which gives optimal order.     

Superconvergence of RT1 mixed finite element approximations for elliptic control problems

In this paper,we investigate the superconvergence property of the numerical solution to a quadratic elliptic control problem by using mixed finite element methods.The state and co-state are approximated by the order k=1 Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions.We prove the superconvergence error estimate of h3/2 in L2-norm between the approximated solution and the average L2 projection of the control.Moreover,by the postprocessing technique,a quadratic superconvergence result of the control is derived.     

Error estimates for mixed finite element methods for nonlinear parabolic problems

We consider a class of mixed finite element methods for nonlinear parabolic problems over a plane domain. The finite element spaces taken are Raviart-Thomas spaces of index k, k ? 0. We obtain optimal order L²- and almost optimal order L^∞-error estimates for the finite element solution and order optimal L²-error estimates for its gradient. We also derive the error estimates for the time derivatives of the solution. Our results extend those previously obtained by Johnson and Thomée for the corresponding linear problems with k ? 1.     

A mixed multiscale finite element method for elliptic problems with oscillating coefficients

The recently introduced multiscale finite element method for solving elliptic equations with oscillating coefficients is designed to capture the large-scale structure of the solutions without resolving all the fine-scale structures. Motivated by the numerical simulation of flow transport in highly heterogeneous porous media, we propose a mixed multiscale finite element method with an over-sampling technique for solving second order elliptic equations with rapidly oscillating coefficients. The multiscale finite element bases are constructed by locally solving Neumann boundary value problems. We provide a detailed convergence analysis of the method under the assumption that the oscillating coefficients are locally periodic. While such a simplifying assumption is not required by our method, it allows us to use homogenization theory to obtain the asymptotic structure of the solutions. Numerical experiments are carried out for flow transport in a porous medium with a random log-normal relative permeability to demonstrate the efficiency and accuracy of the proposed method.     

Discrete and continuous maximum principles for parabolic and elliptic operators

The major qualitative properties of linear parabolic and elliptic operators/PDEs are the different maximum principles (MPs). Another important property is the stabilization property (SP), which connects these two types of operators/PDEs. This means that under some assumptions the solution of the parabolic PDE tends to an equilibrium state when t→∞, which is the solution of the corresponding elliptic PDE. To solve PDEs we need to use some numerical methods, and it is a natural requirement that these qualitative properties are preserved on the discrete level. In this work we investigate this question when a two-level discrete mesh operator is used as the discrete model of the parabolic operator (which is a one-step numerical procedure for solving the parabolic PDE) and a matrix as a discrete elliptic operator (which is a linear algebraic system of equations for solving the elliptic PDE). We clarify the relation between the discrete parabolic maximum principle (DPMP), the discrete elliptic maximum principle (DEMP) and the discrete stabilization property (DSP). The main result is that the DPMP implies the DSP and the DEMP.     

Compactness method in the finite element theory of nonlinear elliptic problems

Summary Finite element approximation of a nonlinear elliptic pseudomonotone second-order boundary value problem in a bounded nonpolygonal domain with mixed Dirichlet-Neumann boundary conditions is studied. In the discretization we approximate the domain by a polygonal one, use linear conforming triangular elements and evaluate integrals by numerical quadratures. We prove the solvability of the discrete problem and on the basis of compactness properties of the corresponding operator (which is not monotone in general) we prove the convergence of approximate solutions to an exact weak solutionuH ¹ ). No additional assumption on the regularity of the exact solution is needed.     

An adaptive least-squares mixed finite element method for nonlinear parabolic problems

A least-squares mixed finite element method for nonlinear parabolic problems is investigated in terms of computational efficiency. An a posteriori error estimator, which is needed in an adaptive refinement algorithm, was composed with the least-squares functional, and a posteriori errors were effectively estimated.