Weak formulation

From HandWiki
(Redirected from Lax–Milgram theorem)

Weak formulations are important tools for the analysis of mathematical equations that permit the transfer of concepts of linear algebra to solve problems in other fields such as partial differential equations. In a weak formulation, equations or conditions are no longer required to hold absolutely (and this is not even well defined) and has instead weak solutions only with respect to certain "test vectors" or "test functions". In a strong formulation, the solution space is constructed such that these equations or conditions are already fulfilled. The Lax–Milgram theorem, named after Peter Lax and Arthur Milgram who proved it in 1954, provides weak formulations for certain systems on Hilbert spaces.

General concept

Let [math]\displaystyle{ V }[/math] be a Banach space, let [math]\displaystyle{ V' }[/math] be the dual space of [math]\displaystyle{ V }[/math], let [math]\displaystyle{ A\colon V \to V' }[/math],[clarification needed] and let [math]\displaystyle{ f \in V' }[/math]. A vector [math]\displaystyle{ u \in V }[/math] is a solution of the equation

[math]\displaystyle{ Au = f }[/math]

if and only if for all [math]\displaystyle{ v \in V }[/math],

[math]\displaystyle{ [Au](v) = f(v). }[/math]

Here, [math]\displaystyle{ v }[/math] is called a test vector (in general) or a test function (if [math]\displaystyle{ V }[/math] is a function space).

To bring this into the generic form of a weak formulation, find [math]\displaystyle{ u\in V }[/math] such that

[math]\displaystyle{ a(u,v) = f(v) \quad \forall v \in V, }[/math]

by defining the bilinear form

[math]\displaystyle{ a(u,v) := [Au](v). }[/math]

Example 1: linear system of equations

Now, let [math]\displaystyle{ V = \mathbb R^n }[/math] and [math]\displaystyle{ A:V \to V }[/math] be a linear mapping. Then, the weak formulation of the equation

[math]\displaystyle{ Au = f }[/math]

involves finding [math]\displaystyle{ u\in V }[/math] such that for all [math]\displaystyle{ v \in V }[/math] the following equation holds:

[math]\displaystyle{ \langle Au,v \rangle = \langle f,v \rangle, }[/math]

where [math]\displaystyle{ \langle \cdot,\cdot \rangle }[/math] denotes an inner product.

Since [math]\displaystyle{ A }[/math] is a linear mapping, it is sufficient to test with basis vectors, and we get

[math]\displaystyle{ \langle Au,e_i\rangle = \langle f,e_i\rangle, \quad i=1,\ldots,n. }[/math]

Actually, expanding [math]\displaystyle{ u = \sum_{j=1}^n u_je_j }[/math], we obtain the matrix form of the equation

[math]\displaystyle{ \mathbf{A}\mathbf{u} = \mathbf{f}, }[/math]

where [math]\displaystyle{ a_{ij} = \langle Ae_j, e_i\rangle }[/math] and [math]\displaystyle{ f_i = \langle f,e_i \rangle }[/math].

The bilinear form associated to this weak formulation is

[math]\displaystyle{ a(u,v) = \mathbf{v}^T\mathbf{A} \mathbf{u}. }[/math]

Example 2: Poisson's equation

To solve Poisson's equation

[math]\displaystyle{ -\nabla^2 u = f, }[/math]

on a domain [math]\displaystyle{ \Omega\subset \mathbb R^d }[/math] with [math]\displaystyle{ u=0 }[/math] on its boundary, and to specify the solution space [math]\displaystyle{ V }[/math] later, one can use the [math]\displaystyle{ L^2 }[/math]-scalar product

[math]\displaystyle{ \langle u,v\rangle = \int_\Omega uv\,dx }[/math]

to derive the weak formulation. Then, testing with differentiable functions [math]\displaystyle{ v }[/math] yields

[math]\displaystyle{ -\int_\Omega ( \nabla^2 u ) v \,dx = \int_\Omega fv \,dx. }[/math]

The left side of this equation can be made more symmetric by integration by parts using Green's identity and assuming that [math]\displaystyle{ v=0 }[/math] on [math]\displaystyle{ \partial\Omega }[/math]:

[math]\displaystyle{ \int_\Omega \nabla u \cdot \nabla v \,dx = \int_\Omega f v \,dx. }[/math]

This is what is usually called the weak formulation of Poisson's equation. Functions in the solution space [math]\displaystyle{ V }[/math] must be zero on the boundary, and have square-integrable derivatives. The appropriate space to satisfy these requirements is the Sobolev space [math]\displaystyle{ H^1_0(\Omega) }[/math] of functions with weak derivatives in [math]\displaystyle{ L^2(\Omega) }[/math] and with zero boundary conditions, so [math]\displaystyle{ V = H^1_0(\Omega) }[/math].

The generic form is obtained by assigning

[math]\displaystyle{ a(u,v) = \int_\Omega \nabla u \cdot \nabla v \,dx }[/math]

and

[math]\displaystyle{ f(v) = \int_\Omega f v \,dx. }[/math]

The Lax–Milgram theorem

This is a formulation of the Lax–Milgram theorem which relies on properties of the symmetric part of the bilinear form. It is not the most general form.

Let [math]\displaystyle{ V }[/math] be a Hilbert space and [math]\displaystyle{ a( \cdot ,\cdot ) }[/math] a bilinear form on [math]\displaystyle{ V }[/math], which is

  1. bounded: [math]\displaystyle{ |a(u,v)| \le C \|u\| \|v\|\,; }[/math] and
  2. coercive: [math]\displaystyle{ a(u,u) \ge c \|u\|^2\,. }[/math]

Then, for any [math]\displaystyle{ f\in V' }[/math], there is a unique solution [math]\displaystyle{ u\in V }[/math] to the equation

[math]\displaystyle{ a(u,v) = f(v) \quad \forall v \in V }[/math]

and it holds

[math]\displaystyle{ \|u\| \le \frac1c \|f\|_{V'}\,. }[/math]

Application to example 1

Here, application of the Lax–Milgram theorem is a stronger result than is needed.

  • Boundedness: all bilinear forms on [math]\displaystyle{ \R^n }[/math] are bounded. In particular, we have [math]\displaystyle{ |a(u,v)| \le \|A\|\,\|u\|\,\|v\| }[/math]
  • Coercivity: this actually means that the real parts of the eigenvalues of [math]\displaystyle{ A }[/math] are not smaller than [math]\displaystyle{ c }[/math]. Since this implies in particular that no eigenvalue is zero, the system is solvable.

Additionally, this yields the estimate [math]\displaystyle{ \|u\| \le \frac1c \|f\|, }[/math] where [math]\displaystyle{ c }[/math] is the minimal real part of an eigenvalue of [math]\displaystyle{ A }[/math].

Application to example 2

Here, choose [math]\displaystyle{ V = H^1_0(\Omega) }[/math] with the norm [math]\displaystyle{ \|v\|_V := \|\nabla v\|, }[/math]

where the norm on the right is the [math]\displaystyle{ L^2 }[/math]-norm on [math]\displaystyle{ \Omega }[/math] (this provides a true norm on [math]\displaystyle{ V }[/math] by the Poincaré inequality). But, we see that [math]\displaystyle{ |a(u,u)| = \|\nabla u\|^2 }[/math] and by the Cauchy–Schwarz inequality, [math]\displaystyle{ |a(u,v)| \le \|\nabla u\|\,\|\nabla v\| }[/math].

Therefore, for any [math]\displaystyle{ f \in [H^1_0(\Omega)]' }[/math], there is a unique solution [math]\displaystyle{ u\in V }[/math] of Poisson's equation and we have the estimate

[math]\displaystyle{ \|\nabla u\| \le \|f\|_{[H^1_0(\Omega)]'}. }[/math]

See also

References

  • Lax, Peter D.; Milgram, Arthur N. (1954), "Parabolic equations", Contributions to the theory of partial differential equations, Annals of Mathematics Studies, 33, Princeton, N. J.: Princeton University Press, pp. 167–190, doi:10.1515/9781400882182-010, ISBN 9781400882182 

External links