An equivalence of moment closure and nonlinear variational approximation of the Fokker-Planck equation for dilute polymeric flow

We establish rigorously the equivalence between classical moment closure and a nonlinear variational approximation of the Fokker-Planck equation for dilute polymeric flow in the linearized Hookean spring chain setting. The variational formulation is based on the Dirac-Frankel principle applied to a Gaussian approximation manifold endowed with the Fisher-Rao information metric. We show that the invariance of this manifold under the linear configurational dynamics yields an exact evolution for the macroscopic conformation tensor, recovering the classical diffusive Oldroyd-B closure. While the equivalence only holds in the linearized setting, the associated variational framework provides an abstract error representation and a starting point for the systematic construction of reduced approximation schemes for polymeric flows with nonlinear forcing laws.

💡 Research Summary

The paper rigorously establishes that, for dilute polymeric flows modeled by a Hookean spring chain, the classical second‑moment closure of the Fokker‑Planck (FP) equation and a nonlinear variational approximation on a Gaussian manifold are mathematically equivalent. The authors begin by recalling the coupled Navier‑Stokes–Fokker‑Planck system: the macroscopic velocity and pressure satisfy incompressible Navier‑Stokes equations, while the microscopic configuration of a polymer chain is described by a probability density ψ(t,x,q) that evolves according to a linear FP equation. The extra‑stress tensor τ entering the momentum balance is the second moment of ψ, so any reduction of ψ to a finite set of moments yields a closed macroscopic model.

In the linearized Hookean setting the configurational operator can be diagonalized using the eigenvectors of the Rouse matrix R. This leads to a block‑diagonal structure where each mode evolves independently. By testing the FP equation with the tensor qqᵀ (or equivalently by a Hermite spectral truncation at second order) one obtains an exact evolution equation for the conformation tensor Cₙ: \