Rename subsection and function

Co-authored-by: hepap <[email protected]>
Co-authored-by: Bruno Blais <[email protected]>

applications_tests/lethe-fluid/heat_transfer_phase_change_constrain_solid_domain.prm

@@ -119,10 +119,10 @@ subsection boundary conditions heat transfer
 end
 
 #---------------------------------------------------
-# Constrain Solid Domain
+# Constrain Stasis
 #---------------------------------------------------
 
-subsection constrain solid domain
+subsection constrain stasis
   set enable                = true
   set number of constraints = 1
   subsection constraint 0

applications_tests/lethe-fluid/heat_transfer_vof_phase_change_constrain_solid_domain.prm

@@ -136,10 +136,10 @@ subsection boundary conditions heat transfer
 end
 
 #---------------------------------------------------
-# Constrain Solid Domain
+# Constrain Stasis
 #---------------------------------------------------
 
-subsection constrain solid domain
+subsection constrain stasis
   set enable                = true
   set number of constraints = 2
   subsection constraint 0

doc/source/parameters/cfd/cfd.rst

@@ -12,7 +12,7 @@ General, CFD and Multiphysics
    boundary_conditions_multiphysics
    box_refinement
    cahn_hilliard
-   constrain_solid_domain
+   constrain_stasis
    dimensionality
    dynamic_flow_control
    fem

doc/source/parameters/cfd/constrain_solid_domain.rst → doc/source/parameters/cfd/constrain_stasis.rst

@@ -1,6 +1,6 @@
-=======================
-Constrain Solid Domain
-=======================
+=================
+Constrain Stasis
+=================
 
 This subsection is used to define temperature-dependent solid domains within a defined fluid.
 Homogenous constraints are applied on velocity and pressure degrees of freedom of cells found within the prescribed temperature range to mimic a solid.
@@ -12,7 +12,7 @@ The subsection with default parameters goes as follows:
 
 .. code-block:: text
 
-    subsection constrain solid domain
+    subsection constrain stasis
       set enable                = false
       set number of constraints = 0
       subsection constraint 0

include/solvers/navier_stokes_base.h

@@ -72,12 +72,13 @@ using namespace dealii;
  * @param[in] filtered_phase_fraction_tolerance Tolerance applied on filtered
  * phase fraction.
  */
-struct SolidDomainConstraint
+struct StasisConstraintWithTemperature
 {
-  SolidDomainConstraint(const unsigned int fluid_id,
-                        const double       min_solid_temperature,
-                        const double       max_solid_temperature,
-                        const double       filtered_phase_fraction_tolerance)
+  StasisConstraintWithTemperature(
+    const unsigned int fluid_id,
+    const double       min_solid_temperature,
+    const double       max_solid_temperature,
+    const double       filtered_phase_fraction_tolerance)
     : fluid_id(fluid_id)
     , min_solid_temperature(min_solid_temperature)
     , max_solid_temperature(max_solid_temperature)
@@ -478,33 +479,33 @@ class NavierStokesBase : public PhysicsSolver<VectorType>
    * @param[in] local_temperature_values Cell's local temperature values at
    * quadrature points.
    *
-   * @param[in,out] solid_domain_constraint_struct Struct containing flagged DOF
+   * @param[in,out] stasis_constraint_struct Struct containing flagged DOF
    * containers, temperature range information and fluid id.
    */
   void
   add_flags_and_constrain_velocity(
     const std::vector<types::global_dof_index> &local_dof_indices,
     const std::vector<double>                  &local_temperature_values,
-    SolidDomainConstraint                      &solid_domain_constraint_struct);
+    StasisConstraintWithTemperature            &stasis_constraint_struct);
 
   /**
    * @brief Check if solid cells are connected to fluid ones and constrain null
    * pressure DOFs if they are not.
    *
    * The check is done by looking if the global DOF indices are located in the
    * flag containers (@p dofs_are_in_solid and @p dofs_are_connected_to_fluid)
-   * of @p solid_domain_constraint_struct.
+   * of @p stasis_constraint_struct.
    *
-   * @param[in] solid_domain_constraint_struct Struct containing flagged DOF
+   * @param[in] stasis_constraint_struct Struct containing flagged DOF
    * containers, temperature range information and fluid id.
    *
    * @param[in,out] local_dof_indices Vector for storing a cell's local DOF
    * indices.
    */
   void
   check_and_constrain_pressure(
-    const SolidDomainConstraint          &solid_domain_constraint_struct,
-    std::vector<types::global_dof_index> &local_dof_indices);
+    const StasisConstraintWithTemperature &stasis_constraint_struct,
+    std::vector<types::global_dof_index>  &local_dof_indices);
 
   /**
    * @brief Constrain velocity DOFs of a solid cell.
@@ -717,7 +718,7 @@ class NavierStokesBase : public PhysicsSolver<VectorType>
    * physic.
    */
   void
