Real analytic solutions to the divergence equation

Abstract: In this paper, we develop a differential-topological method to yield explicit real analytic solutions $v$ to the divergence equation $div_{\mathbb{R}^n} v = f$ on any annali $A(R_1 ,R_2) = { x \in \mathbb{R}ⁿ : R_1 < |x| < R_2}$, with $n \geq 2$, and $0 < R_1 < R_2 < \infty$. The prescribed source term $f$ is supposed to be real analytic on $$\overline{A(R_1 , R_2)} = { x \in \mathbb{R}ⁿ : R_1 \leq |x| \leq R_2}$$ satisfying the zero integral condition on $A(R_1, R_2)$. The resulting solution $v$ is a real analytic vector field on $\overline{A(R_1 , R_2)}$, which vanishes on $\partial \big( A(R_1, R_2 ) \big )$. The method which we develop here is different from the standard Bogovski approach and the Kapitanskii-Pileckas approach. The first main step our method is a clever differential-topological argument, which we develop under the inspiration and guidance of the standard proof of the cohomological statement $$H_cⁿ \big ( \mathbb{R}^n\big ) = \mathbb{R}$$ in Spviak book A Comprehensive Introduction to Differential Geometry, Vol I. This allows us to reduce the problem to that of solving a linear algebra problem.

• The paper constructs explicit real analytic vector fields U on spherical annuli satisfying div U=f and U=0 on both boundary spheres whenever analytic f has zero integral.
• The method combines compactly supported de Rham cohomology, a Poisson equation on the sphere, analytic cutoffs, and finite-dimensional Hermite interpolation to enforce boundary conditions.
• The result extends divergence-equation solvability into the real analytic category, while leaving quantitative estimates, general domains, and higher-order boundary interpolation as open problems.

This paper constructs explicit real analytic solutions to the divergence equation $\dv_{\mathbb{R}^n} U = f$ on annuli $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$, with $n \geq 2$ and $0 < R_1 < R_2 < \infty$, subject to the homogeneous Dirichlet condition $U|_{\partial A(R_1,R_2)} = 0$. The source term $f$ is assumed real analytic on the closed annulus with vanishing total integral. The construction proceeds via a differential-topological method inspired by Spivak's proof of the cohomological identity $H_c^n(\mathbb{R}^n) = \mathbb{R}$, rather than by the classical Bogovskiĭ or Kapitanskii–Pileckas approaches (2602.21925).

Main result

The principal theorem states that for any $f$ real analytic on $\overline{A(R_1,R_2)}$ with $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$, there exists a vector field $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$0 that is real analytic on the closed annulus, satisfies $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$1 pointwise there, and vanishes on the boundary spheres. Real analyticity on the closed annulus is defined via extension to a slightly larger annulus $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$2.

The theorem is stated in two forms: a "loose" existence statement, and a precise version giving an explicit formula for $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$3. The formula involves three ingredients:

• a radial primitive term $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$4, where $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$5 is a corrected version of $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$6;
• a correction term built from the Hodge star and the exterior derivative of a pullback of the codifferential of $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$7, where $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$8 solves a surface Poisson equation $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$9 on the sphere;
• an explicit radial vector field $n \geq 2$0, where $n \geq 2$1 is given by a closed-form expression involving rational functions $n \geq 2$2 of the aspect ratio $n \geq 2$3.

The rational functions $n \geq 2$4 are written out in full for $n \geq 2$5, and their role is to enforce both the vanishing integral condition and the boundary condition $n \geq 2$6. A consequence worth noting: the construction is fully explicit, so it yields a specific analytic solution rather than a mere existence assertion — the paper stresses that solutions to the divergence equation are highly non-unique, so a canonical explicit choice has analytical value.

Method and relation to prior approaches

