arrow
Volume 15, Issue 3
A Geometry-Preserving Finite Volume Method for Compressible Fluids on Schwarzschild Spacetime

Philippe G. LeFloch & Hasan Makhlof

Commun. Comput. Phys., 15 (2014), pp. 827-852.

Published online: 2014-03

Export citation
  • Abstract

We consider the relativistic Euler equations governing spherically symmetric, perfect fluid flows on the outer domain of communication of Schwarzschild spacetime, and we introduce a version of the finite volume method which is formulated from the geometric formulation (and thus takes the geometry into account at the discretization level) and is well-balanced, in the sense that it preserves steady solutions to the Euler equations on the curved geometry under consideration. In order to formulate our method, we first derive a closed formula describing all steady and spherically symmetric solutions to the Euler equations posed on Schwarzschild spacetime. Second, we describe a geometry-preserving, finite volume method which is based on the family of steady solutions to the Euler system. Our scheme is second-order accurate and, as required, preserves the family of steady solutions at the discrete level. Numerical experiments are presented which demonstrate the efficiency and robustness of the proposed method even for solutions containing shock waves and nonlinear interacting wave patterns. As an application, we investigate the late-time asymptotics of perturbed steady solutions and demonstrate its convergence for late time toward another steady solution, taking the overall effect of the perturbation into account.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-15-827, author = {Philippe G. LeFloch and Hasan Makhlof}, title = {A Geometry-Preserving Finite Volume Method for Compressible Fluids on Schwarzschild Spacetime}, journal = {Communications in Computational Physics}, year = {2014}, volume = {15}, number = {3}, pages = {827--852}, abstract = {

We consider the relativistic Euler equations governing spherically symmetric, perfect fluid flows on the outer domain of communication of Schwarzschild spacetime, and we introduce a version of the finite volume method which is formulated from the geometric formulation (and thus takes the geometry into account at the discretization level) and is well-balanced, in the sense that it preserves steady solutions to the Euler equations on the curved geometry under consideration. In order to formulate our method, we first derive a closed formula describing all steady and spherically symmetric solutions to the Euler equations posed on Schwarzschild spacetime. Second, we describe a geometry-preserving, finite volume method which is based on the family of steady solutions to the Euler system. Our scheme is second-order accurate and, as required, preserves the family of steady solutions at the discrete level. Numerical experiments are presented which demonstrate the efficiency and robustness of the proposed method even for solutions containing shock waves and nonlinear interacting wave patterns. As an application, we investigate the late-time asymptotics of perturbed steady solutions and demonstrate its convergence for late time toward another steady solution, taking the overall effect of the perturbation into account.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.291212.160913a}, url = {http://global-sci.org/intro/article_detail/cicp/7117.html} }
TY - JOUR T1 - A Geometry-Preserving Finite Volume Method for Compressible Fluids on Schwarzschild Spacetime AU - Philippe G. LeFloch & Hasan Makhlof JO - Communications in Computational Physics VL - 3 SP - 827 EP - 852 PY - 2014 DA - 2014/03 SN - 15 DO - http://doi.org/10.4208/cicp.291212.160913a UR - https://global-sci.org/intro/article_detail/cicp/7117.html KW - AB -

We consider the relativistic Euler equations governing spherically symmetric, perfect fluid flows on the outer domain of communication of Schwarzschild spacetime, and we introduce a version of the finite volume method which is formulated from the geometric formulation (and thus takes the geometry into account at the discretization level) and is well-balanced, in the sense that it preserves steady solutions to the Euler equations on the curved geometry under consideration. In order to formulate our method, we first derive a closed formula describing all steady and spherically symmetric solutions to the Euler equations posed on Schwarzschild spacetime. Second, we describe a geometry-preserving, finite volume method which is based on the family of steady solutions to the Euler system. Our scheme is second-order accurate and, as required, preserves the family of steady solutions at the discrete level. Numerical experiments are presented which demonstrate the efficiency and robustness of the proposed method even for solutions containing shock waves and nonlinear interacting wave patterns. As an application, we investigate the late-time asymptotics of perturbed steady solutions and demonstrate its convergence for late time toward another steady solution, taking the overall effect of the perturbation into account.

Philippe G. LeFloch and Hasan Makhlof. (2014). A Geometry-Preserving Finite Volume Method for Compressible Fluids on Schwarzschild Spacetime. Communications in Computational Physics. 15 (3). 827-852. doi:10.4208/cicp.291212.160913a
Copy to clipboard
The citation has been copied to your clipboard