Skip to main content

Application of a finite element method for variational inequalities

Abstract

In this paper we explore the application of a finite element method (FEM) to the inequality and Laplacian constrained variational optimization problems. First, we illustrate the connection between the optimization problem and elliptic variational inequalities; secondly, we prove the existence of the solution via the augmented Lagrangian multipliers method. A triangular type FEM is employed in the numerical calculations. Computational results indicate that the present finite element method is a highly efficient technique in these sorts of variational problems involving inequalities.

AMS Subject Classification: 35J86, 26D10.

1 Introduction

Nonlinear variational inequalities (VI) and optimal control problems are a significantly important class of problems having a broad range of applications in mathematics, engineering, mechanics and physics. Elasticity, computational fluid mechanics, computational mechanics are some of the particular applications of these types of problems.

In the present day, nonlinear variational problems or partial differential equation based optimal control problems are not only important but also very challenging problems which usually involve high storage requirements, hard optimization and solution techniques. Therefore, finding reliable and efficient computational and numerical techniques along with fast implementation methods for the solution of these types of problems is quite useful and one of active research areas in the subject. In this paper we explore the application of a finite element method (FEM) [1] to the inequality and Laplacian constrained variational optimization problems [2].

Structure of this paper can be summarized as follows. Section 2 overviews some related methods briefly. In Section 3 we express the model problem as an inequality and Laplacian constrained variational optimization problem. In the same section, we present a theorem which shows the connection between the optimization problem and the second kind variational inequality problems. We discretize the VI problems using a finite element method. We complete the paper with a conclusion section where we discuss the method presented in this paper and point some possible future extensions.

2 Related studies

In this section our goal is to present some related work where some type of FEM is employed as a numerical solution technique. Although there are a vast amount of related studies, we will try to mention only a few of those which have some particular and significant application areas and consider the ones which set up the connection between variational inequalities and optimal control problems besides employing FEM in numerical calculations. Before we start mentioning some of those related studies, we first illustrate an application of the present FEM to a specific variational inequality problem. Let V= H 0 1 (Ω), and

a(u,v)= ∫ Ω ∇u⋅∇vdΩ,

where

∇u⋅∇v= ∂ u ∂ x 1 ∂ v ∂ x 1 + ∂ u ∂ x 2 ∂ v ∂ x 2 ,

and L(v)=〈f,v〉 for f∈ V ∗ = H − 1 (Ω) and v∈V. Let Ψ∈ H 1 (Ω)∩ C 0 ( Ω ¯ ) and Ψ | Γ ≤0. Define K={v∈ H 0 1 (Ω):v≥ψ a.e. on Î©}. Then the variational inequality problem to which to present FEM is applied is defined as finding u such that

a(u,v−u)≥L(v−u)for all v∈K,u∈K.

Detailed information about this variational inequality problem and the application of the present FEM to this problem might be obtained from [3].

In [4] Barret et al. study the existence of a solution to an open problem about a quasi-variational inequality problem arising in a model for sand surface evolution. The authors first introduce a regularized mixed formulation involving both the primal and dual variables, and secondly study two methods for numerical solution of the problem. Detailed information about this study might be obtained in [4].

In [5] Burke et al. pose the Francfort-Marigo model of brittle fracture in terms of the minimization of a quite irregular energy functional. Employing a generalized Ambrosio-Tortorelli functional, the authors obtain a bound-constrained minimization problem, which can be expressed in terms of a variational inequality. By proposing a new FEM, the authors compute the local minimizer of the generalized functional.

In [6] Elliott et al. introduce a finite element approximation of a variational formulation of Bean’s model for the physical configuration of an infinitely long cylindrical superconductor subject to a transverse magnetic field. The authors prove an error between the exact solution and the approximate solution for the current density and the magnetic field in appropriate norms. In the paper some numerical simulations for a variety of applied magnetic fields are also presented.

In [7] Bank et al. consider the application of primal-dual interior methods to the optimization of systems arising in the finite-element discretization of a class of elliptic variational inequalities. These problems lead to very large optimization problems with upper and lower bound constraints. When interior methods are applied to the discretized problem, the resulting linear systems have the same zero/nonzero structure as the finite-element equations solved for the unconstrained case.

