Use puncture initial data

sxs-collaboration · Jan 14, 2025 · a2b935d · a2b935d
1 parent f4ff7f0
commit a2b935d
Show file tree

Hide file tree

Showing 3 changed files with 11 additions and 14 deletions.
diff --git a/src/Evolution/Executables/CurvedScalarWave/EvolveWorldtubeCurvedScalarWave.cpp b/src/Evolution/Executables/CurvedScalarWave/EvolveWorldtubeCurvedScalarWave.cpp
@@ -9,11 +9,13 @@
 #include "Domain/Creators/TimeDependence/RegisterDerivedWithCharm.hpp"
 #include "Domain/FunctionsOfTime/RegisterDerivedWithCharm.hpp"
 #include "Evolution/Systems/CurvedScalarWave/BoundaryCorrections/RegisterDerived.hpp"
+#include "Evolution/Systems/CurvedScalarWave/Worldtube/InitialData/ZerothOrderPuncture.hpp"
 #include "Parallel/CharmMain.tpp"
 #include "Utilities/Serialization/RegisterDerivedClassesWithCharm.hpp"
 
-using metavariables = EvolutionMetavars<gr::Solutions::KerrSchild,
-                                        ScalarWave::Solutions::PlaneWave<3>>;
+using metavariables =
+    EvolutionMetavars<gr::Solutions::KerrSchild,
+                      CurvedScalarWave::AnalyticData::ZerothOrderPuncture>;
 
 extern "C" void CkRegisterMainModule() {
   Parallel::charmxx::register_main_module<metavariables>();

diff --git a/src/Evolution/Executables/CurvedScalarWave/EvolveWorldtubeCurvedScalarWave.hpp b/src/Evolution/Executables/CurvedScalarWave/EvolveWorldtubeCurvedScalarWave.hpp
@@ -227,7 +227,6 @@ struct EvolutionMetavars {
                    tmpl::list<PhaseControl::VisitAndReturn<
                                   Parallel::Phase::LoadBalancing>,
                               PhaseControl::CheckpointAndExitAfterWallclock>>,
-
         tmpl::pair<StepChooser<StepChooserUse::Slab>,
                    tmpl::push_back<StepChoosers::standard_slab_choosers<
                                        system, local_time_stepping>,

diff --git a/tests/InputFiles/CurvedScalarWave/WorldtubeKerrSchild.yaml b/tests/InputFiles/CurvedScalarWave/WorldtubeKerrSchild.yaml
@@ -12,14 +12,10 @@ Parallelization:
   ElementDistribution: NumGridPoints
 
 AnalyticData:
-  PlaneWave:
-    WaveVector: [0., 0., 0.]
-    Center: [1., 0., 0.]
-    Profile:
-      Sinusoid:
-        Amplitude: 0.
-        Wavenumber: 1.
-        Phase: 1.
+  ZerothOrderPuncture:
+    ParticlePosition: [&XCoordWorldtube 10., 0., 0.]
+    ParticleVelocity: [0., 0.01, 0.]
+    ParticleCharge: &Charge 0.5
 
 BackgroundSpacetime:
   KerrSchild:
@@ -55,7 +51,7 @@ DomainCreator:
     ObjectA:
       InnerRadius: 1.6
       OuterRadius: 2.
-      XCoord: 5.
+      XCoord: *XCoordWorldtube
       Interior:
         ExciseWithBoundaryCondition: Worldtube
       UseLogarithmicMap: false
@@ -76,7 +72,7 @@ DomainCreator:
         AsymptoticVelocityOuterBoundary: 0.
         DecayTimescaleOuterBoundaryVelocity: 1.
       RotationMap:
-        InitialAngularVelocity: [0., 0., 0.08944271909999159]
+        InitialAngularVelocity: [0., 0., 0.001]
       TranslationMap:
         InitialValues: [[0.0, 0.0, 0.0], [0.0, 0.0, 0.0], [0.0, 0.0, 0.0]]
       ShapeMapA:
@@ -141,7 +137,7 @@ EventsAndTriggersAtSlabs:
 Worldtube:
   ExcisionSphere: ExcisionSphereA
   ExpansionOrder: 1
-  Charge: 0.5
+  Charge: *Charge
   SelfForceOptions: None
   ObserveCoefficientsTrigger:
     Slabs: