From e9ab0be29d61e82fd4b42ff5ebf64ccfbd6d793b Mon Sep 17 00:00:00 2001
From: Zhi Chen <chenzhi09@me.com>
Date: Sun, 6 Nov 2022 17:24:08 +0800
Subject: [PATCH] Fix bug in correctThermo on boundary patches

diffusivity coefficient missing rho
---
 src/dfChemistryModel/dfChemistryModel.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/dfChemistryModel/dfChemistryModel.C b/src/dfChemistryModel/dfChemistryModel.C
index 332c06a3..575ac542 100644
--- a/src/dfChemistryModel/dfChemistryModel.C
+++ b/src/dfChemistryModel/dfChemistryModel.C
@@ -672,7 +672,7 @@ void Foam::dfChemistryModel<ThermoType>::correctThermo()
                     const scalar RT = constant::physicoChemical::R.value()*1e3*pT[facei];
                     forAll(rhoD_, i)
                     {
-                        rhoD_[i].boundaryFieldRef()[patchi][facei] = dTemp_[i];
+                        rhoD_[i].boundaryFieldRef()[patchi][facei] = prho[facei]*dTemp_[i];
 
                         hai_[i].boundaryFieldRef()[patchi][facei] = hrtTemp_[i]*RT/CanteraGas_->molecularWeight(i);
                     }