07.08.2013 Views

skade eaction ffusion ne-dimensional System - ZIB

skade eaction ffusion ne-dimensional System - ZIB

skade eaction ffusion ne-dimensional System - ZIB

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.

The nom whih KARDOS uses s a weighted L-no, giv by<br />

whe<br />

it<br />

uAbsi<br />

uRelMxi<br />

\\ u<br />

r-<br />

||«<br />

||^||o < f ^ i<br />

uAfe, < ||u||o < uRelMxi<br />

uRelM ||u<br />

ReM RTOL o ,<br />

uA6s ATOL * (O) 1 /<br />

:i-7)<br />

Scaling can be tur<strong>ne</strong>d on or off. If the scaled error norm is used, the tolerances<br />

RTOLi and ATOLi have to be seleced very carefully to reflect accurately<br />

the scale of the problem. The tolra<strong>ne</strong> ATOLi should indicate the absolute<br />

value at which the i-th compo<strong>ne</strong>nt is essentially insignifican. O the other<br />

hand, the value RTOLi affects the relative accuracy of the i-th compo<strong>ne</strong>nt<br />

with respect to its maximal value in time. This control turns out to be quite<br />

efficient and robust for a wid lass of problems. Howeer, it is cear that i<br />

is not a univeral mthod.<br />

To implement o<strong>ne</strong> time step, o<strong>ne</strong> elliptic proble has to be solved at each<br />

tage of the Rosenbrock method. This is do<strong>ne</strong> by means of an adaptiv<br />

ultilevel finite element method elaborated in recent years, see DEUFLHARD<br />

et al. []. This method is an excellent tool to adapt the space discretization<br />

for the current solution in uch a way, that a prescribed tolerance spaceTol<br />

is achieved. KARDOS uses standard li<strong>ne</strong>ar elements con<strong>ne</strong>cted with local<br />

error estimator of BABUSKA-HEINBOLDT [1] yp. We get at each tage,<br />

discretizations of the form<br />

B T v ) = n(v) e SnVvn e S i = 3 1.8)<br />

where Sn denotes the space of al ntinuous, piecewise li<strong>ne</strong>ar functions on<br />

partition An of 0, satisfying the Dirichlet boundary conditions. In S% hey<br />

vanish on the Dirichlet boundary. The equations only differ in the respective<br />

righthand id and in the pecial boundary conditions taken in onsider<br />

ation.

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

Saved successfully!

Ooh no, something went wrong!