[MATLAB] Added massDensity property to ThermoPhase

Modified several samples to include this change

Author: ssun30

interfaces/matlab_experimental/Base/ThermoPhase.m

@@ -66,6 +66,8 @@
         %     Scalar double mean molecular weight. Units: kg/kmol.
         meanMolecularWeight
 
+        massDensity % Mass basis density. Units: kg/m^3.
+
         molarDensity % Molar basis density. Units: kmol/m^3.
 
         molecularWeights % Molecular weights of the species. Units: kg/kmol.
@@ -903,6 +905,10 @@ function display(tp)
             mmw = ctFunc('thermo_meanMolecularWeight', tp.tpID);
         end
 
+        function density = get.massDensity(tp)
+            density = ctFunc('thermo_density', tp.tpID);
+        end
+
         function density = get.molarDensity(tp)
             density = ctFunc('thermo_molarDensity', tp.tpID);
         end

samples/matlab_experimental/conhp.m

@@ -6,14 +6,14 @@
     % It assumes that the ``gas`` object represents a reacting ideal gas mixture.
 
     % Set the state of the gas, based on the current solution vector.
+    gas.basis = 'mass';
     gas.Y = y(2:end);
     gas.TP = {y(1), gas.P};
     nsp = gas.nSpecies;
 
     % energy equation
     wdot = gas.netProdRates;
     H = gas.partialMolarEnthalpies';
-    gas.basis = 'mass';
     tdot =- 1 / (gas.D * gas.cp) .* wdot * H;
 
     % set up column vector for dydt

samples/matlab_experimental/conuv.m

@@ -6,14 +6,14 @@
     % It assumes that the ``gas`` object represents a reacting ideal gas mixture.
 
     % Set the state of the gas, based on the current solution vector.
+    gas.basis = 'mass';
     gas.Y = y(2:end);
     gas.TD = {y(1), gas.D};
     nsp = gas.nSpecies;
 
     % energy equation
     wdot = gas.netProdRates;
     U = gas.partialMolarIntEnergies';
-    gas.basis = 'mass';
     tdot =- 1 / (gas.D * gas.cv) .* wdot * U;
 
     % set up column vector for dydt

samples/matlab_experimental/diamond_cvd.m

@@ -69,7 +69,7 @@
     r = surf_phase.netProdRates;
     carbon_dot = r(iC);
     mdot = mw * carbon_dot;
-    rate = mdot / dbulk.D;
+    rate = mdot / dbulk.massDensity;
     xx = [xx; x(ih)];
     rr = [rr; rate * 1.0e6 * 3600.0];
     cov = [cov; surf_phase.coverages];

samples/matlab_experimental/flame.m

@@ -32,7 +32,7 @@
     f = Sim1D({left flow right});
 
     % set default initial profiles.
-    rho0 = gas.D;
+    rho0 = gas.massDensity;
 
     % find the adiabatic flame temperature and corresponding
     % equilibrium composition

samples/matlab_experimental/ignite.m

@@ -25,7 +25,7 @@
     gas.TPX = {1001.0, OneAtm, 'H2:2,O2:1,N2:4'};
     gas.basis = 'mass';
     y0 = [gas.U
-          1.0 / gas.D
+          1.0 / gas.massDensity
           gas.Y'];
 
     time_interval = [0 0.001];

samples/matlab_experimental/periodic_cstr.m

@@ -72,7 +72,7 @@
     % controller (constant mdot). Set the mass flow rate to 1.25 sccm.
     sccm = 1.25;
     vdot = sccm * 1.0e-6/60.0 * ((OneAtm / gas.P) * (gas.T / 273.15)); % m^3/s
-    mdot = gas.D * vdot; % kg/s
+    mdot = gas.massDensity * vdot; % kg/s
     mfc = MassFlowController(upstream, cstr);
     mfc.massFlowRate = mdot;
 

samples/matlab_experimental/plug_flow_reactor.m

@@ -99,7 +99,7 @@
 
 T_calc(1) = gas_calc.T;
 Y_calc(1, :) = gas_calc.Y;
-rho_calc(1) = gas_calc.D;
+rho_calc(1) = gas_calc.massDensity;
 
 for i = 2:length(x_calc)