[unittests] Cover legacy StFlow class

Author: ischoegl

test/oneD/test_oneD.cpp

@@ -5,6 +5,7 @@
 #include "cantera/core.h"
 #include "cantera/onedim.h"
 #include "cantera/oneD/DomainFactory.h"
+#include "cantera/oneD/StFlow.h"
 #include "cantera/oneD/IonFlow.h"
 
 using namespace Cantera;
@@ -102,6 +103,97 @@ TEST(onedim, freeflame)
     }
 }
 
+TEST(onedim, legacy)
+{
+    auto sol = newSolution("h2o2.yaml", "ohmech", "mixture-averaged");
+    auto gas = sol->thermo();
+    size_t nsp = gas->nSpecies();
+
+    // reactants
+    double uin = .3;
+    double T = 300;
+    double P = 101325;
+    string X = "H2:0.65, O2:0.5, AR:2";
+    gas->setState_TPX(T, P, X);
+    double rho_in = gas->density();
+    vector<double> yin(nsp);
+    gas->getMassFractions(&yin[0]);
+
+    // product estimate
+    gas->equilibrate("HP");
+    vector<double> yout(nsp);
+    gas->getMassFractions(&yout[0]);
+    double rho_out = gas->density();
+    double Tad = gas->temperature();
+
+    // flow
+    suppress_deprecation_warnings();
+    auto flow = newDomain<StFlow>("legacy-flow", sol);
+    make_deprecation_warnings_fatal();
+    flow->setID("flow");
+    flow->setFreeFlow();
+
+    // grid
+    int nz = 21;
+    double lz = 0.02;
+    vector<double> z(nz);
+    double dz = lz;
+    dz /= (double)(nz - 1);
+    for (int iz = 0; iz < nz; iz++) {
+        z[iz] = iz * dz;
+    }
+    flow->setupGrid(nz, &z[0]);
+
+    // inlet
+    auto inlet = newDomain<Inlet1D>("inlet", sol);
+    inlet->setMoleFractions(X);
+    inlet->setMdot(uin * rho_in);
+    inlet->setTemperature(T);
+
+    // outlet
+    auto outlet = newDomain<Outlet1D>("outlet", sol);
+    double uout = inlet->mdot() / rho_out;
+
+    // set up simulation
+    vector<shared_ptr<Domain1D>> domains { inlet, flow, outlet };
+    Sim1D flame(domains);
+    int dom = static_cast<int>(flame.domainIndex("flow"));
+    ASSERT_EQ(dom, 1);
+
+    // set up initial guess
+    vector<double> locs{0.0, 0.3, 0.7, 1.0};
+    vector<double> value{uin, uin, uout, uout};
+    flame.setInitialGuess("velocity", locs, value);
+    value = {T, T, Tad, Tad};
+    flame.setInitialGuess("T", locs, value);
+    for (size_t i = 0; i < nsp; i++) {
+        value = {yin[i], yin[i], yout[i], yout[i]};
+        flame.setInitialGuess(gas->speciesName(i), locs, value);
+    }
+
+    // simulation settings
+    double ratio = 15.0;
+    double slope = 0.3;
+    double curve = 0.5;
+    flame.setRefineCriteria(dom, ratio, slope, curve);
+    flame.setFixedTemperature(0.85 * T + .15 * Tad);
+
+    // solve
+    flow->solveEnergyEqn();
+    bool refine_grid = false;
+    int loglevel = 0;
+    flame.solve(loglevel, refine_grid);
+
+    ASSERT_EQ(flow->nPoints(), static_cast<size_t>(nz + 1));
+    size_t comp = flow->componentIndex("T");
+    double Tprev = flame.value(dom, comp, 0);
+    for (size_t n = 0; n < flow->nPoints(); n++) {
+        T = flame.value(dom, comp, n);
+        ASSERT_GE(T, Tprev);
+        Tprev = T;
+    }
+}
+
 TEST(onedim, flame_types)
 {
     auto sol = newSolution("h2o2.yaml", "ohmech", "mixture-averaged");