[Input] Create ConstDensityThermo objects from YAML definitions

include/cantera/thermo/ConstDensityThermo.h

@@ -136,6 +136,7 @@ public:
         n = 1;
     }
 
+    virtual void initThermo();
     virtual void setParametersFromXML(const XML_Node& eosdata);
 
 protected:

src/thermo/ConstDensityThermo.cpp

@@ -105,6 +105,14 @@ void ConstDensityThermo::_updateThermo() const
     }
 }
 
+void ConstDensityThermo::initThermo()
+{
+    if (m_input.hasKey("density")) {
+        setDensity(m_input.convert("density", "kg/m^3"));
+    }
+    ThermoPhase::initThermo();
+}
+
 void ConstDensityThermo::setParametersFromXML(const XML_Node& eosdata)
 {
     eosdata._require("model","Incompressible");

src/thermo/ThermoFactory.cpp

@@ -52,6 +52,7 @@ ThermoFactory::ThermoFactory()
     reg("IdealGas", []() { return new IdealGasPhase(); });
     m_synonyms["ideal-gas"] = "IdealGas";
     reg("Incompressible", []() { return new ConstDensityThermo(); });
+    m_synonyms["constant-density"] = "Incompressible";
     reg("Surface", []() { return new SurfPhase(); });
     m_synonyms["ideal-surface"] = "Surface";
     reg("Edge", []() { return new EdgePhase(); });

test/data/thermo-models.yaml

@@ -179,6 +179,13 @@ phases:
   thermo: pure-fluid
   pure-fluid-name: nitrogen
 
+- name: const-density
+  thermo: constant-density
+  density: 0.7 g/cm^3
+  species: [NaCl(s), KCl(s)]
+  state: {T: 320, P: 1 atm}
+
+
 species:
 - name: NaCl(s)
   composition: {Na: 1, Cl: 1}

test/thermo/thermoFromYaml.cpp

@@ -332,3 +332,9 @@ TEST(ThermoFromYaml, PureFluid_nitrogen)
     EXPECT_NEAR(thermo->density(), 841.0420151, 1e-6);
     EXPECT_NEAR(thermo->gibbs_mole(), -17654452.8821914, 1e-6);
 }
+
+TEST(ThermoFromYaml, ConstDensityThermo)
+{
+    auto thermo = newThermo("thermo-models.yaml", "const-density");
+    EXPECT_DOUBLE_EQ(thermo->density(), 700.0);
+}