Research reports

Essentially optimal explicit Runge-Kutta methods with application to hyperbolic-parabolic equations

by M. Torrilhon and R. Jeltsch

(Report number 2008-21)

Abstract
Optimal explicit Runge-Kutta methods consider more stages in order to include a particular spectrum in their stability domain and thus reduce time step restrictions. This idea, so far used mostly for real line spectra, is generalized to more general spectra in the form of a thin region. In thin regions the eigenvalues may extend away from the real axis into the imaginary plane. We give a direct characterization of optimal stability polynomials containing a maximal thin region and calculate these polynomials for various cases. Semi-discretizations of hyperbolic-parabolic equations are a relevant application which exhibit a thin region spectrum. As a model, linear, scalar advection-diffusion is investigated. The second order stabilized explicit Runge-Kutta methods derived from the stability polynomials are applied to advection-diffusion and compressible, viscous fluid dynamics in numerical experiments. Due to the stabilization the time step can be controlled solely from the hyperbolic CFL condition even in the presence of viscous fluxes.

Keywords:

BibTeX
@Techreport{TJ08_385,
  author = {M. Torrilhon and R. Jeltsch},
  title = {Essentially optimal explicit Runge-Kutta methods with application to hyperbolic-parabolic equations},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2008-21},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2008/2008-21.pdf },
  year = {2008}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).

JavaScript has been disabled in your browser