28.01.2013 Views

Implicit-Explicit Runge-Kutta schemes for hyperbolic systems ... - utenti

Implicit-Explicit Runge-Kutta schemes for hyperbolic systems ... - utenti

Implicit-Explicit Runge-Kutta schemes for hyperbolic systems ... - utenti

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.

Stability analysis<br />

When studying the A-stability of a RK scheme, one considers a scalar equation of the<br />

<strong>for</strong>m<br />

y ′ = λy, y(0) = 1<br />

with y : R + → R, and λ ∈ R, ℜλ ≤ 0. Such test problem is sufficient to characterize the<br />

stability property of a <strong>Runge</strong>-<strong>Kutta</strong> scheme when applied to a linear system of the <strong>for</strong>m<br />

with y ∈ R m , and B ∈ R m×m .<br />

y ′ = By, y(0) = y0<br />

In the case of IMEX-RK, and in general <strong>for</strong> additive <strong>Runge</strong>-<strong>Kutta</strong> methods, such a stability<br />

analysis has a limited validity. In fact, consider a generic linear system of the <strong>for</strong>m<br />

y ′ = B1 y + B2 y<br />

with y ∈ R m , and B1, B2 ∈ R m×m , and apply an IMEX scheme which is explicit in B1y and<br />

implicit in B2y. The stability of the numerical solution depends on the two matrices B1<br />

and B2, and not only on their eigenvalues since in general the two matrices do not share<br />

the same eigenvectors, and there<strong>for</strong>e they can not be diagonalized simultaneously.<br />

18

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

Saved successfully!

Ooh no, something went wrong!