Existence results for solutions of mixed tensor variational inequalities

By employing the notion of exceptional family of elements, we establish existence results for the mixed tensor variational inequalities. We show that the nonexistence of an exceptional family of elements is a sufficient condition for the solvability of mixed tensor variational inequality. For positive semidefinite mixed tensor variational inequalities, the nonexistence of an exceptional family of elements is proved to be an equivalent characterization of the nonemptiness of the solution sets. We derive several sufficient conditions of the nonemptiness and compactness of the solution sets for the mixed tensor variational inequalities with some special structured tensors. Finally, we show that the mixed tensor variational inequalities can be defined as a class of convex optimization problems.

     

A conforming mixed finite-element method for the coupling of fluid flow with porous media flow

Salim Meddahi ^{We consider a porous medium entirely enclosed within a fluidregion and present a well-posed conforming mixed finite-elementmethod for the corresponding coupled problem. The interfaceconditions refer to mass conservation, balance of normal forcesand the Beavers–Joseph–Saffman law, which yieldsthe introduction of the trace of the porous medium pressureas a suitable Lagrange multiplier. The finite-element subspacesdefining the discrete formulation employ Bernardi–Raugeland Raviart–Thomas elements for the velocities, piecewiseconstants for the pressures and continuous piecewise-linearelements for the Lagrange multiplier. We show stability, convergenceand a priori error estimates for the associated Galerkin scheme.Finally, we provide several numerical results illustrating thegood performance of the method and confirming the theoreticalrates of convergence.}     

A mixed nonconforming finite element for linear elasticity

This article considers a mixed finite element method for linear elasticity. It is based on a modified mixed formulation that enforces the continuity of the stress weakly by adding a jump term of the approximated stress on interior edges. The symmetric stress are approximated by nonconforming linear elements and the displacement by piecewise constants. We establish ??(h) error bound in the (broken) L² norm for the divergence of the stress and ??(h) error bound in the L² norm for both the displacement and the stress tensor. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

An expanded mixed covolume method for elliptic problems

We consider the mixed covolume method combining with the expanded mixed element for a system of first‐order partial differential equations resulting from the mixed formulation of a general self‐adjoint elliptic problem with a full diffusion tensor. The system can be used to model the transport of a contaminant carried by a flow in porous media. We use the lowest order Raviart‐Thomas mixed element space. We show the first‐order error estimate for the approximate solution in L² norm. We show the superconvergence both for pressure and velocity in certain discrete norms. We also get a finite difference scheme by using proper approximate integration formulas. Finally we give some numerical examples. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A family of mixed finite elements for the elasticity problem

Summary A new mixed finite element formulation for the equations of linear elasticity is considered. In the formulation the variables approximated are the displacement, the unsymmetric stress tensor and the rotation. The rotation act as a Lagrange multiplier introduced in order to enforce the symmetry of the stress tensor. Based on this formulation a new family of both two-and three-dimensional mixed methods is defined. Optimal error estimates, which are valid uniformly with respect to the Poisson ratio, are derived. Finally, a new postprocessing scheme for improving the displacement is introduced and analyzed.     

Dual mixed finite element methods for the elasticity problem with Lagrange multipliers

We study a dual mixed formulation of the elasticity system in a polygonal domain of the plane with mixed boundary conditions and its numerical approximation. The (essential) Neumann boundary conditions (or traction boundary condition) are imposed using a discontinuous Lagrange multiplier corresponding to the trace of the displacement field. Moreover, a strain tensor is introduced as a new unknown and its symmetry is relaxed, also by the use of a Lagrange multiplier (the rotation). The singular behaviour of the solution requires us to use refined meshes to restore optimal rates of convergence. Uniform error estimates in the Lamé coefficient λ$\lambda$ are obtained for large λ$\lambda$. The hybridization of the problem is performed and numerical tests are presented confirming our theoretical results.     

Numerical analysis of a mixed finite element method for a flow-transport problem

We consider a generic flow-transport system of partial differential equations, which has wide application in the waste disposal industry. The approximation of this system, using a finite element method for the brine, radionuclides, and heat combined with a mixed finite element method for the pressure and velocity, is analyzed. Optimal order error estimates in H¹ and L² are derived. The error analysis is given with no restriction on the diffusion tensor. That is, we have included the effects of molecular diffusion and dispersion. © 1996 John Wiley & Sons, Inc.     

A computationally efficient modification of mixed finite element methods for flow problems with full transmissivity tensors

The mixed finite element method for approximately solving flow equations in porous media has received a good deal of attention in the literature. The main idea is to solve for the head/pressure and fluid velocity (Darcy velocity) simultaneously to obtain a higher order approximation of the fluid velocity. In the case of a diagonal transmissivity tensor the algebraic equations resulting from the discretization can be reduced to a system of algebraic equations for the head/pressure variable alone. This reduction results in a smaller number of unknows to be solved for in an iterative method such as preconditioned conjugate gradient method. The fluid velocity is then obtained from an algebraic relationship. In the case of full transmissivity tensor, the algebraic reduction is more difficult. This paper investigates some algorithms resulting from the modification of the mixed finite element that take advantage of the mixed finite element method for the diagonal tensor case. The resulting schemes are more efficient implementations that maintain the same order of accuracy as the original schemes. © 1993 John Wiley & Sons, Inc.     

