26.12.2013 Views

Adaptivity with moving grids

Adaptivity with moving grids

Adaptivity with moving grids

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!