Merge pull request #61 from jrgparkinson/min-dt

Allow configuring the min dt for convergence
jrgparkinson · Aug 10, 2023 · dea44db · dea44db
2 parents 696c38d + d95db29
commit dea44db
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 1 deletion.
diff --git a/src/mushyLayerOpt.h b/src/mushyLayerOpt.h
@@ -423,6 +423,9 @@ struct MushyLayerOptions
   /// Minimum dt allowed
   Real minDt;
 
+  // dt below which we don't consider the simulation to have converged
+  Real min_dt_convergence;
+
   /// Maximum dt allowed
   Real max_dt;
 

diff --git a/srcSubcycle/AMRLevelMushyLayer.cpp b/srcSubcycle/AMRLevelMushyLayer.cpp
@@ -218,7 +218,7 @@ bool AMRLevelMushyLayer::convergedToSteadyState()
   //    m_doAutomaticRestart = false;
   //  }
 
-  if (hasConverged && m_dt < 1e-10)
+  if (hasConverged && m_dt < m_opt.min_dt_convergence)
   {
     LOG_INFO("All fields converged but dt < 1e-10 so keep solving");
     return false;

diff --git a/srcSubcycle/MushyLayerSubcycleUtils.cpp b/srcSubcycle/MushyLayerSubcycleUtils.cpp
@@ -908,6 +908,8 @@ void getAMRFactory(RefCountedPtr<AMRLevelMushyLayerFactory> &a_fact)
   //  opt.maxEta = -1;
   //  ppMain.query("max_eta", opt.maxEta); // let user specify different max eta
   //  if they want
+  opt.min_dt_convergence = 1e-10;
+  ppMain.query("min_dt_convergence", opt.min_dt_convergence);
 
   opt.minDt = 1e-7;
   ppMain.query("min_dt", opt.minDt);