The standard tools for this problem operate in weaker topologies. The Bogovskiĭ operator maps zero-mean $n \geq 2$7 data ($n \geq 2$8) on a bounded Lipschitz domain to a $n \geq 2$9 solution with the scale-invariant estimate $0 < R_1 < R_2 < \infty$0; the Kapitanskii–Pileckas approach handles Hölder data $0 < R_1 < R_2 < \infty$1 with $0 < R_1 < R_2 < \infty$2 boundary, producing $0 < R_1 < R_2 < \infty$3. Neither framework reaches the real analytic category, and the authors state that the highest regularity they could find in the existing literature is $0 < R_1 < R_2 < \infty$4; analytic data with analytic solutions under a Dirichlet condition appear to be essentially untreated.

The method here instead works with differential forms and de Rham cohomology with compact support. The key lemma constructs a compactly supported $0 < R_1 < R_2 < \infty$5-form $0 < R_1 < R_2 < \infty$6 with $0 < R_1 < R_2 < \infty$7, supported in $0 < R_1 < R_2 < \infty$8, and real analytic on the closed annulus. Applying the Hodge star and the musical isomorphism converts $0 < R_1 < R_2 < \infty$9 into the desired vector field. The construction of $U|_{\partial A(R_1,R_2)} = 0$0 follows Spivak's proof of $U|_{\partial A(R_1,R_2)} = 0$1 (Theorem 9, Chapter 8 of A Comprehensive Introduction to Differential Geometry, Vol. I), with two substantive modifications:

• Canonical primitive on the sphere: the arbitrary smooth $U|_{\partial A(R_1,R_2)} = 0$2-form $U|_{\partial A(R_1,R_2)} = 0$3 with $U|_{\partial A(R_1,R_2)} = 0$4 in Spivak's argument is replaced by the explicit form $U|_{\partial A(R_1,R_2)} = 0$5, where $U|_{\partial A(R_1,R_2)} = 0$6 is the unique zero-mean solution of the Poisson equation $U|_{\partial A(R_1,R_2)} = 0$7. Elliptic regularity transfers the analyticity of $U|_{\partial A(R_1,R_2)} = 0$8 (which follows from that of $U|_{\partial A(R_1,R_2)} = 0$9) to $f$0 and hence to the whole correction form.
• Analytic cutoff: Spivak's smooth cutoff is replaced by the function $f$1 on the annulus, which is real analytic there and flat-compatible with the constant values $f$2 and $f$3 outside.

The case $f$4 requires a separate argument, since the codifferential route to an exact primitive on $f$5 is replaced by a direct integration: $f$6 is exact on $f$7 by the de Rham isomorphism $f$8, with primitive $f$9, which is real analytic whenever $H_c^n(\mathbb{R}^n) = \mathbb{R}$0 is.

Auxiliary lemmas

Three technical results support the main construction.

Division of analytic functions. A standard lemma shows that if $H_c^n(\mathbb{R}^n) = \mathbb{R}$1 is real analytic on $H_c^n(\mathbb{R}^n) = \mathbb{R}$2 and vanishes identically on the slice $H_c^n(\mathbb{R}^n) = \mathbb{R}$3, then $H_c^n(\mathbb{R}^n) = \mathbb{R}$4 with $H_c^n(\mathbb{R}^n) = \mathbb{R}$5 real analytic and unique. A derived lemma uses this repeatedly to show that a function $H_c^n(\mathbb{R}^n) = \mathbb{R}$6 that is $H_c^n(\mathbb{R}^n) = \mathbb{R}$7, analytic in a slab $H_c^n(\mathbb{R}^n) = \mathbb{R}$8, and vanishes identically on one side of the slab is actually $H_c^n(\mathbb{R}^n) = \mathbb{R}$9 across the transition — this is what upgrades regularity of the radial primitive $f$0 at the inner sphere where it ceases to vanish.

