RAM

A relativistic adaptive mesh refinement hydrodynamics code

Weiqun Zhang, Andrew I. MacFadyen

    Research output: Contribution to journalArticle

    Abstract

    We have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite-difference, weighted essentially nonoscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration we use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. We have also implemented fourth- and fifth-order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM is modular and includes the capability to easily swap hydrodynamics solvers, reconstruction methods, and physics modules. In addition to WENO, we have implemented a finite-volume module with the piecewise parabolic method for reconstruction and the modified Marquina approximate Riemann solver to work with TVD Runge-Kutta time integration. We examine the difficulty of accurately simulating shear flows in numerical relativistic hydrodynamics codes. We show that underresolved simulations of simple test problems with transverse velocity components produce incorrect results and demonstrate the ability of RAM to correctly solve these problems. RAM has been tested in one, two, and three dimensions and in Cartesian, cylindrical, and spherical coordinates. We have demonstrated fifth-order accuracy for WENO in one and two dimensions and performed detailed comparison with other schemes for which we show significantly lower convergence rates. Extensive testing is presented demonstrating the ability of RAM to address challenging open questions in relativistic astrophysics.

    Original languageEnglish (US)
    Pages (from-to)255-279
    Number of pages25
    JournalAstrophysical Journal, Supplement Series
    Volume164
    Issue number1
    DOIs
    StatePublished - May 2006

    Fingerprint

    hydrodynamics
    modules
    essentially non-oscillatory schemes
    parallel computers
    cylindrical coordinates
    Cartesian coordinates
    spherical coordinates
    hydrodynamic equations
    shear flow
    astrophysics
    computer programs
    decomposition
    physics
    code
    simulation
    method
    comparison

    Keywords

    • Hydrodynamics
    • Methods: numerical
    • Relativity

    ASJC Scopus subject areas

    • Space and Planetary Science

    Cite this

    RAM : A relativistic adaptive mesh refinement hydrodynamics code. / Zhang, Weiqun; MacFadyen, Andrew I.

    In: Astrophysical Journal, Supplement Series, Vol. 164, No. 1, 05.2006, p. 255-279.

    Research output: Contribution to journalArticle

    Zhang, Weiqun ; MacFadyen, Andrew I. / RAM : A relativistic adaptive mesh refinement hydrodynamics code. In: Astrophysical Journal, Supplement Series. 2006 ; Vol. 164, No. 1. pp. 255-279.
    @article{79ffbd5efa864ea38a5bde8446e57ef8,
    title = "RAM: A relativistic adaptive mesh refinement hydrodynamics code",
    abstract = "We have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite-difference, weighted essentially nonoscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration we use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. We have also implemented fourth- and fifth-order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM is modular and includes the capability to easily swap hydrodynamics solvers, reconstruction methods, and physics modules. In addition to WENO, we have implemented a finite-volume module with the piecewise parabolic method for reconstruction and the modified Marquina approximate Riemann solver to work with TVD Runge-Kutta time integration. We examine the difficulty of accurately simulating shear flows in numerical relativistic hydrodynamics codes. We show that underresolved simulations of simple test problems with transverse velocity components produce incorrect results and demonstrate the ability of RAM to correctly solve these problems. RAM has been tested in one, two, and three dimensions and in Cartesian, cylindrical, and spherical coordinates. We have demonstrated fifth-order accuracy for WENO in one and two dimensions and performed detailed comparison with other schemes for which we show significantly lower convergence rates. Extensive testing is presented demonstrating the ability of RAM to address challenging open questions in relativistic astrophysics.",
    keywords = "Hydrodynamics, Methods: numerical, Relativity",
    author = "Weiqun Zhang and MacFadyen, {Andrew I.}",
    year = "2006",
    month = "5",
    doi = "10.1086/500792",
    language = "English (US)",
    volume = "164",
    pages = "255--279",
    journal = "Astrophysical Journal, Supplement Series",
    issn = "0067-0049",
    publisher = "IOP Publishing Ltd.",
    number = "1",

    }

    TY - JOUR

    T1 - RAM

    T2 - A relativistic adaptive mesh refinement hydrodynamics code

    AU - Zhang, Weiqun

    AU - MacFadyen, Andrew I.

    PY - 2006/5

    Y1 - 2006/5

    N2 - We have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite-difference, weighted essentially nonoscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration we use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. We have also implemented fourth- and fifth-order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM is modular and includes the capability to easily swap hydrodynamics solvers, reconstruction methods, and physics modules. In addition to WENO, we have implemented a finite-volume module with the piecewise parabolic method for reconstruction and the modified Marquina approximate Riemann solver to work with TVD Runge-Kutta time integration. We examine the difficulty of accurately simulating shear flows in numerical relativistic hydrodynamics codes. We show that underresolved simulations of simple test problems with transverse velocity components produce incorrect results and demonstrate the ability of RAM to correctly solve these problems. RAM has been tested in one, two, and three dimensions and in Cartesian, cylindrical, and spherical coordinates. We have demonstrated fifth-order accuracy for WENO in one and two dimensions and performed detailed comparison with other schemes for which we show significantly lower convergence rates. Extensive testing is presented demonstrating the ability of RAM to address challenging open questions in relativistic astrophysics.

    AB - We have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite-difference, weighted essentially nonoscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration we use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. We have also implemented fourth- and fifth-order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM is modular and includes the capability to easily swap hydrodynamics solvers, reconstruction methods, and physics modules. In addition to WENO, we have implemented a finite-volume module with the piecewise parabolic method for reconstruction and the modified Marquina approximate Riemann solver to work with TVD Runge-Kutta time integration. We examine the difficulty of accurately simulating shear flows in numerical relativistic hydrodynamics codes. We show that underresolved simulations of simple test problems with transverse velocity components produce incorrect results and demonstrate the ability of RAM to correctly solve these problems. RAM has been tested in one, two, and three dimensions and in Cartesian, cylindrical, and spherical coordinates. We have demonstrated fifth-order accuracy for WENO in one and two dimensions and performed detailed comparison with other schemes for which we show significantly lower convergence rates. Extensive testing is presented demonstrating the ability of RAM to address challenging open questions in relativistic astrophysics.

    KW - Hydrodynamics

    KW - Methods: numerical

    KW - Relativity

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

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

    U2 - 10.1086/500792

    DO - 10.1086/500792

    M3 - Article

    VL - 164

    SP - 255

    EP - 279

    JO - Astrophysical Journal, Supplement Series

    JF - Astrophysical Journal, Supplement Series

    SN - 0067-0049

    IS - 1

    ER -