|
|
using | options = tmpl::list< TypeOptionTag > |
| |
|
using | dg_interior_evolved_variables_tags = tmpl::list< gr::Tags::SpacetimeMetric< Dim, Frame::Inertial, DataVector >, Tags::Pi< Dim, Frame::Inertial >, Tags::Phi< Dim, Frame::Inertial > > |
| |
|
using | dg_interior_temporary_tags = tmpl::list< domain::Tags::Coordinates< Dim, Frame::Inertial >, ConstraintDamping::Tags::ConstraintGamma1, ConstraintDamping::Tags::ConstraintGamma2, gr::Tags::Lapse< DataVector >, gr::Tags::Shift< Dim, Frame::Inertial, DataVector >, gr::Tags::InverseSpacetimeMetric< Dim, Frame::Inertial, DataVector >, gr::Tags::SpacetimeNormalVector< Dim, Frame::Inertial, DataVector >, gr::Tags::SpacetimeNormalOneForm< Dim, Frame::Inertial, DataVector >, Tags::ThreeIndexConstraint< Dim, Frame::Inertial >, Tags::GaugeH< Dim, Frame::Inertial >, Tags::SpacetimeDerivGaugeH< Dim, Frame::Inertial > > |
| |
|
using | dg_interior_dt_vars_tags = tmpl::list< ::Tags::dt< gr::Tags::SpacetimeMetric< Dim, Frame::Inertial, DataVector > >, ::Tags::dt< Tags::Pi< Dim, Frame::Inertial > >, ::Tags::dt< Tags::Phi< Dim, Frame::Inertial > > > |
| |
|
using | dg_interior_deriv_vars_tags = tmpl::list< ::Tags::deriv< gr::Tags::SpacetimeMetric< Dim, Frame::Inertial, DataVector >, tmpl::size_t< Dim >, Frame::Inertial >, ::Tags::deriv< Tags::Pi< Dim, Frame::Inertial >, tmpl::size_t< Dim >, Frame::Inertial >, ::Tags::deriv< Tags::Phi< Dim, Frame::Inertial >, tmpl::size_t< Dim >, Frame::Inertial > > |
| |
|
using | dg_gridless_tags = tmpl::list<> |
| |
|
using | creatable_classes = tmpl::list< ConstraintPreservingBjorhus< Dim >, DirichletAnalytic< Dim >, Outflow< Dim >, domain::BoundaryConditions::Periodic< BoundaryCondition< Dim > > > |
| |
|
|
| ConstraintPreservingBjorhus (detail::ConstraintPreservingBjorhusType type) noexcept |
| |
|
| ConstraintPreservingBjorhus (CkMigrateMessage *msg) noexcept |
| |
|
| WRAPPED_PUPable_decl_base_template (domain::BoundaryConditions::BoundaryCondition, ConstraintPreservingBjorhus) |
| |
| auto | get_clone () const noexcept -> std::unique_ptr< domain::BoundaryConditions::BoundaryCondition > override |
| |
|
void | pup (PUP::er &p) override |
| |
|
std::optional< std::string > | dg_time_derivative (gsl::not_null< tnsr::aa< DataVector, Dim, Frame::Inertial > * > dt_spacetime_metric_correction, gsl::not_null< tnsr::aa< DataVector, Dim, Frame::Inertial > * > dt_pi_correction, gsl::not_null< tnsr::iaa< DataVector, Dim, Frame::Inertial > * > dt_phi_correction, const std::optional< tnsr::I< DataVector, Dim, Frame::Inertial > > &face_mesh_velocity, const tnsr::i< DataVector, Dim, Frame::Inertial > &normal_covector, const tnsr::I< DataVector, Dim, Frame::Inertial > &, const tnsr::aa< DataVector, Dim, Frame::Inertial > &spacetime_metric, const tnsr::aa< DataVector, Dim, Frame::Inertial > &pi, const tnsr::iaa< DataVector, Dim, Frame::Inertial > &phi, const tnsr::I< DataVector, Dim, Frame::Inertial > &coords, const Scalar< DataVector > &gamma1, const Scalar< DataVector > &gamma2, const Scalar< DataVector > &lapse, const tnsr::I< DataVector, Dim, Frame::Inertial > &shift, const tnsr::AA< DataVector, Dim, Frame::Inertial > &inverse_spacetime_metric, const tnsr::A< DataVector, Dim, Frame::Inertial > &spacetime_unit_normal_vector, const tnsr::a< DataVector, Dim, Frame::Inertial > &spacetime_unit_normal_one_form, const tnsr::iaa< DataVector, Dim, Frame::Inertial > &three_index_constraint, const tnsr::a< DataVector, Dim, Frame::Inertial > &gauge_source, const tnsr::ab< DataVector, Dim, Frame::Inertial > &spacetime_deriv_gauge_source, const tnsr::aa< DataVector, Dim, Frame::Inertial > &dt_spacetime_metric, const tnsr::aa< DataVector, Dim, Frame::Inertial > &dt_pi, const tnsr::iaa< DataVector, Dim, Frame::Inertial > &dt_phi, const tnsr::iaa< DataVector, Dim, Frame::Inertial > &d_spacetime_metric, const tnsr::iaa< DataVector, Dim, Frame::Inertial > &d_pi, const tnsr::ijaa< DataVector, Dim, Frame::Inertial > &d_phi) const noexcept |
| |
|
| BoundaryCondition (BoundaryCondition &&) noexcept=default |
| |
|
BoundaryCondition & | operator= (BoundaryCondition &&) noexcept=default |
| |
|
| BoundaryCondition (const BoundaryCondition &)=default |
| |
|
BoundaryCondition & | operator= (const BoundaryCondition &)=default |
| |
|
| BoundaryCondition (CkMigrateMessage *msg) noexcept |
| |
|
void | pup (PUP::er &p) override |
| |
 Public Member Functions inherited from domain::BoundaryConditions::BoundaryCondition |
