A Simflowny-based high-performance 3D code for the generalized induction equation

Please use this identifier to cite or link to this item: http://hdl.handle.net/10045/88652
Información del item - Informació de l'item - Item information
Title: A Simflowny-based high-performance 3D code for the generalized induction equation
Authors: Viganò, Daniele | Martínez-Gómez, David | Pons, José A. | Palenzuela, Carlos | Carrasco, Federico | Miñano, Borja | Arbona, Antoni | Bona, Carles | Massó, Joan
Research Group/s: Astrofísica Relativista
Center, Department or Service: Universidad de Alicante. Departamento de Física Aplicada
Keywords: MHD | Neutron stars | Simflowny | Magnetic fields | High resolution shock capturing | AMR
Knowledge Area: Astronomía y Astrofísica
Issue Date: Apr-2019
Publisher: Elsevier
Citation: Computer Physics Communications. 2019, 237: 168-183. doi:10.1016/j.cpc.2018.11.022
Abstract: In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage and a consensus on the general dynamics at least for some simple initial data, only few attempts have been performed in 3D, due to the computational costs and the need for a proper numerical treatment of the intrinsic non-linearity of the equation. Here, we carefully analyze the general induction equation, studying its characteristic structure, and we present a new Cartesian 3D code, generated by the user-friendly, publicly available Simflowny platform. The code uses high-order numerical schemes for the time and spatial discretization, and relies on the highly-scalable SAMRAI architecture for the adaptive mesh refinement. We present the application of the code to several benchmark tests, showing the high order of convergence and accuracy achieved and the capabilities in terms of magnetic shock resolution and three-dimensionality. This paper paves the way for the applications to a realistic, 3D long-term evolution of neutron stars interior and, possibly, of other astrophysical sources.
Sponsor: We acknowledge support from the Spanish Ministry of Economy, Industry and Competitiveness grants AYA2016-80289-P and AYA2017-82089-ERC (AEI/FEDER, UE). CP also acknowledges support from the Spanish Ministry of Education and Science through a Ramon y Cajal grant. JAP acknowledges the Spanish MINECO/FEDER grant AYA2015-66899-C2-2-P, the grant of Generalitat Valenciana, Spain PROMETEOII-2014-069.DMacknowledges support from Vicepresidència i Conselleria d’Innovació, Recerca i Turisme del Govern de les Illes Balears, Spain. The work has been done within the PHAROS COST action CA16214.
URI: http://hdl.handle.net/10045/88652
ISSN: 0010-4655 (Print) | 1879-2944 (Online)
DOI: 10.1016/j.cpc.2018.11.022
Language: eng
Type: info:eu-repo/semantics/article
Rights: © 2018 Elsevier B.V.
Peer Review: si
Publisher version: https://doi.org/10.1016/j.cpc.2018.11.022
Appears in Collections:INV - Astrofísica Relativista - Artículos de Revistas

Files in This Item:
Files in This Item:
File Description SizeFormat 
Thumbnail2019_Vigano_etal_CompPhysComm_final.pdfVersión final (acceso restringido)4,99 MBAdobe PDFOpen    Request a copy
Thumbnail2019_Vigano_etal_CompPhysComm_preprint.pdfPreprint (acceso abierto)6,04 MBAdobe PDFOpen Preview

Items in RUA are protected by copyright, with all rights reserved, unless otherwise indicated.