26.12.2013 Views

Adaptivity with moving grids

Adaptivity with moving grids

Adaptivity with moving grids

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

50 C. J. Budd, W. Huang and R. D. Russell<br />

R ≪ 1andR(a) =R(b) = 0. To leading order, R satisfies the equation<br />

ɛ(Ṙ − γR ξξ) =(MR ξ ) ξ +(M x x ξ R) ξ =(MR ξ + M ξ R) ξ =(MR) ξξ .<br />

Therefore<br />

ɛṘ =(1− γ∂2 ξ )−1 (MR) ξξ ≡ G(MR) ξξ . (3.14)<br />

Here G is a positive compact operator and, as M>0, ER ≡ (MR ξξ )isa<br />

uniformly elliptic operator <strong>with</strong> a negative real spectrum. It follows that R<br />

must decay to zero. Hence the equidistributed solution is locally stable.<br />

(ii) To prove this result, consider a slowly varying monitor function M(x, t)<br />

and an exact solution ˆx of the equidistribution equation (M ˆx ξ ) ξ =0. If<br />

x =ˆx + ɛR is a solution of (3.9) then, extending the calculation in (3.14),<br />

we see that R satisfies the equation<br />

ɛ(ẋ − γẋ ξξ )+O(ɛ 2 )=ɛ(MR) ξξ + O(ɛ 2 ).<br />

Hence, we have<br />

(MR) ξξ = ˙ˆx − γ ˙ˆx ξξ + O(ɛ).<br />

But as (M ˆx ξ ) ξ =0,wehave<br />

(M ˙ˆx ξ ) ξ = −(M tˆx ξ ) ξ .<br />

As the operator Eφ ≡ (Mφ) ξξ is uniformly elliptic, it follows that provided<br />

Ṁ is of order one, then ˙ˆx and hence R and its derivatives are also of order<br />

one. Consequently, the solution x of (3.9) stays ɛ-close to the solution of<br />

the equidistribution equation.<br />

(iii) To show this we need to show that x ξ cannot vanish. This result is a<br />

consequence of the maximum principle. In the case when γ =0,wehave,<br />

on differentiating (3.9), that<br />

ẋ ξ = M ξξ x ξ +2M ξ x ξξ + Mx ξξξ .<br />

Suppose that x ξ is initially positive everywhere, and as x evolves it vanishes<br />

for a first time at (<strong>with</strong>out loss of generality) the point ξ = 0. Then locally<br />

close to this point we have x ξ = aξ 2 + O(ξ 3 ) for some a>0. Hence<br />

ẋ ξ = aξ 2 M ξξ +2aξM ξ + aM + O(ξ).<br />

Thus ẋ ξ > 0 at this point and time. Hence x ξ must remain positive. The<br />

more general result follows from the positivity of the compact operator G.<br />

3.1.1. Coupling a one-dimensional MMPDE method to a PDE<br />

In one dimension, MMPDE methods can be very effectively coupled to an<br />

underlying PDE system by using a variety of different methods, including<br />

finite difference, finite element, collocation and spectral methods. The

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

Saved successfully!

Ooh no, something went wrong!