|
| BoundaryCondition (BoundaryCondition &&) noexcept=default |
| |
|
BoundaryCondition & | operator= (BoundaryCondition &&) noexcept=default |
| |
|
| BoundaryCondition (const BoundaryCondition &)=default |
| |
|
BoundaryCondition & | operator= (const BoundaryCondition &)=default |
| |
|
| BoundaryCondition (CkMigrateMessage *const msg) noexcept |
| |
|
| WRAPPED_PUPable_abstract (BoundaryCondition) |
| |
template<size_t Dim>
class GeneralizedHarmonic::BoundaryConditions::ConstraintPreservingBjorhus< Dim >
Sets constraint preserving boundary conditions using the Bjorhus method.
Details
Boundary conditions for the generalized harmonic evolution system can be divided in to three parts, constraint-preserving, physical and gauge boundary conditions.
The generalized harmonic (GH) evolution system is a first-order reduction of Einstein equations brought about by the imposition of GH gauge. This introduces constraints on the free (evolved) variables in addition to the standard Hamiltonian and momentum constraints. The constraint-preserving portion of the boundary conditions is designed to prevent the influx of constraint violations from external faces of the evolution domain, by damping them away on a controlled and short time-scale. These conditions are imposed as corrections to the projections of the right-hand-sides of the GH evolution equations (i.e. using Bjorhus' method [15]), and are written down in Eq. (63) - (65) of [72] . The gauge degrees of freedom are controlled by imposing a Sommerfeld-type condition ( \(L=0\) member of the hierarchy derived in [11]) that allow gauge perturbations to pass through the boundary without strong reflections. These assume a spherical outer boundary, and can be written down as in Eq. (25) of [93] . Finally, the physical boundary conditions control the influx of inward propagating gravitational-wave solutions from the external boundaries. These are derived by considering the evolution system of the Weyl curvature tensor, and controlling the inward propagating characteristics of the system that are proportional to the Newman-Penrose curvature spinor components \(\Psi_4\) and \(\Psi_0\). Here we use Eq. (68) of [72] to disallow any incoming waves. It is to be noted that all the above conditions are also imposed on characteristic modes with speeds exactly zero.
This class provides two choices of combinations of the above corrections:
ConstraintPreserving : this imposes the constraint-preserving and gauge-controlling corrections;ConstraintPreservingPhysical : this additionally restricts the influx of any physical gravitational waves from the outer boundary, in addition to preventing the influx of constraint violations and gauge perturbations.
We refer to Bjorhus::constraint_preserving_bjorhus_corrections_dt_v_psi(), Bjorhus::constraint_preserving_bjorhus_corrections_dt_v_zero(), Bjorhus::constraint_preserving_bjorhus_corrections_dt_v_minus(), and Bjorhus::constraint_preserving_physical_bjorhus_corrections_dt_v_minus() for the further details on implementation.
- Note
- These boundary conditions assume a spherical outer boundary. Also, we do not yet have an option to inject incoming gravitational waves at the outer boundary.