From 6b3e8c972b01a29e5efd624f1b17ea7b3a5dae8e Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Mon, 30 Apr 2018 05:09:10 +0900
Subject: [PATCH] Qc and Qc2 field calculation dimension correctly set

---
 .../basicSpecieMixture/basicSpecieMixture.C      | 16 ++++++++++++----
 1 file changed, 12 insertions(+), 4 deletions(-)

diff --git a/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C b/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C
index 0ff307384..64af9f76c 100644
--- a/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C
+++ b/src/thermophysicalModels/reactionThermo/mixtures/basicSpecieMixture/basicSpecieMixture.C
@@ -96,15 +96,19 @@ Foam::tmp<Foam::volScalarField> Foam::basicSpecieMixture::Qc() const
                 Y[0].mesh()
             ),
             Y[0].mesh(),
-            dimensionedScalar("zero", dimless, 0)
+            dimensionedScalar("zero", dimCurrent*dimTime/dimMass, 0)
         )
     );
 
     volScalarField& Qc = tQc.ref();
 
+    dimensionedScalar Qci("Qci", dimCurrent*dimTime/dimMass, 0);
+
     forAll(Y, i)
     {
-        Qc += Y[i]*this->Qc(i);
+        Qci.value() = this->Qc(i);
+        Qc += Y[i]*Qci;
+        // Qc += Y[i]*this->Qc(i);
     }
 
     return tQc;
@@ -126,15 +130,19 @@ Foam::tmp<Foam::volScalarField> Foam::basicSpecieMixture::Qc2() const
                 Y[0].mesh()
             ),
             Y[0].mesh(),
-            dimensionedScalar("zero", dimless, 0)
+            dimensionedScalar("zero", sqr(dimCurrent*dimTime)/dimMass, 0)
         )
     );
 
     volScalarField& Qc = tQc.ref();
 
+    dimensionedScalar Qc2i("Qc2i", sqr(dimCurrent*dimTime)/dimMass, 0);
+
     forAll(Y, i)
     {
-        Qc += Y[i]*this->Qc2(i);
+        Qc2i.value() = this->Qc2(i);
+        Qc += Y[i]*Qc2i;
+        //Qc += Y[i]*this->Qc2(i);
     }
 
     return tQc;