skade eaction ffusion ne-dimensional System - ZIB
skade eaction ffusion ne-dimensional System - ZIB
skade eaction ffusion ne-dimensional System - ZIB
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.