A Finite Volume Method for the Relativistic Burgers Equation on a FLRW Background Spacetime


Creative Commons License

Ceylan T., Lefloch P. G., Okutmuştur B.

COMMUNICATIONS IN COMPUTATIONAL PHYSICS, cilt.23, ss.500-519, 2018 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 23
  • Basım Tarihi: 2018
  • Doi Numarası: 10.4208/cicp.020415.260717a
  • Dergi Adı: COMMUNICATIONS IN COMPUTATIONAL PHYSICS
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Sayfa Sayıları: ss.500-519
  • Orta Doğu Teknik Üniversitesi Adresli: Evet

Özet

A relativistic generalization of the inviscid Burgers equation was introduced by LeFloch and co-authors and was recently investigated numerically on a Schwarzschild background. We extend this analysis to a Friedmann-Lemaitre-Robertson-Walker (FLRW) background, which is more challenging due to the existence of time-dependent, spatially homogeneous solutions. We present a derivation of the model of interest and we study its basic properties, including the class of spatially homogeneous solutions. Then, we design a second-order accurate scheme based on the finite volume methodology, which provides us with a tool for investigating the properties of solutions. Computational experiments demonstrate the efficiency of the proposed scheme for numerically capturing weak solutions.

A relativistic generalization of the inviscid Burgers equation was introduced by LeFloch and co-authors and was recently investigated numerically on a Schwarzschild background. We extend this analysis to a Friedmann-Lemaitre-Robertson-Walker (FLRW) background, which is more challenging due to the existence of time-dependent, spatially homogeneous solutions. We present a derivation of the model of interest and we study its basic properties, including the class of spatially homogeneous solutions. Then, we design a second-order accurate scheme based on the finite volume methodology, which provides us with a tool for investigating the properties of solutions. Computational experiments demonstrate the efficiency of the proposed scheme for numerically capturing weak solutions