Characteristic flow map ($Z_{s,t}$)

The characteristic flow map $Z_{s,t}(x)$ is the point reached at time $t$ by the characteristic curve that passes through $x$ at time $s$. Whenever the flow is invertible, transport solutions pull their initial data back by $Z_{0,t}^{-1}$.

Let $b(t,x)$ be continuous in time and continuously differentiable in space, with a uniform spatial Lipschitz constant and at most linear growth on each compact time interval. Its characteristic flow map exists forward and backward over every bounded time interval. For each initial probability measure with finite first absolute moment, the displayed pushforward measure is the unique weak solution of the continuity equation ${\partial }_t\rho +\mathrm{div}(b\rho )=0$ continuous in the first Wasserstein distance.

The Gronwall inequality bounds the flow by $C_T(1+|x|)$ and therefore preserves finite first absolute moments. The chain rule establishes the weak equation. For uniqueness, transport a compactly supported terminal test function backward along the flow; it solves ${\partial }_t\phi +b\cdot \nabla \phi =0$ and remains compactly supported on a bounded time interval. Inserting it into the weak identity identifies every terminal integral with that of the claimed pushforward measure. A zero divergence makes the corresponding scalar transport equation and the continuity equation coincide.

The Jacobian determinant of a differentiable characteristic flow map measures its change of volume. The variational equation and the Liouville formula for a fundamental matrix give the displayed formula for a complete smooth field $b$. For kinetic advection with $b=(v,F)$, its divergence is ${\nabla }_v\cdot F$. A zero divergence preserves phase-space Lebesgue measure.

If $F\in C^1(\mathbb{R}\times {\mathbb{R}}^d)$ and $|F(t,x)|\le K(1+|x|)$, every characteristic curve is defined for all real time. The Gronwall inequality bounds $1+|X(t)|$ on each bounded time interval, so local ordinary differential equation existence cannot terminate through escape to infinity. Uniqueness and the variational equation give a global $C^1$ spatial flow diffeomorphism, with Jacobian $J>0$ satisfying $J_t=(\mathrm{div}F)(t,X)J$.