Analysis of augmented three-field macro-hybrid mixed finite element schemes

On the basis of composition duality principles, augmented three-field macro- hybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladyenskaja-Babuka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decomp...     

Composition Duality Methods for Mixed Variational Inclusions

Composition duality methods for mixed variational inclusions are studied in a functional framework of reflexive Banach spaces. On the basis of duality principles, the solvability of maximal monotone and subdifferential mixed variational inclusions is established. For computational purposes, mass-preconditioned augmented formulations are introduced for regularization, as well as three-field and macro-hybrid variational versions. At a finite-dimensional level, corresponding discrete mixed and macro-hybrid internal approximations are discussed, as well as proximal-point iterative algorithms. Primal and dual mixed variational inclusions from contact mechanics illustrate the theory.     

Proximal penalty-duality algorithms for mixed optimality conditions

Preconditioned proximal penalty-duality two- and three-field algorithms for mixed optimality conditions, of evolution mixed constrained optimal control problems, are considered. Fixed point existence analysis is performed for corresponding evolution mixed governing variational state systems, in reflexive Banach spaces. Further, convergence analysis of the proximal penalty-duality algorithms is established via fixed point characterizations. In both analysis, a resolvent fixed point variational strategy is applied.     

Augmented Macro-Hybrid Mixed Finite Element Schemes for Elastic Contact Problems

On a setting of subdifferential models, variational augmented macro-hybrid mixed finite element schemes are formulated and analyzed for elastic unilateral contact problems with prescribed friction. Composition duality principles determine primal and dual mixed solvability, adopting coupling surjectivity for dualization. Macro-hybridization corresponds to nonoverlapping decompositions of elastic solid body systems, with displacement continuity and traction equilibrium transmission conditions dualized. In general, traction and displacement multipliers synchronize sub-bodies through nonmatching finite element interfaces. Three-field formulations give the basis for variational augmentation, in a sense of exact penalization, allowing speed-up of rates of convergence as well as proximation procedures of parallel numerical resolution algorithms.     

Composition duality methods for quasistatic evolution elastoviscoplastic variational problems

Composition duality methods for dual quasistatic evolution elastoviscoplastic variational problems are studied. Dual evolution mixed analysis is performed, as well as corresponding primal static mixed analysis. For multi-constitutive modeling and parallel computing, macro-hybrid variational formulations are further considered at the continuous level.     

On Gabay’s algorithms for mixed variational inequalities

Iterative numerical algorithms for variational inequalities are systematically constructed from fixed-point problem characterizations in terms of resolvent operators. The applied method is the one introduced by Gabay in [Ga], used here in the context of discrete variational inequalities, and with the emphasis on mixed finite element models. The algorithms apply to nonnecessarily potential problems, generalizing primal and mixed Uzawa and augmented Lagrangian-type algorithms. They are also identified with Euler and operator splitting methods for the time discretization of evolution first-order problems.     

Finite-element formulations of shear localization in granular materials

Displacement and mixed finite element formulations of shear localization in granular materials are presented. The formulations are based on hypoplastic constitutive laws for soils and the mixed-enhanced treatment involving displacement, strain and stress rates as independently varied fields. Included in these formulations are the standard displacement method, the three-field mixed formulation, the method of incompatible modes, the enhanced assumed strain method and the mixed enhanced strain method. Several numerical examples demonstrating the capability and performance of the different finite element formulations are presented. The numerical results are compared with available experimental data of Hostun RF sand and numerical results of Karlsruhe sand on biaxial tests. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Steady filtration problems with seawater intrusion: macro-hybrid primal mixed variational analysis

Macro-hybrid penalized, primal mixed continuous and discrete variational formulations, for steady filtration problems, with seawater intrusion, are studied. This is a mixed version of our previous paper on macro-hybrid penalized approximations (G. Alduncin, J. Esquivel-Avila, and N. Vera-Guzman, Steady filtration problems with seawater intrusion: Macro-hybrid penalized finite element approximations, Int. J. Numer. Meth. Fluids (2005) 49, pp. 935–957). Penalized pressure–velocity mixed variational formulations are introduced for non-overlapping domain decompositions, with vertical interfaces, of sections of coastal aquifers. Well-posedness and stability conditions are established at the continuous and discrete levels. Macro-hybrid primal mixed internal approximations are defined on independent subdomain grids, with transmission conditions imposed in a dual variational sense. Parallel relaxation penalty-duality algorithms are discussed from fixed-point characterizations, for iterative numerical resolution.     

