Move parameters to input file

Author: mrp089

Code/Source/svFSI/ComMod.h

@@ -383,6 +383,45 @@ class stModelType
     fibStrsType Tf;
 };
 
+class grModelType
+{
+  public:
+    int n_t_pre = 0;
+    int n_t_end = 0;
+    int example = 0;
+    double KsKi = 0.0;
+    double curve = 0.0;
+    double mult = 0.0;
+    double rIo = 0.0;
+    double hwaves = 0.0;
+    double lo = 0.0;
+    double phieo = 0.0;
+    double phimo = 0.0;
+    double phico = 0.0;
+    double eta = 0.0;
+    double mu = 0.0;
+    double Get = 0.0;
+    double Gez = 0.0;
+    double alpha = 0.0;
+    double cm = 0.0;
+    double dm = 0.0;
+    double Gm = 0.0;
+    double cc = 0.0;
+    double dc = 0.0;
+    double Gc = 0.0;
+    double betat = 0.0;
+    double betaz = 0.0;
+    double betad = 0.0;
+    double Tmax = 0.0;
+    double lamM = 0.0;
+    double lam0 = 0.0;
+    double KfKi = 0.0;
+    double inflam = 0.0;
+    double aexp = 0.0;
+    double delta = 0.0;
+
+};
+
 /// @brief Fluid viscosity model type
 //
 class viscModelType
@@ -439,6 +478,9 @@ class dmnType
 
     // Viscosity model for fluids
     viscModelType visc;
+
+    // G&R model for growth and remodeling
+    grModelType grM;
 };
 
 /// @brief Mesh adjacency (neighboring element for each element)

Code/Source/svFSI/Parameters.cpp

@@ -730,11 +730,58 @@ void StVenantKirchhoffParameters::print_parameters()
 /// @brief Parameters associated with GR_equilibrated material model
 GREquilibratedParameters::GREquilibratedParameters()
 {
-  value_set = true;
-}
-
-void GREquilibratedParameters::set_values(tinyxml2::XMLElement* con_params)
+  // A parameter that must be defined.
+  const bool required = true;
+
+  set_parameter("n_t_pre", 0, required, n_t_pre);
+  set_parameter("n_t_end", 0, required, n_t_end);
+  set_parameter("example", 0, required, example);
+  set_parameter("KsKi", 0.0, required, KsKi);
+  set_parameter("curve", 0.0, required, curve);
+  set_parameter("mult", 0.0, required, mult);
+  set_parameter("rIo", 0.0, required, rIo);
+  set_parameter("hwaves", 0.0, required, hwaves);
+  set_parameter("lo", 0.0, required, lo);
+  set_parameter("phieo", 0.0, required, phieo);
+  set_parameter("phimo", 0.0, required, phimo);
+  set_parameter("phico", 0.0, required, phico);
+  set_parameter("eta", 0.0, required, eta);
+  set_parameter("mu", 0.0, required, mu);
+  set_parameter("Get", 0.0, required, Get);
+  set_parameter("Gez", 0.0, required, Gez);
+  set_parameter("alpha", 0.0, required, alpha);
+  set_parameter("cm", 0.0, required, cm);
+  set_parameter("dm", 0.0, required, dm);
+  set_parameter("Gm", 0.0, required, Gm);
+  set_parameter("cc", 0.0, required, cc);
+  set_parameter("dc", 0.0, required, dc);
+  set_parameter("Gc", 0.0, required, Gc);
+  set_parameter("betat", 0.0, required, betat);
+  set_parameter("betaz", 0.0, required, betaz);
+  set_parameter("betad", 0.0, required, betad);
+  set_parameter("Tmax", 0.0, required, Tmax);
+  set_parameter("lamM", 0.0, required, lamM);
+  set_parameter("lam0", 0.0, required, lam0);
+  set_parameter("KfKi", 0.0, required, KfKi);
+  set_parameter("inflam", 0.0, required, inflam);
+  set_parameter("aexp", 0.0, required, aexp);
+  set_parameter("delta", 0.0, required, delta);
+
+  set_xml_element_name("Constitutive_model type=GR_equilibrated");
+}
+
+void GREquilibratedParameters::set_values(tinyxml2::XMLElement* xml_elem)
 {
+  
+  std::string error_msg = "Unknown Constitutive_model type=GR_equilibrated XML element '";
+
+  using std::placeholders::_1;
+  using std::placeholders::_2;
+  std::function<void(const std::string&, const std::string&)> ftpr =
+      std::bind( &LeeSacksParameters::set_parameter_value, *this, _1, _2);
+
+  xml_util_set_parameters(ftpr, xml_elem, error_msg);
+
   value_set = true;
 }
 

Code/Source/svFSI/Parameters.h

@@ -510,6 +510,40 @@ class GREquilibratedParameters : public ParameterLists
     void set_values(tinyxml2::XMLElement* modl_params);
     void print_parameters();
     bool value_set = false;
