fix cubic profile and add user input for damping strength

RevathiJambunathan · Oct 6, 2023 · aeb088f · aeb088f
1 parent 561a34a
commit aeb088f
Show file tree

Hide file tree

Showing 3 changed files with 19 additions and 0 deletions.
diff --git a/Source/BoundaryConditions/PML.cpp b/Source/BoundaryConditions/PML.cpp
@@ -212,9 +212,23 @@ SigmaBox::SigmaBox (const Box& box, const BoxArray& grids, const Real* dx, const
         sigma_star_cumsum_fac[idim].m_hi = hi[idim]+1;
     }
 
+#ifdef PULSAR
+    int pml_cubic_sigma = Pulsar::m_pml_cubic_sigma;
+    amrex::Real damping_strength = Pulsar::m_damping_strength;
+#endif
     Array<Real,AMREX_SPACEDIM> fac;
     for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
+#ifdef PULSAR
+        if (pml_cubic_sigma == 1) {
+            fac[idim] = damping_strength * PhysConst::c
+                        / ( dx[idim]*static_cast<Real>(delta[idim]*delta[idim]*delta[idim]) );
+        } else {
+            fac[idim] = damping_strength * PhysConst::c
+                        / ( dx[idim]*static_cast<Real>(delta[idim]*delta[idim]) );
+        }
+#else
         fac[idim] = 4.0_rt*PhysConst::c/(dx[idim]*static_cast<Real>(delta[idim]*delta[idim]));
+#endif
     }
 
     if (regdomain.ok()) { // The union of the regular grids is a single box

diff --git a/Source/Particles/PulsarParameters.H b/Source/Particles/PulsarParameters.H
@@ -704,6 +704,7 @@ public:
     static int m_pml_cubic_sigma;
     static amrex::Real m_gammarad_real;
     static amrex::Real m_gammarad_scaled;
+    static amrex::Real m_damping_strength;
 private:
 };
 

diff --git a/Source/Particles/PulsarParameters.cpp b/Source/Particles/PulsarParameters.cpp
@@ -145,6 +145,7 @@ amrex::Real Pulsar::m_depositbuffer_max = 0.;
 int Pulsar::m_pml_cubic_sigma;
 amrex::Real Pulsar::m_gammarad_real = 1.;
 amrex::Real Pulsar::m_gammarad_scaled = 1.;
+amrex::Real Pulsar::m_damping_strength = 4.;
 
 Pulsar::Pulsar ()
 {
@@ -378,6 +379,9 @@ Pulsar::ReadParameters () {
     pp.get("pml_cubic_sigma", m_pml_cubic_sigma);
     pp.query("gammarad_real", m_gammarad_real);
     pp.query("gammarad_scaled", m_gammarad_scaled);
+    pp.get("damping_strength", m_damping_strength);
+    amrex::Print() << " damping strength in PML " << m_damping_strength << "\n";
+    amrex::Print() << " cubic pml flag : " << m_pml_cubic_sigma  << "\n";
 }