For a matrix field g=(gij)g=(g_{ij}) with entries in L1(Rn)L^1(\mathbb R^n), define its canonical pressure π∗\pi_* by

π∗^(ξ)=−∑i,jξiξj∣ξ∣2gij^(ξ)(ξ≠0).\widehat{\pi_*}(\xi)=-\sum_{i,j}\frac{\xi_i\xi_j}{|\xi|^2}\widehat{g_{ij}}(\xi)\quad(\xi\ne0).

The right side is a bounded measurable function. Its inverse Fourier transform is in H−sH^{-s} for every s>n/2s>n/2, by the .

Pressure equation

The definition gives

−Δπ∗=∑i,j∂i∂jgij,∥π∗∥H−s≤Cn,s∑i,j∥gij∥1.-\Delta\pi_*=\sum_{i,j}\partial_i\partial_jg_{ij},\qquad \|\pi_*\|_{H^{-s}}\le C_{n,s}\sum_{i,j}\|g_{ij}\|_1.

For inputs also in L2L^2, this agrees with ∑RiRjgij\sum R_iR_jg_{ij}. For an incompressible velocity, taking divergence of the momentum equation produces this pressure source with g=u⊗ug=u\otimes u, when the other terms have zero divergence. A difference of quadratic tensors gives a difference-pressure source.

Actual pressures

Another distributional solution differs by a harmonic distribution. Identifying its gradient with ∇π∗\nabla\pi_* requires a global growth or Sobolev bound on that gradient; the Poisson equation alone does not supply it.