-  constrain_solid_domain(const DoFHandler<dim> *dof_handler_ht);
+  constrain_stasis_with_temperature(const DoFHandler<dim> *dof_handler_ht);
 
   /**
    * @brief Constrain fluids' subdomains according to the temperature field to
@@ -730,8 +731,8 @@ class NavierStokesBase : public PhysicsSolver<VectorType>
    * physic.
    */
   void
-  constrain_solid_domain_vof(const DoFHandler<dim> *dof_handler_vof,
-                             const DoFHandler<dim> *dof_handler_ht);
+  constrain_stasis_with_temperature_vof(const DoFHandler<dim> *dof_handler_vof,
+                                        const DoFHandler<dim> *dof_handler_ht);
 
   /**
    * @brief write_checkpoint
@@ -870,7 +871,7 @@ class NavierStokesBase : public PhysicsSolver<VectorType>
   std::shared_ptr<FEValues<dim>> fe_values_vof;
   /// Vector containing solid domain constraint structs for
   /// temperature-dependent solid domain constraints in VOF simulations
-  std::vector<SolidDomainConstraint> solid_domain_constraint_structs;
+  std::vector<StasisConstraintWithTemperature> stasis_constraint_structs;
   /// Dynamic homogeneous constraints used for temperature-dependent solid
   /// domain constraints
   AffineConstraints<double> dynamic_zero_constraints;

source/core/parameters.cc

@@ -504,7 +504,7 @@ namespace Parameters
     dealii::ParameterHandler &prm,
     const unsigned int        number_of_constraints)
   {
-    prm.enter_subsection("constrain solid domain");
+    prm.enter_subsection("constrain stasis");
     {
       prm.declare_entry(
         "enable",
@@ -564,7 +564,7 @@ namespace Parameters
   void
   ConstrainSolidDomain::parse_parameters(dealii::ParameterHandler &prm)
   {
-    prm.enter_subsection("constrain solid domain");
+    prm.enter_subsection("constrain stasis");
     {
       this->enable                = prm.get_bool("enable");
       this->number_of_constraints = prm.get_integer("number of constraints");

source/solvers/gls_navier_stokes.cc

@@ -399,12 +399,13 @@ GLSNavierStokesSolver<dim>::define_dynamic_zero_constraints()
     this->multiphysics->get_dof_handler(PhysicsID::heat_transfer);
 
   if (!this->simulation_parameters.multiphysics.VOF)
-    this->constrain_solid_domain(dof_handler_ht);
+    this->constrain_stasis_with_temperature(dof_handler_ht);
   else
     {
       const DoFHandler<dim> *dof_handler_vof =
         this->multiphysics->get_dof_handler(PhysicsID::VOF);
-      this->constrain_solid_domain_vof(dof_handler_vof, dof_handler_ht);
+      this->constrain_stasis_with_temperature_vof(dof_handler_vof,
+                                                  dof_handler_ht);
     }
 
   this->dynamic_zero_constraints.merge(
@@ -414,11 +415,11 @@ GLSNavierStokesSolver<dim>::define_dynamic_zero_constraints()
   this->dynamic_zero_constraints.close();
 
   // Clear sets for next time step
-  for (SolidDomainConstraint &solid_domain_constraint_struct :
-       this->solid_domain_constraint_structs)
+  for (StasisConstraintWithTemperature &stasis_constraint_struct :
+       this->stasis_constraint_structs)
     {
-      solid_domain_constraint_struct.dofs_are_in_solid.clear();
-      solid_domain_constraint_struct.dofs_are_connected_to_fluid.clear();
+      stasis_constraint_struct.dofs_are_in_solid.clear();
+      stasis_constraint_struct.dofs_are_connected_to_fluid.clear();
     }
 }
 

source/solvers/navier_stokes_base.cc

@@ -623,22 +623,21 @@ NavierStokesBase<dim, VectorType, DofsType>::
   if (!this->simulation_parameters.constrain_solid_domain.enable)
     return;
 
-  // Initialize SolidDomainConstraint structs for each constraint
+  // Initialize StasisConstraintWithTemperature structs for each constraint
   for (unsigned int c_id = 0;
        c_id <
        this->simulation_parameters.constrain_solid_domain.number_of_constraints;
        c_id++)
     {
-      SolidDomainConstraint solid_domain_constraint_struct(
+      StasisConstraintWithTemperature stasis_constraint_struct(
         this->simulation_parameters.constrain_solid_domain.fluid_ids[c_id],
         this->simulation_parameters.constrain_solid_domain
           .temperature_min_values[c_id],
         this->simulation_parameters.constrain_solid_domain
           .temperature_max_values[c_id],
         this->simulation_parameters.constrain_solid_domain
           .filtered_phase_fraction_tolerance[c_id]);
-      this->solid_domain_constraint_structs.push_back(
-        solid_domain_constraint_struct);
+      this->stasis_constraint_structs.push_back(stasis_constraint_struct);
     }
 
   // For temperature-dependent constraints
@@ -1821,16 +1820,16 @@ NavierStokesBase<dim, VectorType, DofsType>::establish_solid_domain(
 
 template <int dim, typename VectorType, typename DofsType>
 void
-NavierStokesBase<dim, VectorType, DofsType>::constrain_solid_domain(
+NavierStokesBase<dim, VectorType, DofsType>::constrain_stasis_with_temperature(
   const DoFHandler<dim> *dof_handler_ht)
 {
   const unsigned int                   dofs_per_cell = this->fe->dofs_per_cell;
   std::vector<types::global_dof_index> local_dof_indices(dofs_per_cell);
 
   // Get struct containing temperature range information and flag
   // containers for DOFs.
-  SolidDomainConstraint &solid_domain_constraint_struct =
-    this->solid_domain_constraint_structs[0];
+  StasisConstraintWithTemperature &stasis_constraint_struct =
+    this->stasis_constraint_structs[0];
 
   // Get temperature solution
   const auto temperature_solution =
@@ -1848,18 +1847,17 @@ NavierStokesBase<dim, VectorType, DofsType>::constrain_solid_domain(
                                       local_temperature_values);
           add_flags_and_constrain_velocity(local_dof_indices,
                                            local_temperature_values,
-                                           solid_domain_constraint_struct);
+                                           stasis_constraint_struct);
         }
     }
-  check_and_constrain_pressure(solid_domain_constraint_struct,
-                               local_dof_indices);
+  check_and_constrain_pressure(stasis_constraint_struct, local_dof_indices);
 }
 
 template <int dim, typename VectorType, typename DofsType>
 void
-NavierStokesBase<dim, VectorType, DofsType>::constrain_solid_domain_vof(
-  const DoFHandler<dim> *dof_handler_vof,
-  const DoFHandler<dim> *dof_handler_ht)
+NavierStokesBase<dim, VectorType, DofsType>::
+  constrain_stasis_with_temperature_vof(const DoFHandler<dim> *dof_handler_vof,
+                                        const DoFHandler<dim> *dof_handler_ht)
 {
   const unsigned int                   dofs_per_cell = this->fe->dofs_per_cell;
   std::vector<types::global_dof_index> local_dof_indices(dofs_per_cell);
@@ -1877,8 +1875,8 @@ NavierStokesBase<dim, VectorType, DofsType>::constrain_solid_domain_vof(
 
   // Loop over structs containing fluid id, temperature range information and
   // flag containers for DOFs.
-  for (SolidDomainConstraint &solid_domain_constraint_struct :
-       this->solid_domain_constraint_structs)
+  for (StasisConstraintWithTemperature &stasis_constraint_struct :
+       this->stasis_constraint_structs)
     {
       for (const auto &cell : dof_handler.active_cell_iterators())
         {
@@ -1899,9 +1897,9 @@ NavierStokesBase<dim, VectorType, DofsType>::constrain_solid_domain_vof(
               for (const double &filtered_phase_fraction :
                    local_filtered_phase_fraction_values)
                 {
-                  if (abs(solid_domain_constraint_struct.fluid_id -
+                  if (abs(stasis_constraint_struct.fluid_id -
                           filtered_phase_fraction) >=
-                      solid_domain_constraint_struct
+                      stasis_constraint_struct
                         .filtered_phase_fraction_tolerance)
                     {
                       cell_is_in_right_fluid = false;
@@ -1916,7 +1914,7 @@ NavierStokesBase<dim, VectorType, DofsType>::constrain_solid_domain_vof(
                 {
                   flag_dofs_connected_to_fluid(
                     local_dof_indices,
-                    solid_domain_constraint_struct.dofs_are_connected_to_fluid);
+                    stasis_constraint_struct.dofs_are_connected_to_fluid);
                   continue;
                 }
 
@@ -1926,11 +1924,10 @@ NavierStokesBase<dim, VectorType, DofsType>::constrain_solid_domain_vof(
                                           local_temperature_values);
               add_flags_and_constrain_velocity(local_dof_indices,
                                                local_temperature_values,
-                                               solid_domain_constraint_struct);
+                                               stasis_constraint_struct);
             }
         }
-      check_and_constrain_pressure(solid_domain_constraint_struct,
-                                   local_dof_indices);
+      check_and_constrain_pressure(stasis_constraint_struct, local_dof_indices);
     }
 }
 
@@ -1939,18 +1936,18 @@ void
 NavierStokesBase<dim, VectorType, DofsType>::add_flags_and_constrain_velocity(
   const std::vector<types::global_dof_index> &local_dof_indices,
   const std::vector<double>                  &local_temperature_values,
-  SolidDomainConstraint                      &solid_domain_constraint_struct)
+  StasisConstraintWithTemperature            &stasis_constraint_struct)
 {
   bool cell_is_solid = false;
 
   // Flag cell if the cell is solid and constrain its velocity DOFs
   for (const double &temperature : local_temperature_values)
     {
-      if (temperature >= solid_domain_constraint_struct.min_solid_temperature &&
-          temperature <= solid_domain_constraint_struct.max_solid_temperature)
+      if (temperature >= stasis_constraint_struct.min_solid_temperature &&
+          temperature <= stasis_constraint_struct.max_solid_temperature)
         {
           flag_dofs_in_solid(local_dof_indices,
-                             solid_domain_constraint_struct.dofs_are_in_solid);
+                             stasis_constraint_struct.dofs_are_in_solid);
           constrain_solid_cell_velocity_dofs(false,
                                              local_dof_indices,
                                              this->dynamic_zero_constraints);
@@ -1963,28 +1960,27 @@ NavierStokesBase<dim, VectorType, DofsType>::add_flags_and_constrain_velocity(
   // be constrained.
   if (!cell_is_solid)
     flag_dofs_connected_to_fluid(
-      local_dof_indices,
-      solid_domain_constraint_struct.dofs_are_connected_to_fluid);
+      local_dof_indices, stasis_constraint_struct.dofs_are_connected_to_fluid);
 }
 
 template <int dim, typename VectorType, typename DofsType>
 void
 NavierStokesBase<dim, VectorType, DofsType>::check_and_constrain_pressure(
-  const SolidDomainConstraint          &solid_domain_constraint_struct,
-  std::vector<types::global_dof_index> &local_dof_indices)
+  const StasisConstraintWithTemperature &stasis_constraint_struct,
+  std::vector<types::global_dof_index>  &local_dof_indices)
 {
   for (const auto &cell : dof_handler.active_cell_iterators())
     {
       if (cell->is_locally_owned() || cell->is_ghost())
         {
           cell->get_dof_indices(local_dof_indices);
-          bool cell_is_solid = check_cell_is_in_solid(
-            solid_domain_constraint_struct.dofs_are_in_solid,
-            local_dof_indices);
+          bool cell_is_solid =
+            check_cell_is_in_solid(stasis_constraint_struct.dofs_are_in_solid,
+                                   local_dof_indices);
           if (cell_is_solid)
             {
               bool connected_to_fluid = check_cell_is_connected_to_fluid(
-                solid_domain_constraint_struct.dofs_are_connected_to_fluid,
+                stasis_constraint_struct.dofs_are_connected_to_fluid,
                 local_dof_indices);
               if (!connected_to_fluid)
                 {