Add linear solver controls for the solve object. (idaholab#28819)

grmnptr · Oct 10, 2024 · e3622e4 · e3622e4
1 parent 00c0f69
commit e3622e4
Show file tree

Hide file tree

Showing 2 changed files with 23 additions and 2 deletions.
diff --git a/modules/navier_stokes/include/executioners/SIMPLESolve.h b/modules/navier_stokes/include/executioners/SIMPLESolve.h
@@ -11,6 +11,7 @@
 
 #include "SolveObject.h"
 #include "RhieChowMassFlux.h"
+#include "SegregatedSolverBase.h"
 
 class SIMPLESolve : public SolveObject
 {
@@ -30,7 +31,7 @@ class SIMPLESolve : public SolveObject
     mooseError("Cannot set inner solve for SIMPLESolve");
   }
 
-  /// Fetch the Rhie Chow
+  /// Fetch the Rhie Chow user object that
   void linkRhieChowUserObject();
 
   const RhieChowMassFlux & getRCUserObject() { return *_rc_uo; }
@@ -67,4 +68,13 @@ class SIMPLESolve : public SolveObject
 
   /// Pointer to the segregated RhieChow interpolation object
   RhieChowMassFlux * _rc_uo;
+
+  /// Options for the linear solver of the momentum equation
+  SIMPLESolverConfiguration _momentum_linear_control;
+
+  /// Options for the linear solver of the pressure equation
+  SIMPLESolverConfiguration _pressure_linear_control;
+
+  /// Debug parameter which allows printing the coupling and solution vectors/matrices
+  const bool _print_fields;
 };
diff --git a/modules/navier_stokes/src/executioners/SIMPLESolve.C b/modules/navier_stokes/src/executioners/SIMPLESolve.C
@@ -22,7 +22,8 @@ SIMPLESolve::SIMPLESolve(Executioner & ex)
     _momentum_system_names(getParam<std::vector<SolverSystemName>>("momentum_systems")),
     _pressure_system_name(getParam<SolverSystemName>("pressure_system")),
     _pressure_sys_number(_problem.linearSysNum(_pressure_system_name)),
-    _pressure_system(_problem.getLinearSystem(_pressure_sys_number))
+    _pressure_system(_problem.getLinearSystem(_pressure_sys_number)),
+    _print_fields(getParam<bool>("print_fields"))
 {
   // We fetch the system numbers for the momentum components plus add vectors
   // for removing the contribution from the pressure gradient terms.
@@ -31,6 +32,16 @@ SIMPLESolve::SIMPLESolve(Executioner & ex)
     _momentum_system_numbers.push_back(_problem.linearSysNum(_momentum_system_names[system_i]));
     _momentum_systems.push_back(&_problem.getLinearSystem(_momentum_system_numbers[system_i]));
   }
+
+  _momentum_linear_control.real_valued_data["rel_tol"] = getParam<Real>("momentum_l_tol");
+  _momentum_linear_control.real_valued_data["abs_tol"] = getParam<Real>("momentum_l_abs_tol");
+  _momentum_linear_control.int_valued_data["max_its"] =
+      getParam<unsigned int>("momentum_l_max_its");
+
+  _pressure_linear_control.real_valued_data["rel_tol"] = getParam<Real>("pressure_l_tol");
+  _pressure_linear_control.real_valued_data["abs_tol"] = getParam<Real>("pressure_l_abs_tol");
+  _pressure_linear_control.int_valued_data["max_its"] =
+      getParam<unsigned int>("pressure_l_max_its");
 }
 
 bool