Positivity of a polynomial. For the case $f$1 of the interpolation lemma, uniqueness of a linear system reduces to strict positivity of the degree-$f$2 polynomial $f$3 for $f$4. The paper proves $f$5 with all coefficients $f$6, via a decomposition into three polynomials $f$7 whose coefficients are shown positive by direct combinatorial estimates. For $f$8 the analogous polynomial $f$9 is manifestly positive for $\overline{A(R_1,R_2)}$0.

Hermite-type interpolation with vanishing integral. Given analytic boundary data $\overline{A(R_1,R_2)}$1 on $\overline{A(R_1,R_2)}$2 ($\overline{A(R_1,R_2)}$3, $\overline{A(R_1,R_2)}$4 or $\overline{A(R_1,R_2)}$5), there exist unique analytic coefficient functions $\overline{A(R_1,R_2)}$6 on $\overline{A(R_1,R_2)}$7 such that the radial polynomial $\overline{A(R_1,R_2)}$8 matches the prescribed boundary jets at $\overline{A(R_1,R_2)}$9 and satisfies $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$0. The proof is by row-reduction of an explicit $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$1 Vandermonde-type linear system, with invertibility following from the positivity of $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$2 or $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$3. This is the step the authors describe as reducing the analytic problem "to that of solving a linear algebra problem," and it is what allows the correction term $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$4 to fix the boundary values without disturbing the divergence equation or the zero-mean condition.

Structure of the proof of the main theorem

The proof combines the two major steps. Given analytic $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$5 with vanishing integral, define the corrected function $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$6 inside the annulus and zero outside. The interpolation lemma guarantees $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$7 and that $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$8 is continuous, compactly supported in the closed annulus, and analytic there. The cohomological corollary then produces an analytic, compactly supported vector field $\int_{A(R_1,R_2)} f\,\Vol_{\mathbb{R}^n} = 0$9 with $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$00 and $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$01. The residual $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$02 is exactly the divergence of the explicit radial field $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$03, and the interpolation lemma also gives $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$04. Setting $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$05 completes the construction; a direct computation identifies this with the closed-form expression in the theorem statement.

Limitations and open questions

Several restrictions are intrinsic to what is proved. The domain is exclusively a spherical annulus; the differential-topological argument relies on the radial retraction $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$06 and on solving a Poisson equation on $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$07, and the paper offers no extension to general domains. The interpolation step is carried out only for $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$08 and $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$09 boundary jets, and the cutoff exponent $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$10 is tied to these cases; whether the construction extends to arbitrary boundary regularity is not addressed. No quantitative estimates (e.g., analytic norms of $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$11 in terms of those of $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$12, or dependence on the ratio $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$13) are derived, in contrast with the scale-invariant $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$14 estimate of Bogovskiĭ or the $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$15 estimate of Kapitanskii–Pileckas; establishing such bounds would be a natural next question. The paper also does not address uniqueness or minimality of any analytic norm among solutions, and the authors themselves characterize the contribution as a first, "modest" step toward the analytic Dirichlet problem for the divergence equation. Finally, the connection to applications — for instance, analytic solutions of the incompressible Navier–Stokes or Stokes systems where Bogovskiĭ-type operators typically enter — is left implicit rather than developed.

Conclusion

The paper provides an explicit, self-contained construction of real analytic divergence-free-correcting vector fields on annuli with homogeneous Dirichlet boundary data, by transplanting Spivak's compactly supported cohomology argument into the analytic category and supplementing it with an analytic Hermite-type interpolation scheme that reduces the boundary correction to finite-dimensional linear algebra. The result fills a gap in the regularity scale of the divergence equation, which had previously been treated at most in the $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$16 category, and the explicit nature of the solution distinguishes it from the non-constructive existence available through Sobolev- and Hölder-space methods. The main open directions left by the paper are quantitative estimates for the constructed solution and extensions beyond the annular geometry and the $$A(R_1, R_2) = \{x \in \mathbb{R}^n : R_1 < |x| < R_2\}$$17 interpolation range (2602.21925).