As we have pointed above, although it is highly possible to extend this list of related studies, we refer the interested reader to the literature for further related work within the present paper.

3 Main results

Let Ω be a bounded polygonal domain in R 2 . We consider a model problem defined as follows.

For given f∈ L 2 (Ω), find a scalar function ϕ on Ω being a minimizer of the cost functional

J(u):= 1 2 ∫ Ω ΔuΔudΩ+ ∫ Ω |∇u|dΩ− ∫ Ω fudΩ
(1)

over the convex set K defined by

K= { u ∈ H 0 1 ( Ω ) : | ∇ u | ≤ 1  almost everywhere on  Ω } ,

where the Sobolev space H 0 1 (Ω) is defined as follows:

Variational inequalities might be separated into two main groups as elliptic and parabolic variational inequalities. Glowinski studies these sorts of inequalities in [3] in detail. He considers the elliptic variational inequalities (EVI) as the first kind and second kind EVI and defines those in a functional context as follows.

Let V denote a real Hilbert space with the inner product (⋅,⋅) and the associated norm ∥⋅∥. V ∗ is a dual space of V, a(⋅,⋅):V×V→R is a bilinear, continuous and V-elliptic form on V×V, L:V→R continuous linear functional, K is a closed convex nonempty subset of V, j(⋅):V→R∪{∞} is a convex lower semi-continuous functional. The first and second kind EVI are typically defined in the following way.

The first kind EVI: find ϕ∈V such that ϕ is a solution of the problem

a(Ï•,v−ϕ)≥L(v−ϕ),every v∈K,ϕ∈K.

The second kind EVI: find ϕ∈V such that ϕ is a solution of the problem

a(Ï•,v−ϕ)+j(v)−j(Ï•)≥L(v−ϕ),every v∈V,ϕ∈V.

The next lemma sets up the connection between optimization and VI problems.

Lemma 1 [3]

Let b:V×V→R be a symmetric continuous bilinear V-elliptic form. Let L∈ V ∗ and j:V→R∪{∞} be a convex lower semi-continuous proper functional. Let

J(v)= 1 2 b(v,v)+j(v)−L(v).

Then the minimization problem, find Ï• such that

J(Ï•)≤J(v),every v∈V,ϕ∈V,

has a unique solution which is characterized by

b(Ï•,v−ϕ)+j(v)−j(Ï•)≥L(v−ϕ),every v∈V,ϕ∈V.

Proof The proof can be seen on p.7 of [3]. □

It is clear that we can define the following in terms of the notations of Lemma 1:

We will approximate (1) with a finite element method introduced in [3]. Assume that Ω is a polygonal domain of R 2 . Consider a triangulation F h of Ω in the following sense: F h is a finite set of triangles T such that

Here T i ∘ denotes the inner part of the corresponding triangle. Furthermore, for all T 1 , T 2 ∈ F h and T 1 ≠ T 2 , exactly one of the following conditions must hold:

(i) T 1 ∩ T 2 =∅,

(ii) T 1 and T 2 have only one common vertex,

(iii) T 1 and T 2 have only a whole common edge.

h is the length of the largest edge of the triangles in the triangulation. Define P k as a space of polynomials in x 1 and x 2 of degree less than or equal to k, and

∑ h :={P∈ Ω ¯ ,P is a vertex of T∈ F h }.

The space V= H 0 1 (Ω) is approximated by the family of subspaces ( V h k ) h with k=1 or k=2, where

V h k := { v h ∈ C 0 ( Ω ¯ ) , v h | ∂ Ω = 0  and  v h | T ∈ P k , ∀ T ∈ F h } ,k=1,2.

It is obvious that the V h k are finite dimensional. Then the space K is approximated by

K h k = { v h ∈ V h k , v h ( P ) ≥ ψ ( P ) , ∀ P ∈ ∑ h k } ,k=1,2.

Notice that K h k for k=1,2 are closed convex nonempty subsets of V h k .

With these settings, the solution ϕ∈K is approximated by

b( Ï• h ,u− Ï• h )+j(u)−j( Ï• h )≥L(u− Ï• h ),every u∈ K h ,ϕ∈ K h ,
(2)

or equivalently,

