phi calculation due to turbulence.phi() returning volume flux

applications/solvers/combustion/postSootFoam/Make/options

@@ -30,6 +30,7 @@ EXE_LIBS = \
     -lchemistryModel \
     -lcombustionModels \
     -lradiationModels \
+    -lradiationSootModels \
     -llagrangian \
     -llagrangianIntermediate \
     -llagrangianTurbulence \

src/thermophysicalModels/radiationSootModel/MossBrookesSoot/MossBrookesSoot.C

@@ -351,16 +351,27 @@ void Foam::radiation::MossBrookesSoot::correct()
 
     calcSource();
 
-    const volScalarField rho_ = thermo_.rho();
+    const volScalarField &rho_ = thermo_.rho();
 
-    const volScalarField &nut_ = turbulence().nut();
+    const volVectorField &U_ = turbulence().U();
 
-    const surfaceScalarField phi_ (turbulence().phi());
+    volScalarField Dsoot (turbulence().mut()/PrtSoot);
+
+    const surfaceScalarField phi_
+    (
+        IOobject
+        (
+            "phiSoot",
+            mesh().time().timeName(),
+            mesh()
+        ),
+        linearInterpolate(rho_*U_) & mesh().Sf()
+    );
 
     fvScalarMatrix SootEqn
     (
         fvm::div(phi_, soot_)
-      - fvm::laplacian((rho_*nut_/PrtSoot), soot_)
+      - fvm::laplacian(Dsoot, soot_)
      ==
         Snet_
     );
@@ -372,7 +383,7 @@ void Foam::radiation::MossBrookesSoot::correct()
     fvScalarMatrix NucConcEqn
     (
         fvm::div(phi_, NucConc_)
-      - fvm::laplacian((rho_*nut_/PrtSoot), NucConc_)
+      - fvm::laplacian(Dsoot, NucConc_)
      ==
         NSnet_
     );
@@ -393,9 +404,9 @@ void Foam::radiation::MossBrookesSoot::correct()
 
 void Foam::radiation::MossBrookesSoot::calcSource()
 {
-    const volScalarField rho_ = thermo_.rho();
+    const volScalarField &rho_ = thermo_.rho();
-    const volScalarField T_ = thermo_.T();
+    const volScalarField &T_ = thermo_.T();
-    const volScalarField p_ = thermo_.p();
+    const volScalarField &p_ = thermo_.p();
 
     const scalar precI = mixture_.species()[coeffsDict_.lookup("Precursor")];
     const scalar growI = mixture_.species()[coeffsDict_.lookup("SurfaceGrowSpecie")];

src/thermophysicalModels/radiationSootModel/MossBrookesSoot/MossBrookesSoot.H

@@ -153,15 +153,15 @@ class MossBrookesSoot
         //- Soot particle mean diameter
         volScalarField sootMeanDiameter_;
 
+        //- soot number density normalization factor
+        const dimensionedScalar Nnorm;
+
         //- Soot turbulent Prandtl number
         scalar PrtSoot;
 
         //- Density of soot particle
         scalar rhoSoot;
 
-        //- soot number density normalization factor
-        const dimensionedScalar Nnorm;
-
         //- Model constant for soot inception rate
         scalar Calpha;
 

src/thermophysicalModels/radiationSootModel/khanGreeveSoot/khanGreeveSoot.C

@@ -271,15 +271,24 @@ void Foam::radiation::khanGreeveSoot::correct()
 
     const volScalarField &rho_ = thermo_.rho();
 
-    const volScalarField &nut_ = turbulence().nut();
+    const volVectorField &U_ = turbulence().U();
 
-    const surfaceScalarField phi_ (turbulence().phi());
+    const surfaceScalarField phi_
+    (
+        IOobject
+        (
+            "phiSoot",
+            mesh().time().timeName(),
+            mesh()
+        ),
+        linearInterpolate(rho_*U_) & mesh().Sf()
+    );
 
     fvScalarMatrix SootEqn
     (
         fvm::ddt(rho_, soot_)
       + fvm::div(phi_, soot_)
-      - fvm::laplacian((rho_*nut_/PrtSoot), soot_)
+      - fvm::laplacian((turbulence().mut()/PrtSoot), soot_)
      ==
         Snet_
     );