Wave properties of fourth-order fully implicit Runge-Kutta time integration schemes

Swagata Bhaumik, Soumyo Sengupta, Aditi Sengupta

Research output: Contribution to journalArticle

4 Citations (Scopus)

Abstract

A global spectral analysis of fully implicit, two-stage, fourth-order Runge-Kutta method (IRK2) is presented here. The essential numerical properties of the IRK2 time-integration scheme, viz., nodal amplification factor, phase speed and group velocity have been quantified and plotted at all the nodal points of the computational domain as functions of the CFL (Courant-Friedrichs-Lewy number) and wavenumber for various explicit and implicit spatial discretization schemes. Results show that IRK2 scheme produces neutral stability at all interior nodes for all CFL numbers. This has been shown that while using IRK2 scheme one has to use refined grid at high CFL numbers for dispersion relation preservation property.

Original languageEnglish
Pages (from-to)110-121
Number of pages12
JournalComputers and Fluids
Volume81
DOIs
Publication statusPublished - 20-07-2013
Externally publishedYes

Fingerprint

Runge Kutta methods
Spectrum analysis
Amplification

All Science Journal Classification (ASJC) codes

  • Computer Science(all)
  • Engineering(all)

Cite this

Bhaumik, Swagata ; Sengupta, Soumyo ; Sengupta, Aditi. / Wave properties of fourth-order fully implicit Runge-Kutta time integration schemes. In: Computers and Fluids. 2013 ; Vol. 81. pp. 110-121.
@article{ee0e487b068f4a2297dde800f8be2542,
title = "Wave properties of fourth-order fully implicit Runge-Kutta time integration schemes",
abstract = "A global spectral analysis of fully implicit, two-stage, fourth-order Runge-Kutta method (IRK2) is presented here. The essential numerical properties of the IRK2 time-integration scheme, viz., nodal amplification factor, phase speed and group velocity have been quantified and plotted at all the nodal points of the computational domain as functions of the CFL (Courant-Friedrichs-Lewy number) and wavenumber for various explicit and implicit spatial discretization schemes. Results show that IRK2 scheme produces neutral stability at all interior nodes for all CFL numbers. This has been shown that while using IRK2 scheme one has to use refined grid at high CFL numbers for dispersion relation preservation property.",
author = "Swagata Bhaumik and Soumyo Sengupta and Aditi Sengupta",
year = "2013",
month = "7",
day = "20",
doi = "10.1016/j.compfluid.2013.04.007",
language = "English",
volume = "81",
pages = "110--121",
journal = "Computers and Fluids",
issn = "0045-7930",
publisher = "Elsevier Limited",

}

Wave properties of fourth-order fully implicit Runge-Kutta time integration schemes. / Bhaumik, Swagata; Sengupta, Soumyo; Sengupta, Aditi.

In: Computers and Fluids, Vol. 81, 20.07.2013, p. 110-121.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Wave properties of fourth-order fully implicit Runge-Kutta time integration schemes

AU - Bhaumik, Swagata

AU - Sengupta, Soumyo

AU - Sengupta, Aditi

PY - 2013/7/20

Y1 - 2013/7/20

N2 - A global spectral analysis of fully implicit, two-stage, fourth-order Runge-Kutta method (IRK2) is presented here. The essential numerical properties of the IRK2 time-integration scheme, viz., nodal amplification factor, phase speed and group velocity have been quantified and plotted at all the nodal points of the computational domain as functions of the CFL (Courant-Friedrichs-Lewy number) and wavenumber for various explicit and implicit spatial discretization schemes. Results show that IRK2 scheme produces neutral stability at all interior nodes for all CFL numbers. This has been shown that while using IRK2 scheme one has to use refined grid at high CFL numbers for dispersion relation preservation property.

AB - A global spectral analysis of fully implicit, two-stage, fourth-order Runge-Kutta method (IRK2) is presented here. The essential numerical properties of the IRK2 time-integration scheme, viz., nodal amplification factor, phase speed and group velocity have been quantified and plotted at all the nodal points of the computational domain as functions of the CFL (Courant-Friedrichs-Lewy number) and wavenumber for various explicit and implicit spatial discretization schemes. Results show that IRK2 scheme produces neutral stability at all interior nodes for all CFL numbers. This has been shown that while using IRK2 scheme one has to use refined grid at high CFL numbers for dispersion relation preservation property.

UR - http://www.scopus.com/inward/record.url?scp=84878165094&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84878165094&partnerID=8YFLogxK

U2 - 10.1016/j.compfluid.2013.04.007

DO - 10.1016/j.compfluid.2013.04.007

M3 - Article

VL - 81

SP - 110

EP - 121

JO - Computers and Fluids

JF - Computers and Fluids

SN - 0045-7930

ER -