From 4354360884cae6cafb655e05f722e3bb0e7964cb Mon Sep 17 00:00:00 2001
From: Yeongdo Park <ignis@postech.ac.kr>
Date: Wed, 2 Jan 2019 12:37:43 -0500
Subject: [PATCH] fix mixture density to pure composition density

---
 diffusivityModel/Particle/ParticleI.H | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/diffusivityModel/Particle/ParticleI.H b/diffusivityModel/Particle/ParticleI.H
index c8c137c..76b69f8 100644
--- a/diffusivityModel/Particle/ParticleI.H
+++ b/diffusivityModel/Particle/ParticleI.H
@@ -146,7 +146,7 @@ inline Foam::scalar Foam::Particle::fTrans(const scalar rSc)
 inline Foam::scalar Foam::Particle::lambda(const scalar mu, const scalar D, const scalar Zrot, const GasState &state)
     const
 {
-    const scalar fVib = state.rho() * D / mu;
+    const scalar fVib = thermo_.rho(state.p(),state.T()) * D / mu;
 
     const scalar A = 5./2. - fVib;
     const scalar B = Zrot + (2./pi) * ((5./3.)*CvRot()/R() + fVib);