EVOLUTION OF SELF-SIMILARITY, AND OTHER PROPERTIES OF WAITING-TIME SOLUTIONS OF THE POROUS-MEDIUM EQUATION - THE CASE OF VISCOUS GRAVITY CURRENTS

The one-dimensional porous medium equation h(l) = (h(m)h(x))(x) (m > 0 ) admits waiting-time solutions, whose front remains motionless during a finite time interval t(w) before starting to move. We consider a fa mily of initial value problems, and investigate the asymptotics, close to the front and near start-up, which we expect to be self-similar. W e obtain numerical solutions for viscous gravity currents (m = 3) and power-law initial conditions (h proportional to x(p), h is proportiona l to the thickness of the fluid, x is the distance to the front). We f ind that: (a) if p < 2/3 the front starts moving immediately, (b) if p = 2/3 the front remains motionless during a finite time, (c) if p > 2 /3 one obtains waiting-time solutions in which a moving corner layer ( a small interval Delta x in which h(x) varies strongly) appears behind the front; the front starts moving when it is overrun by the corner l ayer. The corner layer strengthens (Delta x reduces and the variation of h(x) increases) as it approaches the front. Our initial conditions produce waiting-time solutions whose front starts moving with nonzero velocity. We determine t(w)(p) and study the motion of the corner laye r and the front, as well as other properties of the solutions. We comp are the results with the theoretical upper and lower bounds of t(w). W e investigate the asymptotics of the numerical solutions for p > 2/3, close to the corner layer and the front, and near start-up. To represe nt this asymptotics various kinds of similarity solutions are availabl e, that can be classified according to the self-similarity exponent de lta. We find that only two types (called L and A) are relevant. The L solutions correspond to 1 < delta < 13/10, and have an infinite series of corner layers that accumulate at the front. The part of these solu tions behind the first corner layer of the series represents the asymp totics of the numerical solutions in a domain that excludes the region between the corner layer and the front, for a time interval excluding the neighbourhood of start-up. The A solutions have delta less than o r equal to 1, and represent the evolution of the strong corner layer t hat is arriving at the front. The numerical evidence shows that the co nstant front velocity solution (type A with delta = 1) describes the a symptotics close to, and including start-up, so that the motion of the corner layer joins smoothly with that of the front.