Dual finite element analysis for some elliptic variational equations and inequalities

In some boundary-value problems the gradient or the cogradient of the solution is more important than the solution itself. Dual variational formulation of elliptic problems is utilized to define finiteelement approximations of the cogradient. A priori error estimates are presented for a class of second-order elliptic problems, including problems of elastostatics. If the boundary conditions are classical (i.e., of Dirichlet, Neumann. Newton, or mixed type), the primal and dual formulations are equivalent with variational equations, whereas the unilateral boundary conditions lead to variational inequalities. The paper has a surveyable character.     

Lagrange Multipliers and mixed formulations for some inequality constrained variational inequalities and some nonsmooth unilateral problems

Abstract

This paper addresses a class of inequality constrained variational inequalities and nonsmooth unilateral variational problems. We present mixed formulations arising from Lagrange multipliers. First we treat in a reflexive Banach space setting the canonical case of a variational inequality that has as essential ingredients a bilinear form and a non-differentiable sublinear, hence convex functional and linear inequality constraints defined by a convex cone. We extend the famous Brezzi splitting theorem that originally covers saddle point problems with equality constraints, only, to these nonsmooth problems and obtain independent Lagrange multipliers in the subdifferential of the convex functional and in the ordering cone of the inequality constraints. For illustration of the theory we provide and investigate an example of a scalar nonsmooth boundary value problem that models frictional unilateral contact problems in linear elastostatics. Finally we discuss how this approach to mixed formulations can be further extended to variational problems with nonlinear operators and equilibrium problems, and moreover, to hemivariational inequalities.     

A NEW SELF-ADAPTIVE ITERATIVE METHOD FOR GENERAL MIXED QUASI VARIATIONAL INEQUALITIES

The general mixed quasi variational inequality containing a nonlinear term φ is a useful and an important generalization of variational inequalities. The projection method can not be applied to solve this problem due to the presence of nonlinear term. It is well known that the variational inequalities involving the nonlinear term φ are equivalent to the fixed point problems and resolvent equations. In this article, the authors use these alternative equivalent formulations to suggest and analyze a new self-adaptive iterative method for solving general mixed quasi variational inequalities. Global convergence of the new method is proved. An example is given to illustrate the efficiency of the proposed method.     

Advanced finite element formulations for modeling thin piezoelectric structures

This contribution is concerned with mixed finite element formulations for modeling piezoelectric beam and shell structures. Due to the electromechanical coupling, specific deformation modes are joined with electric field components. In bending dominated problems incompatible approximation functions of these fields cause incorrect results. These effects occur in standard finite element formulations, where interpolation functions of lowest order are used. A mixed variational approach is introduced to overcome these problems. The mixed formulation allows for a consistent approximation of the electromechanical coupled problem. It utilizes six independent fields and could be derived from a Hu-Washizu variational principle. Displacements, rotations and the electric potential are employed as nodal degrees of freedom. According to the Timoshenko theory (beam) and the Reissner-Mindlin theory (shell), the formulations account for constant transversal shear strains. To incorporate three dimensional constitutive relations all transversal components of the electric field and the strain field are enriched by mixed finite element interpolations. Thus the complete piezoelectric coupling is appropriately captured. The common assumption of vanishing transversal stress and dielectric displacement components is enforced in an integral sense. Some numerical examples will demonstrate the capability of the presented finite element formulation. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Variational formulations of nonlinear constrained boundary value problems

Variational formulations of nonlinear constrained boundary value problems in reflexive Banach spaces are discussed from a compositional duality approach. The mixed variational compatibility conditions of the theory correspond to the surjectivity of the primal coupling boundary and interior operators.     

OPTIMAL MIXED H-P FINITE ELEMENT METHODS FOR STOKES AND NON-NEWTONIAN FLOW

1. IntroductionMOtivated,,by some advalltages of h -- p FEM over the classic FEM uncovered by recelltcomputation works(see-[19]), Schwab and Sari [11] have considered the ~d h--p lhate elementmethod for Non-Newtonian flow based upon a three-field Stokes formulation emallating fromMarization of some dmerellt models of Non-Newtoulan flow, in which, stress, velocity andPressure are coupled. Theoretical analysis and tailored numerical expert~8 show that theabed h -- p finite elemellt method …     

