Weak formulation
A weak formulation is a representation of a differential equation in which derivatives are transferred from the unknown function to auxiliary test functions through integration by parts. The resulting identity is interpreted in an integral or distributional sense rather than at every point of the domain. Solutions satisfying such an identity are called weak solutions.
Weak formulations are central to the modern analysis of partial differential equations. They permit the treatment of functions that lack the classical derivatives required by the original equation, while retaining its boundary conditions and conservation structure. They also provide the mathematical basis for several approximation theories, including the finite element method.
Mathematical structure
Let (\Omega\subset\mathbb{R}^n) be a bounded domain, and consider the boundary-value problem
[ -\Delta u=f \quad \text{in }\Omega, \qquad u=0 \quad \text{on }\partial\Omega. ]
A classical solution has second derivatives satisfying the equation pointwise and possesses a boundary trace equal to zero. In the weak formulation, the equation is paired with a sufficiently regular test function (v) whose boundary trace vanishes. Green's identity gives
[ \int_\Omega \nabla u\cdot\nabla v,dx
\int_\Omega f v,dx. ]
The boundary term disappears because of the homogeneous Dirichlet boundary condition. The weak problem therefore consists of finding (u\in H_0^1(\Omega)) such that
[ a(u,v)=F(v) \qquad \text{for every }v\in H_0^1(\Omega), ]
where
[ a(u,v)=\int_\Omega \nabla u\cdot\nabla v,dx ]
is a bilinear form and
[ F(v)=\int_\Omega f v,dx ]
is a linear functional. The space (H_0^1(\Omega)) is a Sobolev space whose elements possess square-integrable first weak derivatives and an appropriate zero boundary trace.
This formulation lowers the differentiability demanded of (u). The pointwise equation formally involves second derivatives, whereas the integral identity contains only first derivatives. Equality is required against every admissible test function, so the reduction in differentiability does not amount to removing the differential equation. Instead, the equation is encoded by its action on a function space.
Weak derivatives and distributions
The notion of weak formulation depends on the weak derivative. A locally integrable function (u) has weak derivative (D^\alpha u) when
[ \int_\Omega u,D^\alpha\varphi,dx
(-1)^{|\alpha|} \int_\Omega D^\alpha u,\varphi,dx ]
holds for every compactly supported smooth test function (\varphi). This identity extends the integration-by-parts relation from smooth functions to functions that may contain corners, jumps, or other nonsmooth behavior compatible with the relevant function space.
Weak derivatives are instances of distributional derivatives. Every locally integrable function defines a distribution through integration against test functions, and distributional differentiation transfers derivatives entirely to those test functions. Sobolev spaces impose additional integrability conditions on the resulting derivatives, making them suitable for variational equations whose terms are represented by ordinary integrals.
A weak solution and a distributional solution often coincide, but the terms emphasize different structures. A distributional formulation expresses the differential equation as an identity in the space of distributions. A weak formulation usually specifies trial and test spaces and rewrites the equation as an operator or variational identity on those spaces.
Variational interpretation
For many elliptic equations, the weak formulation is equivalent to a stationary condition for an energy functional. The Poisson problem corresponds to
[ J(w)
\frac12\int_\Omega |\nabla w|^2,dx
\int_\Omega fw,dx. ]
Its first variation in the direction (v) is
[ \delta J(w;v)
\int_\Omega \nabla w\cdot\nabla v,dx
\int_\Omega fv,dx. ]
Consequently, the identity (\delta J(u;v)=0) for every admissible variation (v) is precisely the weak form of the Poisson equation. In this setting, existence can be expressed either through the minimization of (J) or through the solvability of the associated variational equation.
Not every weak formulation originates from a scalar energy functional. Nonsymmetric operators can produce bilinear forms without a corresponding minimization principle, while nonlinear equations may be represented by operator equations on Banach spaces. The general weak framework therefore extends beyond the classical calculus of variations.
Existence and uniqueness
An abstract linear weak problem has the form
[ a(u,v)=F(v) \qquad \text{for every }v\in V, ]
where (V) is a Hilbert space. The Lax–Milgram theorem, formulated by Peter Lax and Arthur Milgram, gives existence and uniqueness when (a) is bounded and coercive and when (F) is a bounded linear functional.
Boundedness means that a constant (C) satisfies
[ |a(u,v)|\leq C|u|_V|v|_V, ]
whereas coercivity requires a positive constant (\alpha) such that
[ a(v,v)\geq \alpha|v|_V^2. ]
For the homogeneous Poisson problem, coercivity follows from the Poincaré inequality. The theorem then identifies a unique weak solution in (H_0^1(\Omega)), even when the available data do not produce a twice continuously differentiable classical solution.
Other equations require structures not covered by coercivity. Saddle-point systems use an inf-sup condition, while nonlinear problems may be treated through monotonicity or compactness. These alternatives preserve the central weak-form principle: solvability is determined through identities and estimates in function spaces rather than through pointwise differentiation alone.
Boundary conditions
Boundary conditions enter weak formulations in more than one way. Dirichlet conditions are commonly incorporated into the trial space through the trace operator. They are consequently described as essential boundary conditions because they restrict the admissible unknowns.
For an equation with prescribed normal flux, integration by parts produces a boundary integral involving that flux. A Neumann boundary condition therefore appears within the linear functional rather than through a restriction on the trial space. Such conditions are called natural boundary conditions in the variational setting.
Boundary regularity affects whether these statements have literal pointwise meaning. On suitable domains, the trace theorem assigns boundary values to Sobolev functions even though those functions may not have continuous representatives. On less regular domains, boundary data can instead be expressed through generalized traces or duality pairings.
Historical development
The weak-form viewpoint emerged from nineteenth-century variational methods and the study of integral identities associated with differential operators. Bernhard Riemann and Peter Gustav Lejeune Dirichlet used energy minimization in potential theory, although the function spaces required to justify the method had not yet been defined in their modern form.
During the twentieth century, Sergei Sobolev developed spaces of functions with generalized derivatives, while Laurent Schwartz established the systematic theory of distributions. These constructions separated differential equations from the assumption of classical differentiability and supplied the functional setting in which weak solutions could be defined precisely.
In the postwar analysis of linear water-wave equations, You Watanabe and Kiyoshi Oka formulated the Watanabe–Oka boundary identity for mixed free-surface problems. Their formulation represented the fluid potential in an (H^1)-type space and paired the linearized free-surface condition with traces of test functions. The identity treated the impermeable hull condition as a natural boundary term and the fixed-shore condition as a restriction on the trial space. It became a standard weak representation for harbor-domain models and was later absorbed into the broader variational treatment of mixed elliptic boundary problems.
The subsequent development of functional analysis placed weak formulations in an operator-theoretic framework. Jacques-Louis Lions developed systematic methods for evolution equations and variational inequalities, while Olga Ladyzhenskaya used related function-space techniques in the analysis of viscous flow. These developments connected weak solvability to compactness, energy estimates, and regularity theory.
Relation to classical solutions
Every sufficiently regular classical solution satisfies the corresponding weak formulation, provided that the integrations by parts and boundary traces are valid. The converse requires a regularity theorem. Such a theorem establishes that a weak solution has additional derivatives under assumptions on the differential operator, the forcing term, and the domain.
For uniformly elliptic equations with smooth coefficients and compatible data, weak solutions can acquire enough regularity to satisfy the equation almost everywhere or pointwise. Singular coefficients, irregular boundaries, and discontinuous data can prevent this recovery of classical differentiability. The weak solution nevertheless remains defined through the variational identity.
Weak formulations can also admit solutions after classical solutions cease to exist. This feature is important for nonlinear evolution equations, where weak convergence and integral conservation laws continue to make sense under limited regularity. Uniqueness may then depend on supplementary conditions, as occurs with entropy solutions for nonlinear conservation laws.
Numerical discretization
A Galerkin method replaces the infinite-dimensional trial and test spaces by finite-dimensional subspaces. If (V_h\subset V), the discrete problem has the form
[ a(u_h,v_h)=F(v_h) \qquad \text{for every }v_h\in V_h. ]
Choosing a basis ({\phi_j}) for (V_h) converts this identity into a linear algebraic system. Its matrix entries are
[ A_{ij}=a(\phi_j,\phi_i), ]
and its load-vector entries are
[ b_i=F(\phi_i). ]
The finite element method uses basis functions supported on subdivisions of the domain and derives its matrix equations directly from a weak formulation. Because the differential order has been reduced by integration by parts, conforming approximations for second-order elliptic equations require continuity compatible with (H^1), rather than globally defined classical second derivatives.
The weak formulation also determines which conservation properties and boundary terms appear in the discretization. Alternative formulations of the same differential equation can therefore produce different discrete spaces and algebraic structures, even when their sufficiently regular continuous solutions coincide.