with A. Bertozzi, Nonlinearity, 7(1994)1535-1564.
For simplicity we consider periodic boundary conditions which has the physical interpretation of modeling a periodic array of droplets. In a previous work we studied the above equation without the second order `porous media' term. In particular we showed the existence of nonnegative weak solutions with increasing support for 0<n<3 but the techniques failed for . This is consistent with the fact that, in this case, nonnegative self-similar source-type solutions do not exist for .
In this work, we discuss a physical justification for the `porous media' term when n=3 and 1<m<2. We propose such behavior as a cut off of the singular `disjoining pressure' modeling long range Van der Waals interactions.
For all n>0 and 1<m<2 we discuss possible behavior at the edge of the support of the solution via leading order asymptotic analysis of traveling wave solutions. This analysis predicts a certain `competition' between the second and fourth order terms. We present rigorous weak existence theory for equation (1) for all n>0 and 1<m<2. In particular, the presence of a second order `porous media' term in equation (1) yields nonnegative weak solutions that converge to their mean as and that have additional regularity. Moreover, we show that there exists a time T* after which the weak solution is a positive strong solution. For n > 3/2, we show that the regularity of the weak solutions is in exact agreement with that predicted by the asymptotics.
Finally, we present several numerical computations of solutions.
The simulations
use a weighted implicit-explicit scheme on a dynamically adaptive mesh.
The numerics suggest that the weak solution
described by our existence theory has
compact support with a finite speed of propagation.
The data confirms the local `power law'
behavior at the edge of the support predicted by
asymptotics.
click here for the *.ps file of
the article.
click here for the *.pdf file of
the article.
Some of the computations presented in this paper were performed at the
Advanced Computing Research Facility, Mathematics and Computer Science
Division, Argonne National Laboratory.
Both authors are supported by NSF postdoctoral fellowships. AB is
also partially supported by the Materials Research Laboratory at the
Univ. of Chicago and the DOE and MP by NSF grant DMS-9305996.