Variational inequality explained

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

Definition

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

\boldsymbol{E}

, a subset

\boldsymbol{K}

of

\boldsymbol{E}

, and a functional

F\colon\boldsymbol{K}\to\boldsymbol{E}\ast

from

\boldsymbol{K}

to the dual space

\boldsymbol{E}\ast

of the space

\boldsymbol{E}

, the variational inequality problem is the problem of solving for the variable

x

belonging to

\boldsymbol{K}

the following inequality:

\langleF(x),y-x\rangle\geq0    \forally\in\boldsymbol{K}

where

\langle,\rangle\colon\boldsymbol{E}\ast x \boldsymbol{E}\to 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

f

over a closed interval

I=[a,b]

. Let

x\ast

be a point in

I

where the minimum occurs. Three cases can occur:
  1. if

a<x\ast<b,

then

f\prime(x\ast)=0;

  1. if

x\ast=a,

then

f\prime(x\ast)\ge0;

  1. if

x\ast=b,

then

f\prime(x\ast)\le0.

These necessary conditions can be summarized as the problem of finding

x\ast\inI

such that

f\prime(x\ast)(y-x\ast)\geq0  

for

\forally\inI.

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

Rn

is the following: given a subset

K

of

Rn

and a mapping

F\colonK\toRn

, the finite-dimensional variational inequality problem associated with

K

consist of finding a

n

-dimensional
vector

x

belonging to

K

such that

\langleF(x),y-x\rangle\geq0    \forally\inK

where

\langle,\rangle\colonRn x Rn\toR

is the standard inner product on the vector space

Rn

.

The variational inequality for the Signorini problem

\boldsymbol{u}(\boldsymbol{x}) =\left(u1(\boldsymbol{x}),u2(\boldsymbol{x}),u3(\boldsymbol{x})\right)

of an anisotropic non-homogeneous elastic body that lies in a subset

A

of the three-dimensional euclidean space whose boundary is

\partialA

, resting on a rigid frictionless surface and subject only to its mass forces. The solution

u

of the problem exists and is unique (under precise assumptions) in the set of admissible displacements

l{U}\Sigma

i.e. the set of displacement vectors satisfying the system of ambiguous boundary conditions if and only if

B(\boldsymbol{u},\boldsymbol{v}-\boldsymbol{u})-F(\boldsymbol{v}-\boldsymbol{u})\geq0    \forall\boldsymbol{v}\inl{U}\Sigma

where

B(\boldsymbol{u},\boldsymbol{v})

and

F(\boldsymbol{v})

are the following functionals, written using the Einstein notation

B(\boldsymbol{u},\boldsymbol{v})=-\intA\sigmaik(\boldsymbol{u})\varepsilonik(\boldsymbol{v})dx

,   

F(\boldsymbol{v})=\intAvifidx+\int\partialvigid\sigma

,   

\boldsymbol{u},\boldsymbol{v}\inl{U}\Sigma

where, for all

\boldsymbol{x}\inA

,

\Sigma

is the contact surface (or more generally a contact set),

\boldsymbol{f}(\boldsymbol{x})=\left(f1(\boldsymbol{x}),f2(\boldsymbol{x}),f3(\boldsymbol{x})\right)

is the body force applied to the body,

\boldsymbol{g}(\boldsymbol{x})=\left(g1(\boldsymbol{x}),g2(\boldsymbol{x}),g3(\boldsymbol{x})\right)

is the surface force applied to

\partialA\setminus\Sigma

,

\boldsymbol{\varepsilon}=\boldsymbol{\varepsilon}(\boldsymbol{u})=\left(\varepsilonik(\boldsymbol{u})\right)=\left(

1
2

\left(

\partialui
\partialxk

+

\partialuk
\partialxi

\right)\right)

is the infinitesimal strain tensor,

\boldsymbol{\sigma}=\left(\sigmaik\right)

is the Cauchy stress tensor, defined as

\sigmaik=-

\partialW
\partial\varepsilonik

   \foralli,k=1,2,3

where

W(\boldsymbol{\varepsilon})=aikjh(\boldsymbol{x})\varepsilonik\varepsilonjh

is the elastic potential energy and

\boldsymbol{a}(\boldsymbol{x})=\left(aikjh(\boldsymbol{x})\right)

is the elasticity tensor.

See also

References

Historical references

Scientific works

External links