SpECTRE
v2025.01.30
|
Initializes the constraint damping parameters \(\gamma_1\) and \(\gamma_2\). More...
#include <InitializeConstraintGammas.hpp>
Public Types | |
using | return_tags = tmpl::list< CurvedScalarWave::Tags::ConstraintGamma1, CurvedScalarWave::Tags::ConstraintGamma2 > |
using | argument_tags = tmpl::list< domain::Tags::Coordinates< Dim, Frame::Inertial > > |
using | simple_tags = return_tags |
using | compute_tags = tmpl::list<> |
using | simple_tags_from_options = tmpl::list<> |
using | const_global_cache_tags = tmpl::list<> |
using | mutable_global_cache_tags = tmpl::list<> |
Static Public Member Functions | |
static void | apply (const gsl::not_null< Scalar< DataVector > * > gamma1, const gsl::not_null< Scalar< DataVector > * > gamma2, const tnsr::I< DataVector, Dim, Frame::Inertial > &inertial_coords) |
Initializes the constraint damping parameters \(\gamma_1\) and \(\gamma_2\).
The constraint damping parameters of the CurvedScalarWave system are initialized according to a Gaussian
\begin{align*} \gamma_1 &= 0 \\ \gamma_2 &= A e^{-(\sigma r)^2 } + c, \end{align*}
where \(r = \sqrt{\delta_{ij} x^i x^j}\) and with \(A = 10\), \(\sigma = 10^{-1}\), \(c = 10^{-4}\). These values were experimentally found to ensure a stable evolution when using the worldtube scheme.
DataBox changes: