We consider a one-dimensional aggregation equation for a non-negative density (x, t) associated with a quartic potential W(x) = x 2 + x 4 ( > 0, ∈ R). We show that for the case of symmetric initial data [ (x, 0) ≡ (−x, 0)], the solution of the aggregation equation can be expressed in terms of an explicit function of x, L(t), and Φ(t), where the functions L(t) and Φ(t) are determined by an ordinary differential equation initial value problem, the numerical solution of which is relatively straightforward. The function L(t), which can be interpreted as an upper bound for the radius of the support of (x, t), is finite for all t > 0, even if the support of (x, 0) is unbounded, while Φ(t) is a potential related to the time integral of the second spatial moment (t) of (x, t). We develop various bounds on L(t) and (t) and a number of results concerning convexity and monotonicity that require no knowledge of (x, 0) other than its symmetry. We show that many, but not all, of these results persist if the assumption of symmetry is relaxed. Our general results are tested against numerical solutions for two examples of symmetric (x, 0). We also exhibit some counterintuitive behaviour of the model, byfinding conditions under which xx (0, t) → ∞ as t → ∞, even if we start with xx (0, 0) < 0 and take ≥ 0, which ensures ultimate collapse to a single Dirac component at the origin.
KEYWORDSaggregation models, genetics and population dynamics, initial value problems for non-linear first-order systems, integro-partial differential equations, method of characteristics, numerical solutions
MSC CLASSIFICATION35R09; 35F55; 35L45; 45K05; 65M25