Adaptivity with moving grids
Adaptivity with moving grids
Adaptivity with moving grids
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
94 C. J. Budd, W. Huang and R. D. Russell<br />
Whilst the PDE may be invariant under the action of a symmetry group,<br />
not all of the solutions have this property, although any one solution can<br />
be mapped into another by the action of the symmetry operator. Those<br />
solutions which are themselves invariant under the action of the symmetry<br />
operator are termed the self-similar solutions. A self-similar solution<br />
satisfies the functional equation<br />
u(λt, λ α x)=λ β u(t, x). (5.4)<br />
Such a solution can be described in terms of a new set of coordinates,<br />
typically of the form<br />
u(t, x) =t β v(z), z = x/t β . (5.5)<br />
In the case of the porous medium equation, the self-similar solution <strong>with</strong><br />
constant mass takes the form<br />
u(t, x) =t −1/3 v(x/t 1/3 ).<br />
The function v(z) generally satisfies an ordinary differential equation, which<br />
is much simpler than the original PDE. The solutions of this ODE which<br />
correspond to solutions of the PDE are generally those <strong>with</strong> certain decay or<br />
boundedness conditions as |z| →∞. This allows the construction of exact<br />
solutions to PDEs in many cases. In the example of the porous medium<br />
equation we have the famous Barenblatt–Pattle solution (Barenblatt 1996),<br />
given by<br />
u(x, t) =t −1/3 (a − x 2 /t 2/3 ) + ,<br />
where a is a constant. Significantly such solutions are globally attracting (in<br />
the sense of L 1 -convergence. An extensive description of such problems is<br />
given in Barenblatt (1996) and Olver (1986).<br />
More generally, solutions of PDEs <strong>with</strong> symmetries may also take the<br />
form<br />
u(t, x) =U(t)v(z), where z = x/L(t). (5.6)<br />
Here U(t) andL(t) are appropriate solution and length scales. In the case<br />
of self-similar solutions these are pure powers of t (or a translation of t) but<br />
they can take more general forms. For example, in the case of the blow-up<br />
equation,<br />
u t = u xx + u 3 ,<br />
<strong>with</strong> u →∞as t → T ,wehave<br />
U(t) =(T − t) −1/2 and L(t) =(T − t) 1/2 | log(T − t)|.<br />
Such solutions are called approximately self-similar solutions (Samarskii,<br />
Galaktionov, Kurdyumov and Mikhailov 1973). Examples of this type of