In mathematics, a variational inequality is an inequality involving a functional, which has to be solved for all possible values of a given variable, belonging usually to a convex set. The mathematical theory of variational inequalities was initially developed to deal with equilibrium problems, precisely the Signorini problem: in that model problem, the functional involved was obtained as the first variation of the involved potential energy. Therefore, it has a variational origin, recalled by the name of the general abstract problem. The applicability of the theory has since been expanded to include problems from economics, finance, optimization and game theory.

History

The first problem involving a variational inequality was the Signorini problem, posed by Antonio Signorini in 1959 and solved by Gaetano Fichera in 1963, according to the references (Antman 1983, pp. 282–284) and (Fichera 1995): the first papers of the theory were (Fichera 1963) and (Fichera 1964a), (Fichera 1964b). Later on, Guido Stampacchia proved his generalization to the Lax–Milgram theorem in (Stampacchia 1964) in order to study the regularity problem for partial differential equations and coined the name "variational inequality" for all the problems involving inequalities of this kind. Georges Duvaut encouraged his graduate students to study and expand on Fichera's work, after attending a conference in Brixen on 1965 where Fichera presented his study of the Signorini problem, as Antman 1983, p. 283 reports: thus the theory become widely known throughout France. Also in 1965, Stampacchia and Jacques-Louis Lions extended earlier results of (Stampacchia 1964), announcing them in the paper (Lions & Stampacchia 1965): full proofs of their results appeared later in the paper (Lions & Stampacchia 1967).

Definition

Following Antman (1983, p. 283), the definition of a variational inequality is the following one.

Definition 1. Given a Banach space ${\displaystyle {\boldsymbol {E))}$, a subset ${\displaystyle {\boldsymbol {K))}$ of ${\displaystyle {\boldsymbol {E))}$, and a functional ${\displaystyle F\colon {\boldsymbol {K))\to {\boldsymbol {E))^{\ast ))$ from ${\displaystyle {\boldsymbol {K))}$ to the dual space ${\displaystyle {\boldsymbol {E))^{\ast ))$ of the space ${\displaystyle {\boldsymbol {E))}$, the variational inequality problem is the problem of solving for the variable ${\displaystyle x}$ belonging to ${\displaystyle {\boldsymbol {K))}$ the following inequality:

${\displaystyle \langle F(x),y-x\rangle \geq 0\qquad \forall y\in {\boldsymbol {K))}$

where ${\displaystyle \langle \cdot ,\cdot \rangle \colon {\boldsymbol {E))^{\ast }\times {\boldsymbol {E))\to \mathbb {R} }$ is the duality pairing.

In general, the variational inequality problem can be formulated on any finite – or infinite-dimensional Banach space. The three obvious steps in the study of the problem are the following ones:

1. Prove the existence of a solution: this step implies the mathematical correctness of the problem, showing that there is at least a solution.
2. Prove the uniqueness of the given solution: this step implies the physical correctness of the problem, showing that the solution can be used to represent a physical phenomenon. It is a particularly important step since most of the problems modeled by variational inequalities are of physical origin.
3. Find the solution or prove its regularity.

Examples

The problem of finding the minimal value of a real-valued function of real variable

This is a standard example problem, reported by Antman (1983, p. 283): consider the problem of finding the minimal value of a differentiable function ${\displaystyle f}$ over a closed interval ${\displaystyle I=[a,b]}$. Let ${\displaystyle x^{\ast ))$ be a point in ${\displaystyle I}$ where the minimum occurs. Three cases can occur:

1. if ${\displaystyle a<x^{\ast }<b,}$ then ${\displaystyle f^{\prime }(x^{\ast })=0;}$
2. if ${\displaystyle x^{\ast }=a,}$ then ${\displaystyle f^{\prime }(x^{\ast })\geq 0;}$
3. if ${\displaystyle x^{\ast }=b,}$ then ${\displaystyle f^{\prime }(x^{\ast })\leq 0.}$

These necessary conditions can be summarized as the problem of finding ${\displaystyle x^{\ast }\in I}$ such that

${\displaystyle f^{\prime }(x^{\ast })(y-x^{\ast })\geq 0\quad }$ for ${\displaystyle \quad \forall y\in I.}$

The absolute minimum must be searched between the solutions (if more than one) of the preceding inequality: note that the solution is a real number, therefore this is a finite dimensional variational inequality.

The general finite-dimensional variational inequality

