New methods for all species Cp and Cv

src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.C

@@ -98,6 +98,40 @@ Foam::scalar Foam::SpecieMixture<MixtureType>::Qc2
 }
 
 
+template<class MixtureType>
+void Foam::SpecieMixture<MixtureType>::Cp
+(
+    scalarField& Cps,
+    const scalar p,
+    const scalar T
+) const
+{
+    forAll(this->species(), speciei)
+    {
+        Cps[speciei] = this->getLocalThermo(speciei).Cp(p, T);
+    }
+
+    return;
+}
+
+
+template<class MixtureType>
+void Foam::SpecieMixture<MixtureType>::Cv
+(
+    scalarField& Cvs,
+    const scalar p,
+    const scalar T
+) const
+{
+    forAll(this->species(), speciei)
+    {
+        Cvs[speciei] = this->getLocalThermo(speciei).Cv(p, T);
+    }
+
+    return;
+}
+
+
 template<class MixtureType>
 Foam::scalar Foam::SpecieMixture<MixtureType>::Cp
 (

src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.H

@@ -91,6 +91,25 @@ public:
             virtual scalar Qc2(const label specieI) const;
 
 
+        // All species thermo properties
+
+            //- Heat capacity at constant pressure [J/(kg K)]
+            virtual void Cp
+            (
+                scalarField& Cps,
+                const scalar p,
+                const scalar T
+            ) const;
+
+            //- Heat capacity at constant volume [J/(kg K)]
+            virtual void Cv
+            (
+                scalarField& Cvs,
+                const scalar p,
+                const scalar T
+            ) const;
+
+
         // Per specie thermo properties
 
             //- Heat capacity at constant pressure [J/(kg K)]

src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.H

@@ -107,6 +107,25 @@ public:
             tmp<volScalarField> Qc2() const;
 
 
+        // All species thermo properties
+
+            //- Heat capacity at constant pressure [J/(kg K)]
+            virtual void Cp
+            (
+                scalarField& Cps,
+                const scalar p,
+                const scalar T
+            ) const = 0;
+
+            //- Heat capacity at constant volume [J/(kg K)]
+            virtual void Cv
+            (
+                scalarField& Cvs,
+                const scalar p,
+                const scalar T
+            ) const = 0;
+
+
         // Per specie thermo properties
 
             //- Heat capacity at constant pressure [J/(kg K)]