Related papers: Computational Relativistic Astrophysics With Adapt…
Computations have helped elucidate the dynamics of Earth's mantle for several decades already. The numerical methods that underlie these simulations have greatly evolved within this time span, and today include dynamically changing and…
This article deals with the adaptive and approximative computation of the Lam\'e equations. The equations of linear elasticity are considered as boundary integral equations and solved in the setting of the boundary element method (BEM).…
The averaged alternating modified reflections (AAMR) method is a projection algorithm for finding the closest point in the intersection of convex sets to any arbitrary point in a Hilbert space. This method can be seen as an adequate…
We present CubismAMR, a C++ library for distributed simulations with block-structured grids and Adaptive Mesh Refinement. A numerical method to solve the incompressible Navier-Stokes equations is proposed, that comes with a novel approach…
GPUs and other accelerators are increasingly used for scientific computing. In the future, we want to add GPU support to parallel adaptive mesh refinement (AMR) codes written in Fortran. To understand which changes are necessary to obtain…
Computationally solving the equations of elasticity is a key component in many materials science and mechanics simulations. Phenomena such as deformation-induced microstructure evolution, microfracture, and microvoid nucleation are examples…
We describe a method for the numerical solution of the relativistic Euler equations which we have found to be both robust and efficient, and which has enabled us to simulate relativistic jets. The technique employs a solver of the…
Computing the inverse covariance matrix (or precision matrix) of large data vectors is crucial in weak lensing (and multi-probe) analyses of the large scale structure of the universe. Analytically computed covariances are noise-free and…
We present here the result of continuation work, performed to further fulfill the vision we outlined in [Harel,Lekien,P\'eba\"y-2017] for the visualization and analysis of tree-based adaptive mesh refinement (AMR) simulations, using the…
We present preliminar results and tests of a new general relativistic code to simulate the hydrodynamic collapse of a 21 solar masses star. We have assumed spherical symmetry and used the formalism of Misner and Sharp to construct a…
Data analysis of gravitational waves detected by the Ligo-Virgo-Kagra collaboration and future observatories relies on precise modelling of the sources. In order to build, calibrate and validate current models, we resort to expensive…
We discuss refinement criteria for the Berger-Rigoutsos (block-based) refinement algorithm in our numerical relativity code GR-Athena++ in the context of binary black hole merger simulations. We compare three different strategies: the…
Algorithms that promise to leverage resources of quantum computers efficiently to accelerate the finite element method have emerged. However, the finite element method is usually incorporated into a high-level numerical scheme which allows…
Training deep reinforcement learning (RL) agents necessitates overcoming the highly unstable nonconvex stochastic optimization inherent in the trial-and-error mechanism. To tackle this challenge, we propose a physics-inspired optimization…
A new code, named MAP, is written in Fortran language for magnetohydrodynamics (MHD) calculation with the adaptive mesh refinement (AMR) and Message Passing Interface (MPI) parallelization. There are several optional numerical schemes for…
Adaptive Mesh Refinement (AMR) is becoming a prevalent data representation for scientific visualization. Resulting from large fluid mechanics simulations, the data is usually cell centric, imposing a number of challenges for high quality…
The computational modeling of high-speed flows (e.g. hypersonic) and space plasmas is characterized by a plethora of complex physical phenomena, in particular involving strong oblique shocks, bow shocks and/or shock waves boundary layer…
We propose a new adaptive algorithm for the approximation of the Landau-Lifshitz-Gilbert equation via a higher-order tangent plane scheme. We show that the adaptive approximation satisfies an energy inequality and demonstrate numerically,…
Einstein's field equation of General Relativity (GR) has been known for over 100 years, yet it remains challenging to solve analytically in strongly relativistic regimes, particularly where there is a lack of a priori symmetry. Numerical…
We describe an implicit general relativistic hydrodynamics code. The evolution equations are formulated in comoving coordinates. A conservative finite differencing of the Einstein equations is outlined, and artificial viscosity and…