update psi calculation in dfChemistryModel and dfCanteraMixture (#252)

* Update psi in CanteraMixture.H

replace W/RT with rho/p

* Update psi calculation dfChemistryModel.C

call psi function from dfCanteraMixture directly

* Update dfChemistryModel.C

* Update dfChemistryModel.C

* Update dfChemistryModel.C

* Update CanteraMixture.H

Author: wengzf20

src/dfCanteraMixture/CanteraMixture.H

@@ -117,7 +117,8 @@ public:
             const scalar& T
         ) const
         {
-            return CanteraGas_->meanMolecularWeight()/CanteraGas_->RT();
+            // for ideal gas, psi = CanteraGas_->density()/p = 1/RT
+            return CanteraGas_->density()/p;
         }
 
         scalar mu

src/dfChemistryModel/dfChemistryModel.C

@@ -383,8 +383,7 @@ void Foam::dfChemistryModel<ThermoType>::correctThermo()
 
         T_[celli] = CanteraGas_->temperature();
 
-        // meanMolecularWeight() kg/kmol    RT() Joules/kmol
-        psi_[celli] = CanteraGas_->meanMolecularWeight()/CanteraGas_->RT();
+        psi_[celli] = mixture_.psi(p_[celli],T_[celli]);
 
         mu_[celli] = mixture_.CanteraTransport()->viscosity(); // Pa-s
 
@@ -453,7 +452,7 @@ void Foam::dfChemistryModel<ThermoType>::correctThermo()
 
                 ph[facei] = CanteraGas_->enthalpy_mass();
 
-                ppsi[facei] = CanteraGas_->meanMolecularWeight()/CanteraGas_->RT();
+                ppsi[facei] = mixture_.psi(pp[facei],pT[facei]);
 
                 pmu[facei] = mixture_.CanteraTransport()->viscosity();
 
@@ -494,7 +493,7 @@ void Foam::dfChemistryModel<ThermoType>::correctThermo()
 
                 pT[facei] = CanteraGas_->temperature();
 
-                ppsi[facei] = CanteraGas_->meanMolecularWeight()/CanteraGas_->RT();
+                ppsi[facei] = mixture_.psi(pp[facei],pT[facei]);
 
                 pmu[facei] = mixture_.CanteraTransport()->viscosity();