Using the augmented Lagrangian multipliers method, we find a discrete solution of (1) as follows. First, let us introduce the Lagrange functional

L ˆ (u,μ)= 1 2 b(u,u)+ ∫ Ω |∇u|dΩ− ∫ Ω fudΩ+ ∫ Ω μ ( | ∇ u | − 1 ) dΩ.

For r≥0, an augmented Lagrangian L r is defined by

L r (u,μ)= L ˆ (u,μ)+ r 2 ∫ Ω | ∇ u − 1 | 2 dΩ.
(3)

Augmented Lagrangian multipliers methods for VI problems have been introduced by Glowinski and Marrocco (see [8]). Theorem 2.1 on p.168 in [3] guarantees the existence of a solution of this optimization problem.

Let us do notice that the first component Ï• h of (2) is then the solution of the original problem (1). Using the techniques of the variational calculus, we can write that

( ( 1 + μ h ) Δ ϕ h , Δ ( v − ϕ h ) ) +j(v)−j( ϕ h )−L(u− ϕ h )= ( μ h , Δ ( v − ϕ h ) ) ∀ϕ∈ V h ,

or equivalently,

We can describe the solution algorithm as follows:

(i) Choose an initial iterate μ h and λ≥0;

(ii) Solve the linear problem ((1+ μ h )Δ ϕ h ,Δ(v− ϕ h ))+j(v)−j( ϕ h )−L(u− ϕ h )=( μ h ,Δ(v− ϕ h )), ∀ϕ∈ V h ;

(iii) Update μ h α + 1 =max{0, μ h +α(∇ u h α −1)} on each cell;

(iv) Set α=α+1 and go back to (ii).

4 Conclusion

In this note we have applied a finite element method to the inequality and Laplacian constrained variational optimization problems. Firstly, we illustrate the connection between the optimization problem and elliptic variational inequalities; and secondly, we solve the problem using the augmented Lagrangian multipliers method. In a future work, we are planning to apply the present method to some specific application areas including modeling with optimal control problems.

References

  1. Lu Z, Chen Y, Huang Y: A priori error estimates for higher order variational discretization and mixed finite element methods of optimal control problems. J. Inequal. Appl. 2012., 2012: Article ID 95

    Google Scholar 

  2. Xu H: Efficient implementation of a modified and relaxed hybrid steepest-descent method for a type of variational inequality. J. Inequal. Appl. 2012., 2012: Article ID 93

    Google Scholar 

  3. Glowinski R: Numerical Methods for Nonlinear Variational Problems. Springer, Berlin; 2008.

    MATH  Google Scholar 

  4. Barrett JW, Prigozhin L: A quasi-variational inequality problem arising in the modeling of growing sandpiles. ESAIM, Math. Model. Numer. Anal. 2012. (accepted). doi:10.1051/m2an/2012062

    Google Scholar 

  5. Burke, S, Ortner, C, Suli, E: An adaptive finite element approximation of a generalised Ambrosio-Tortorelli functional. http://www2.maths.ox.ac.uk/oxmos/reports/pdfs/oxmos29.pdf

  6. Elliott CM, Kay D, Styles V: Finite element approximation of a variational inequality formulation of Bean’s model for superconductivity. SIAM J. Numer. Anal. 2004, 42(3):1324–1341. 10.1137/S0036142902412324

    Article  MATH  MathSciNet  Google Scholar 

  7. Bank RE, Gill PE, Marcia RF: Interior Methods for a Class of Elliptic Variational Inequalities, Large-Scale PDE-Constrained Optimization. Springer, Berlin; 2003.

    Google Scholar 

  8. Glowinski R, Marrocco A: Numerical solution of two-dimensional magneto-static problems by augmented Lagrangian methods. Comput. Methods Appl. Mech. Eng. 1977, 12: 33–46. 10.1016/0045-7825(77)90049-4

    Article  MATH  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mehmet Ali Akinlar.

Additional information

Competing interests

The author declares that they have no competing interests.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Akinlar, M.A. Application of a finite element method for variational inequalities. J Inequal Appl 2013, 45 (2013). https://doi.org/10.1186/1029-242X-2013-45

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1029-242X-2013-45

Keywords