A formulation of the general problem in ${\displaystyle \mathbb {R} ^{n))$ is the following: given a subset ${\displaystyle K}$ of ${\displaystyle \mathbb {R} ^{n))$ and a mapping ${\displaystyle F\colon K\to \mathbb {R} ^{n))$, the finite-dimensional variational inequality problem associated with ${\displaystyle K}$ consist of finding a ${\displaystyle n}$-dimensional vector ${\displaystyle x}$ belonging to ${\displaystyle K}$ such that

${\displaystyle \langle F(x),y-x\rangle \geq 0\qquad \forall y\in K}$

where ${\displaystyle \langle \cdot ,\cdot \rangle \colon \mathbb {R} ^{n}\times \mathbb {R} ^{n}\to \mathbb {R} }$ is the standard inner product on the vector space ${\displaystyle \mathbb {R} ^{n))$.

The variational inequality for the Signorini problem

In the historical survey (Fichera 1995), Gaetano Fichera describes the genesis of his solution to the Signorini problem: the problem consist in finding the elastic equilibrium configuration ${\displaystyle {\boldsymbol {u))({\boldsymbol {x)))=\left(u_{1}({\boldsymbol {x))),u_{2}({\boldsymbol {x))),u_{3}({\boldsymbol {x)))\right)}$ of an anisotropic non-homogeneous elastic body that lies in a subset ${\displaystyle A}$ of the three-dimensional euclidean space whose boundary is ${\displaystyle \partial A}$, resting on a rigid frictionless surface and subject only to its mass forces. The solution ${\displaystyle u}$ of the problem exists and is unique (under precise assumptions) in the set of admissible displacements ${\displaystyle {\mathcal {U))_{\Sigma ))$ i.e. the set of displacement vectors satisfying the system of ambiguous boundary conditions if and only if

${\displaystyle B({\boldsymbol {u)),{\boldsymbol {v))-{\boldsymbol {u)))-F({\boldsymbol {v))-{\boldsymbol {u)))\geq 0\qquad \forall {\boldsymbol {v))\in {\mathcal {U))_{\Sigma ))$

where ${\displaystyle B({\boldsymbol {u)),{\boldsymbol {v)))}$ and ${\displaystyle F({\boldsymbol {v)))}$ are the following functionals, written using the Einstein notation

${\displaystyle B({\boldsymbol {u)),{\boldsymbol {v)))=-\int _{A}\sigma _{ik}({\boldsymbol {u)))\varepsilon _{ik}({\boldsymbol {v)))\,\mathrm {d} x}$,    ${\displaystyle F({\boldsymbol {v)))=\int _{A}v_{i}f_{i}\,\mathrm {d} x+\int _{\partial A\setminus \Sigma }\!\!\!\!\!v_{i}g_{i}\,\mathrm {d} \sigma }$,    ${\displaystyle {\boldsymbol {u)),{\boldsymbol {v))\in {\mathcal {U))_{\Sigma ))$

where, for all ${\displaystyle {\boldsymbol {x))\in A}$,

• ${\displaystyle \Sigma }$ is the contact surface (or more generally a contact set),
• ${\displaystyle {\boldsymbol {f))({\boldsymbol {x)))=\left(f_{1}({\boldsymbol {x))),f_{2}({\boldsymbol {x))),f_{3}({\boldsymbol {x)))\right)}$ is the body force applied to the body,
• ${\displaystyle {\boldsymbol {g))({\boldsymbol {x)))=\left(g_{1}({\boldsymbol {x))),g_{2}({\boldsymbol {x))),g_{3}({\boldsymbol {x)))\right)}$ is the surface force applied to ${\displaystyle \partial A\!\setminus \!\Sigma }$,
• ${\displaystyle {\boldsymbol {\varepsilon ))={\boldsymbol {\varepsilon ))({\boldsymbol {u)))=\left(\varepsilon _{ik}({\boldsymbol {u)))\right)=\left({\frac {1}{2))\left({\frac {\partial u_{i)){\partial x_{k))}+{\frac {\partial u_{k)){\partial x_{i))}\right)\right)}$ is the infinitesimal strain tensor,
• ${\displaystyle {\boldsymbol {\sigma ))=\left(\sigma _{ik}\right)}$ is the Cauchy stress tensor, defined as

${\displaystyle \sigma _{ik}=-{\frac {\partial W}{\partial \varepsilon _{ik))}\qquad \forall i,k=1,2,3}$

where ${\displaystyle W({\boldsymbol {\varepsilon )))=a_{ikjh}({\boldsymbol {x)))\varepsilon _{ik}\varepsilon _{jh))$ is the elastic potential energy and ${\displaystyle {\boldsymbol {a))({\boldsymbol {x)))=\left(a_{ikjh}({\boldsymbol {x)))\right)}$ is the elasticity tensor.