Abstract
When a numerical simulation has to handle a physics problem with a wide range of time-dependent length scales, dynamically adaptive discretizations can be the method of choice. We present a major upgrade to the numerical relativity code bamps in the form of fully adaptive, physics-agnostic hp refinement. We describe the foundations of mesh refinement in the context of spectral element methods, the precise algorithm used to perform refinement in bamps, as well as several indicator functions used to drive it. Finally, we test the performance, scaling, and the accuracy of the code in treating several 1D and 2D example problems, showing clear improvements over static mesh configurations. In particular, we consider a simple nonlinear wave equation, the evolution of a real scalar field minimally coupled to gravity, as well as nonlinear gravitational waves.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