An augmented mixed finite element method for 3D linear elasticity problems

In this paper we introduce and analyze a new augmented mixed finite element method for linear elasticity problems in 3D. Our approach is an extension of a technique developed recently for plane elasticity, which is based on the introduction of consistent terms of Galerkin least-squares type. We consider non-homogeneous and homogeneous Dirichlet boundary conditions and prove that the resulting augmented variational formulations lead to strongly coercive bilinear forms. In this way, the associated Galerkin schemes become well posed for arbitrary choices of the corresponding finite element subspaces. In particular, Raviart-Thomas spaces of order 0 for the stress tensor, continuous piecewise linear elements for the displacement, and piecewise constants for the rotation can be utilized. Moreover, we show that in this case the number of unknowns behaves approximately as 9.5 times the number of elements (tetrahedrons) of the triangulation, which is cheaper, by a factor of 3, than the classical PEERS in 3D. Several numerical results illustrating the good performance of the augmented schemes are provided.     

Numerical Solution of Viscoplastic Flow Problems by Augmented Lagrangians

This paper describes an application of augmented Lagrangiantechniques to the numerical solution of quasistatic flow problemsin incompressible viscoplasticity, focusing on cases where theinternal viscoplastic dissipation potential is not a differentiablefunction of the material deformation rate. The stresses of elasticorigin are neglected, and the variational formulation of theseproblems is approximated via low-order mixed finite elements,which reduces the original problems to the constrained minimizationof a convex, but possibly not differentiable functional. Convergenceresults are proved or recalled, both for the finite elementapproximation and for the augmented Lagrangian algorithm. Adetailed study of the local minimization problems which occurin the augmented Lagrangian decomposition of the above problemsis also presented, together with several numerical results.     

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.     

Smooth methods of multipliers for complementarity problems

This paper describes several methods for solving nonlinear complementarity problems. A general duality framework for pairs of monotone operators is developed and then applied to the monotone complementarity problem, obtaining primal, dual, and primal-dual formulations. We derive Bregman-function-based generalized proximal algorithms for each of these formulations, generating three classes of complementarity algorithms. The primal class is well-known. The dual class is new and constitutes a general collection of methods of multipliers, or augmented Lagrangian methods, for complementarity problems. In a special case, it corresponds to a class of variational inequality algorithms proposed by Gabay. By appropriate choice of Bregman function, the augmented Lagrangian subproblem in these methods can be made continuously differentiable. The primal-dual class of methods is entirely new and combines the best theoretical features of the primal and dual methods. Some preliminary computation shows that this class of algorithms is effective at solving many of the standard complementarity test problems. Received February 21, 1997 / Revised version received December 11, 1998? Published online May 12, 1999     

Optimal Control of Nonlinear Transport-Flow Mixed Variational Problems

Optimal control of nonlinear transport-flow mixed variational problems are studied qualitatively, and the solvability analysis of the transport and flow mixed state systems is performed on the basis of primal and dual evolution duality principles. Corresponding primal and dual mixed optimality conditions are established by the application of some fundamental perturbation conjugate duality results recently proposed [8 G. Alduncin (2013). Optimal control of evolution mixed variational inclusions. Applied Mathematics and Optimization 68:445–473.[Crossref], [Web of Science ®] , [Google Scholar]]. Further, for computational purposes, two- and three-field proximation penalty-duality algorithms in the resolution of the mixed optimality conditions are finally presented and discussed.     

A variational inequality approach to generalized two-phase Stefan problem in several space variables

Summary The paper offers a study of a broad class of multidimensional two- phase problems of Stefan type by means of variational inequality techniques. The problems for quasilinear equations of alternatively parabolic or mixed parabolicelliptic type, mixed type nonlinear conditions at the fixed lateral boundary, involving free boundary conditions corresponding to phase transitions of both first (latent heat positive) and second kind (latent heat equal to zero) are taken into consideration. Results concerning existence of weak solutions, their uniqueness and stability are established.The preparation of the paper was partially carried out while the author's visiting Istituto di Analisi Numerica del C.N.R., Pavia, due to support of the C.N.R.