+    enum gr_example { none, aneurysm, tortuosity, stenosis };
+    Parameter<int> n_t_pre;
+    Parameter<int> n_t_end;
+    Parameter<int> example;
+    Parameter<double> KsKi;
+    Parameter<double> curve;
+    Parameter<double> mult;
+    Parameter<double> rIo;
+    Parameter<double> hwaves;
+    Parameter<double> lo;
+    Parameter<double> phieo;
+    Parameter<double> phimo;
+    Parameter<double> phico;
+    Parameter<double> eta;
+    Parameter<double> mu;
+    Parameter<double> Get;
+    Parameter<double> Gez;
+    Parameter<double> alpha;
+    Parameter<double> cm;
+    Parameter<double> dm;
+    Parameter<double> Gm;
+    Parameter<double> cc;
+    Parameter<double> dc;
+    Parameter<double> Gc;
+    Parameter<double> betat;
+    Parameter<double> betaz;
+    Parameter<double> betad;
+    Parameter<double> Tmax;
+    Parameter<double> lamM;
+    Parameter<double> lam0;
+    Parameter<double> KfKi;
+    Parameter<double> inflam;
+    Parameter<double> aexp;
+    Parameter<double> delta;
 };
 
 /// @brief The ConstitutiveModelParameters class store parameters

Code/Source/svFSI/distribute.cpp

@@ -1052,6 +1052,9 @@ void dist_eq(ComMod& com_mod, const CmMod& cm_mod, const cmType& cm, const std::
         (dmn.phys == EquationType::phys_gr)) {
       dist_mat_consts(com_mod, cm_mod, cm, dmn.stM);
     } 
+    if (dmn.phys == EquationType::phys_gr) {
+      dist_gr_consts(com_mod, cm_mod, cm, dmn.grM);
+    } 
 
     if ((dmn.phys == EquationType::phys_fluid) || 
         (dmn.phys == EquationType::phys_stokes) || 
@@ -1212,6 +1215,44 @@ void dist_mat_consts(const ComMod& com_mod, const CmMod& cm_mod, const cmType& c
 
 }
 
+/// @brief Distribute G&R properties to all processors.
+void dist_gr_consts(const ComMod& com_mod, const CmMod& cm_mod, const cmType& cm, grModelType& grM)
+{
+  using namespace consts;
+
+  cm.bcast(cm_mod, &grM.KsKi);
+  cm.bcast(cm_mod, &grM.curve);
+  cm.bcast(cm_mod, &grM.mult);
+  cm.bcast(cm_mod, &grM.rIo);
+  cm.bcast(cm_mod, &grM.hwaves);
+  cm.bcast(cm_mod, &grM.lo);
+  cm.bcast(cm_mod, &grM.phieo);
+  cm.bcast(cm_mod, &grM.phimo);
+  cm.bcast(cm_mod, &grM.phico);
+  cm.bcast(cm_mod, &grM.eta);
+  cm.bcast(cm_mod, &grM.mu);
+  cm.bcast(cm_mod, &grM.Get);
+  cm.bcast(cm_mod, &grM.Gez);
+  cm.bcast(cm_mod, &grM.alpha);
+  cm.bcast(cm_mod, &grM.cm);
+  cm.bcast(cm_mod, &grM.dm);
+  cm.bcast(cm_mod, &grM.Gm);
+  cm.bcast(cm_mod, &grM.cc);
+  cm.bcast(cm_mod, &grM.dc);
+  cm.bcast(cm_mod, &grM.Gc);
+  cm.bcast(cm_mod, &grM.betat);
+  cm.bcast(cm_mod, &grM.betaz);
+  cm.bcast(cm_mod, &grM.betad);
+  cm.bcast(cm_mod, &grM.Tmax);
+  cm.bcast(cm_mod, &grM.lamM);
+  cm.bcast(cm_mod, &grM.lam0);
+  cm.bcast(cm_mod, &grM.KfKi);
+  cm.bcast(cm_mod, &grM.inflam);
+  cm.bcast(cm_mod, &grM.aexp);
+  cm.bcast(cm_mod, &grM.delta);
+  cm.bcast(cm_mod, &grM.example);
+}
+
 
 void dist_visc_model(const ComMod& com_mod, const CmMod& cm_mod, const cmType& cm, viscModelType& lVis)
 {

Code/Source/svFSI/distribute.h

@@ -45,6 +45,8 @@ void dist_eq(ComMod& com_mod, const CmMod& cm_mod, const cmType& cm, const std::
 
 void dist_mat_consts(const ComMod& com_mod, const CmMod& cm_mod, const cmType& cm, stModelType& lStM);
 
+void dist_gr_consts(const ComMod& com_mod, const CmMod& cm_mod, const cmType& cm, grModelType& grM);
+
 void dist_visc_model(const ComMod& com_mod, const CmMod& cm_mod, const cmType& cm, viscModelType& lVis);
 
 void part_face(Simulation* simulation, mshType& lM, faceType& lFa, faceType& gFa, Vector<int>& gmtl);