From cf3a6a5cbe46c3143aff9ba0b1c76d3f876489ef Mon Sep 17 00:00:00 2001
From: ekinaci <efekinaci@yahoo.com>
Date: Fri, 24 Feb 2017 09:41:21 +0100
Subject: [PATCH] changed molecular weight

---
 applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H      | 3 +++
 .../solvers/cfdemSolverRhoPimpleChem/createFields.H       | 8 +++-----
 2 files changed, 6 insertions(+), 5 deletions(-)

diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H b/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
index df180fa7..9b643abd 100644
--- a/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
+++ b/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
@@ -50,8 +50,11 @@ tmp<fv::convectionScheme<scalar> > mvConvection
         }
     }
 
+    // calculate total molar number in volume
     N = Yt*rhoeps/W;
 
+    Info << "molecular weight of mixture" << W << endl;
+
     Y[inertIndex] = scalar(1) - Yt;
     Y[inertIndex].max(0.0);
 }
diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H b/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
index 89b7d252..e7a19636 100644
--- a/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
+++ b/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
@@ -18,9 +18,8 @@
     basicSpecieMixture& composition = thermo.composition();
     PtrList<volScalarField>& Y = composition.Y();
 
-    // The molecular weight [kj/mol]  ====== EKI
-    //const volScalarField W(composition.W());
-    const volScalarField& W = composition.W();
+    // The molecular weight [kj/mol]
+    const volScalarField W(composition.W());
 
     const word inertSpecie(thermo.lookup("inertSpecie"));
 
@@ -42,7 +41,7 @@
     }
     fields.add(thermo.he());
 
-    // define total number of moles in volume ====== EKI
+    // Total Mole in Volume
     volScalarField N
     (
         IOobject
@@ -55,7 +54,6 @@
         ),
         mesh
     );
-    // =============================================
 
     volScalarField rho
     (