Second order divergence constraint preserving schemes for two-fluid relativistic plasma flow equations
Fuente:
arXiv
Guardado en:
| Autores principales: | , , , , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2025
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866915213853327360 |
|---|---|
| author | Agnihotri, Jaya Bhoriya, Deepak Kumar, Harish Chandrashekar, Praveen Balsara, Dinshaw S. |
| author_facet | Agnihotri, Jaya Bhoriya, Deepak Kumar, Harish Chandrashekar, Praveen Balsara, Dinshaw S. |
| contents | Two-fluid relativistic plasma flow equations combine the equations of relativistic hydrodynamics with Maxwell's equations for electromagnetic fields, which involve divergence constraints for the magnetic and electric fields. When developing numerical schemes for the model, the divergence constraints are ignored, or Maxwell's equations are reformulated as Perfectly Hyperbolic Maxwell's (PHM) equations by introducing additional equations for correction potentials. In the latter case, the divergence constraints are preserved only as the limiting case.
In this article, we present second-order numerical schemes that preserve the divergence constraint for electric and magnetic fields at the discrete level. The schemes are based on using a multidimensional Riemann solver at the vertices of the cells to define the numerical fluxes on the edges. The second-order accuracy is obtained by reconstructing the electromagnetic fields at the corners using a MinMod limiter. The discretization of Maxwell's equations can be combined with any consistent and stable discretization of the fluid parts. In particular, we consider entropy-stable schemes for the fluid part. The resulting schemes are second-order accurate, entropy stable, and preserve the divergence constraints of the electromagnetic fields. We use explicit and IMEX-based time discretizations. We then test these schemes using several one- and two-dimensional test cases. We also compare the divergence constraint errors of the proposed schemes with schemes having no divergence constraints treatment and schemes based on the PHM-based divergence cleaning. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2503_20372 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Second order divergence constraint preserving schemes for two-fluid relativistic plasma flow equations Agnihotri, Jaya Bhoriya, Deepak Kumar, Harish Chandrashekar, Praveen Balsara, Dinshaw S. Numerical Analysis Two-fluid relativistic plasma flow equations combine the equations of relativistic hydrodynamics with Maxwell's equations for electromagnetic fields, which involve divergence constraints for the magnetic and electric fields. When developing numerical schemes for the model, the divergence constraints are ignored, or Maxwell's equations are reformulated as Perfectly Hyperbolic Maxwell's (PHM) equations by introducing additional equations for correction potentials. In the latter case, the divergence constraints are preserved only as the limiting case. In this article, we present second-order numerical schemes that preserve the divergence constraint for electric and magnetic fields at the discrete level. The schemes are based on using a multidimensional Riemann solver at the vertices of the cells to define the numerical fluxes on the edges. The second-order accuracy is obtained by reconstructing the electromagnetic fields at the corners using a MinMod limiter. The discretization of Maxwell's equations can be combined with any consistent and stable discretization of the fluid parts. In particular, we consider entropy-stable schemes for the fluid part. The resulting schemes are second-order accurate, entropy stable, and preserve the divergence constraints of the electromagnetic fields. We use explicit and IMEX-based time discretizations. We then test these schemes using several one- and two-dimensional test cases. We also compare the divergence constraint errors of the proposed schemes with schemes having no divergence constraints treatment and schemes based on the PHM-based divergence cleaning. |
| title | Second order divergence constraint preserving schemes for two-fluid relativistic plasma flow equations |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2503.20372 |