SpECTRE  v2023.05.16
CurvedScalarWave::Worldtube::Initialization::InitializeConstraintDampingGammas< Dim > Struct Template Reference

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)

Detailed Description

template<size_t Dim>
struct CurvedScalarWave::Worldtube::Initialization::InitializeConstraintDampingGammas< Dim >

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:

The documentation for this struct was generated from the following file: