Accurately simulating anisotropic thermal conduction on a moving mesh
We present a novel implementation of an extremum preserving anisotropic diffusion solver for thermal conduction on the unstructured moving Voronoi mesh of the Arepo code. The method relies on splitting the one-sided facet fluxes into normal and oblique components, with the oblique fluxes being limit...
| Main Authors: | , , |
|---|---|
| Other Authors: | |
| Format: | Article |
| Language: | en_US |
| Published: |
Oxford University Press
2017
|
| Online Access: | http://hdl.handle.net/1721.1/108512 https://orcid.org/0000-0002-3074-2326 https://orcid.org/0000-0003-3816-7028 https://orcid.org/0000-0001-8593-7692 |
| _version_ | 1826188737878425600 |
|---|---|
| author | Kannan, Rahul Marinacci, Federico Vogelsberger, Mark |
| author2 | Massachusetts Institute of Technology. Department of Physics |
| author_facet | Massachusetts Institute of Technology. Department of Physics Kannan, Rahul Marinacci, Federico Vogelsberger, Mark |
| author_sort | Kannan, Rahul |
| collection | MIT |
| description | We present a novel implementation of an extremum preserving anisotropic diffusion solver for thermal conduction on the unstructured moving Voronoi mesh of the Arepo code. The method relies on splitting the one-sided facet fluxes into normal and oblique components, with the oblique fluxes being limited such that the total flux is both locally conservative and extremum preserving. The approach makes use of harmonic averaging points and a simple, robust interpolation scheme that works well for strong heterogeneous and anisotropic diffusion problems. Moreover, the required discretization stencil is small. Efficient fully implicit and semi-implicit time integration schemes are also implemented. We perform several numerical tests that evaluate the stability and accuracy of the scheme, including applications such as point explosions with heat conduction and calculations of convective instabilities in conducting plasmas. The new implementation is suitable for studying important astrophysical phenomena, such as the conductive heat transport in galaxy clusters, the evolution of supernova remnants, or the distribution of heat from black hole-driven jets into the intracluster medium. |
| first_indexed | 2024-09-23T08:04:16Z |
| format | Article |
| id | mit-1721.1/108512 |
| institution | Massachusetts Institute of Technology |
| language | en_US |
| last_indexed | 2024-09-23T08:04:16Z |
| publishDate | 2017 |
| publisher | Oxford University Press |
| record_format | dspace |
| spelling | mit-1721.1/1085122022-09-30T07:17:02Z Accurately simulating anisotropic thermal conduction on a moving mesh Kannan, Rahul Marinacci, Federico Vogelsberger, Mark Massachusetts Institute of Technology. Department of Physics MIT Kavli Institute for Astrophysics and Space Research Kannan, Rahul Marinacci, Federico Vogelsberger, Mark We present a novel implementation of an extremum preserving anisotropic diffusion solver for thermal conduction on the unstructured moving Voronoi mesh of the Arepo code. The method relies on splitting the one-sided facet fluxes into normal and oblique components, with the oblique fluxes being limited such that the total flux is both locally conservative and extremum preserving. The approach makes use of harmonic averaging points and a simple, robust interpolation scheme that works well for strong heterogeneous and anisotropic diffusion problems. Moreover, the required discretization stencil is small. Efficient fully implicit and semi-implicit time integration schemes are also implemented. We perform several numerical tests that evaluate the stability and accuracy of the scheme, including applications such as point explosions with heat conduction and calculations of convective instabilities in conducting plasmas. The new implementation is suitable for studying important astrophysical phenomena, such as the conductive heat transport in galaxy clusters, the evolution of supernova remnants, or the distribution of heat from black hole-driven jets into the intracluster medium. 2017-04-28T20:10:04Z 2017-04-28T20:10:04Z 2016-02 2016-02 Article http://purl.org/eprint/type/JournalArticle 0035-8711 1365-2966 http://hdl.handle.net/1721.1/108512 Kannan, Rahul; Springel, Volker; Pakmor, Rüdiger; Marinacci, Federico and Vogelsberger, Mark. “Accurately Simulating Anisotropic Thermal Conduction on a Moving Mesh.” Monthly Notices of the Royal Astronomical Society 458, no. 1 (February 8, 2016): 410–424. https://orcid.org/0000-0002-3074-2326 https://orcid.org/0000-0003-3816-7028 https://orcid.org/0000-0001-8593-7692 en_US http://dx.doi.org/10.1093/mnras/stw294 Monthly Notices of the Royal Astronomical Society Creative Commons Attribution-Noncommercial-Share Alike http://creativecommons.org/licenses/by-nc-sa/4.0/ application/pdf Oxford University Press arXiv |
| spellingShingle | Kannan, Rahul Marinacci, Federico Vogelsberger, Mark Accurately simulating anisotropic thermal conduction on a moving mesh |
| title | Accurately simulating anisotropic thermal conduction on a moving mesh |
| title_full | Accurately simulating anisotropic thermal conduction on a moving mesh |
| title_fullStr | Accurately simulating anisotropic thermal conduction on a moving mesh |
| title_full_unstemmed | Accurately simulating anisotropic thermal conduction on a moving mesh |
| title_short | Accurately simulating anisotropic thermal conduction on a moving mesh |
| title_sort | accurately simulating anisotropic thermal conduction on a moving mesh |
| url | http://hdl.handle.net/1721.1/108512 https://orcid.org/0000-0002-3074-2326 https://orcid.org/0000-0003-3816-7028 https://orcid.org/0000-0001-8593-7692 |
| work_keys_str_mv | AT kannanrahul accuratelysimulatinganisotropicthermalconductiononamovingmesh AT marinaccifederico accuratelysimulatinganisotropicthermalconductiononamovingmesh AT vogelsbergermark accuratelysimulatinganisotropicthermalconductiononamovingmesh |