Exponentially convergent parallel algorithm for nonlinear eigenvalue problems

A. V. Klimenko, V. L. Makarov ^{A new algorithm for nonlinear eigenvalue problems is proposed.The numerical technique is based on a perturbation of the coefficientsof differential equation combined with the Adomian decompositionmethod for the nonlinear part. The approach provides an exponentialconvergence rate with a base which is inversely proportionalto the index of the eigenvalue under consideration. The eigenpairscan be computed in parallel. Numerical examples are presentedto support the theory. They are in good agreement with the spectralasymptotics obtained by other authors.}     

Symmetric mixed covolume methods for parabolic problems

We present a mixed covolume method for a system of first order partial differential equations resulting from the mixed formulation of the general self‐adjoint parabolic problem with a variable nondiagonal diffusion tensor. The lowest order Raviart‐Thomas mixed element space on rectangles is used. We prove the first order optimal rate of convergence for approximate pressure as well as for approximate velocity. We also prove the second order superconvergence both for approximate velocity and pressure in certain discrete norms. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 561–583, 2002     

Stabilized dual-mixed method for the problem of linear elasticity with mixed boundary conditions

We extend the applicability of the augmented dual-mixed method introduced recently in Gatica (2007), Gatica et al. (2009) to the problem of linear elasticity with mixed boundary conditions. The method is based on the Hellinger–Reissner principle and the symmetry of the stress tensor is imposed in a weak sense. The Neumann boundary condition is prescribed in the finite element space. Then, suitable Galerkin least-squares type terms are added in order to obtain an augmented variational formulation which is coercive in the whole space. This allows to use any finite element subspaces to approximate the displacement, the Cauchy stress tensor and the rotation.     

A priori error estimation for the dual mixed finite element method of the elastodynamic problem in a polygonal domain, I

In this paper we analyze a new dual mixed formulation of the elastodynamic system in polygonal domains. In this formulation the symmetry of the strain tensor is relaxed by the rotation of the displacement. For the time discretization of this new dual mixed formulation, we use an explicit scheme. After the analysis of stability of the fully discrete scheme, L^∞ in time, L² in space a priori error estimates are derived for the approximation of the displacement, the strain, the pressure and the rotation. Numerical experiments confirm our theoretical predictions.     

A cell filtration of mixed tensor space

We construct a cellular basis of the walled Brauer algebra which has similar properties as the Murphy basis of the group algebra of the symmetric group. In particular, the restriction of a cell module to a certain subalgebra can be easily described via this basis. Furthermore, the mixed tensor space possesses a filtration by cell modules—although not by cell modules of the walled Brauer algebra itself, but by cell modules of its image in the algebra of endomorphisms of mixed tensor space.     

Analysis of expanded mixed methods for fourth-order elliptic problems

The recently proposed expanded mixed formulation for numerical solution of second-order elliptic problems is here extended to fourth-order elliptic problems. This expanded formulation for the differential problems under consideration differs from the classical formulation in that three variables are treated, i.e., the displacement, the stress, and the moment tensors. It works for the case where the coefficient of the differential equations is small and does not need to be inverted, or for the case in which the stress tensor of the equations does not need to be symmetric. Based on this new formulation, various mixed finite elements for fourth-order problems are considered; error estimates of quasi-optimal or optimal order depending upon the mixed elements are derived. Implementation techniques for solving the linear system arising from these expanded mixed methods are discussed, and numerical results are presented. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 483–503, 1997     

Adaptive frame methods for elliptic operator equations: the steepest descent approach

Massimo Fornasier^¶ Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate, Università "La Sapienza" in Roma, Via Antonio Scarpa, 16/B, I-00161 Roma, Italy Rob Stevenson^|| Department of Mathematics, Utrecht University, PO Box 80.010, NL-3508 TA Utrecht, The Netherlands ^{This paper is concerned with the development of adaptive numericalmethods for elliptic operator equations. We are particularlyinterested in discretization schemes based on wavelet frames.We show that by using three basic subroutines an implementable,convergent scheme can be derived, which, moreover, has optimalcomputational complexity. The scheme is based on adaptive steepestdescent iterations. We illustrate our findings by numericalresults for the computation of solutions of the Poisson equationwith limited Sobolev smoothness on intervals in 1D and L-shapeddomains in 2D.}     

A new mixed finite element method for the Stokes problem

We propose a new mixed formulation of the Stokes problem where the extra stress tensor is considered. Based on such a formulation, a mixed finite element is constructed and analyzed. This new finite element has properties analogous to the finite volume methods, namely, the local conservation of the momentum and the mass. Optimal error estimates are derived. For the numerical implementation of this finite element, a hybrid form is presented. This work is a first step towards the treatment of viscoelastic fluid flows by mixed finite element methods.