Use the conventional kinematic plume fluxes
Then is the volume flux, the dimensional mass flux, the dimensional momentum flux, and the kinematic buoyancy flux. Calling a mass flux uses a density-normalized convention: the reference density has been suppressed. If density-weighted momentum and buoyancy fluxes are desired, multiply and by the same constant throughout. Here the conventional kinematic definitions make the Boussinesq approximation explicit.
For an upward plume with , the inverse relations are
The volume conservation, vertical momentum conservation, and buoyancy balances from part a become
Here the vertical momentum conservation equation uses in the inertia and retains in the buoyancy force. The ambient has no imposed vertical momentum at entrainment.
To classify the resulting partial differential equations, work in the nonsingular primitive variables with . Subtract times the volume equation from the momentum equation and times the volume equation from the buoyancy equation. This gives
The principal matrix is therefore
All three characteristic speeds are real and equal to , so the model is a hyperbolic system if that term means only reality of the characteristic roots. However, for , requires , leaving only two independent eigenvectors. Hence
This weak hyperbolicity of the unsteady top-hat plume prevents a complete set of independent characteristic fields. For example, freezing in the homogeneous principal part gives and . A Fourier perturbation of wave number therefore produces a term proportional to in . Real characteristic roots alone do not provide uniform control of high-frequency perturbations. The flux-variable transformation above is invertible in the stated region, so changing to does not remove this degeneracy.