BUG: Multiply phase averaged he by rhos

jheylmun · jheylmun · commit d8e89513b8ef · 2026-02-24T10:38:40.000-06:00
diff --git a/src/thermodynamicModels/derived/multiphaseFluidBlastThermo/multiphaseFluidBlastThermo.C b/src/thermodynamicModels/derived/multiphaseFluidBlastThermo/multiphaseFluidBlastThermo.C
@@ -765,10 +765,18 @@ Foam::multiphaseFluidBlastThermo::he
     const volScalarField& T
 ) const
 {
-    tmp<volScalarField> tmpF(volumeFractions_[0]*thermos_[0].he(p, T));
+    tmp<volScalarField> tmpF
+    (
+        volumeFractions_[0]
+       *rhos_[0]
+       *thermos_[0].he(p, T)
+    );
     for (label phasei = 1; phasei < thermos_.size(); phasei++)
     {
-        tmpF.ref() += volumeFractions_[phasei]*thermos_[phasei].he(p, T);
+        tmpF.ref() +=
+            volumeFractions_[phasei]
+           *rhos_[phasei]
+           *thermos_[phasei].he(p, T);
     }
     return massNormalise(tmpF);
 }

Original file line number	Diff line number	Diff line change
`@@ -765,10 +765,18 @@ Foam::multiphaseFluidBlastThermo::he`
`765`	`765`	`const volScalarField& T`
`766`	`766`	`) const`
`767`	`767`	`{`
`768`		`- tmp<volScalarField> tmpF(volumeFractions_[0]*thermos_[0].he(p, T));`
	`768`	`+ tmp<volScalarField> tmpF`
	`769`	`+ (`
	`770`	`+ volumeFractions_[0]`
	`771`	`+ *rhos_[0]`
	`772`	`+ *thermos_[0].he(p, T)`
	`773`	`+ );`
`769`	`774`	`for (label phasei = 1; phasei < thermos_.size(); phasei++)`
`770`	`775`	`{`
`771`		`- tmpF.ref() += volumeFractions_[phasei]*thermos_[phasei].he(p, T);`
	`776`	`+ tmpF.ref() +=`
	`777`	`+ volumeFractions_[phasei]`
	`778`	`+ *rhos_[phasei]`
	`779`	`+ *thermos_[phasei].he(p, T);`
`772`	`780`	`}`
`773`	`781`	`return massNormalise(tmpF);`
`774`	`782`	`}`