  We model the unsteady evolution of turbulent buoyant plumes following
temporal changes to the source conditions. The integral model is derived from
radial integration of the governing equations expressing the conservation of
mass, axial momentum and buoyancy. The non-uniform radial profiles of the axial
velocity and density deficit in the plume are explicitly described by shape
factors in the integral equations; the commonly-assumed top-hat profiles lead
to shape factors equal to unity. The resultant model is hyperbolic when the
momentum shape factor, determined from the radial profile of the mean axial
velocity, differs from unity. The solutions of the model when source conditions
are maintained at constant values retain the form of the well-established
steady plume solutions. We demonstrate that the inclusion of a momentum shape
factor that differs from unity leads to a well-posed integral model. Therefore,
our model does not exhibit the mathematical pathologies that appear in
previously proposed unsteady integral models of turbulent plumes. A stability
threshold for the value of the shape factor is identified, resulting in a range
of its values where the amplitude of small perturbations to the steady
solutions decay with distance from the source. The hyperbolic character of the
system allows the formation of discontinuities in the fields describing the
plume properties during the unsteady evolution. We compute numerical solutions
to illustrate the transient development following an abrupt change in the
source conditions. The adjustment to the new source conditions occurs through
the propagation of a pulse of fluid through the plume. The dynamics of this
pulse are described by a similarity solution and, by constructing this new
similarity solution, we identify three regimes in which the evolution of the
transient pulse following adjustment of the source qualitatively differ.
