Variational inequality

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 Template:Harv and Template:Harv: the first papers of the theory were Template:Harv and Template:Harv, Template:Harv. Later on, Guido Stampacchia proved his generalization to the Lax–Milgram theorem in Template:Harv 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 Template:Harvnb reports: thus the theory become widely known throughout France. Also in 1965, Stampacchia and Jacques-Louis Lions extended earlier results of Template:Harv, announcing them in the paper Template:Harv: full proofs of their results appeared later in the paper Template:Harv.

Definition

Following Template:Harvtxt, the definition of a variational inequality is the following one.

Template:EquationRef Given a Banach space <math>\boldsymbol{E}</math>, a subset <math>\boldsymbol{K}</math> of <math>\boldsymbol{E}</math>, and a functional <math>F\colon \boldsymbol{K}\to \boldsymbol{E}^{\ast}</math> from <math>\boldsymbol{K}</math> to the dual space <math>\boldsymbol{E}^{\ast}</math> of the space <math>\boldsymbol{E}</math>, the variational inequality problem is the problem of solving for the variable <math>x</math> belonging to <math>\boldsymbol{K}</math> the following inequality:

<math>\langle F(x), y-x \rangle \geq 0\qquad\forall y \in \boldsymbol{K}</math>

where <math>\langle\cdot,\cdot\rangle\colon \boldsymbol{E}^{\ast}\times\boldsymbol{E}\to \mathbb{R}</math> 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 Template:Harvtxt: consider the problem of finding the minimal value of a differentiable function <math>f</math> over a closed interval <math>I = [a,b]</math>. Let <math>x^{\ast}</math> be a point in <math>I</math> where the minimum occurs. Three cases can occur:

1. if <math>a<x^{\ast}< b,</math> then <math>f^{\prime}(x^{\ast}) = 0;</math>
2. if <math>x^{\ast}=a,</math> then <math>f^{\prime}(x^{\ast}) \ge 0;</math>
3. if <math>x^{\ast}=b,</math> then <math>f^{\prime}(x^{\ast}) \le 0.</math>

These necessary conditions can be summarized as the problem of finding <math>x^{\ast}\in I</math> such that

<math>f^{\prime}(x^{\ast})(y-x^{\ast}) \geq 0\quad</math> for <math>\quad\forall y \in I.</math>

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 <math>\mathbb{R}^n</math> is the following: given a subset <math>K</math> of <math>\mathbb{R}^{n}</math> and a mapping <math>F\colon K\to\mathbb{R}^{n}</math>, the finite-dimensional variational inequality problem associated with <math>K</math> consist of finding a <math>n</math>-dimensional vector <math>x</math> belonging to <math>K</math> such that

<math>\langle F(x), y-x \rangle \geq 0\qquad\forall y \in K</math>

where <math>\langle\cdot,\cdot\rangle\colon\mathbb{R}^{n}\times\mathbb{R}^{n}\to\mathbb{R}</math> is the standard inner product on the vector space <math>\mathbb{R}^{n}</math>.

The variational inequality for the Signorini problem

File:Classical Signorini problem.svg

The classical Signorini problem: what will be the equilibrium configuration of the orange spherically shaped elastic body resting on the blue rigid frictionless plane?

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

<math>B(\boldsymbol{u},\boldsymbol{v} - \boldsymbol{u}) - F(\boldsymbol{v} - \boldsymbol{u}) \geq 0 \qquad \forall \boldsymbol{v} \in \mathcal{U}_\Sigma </math>

where <math>B(\boldsymbol{u},\boldsymbol{v}) </math> and <math>F(\boldsymbol{v}) </math> are the following functionals, written using the Einstein notation

<math>B(\boldsymbol{u},\boldsymbol{v}) = -\int_A \sigma_{ik}(\boldsymbol{u})\varepsilon_{ik}(\boldsymbol{v})\,\mathrm{d}x</math>,    <math>F(\boldsymbol{v}) = \int_A v_i f_i\,\mathrm{d}x + \int_{\partial A\setminus\Sigma}\!\!\!\!\! v_i g_i \,\mathrm{d}\sigma</math>,    <math>\boldsymbol{u},\boldsymbol{v} \in \mathcal{U}_\Sigma </math>

where, for all <math>\boldsymbol{x}\in A</math>,

• <math>\Sigma</math> is the contact surface (or more generally a contact set),
• <math>\boldsymbol{f}(\boldsymbol{x}) = \left( f_1(\boldsymbol{x}), f_2(\boldsymbol{x}), f_3(\boldsymbol{x}) \right)</math> is the body force applied to the body,
• <math>\boldsymbol{g}(\boldsymbol{x})=\left(g_1(\boldsymbol{x}),g_2(\boldsymbol{x}),g_3(\boldsymbol{x})\right)</math> is the surface force applied to <math>\partial A\!\setminus\!\Sigma</math>,
• <math>\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)</math> is the infinitesimal strain tensor,
• <math>\boldsymbol{\sigma}=\left(\sigma_{ik}\right)</math> is the Cauchy stress tensor, defined as

<math>\sigma_{ik}= - \frac{\partial W}{\partial \varepsilon_{ik}} \qquad\forall i,k=1,2,3</math>

where <math>W(\boldsymbol{\varepsilon})=a_{ikjh}(\boldsymbol{x})\varepsilon_{ik}\varepsilon_{jh}</math> is the elastic potential energy and <math>\boldsymbol{a}(\boldsymbol{x})=\left(a_{ikjh}(\boldsymbol{x})\right)</math> is the elasticity tensor.

See also

• Complementarity theory
• Differential variational inequality
• Extended Mathematical Programming for Equilibrium Problems
• Mathematical programming with equilibrium constraints
• Obstacle problem
• Projected dynamical system
• Signorini problem
• Unilateral contact

References

Historical references

• Template:Citation. An historical paper about the fruitful interaction of elasticity theory and mathematical analysis: the creation of the theory of variational inequalities by Gaetano Fichera is described in §5, pages 282–284.
• Template:Citation. A brief research survey describing the field of variational inequalities, precisely the sub-field of continuum mechanics problems with unilateral constraints.
• Template:Citation. The birth of the theory of variational inequalities remembered thirty years later (English translation of the title) is an historical paper describing the beginning of the theory of variational inequalities from the point of view of its founder.

Scientific works

• Template:Citation
• Template:Citation
• Template:Citation. A short research note announcing and describing (without proofs) the solution of the Signorini problem.
• Template:Citation. The first paper where an existence and uniqueness theorem for the Signorini problem is proved.
• Template:Citation. An English translation of Template:Harv.
• Template:Citation
• Template:Citation.
• Template:Citation, available at Gallica. Announcements of the results of paper Template:Harv.
• Template:Citation. An important paper, describing the abstract approach of the authors to the theory of variational inequalities.
• Template:Citation.
• Template:Citation, available at Gallica. The paper containing Stampacchia's generalization of the Lax–Milgram theorem.

External links

• Template:Springer
• Alessio Figalli, On global homogeneous solutions to the Signorini problem,