-
Simflowny 3: An upgraded platform for scientific modelling and simulation
Authors:
C. Palenzuela,
B. Miñano,
A. Arbona,
C. Bona-Casas,
C. Bona,
J. Massó
Abstract:
Simflowny is an open platform which automatically generates efficient parallel code of scientific dynamical models for different simulation frameworks. Here we present major upgrades on this software to support simultaneously a quite generic family of partial differential equations. These equations can be discretized using: (i) standard finite-difference for systems with derivatives up to any orde…
▽ More
Simflowny is an open platform which automatically generates efficient parallel code of scientific dynamical models for different simulation frameworks. Here we present major upgrades on this software to support simultaneously a quite generic family of partial differential equations. These equations can be discretized using: (i) standard finite-difference for systems with derivatives up to any order, (ii) High-Resolution-Shock-Capturing methods to deal with shocks and discontinuities of balance law equations, and (iii) particle-based methods. We have improved the adaptive-mesh-refinement algorithms to preserve the convergence order of the numerical methods, which is a requirement for improving scalability. Finally, we have also extended our graphical user interface (GUI) to accommodate these and future families of equations. This paper summarizes the formal representation and implementation of these new families, providing several validation results.
△ Less
Submitted 2 October, 2020;
originally announced October 2020.
-
A Simflowny-based high-performance 3D code for the generalized induction equation
Authors:
Daniele Viganò,
David Martínez-Gómez,
José A. Pons,
Carlos Palenzuela,
Federico Carrasco,
Borja Miñano,
Antoni Arbona,
Carles Bona,
Joan Massó
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…
▽ More
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 {\em Simflowny} platform. The code uses high-order numerical schemes for the time and spatial discretization, and relies on the highly-scalable {\em 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.
△ Less
Submitted 11 December, 2018; v1 submitted 20 November, 2018;
originally announced November 2018.
-
A Simflowny-based finite-difference code for high-performance computing in Relativity
Authors:
Carlos Palenzuela,
Borja Miñano,
Daniele Viganò,
Antoni Arbona,
Carles Bona-Casas,
Andreu Rigo,
Miguel Bezares,
Carles Bona,
Joan Massó
Abstract:
The tremendous challenge of comparing our theoretical models with the gravitational-wave observations in the new era of multimessenger astronomy requires accurate and fast numerical simulations of complicated physical systems described by the Einstein and the matter equations. These requirements can only be satisfied if the simulations can be parallelized efficiently on a large number of processor…
▽ More
The tremendous challenge of comparing our theoretical models with the gravitational-wave observations in the new era of multimessenger astronomy requires accurate and fast numerical simulations of complicated physical systems described by the Einstein and the matter equations. These requirements can only be satisfied if the simulations can be parallelized efficiently on a large number of processors and advanced computational strategies are adopted. To achieve this goal we have developed Simflowny, an open platform for scientific dynamical models which automatically generates parallel code for different simulation frameworks, allowing the use of HPC infrastructures to non-specialist scientists. One of these frameworks is SAMRAI, a mature patch-based structured adaptive mesh refinement infrastructure, capable of reaching exascale in some specific problems. Here we present the numerical techniques that we have implemented on this framework by using Simflowny in order to perform fast, efficient, accurate and highly-scalable simulations. These techniques involve high-order schemes for smooth and non-smooth solutions, Adaptive Mesh Refinement with arbitrary resolution ratios and an optimal strategy for the sub-cycling in time. We validate the automatically generated codes for the SAMRAI infrastructure with some simple test examples (i.e., wave equation and Newtonian MHD) and finally with the Einstein equations.
△ Less
Submitted 11 June, 2018;
originally announced June 2018.
-
Linear high-resolution schemes for hyperbolic conservation laws: TVB numerical evidence
Authors:
C. Bona,
C. Bona-Casas,
J. Terradas
Abstract:
The Osher-Chakrabarthy family of linear flux-modification schemes is considered. Improved lower bounds on the compression factors are provided, which suggest the viability of using the unlimited version. The LLF flux formula is combined with these schemes in order to obtain efficient finite-difference algorithms. The resulting schemes are applied to a battery of numerical tests, going from advec…
▽ More
The Osher-Chakrabarthy family of linear flux-modification schemes is considered. Improved lower bounds on the compression factors are provided, which suggest the viability of using the unlimited version. The LLF flux formula is combined with these schemes in order to obtain efficient finite-difference algorithms. The resulting schemes are applied to a battery of numerical tests, going from advection and Burgers equations to Euler and MHD equations, including the double Mach reflection and the Orszag-Tang 2D vortex problem. Total-variation-bounded behavior is evident in all cases, even with time-independent upper bounds. The proposed schemes, however, do not deal properly with compound shocks, arising from non-convex fluxes, as shown by Buckley-Leverett test simulations.
△ Less
Submitted 12 December, 2008; v1 submitted 13 October, 2008;
originally announced October 2008.