turn on radaition model

Author: chenzhenyi-123

Allwclean

@@ -21,6 +21,7 @@ wclean ./applications/solvers/dfSprayFoam
 wclean ./applications/solvers/df0DFoam
 wclean ./applications/solvers/dfLowMachFoam
 wclean ./applications/solvers/dfHighSpeedFoam
+wclean ./applications/solvers/dfBuoyancyFoam
 
 rm -rf src_orig/
 rm -rf bin/

applications/solvers/dfBuoyancyFoam/Make/options

@@ -32,6 +32,7 @@ EXE_INC = -std=c++14 \
     -I$(DF_SRC)/dfCanteraMixture/lnInclude \
     -I$(DF_SRC)/dfChemistryModel/lnInclude \
     -I$(DF_SRC)/dfCombustionModels/lnInclude \
+    -I$(DF_SRC)/radiationModels/lnInclude \
     -I$(CANTERA_ROOT)/include \
     $(if $(LIBTORCH_ROOT),-I$(LIBTORCH_ROOT)/include,) \
     $(if $(LIBTORCH_ROOT),-I$(LIBTORCH_ROOT)/include/torch/csrc/api/include,) \
@@ -62,6 +63,7 @@ EXE_LIBS = \
     -ldfCanteraMixture \
     -ldfChemistryModel \
     -ldfCombustionModels \
+    -ldfRadiationModels \
     $(CANTERA_ROOT)/lib/libcantera.so \
     $(if $(LIBTORCH_ROOT),$(LIBTORCH_ROOT)/lib/libtorch.so,) \
     $(if $(LIBTORCH_ROOT),$(LIBTORCH_ROOT)/lib/libc10.so,) \

applications/solvers/dfBuoyancyFoam/YEEqn.H

@@ -9,7 +9,7 @@ tmp<fv::convectionScheme<scalar>> mvConvection
     )
 );
 {
-    // radiation->correct();
+    radiation->correct();
     combustion->correct();
 
     HRR = combustion->Qdot();
@@ -68,8 +68,8 @@ tmp<fv::convectionScheme<scalar>> mvConvection
       - fvm::laplacian(turbulence->alphaEff(), he)
      ==
         // combustion->Qdot() // ha Eqn has no reaction source
-    //   + radiation->Sh(thermo, he)
-         parcels.Sh(he)
+        radiation->Sh(thermo, he)
+      + parcels.Sh(he)
     //   + surfaceFilm.Sh()
     //   + fvOptions(rho, he)
     );

applications/solvers/dfBuoyancyFoam/createFields.H

@@ -164,5 +164,5 @@ volScalarField K("K", 0.5*magSqr(U));
 #include "createClouds.H"
 #include "createSurfaceFilmModel.H"
 // #include "createPyrolysisModel.H"
-// #include "createRadiationModel.H"
+#include "createRadiationModel.H"
 // #include "createFvOptions.H"

applications/solvers/dfBuoyancyFoam/dfBuoyancyFoam.C

@@ -42,7 +42,7 @@ Description
 #include "basicReactingCloud.H"
 #include "surfaceFilmModel.H"
 #include "pyrolysisModelCollection.H"
-// #include "radiationModel.H"
+#include "radiationModel.H"
 #include "SLGThermo.H"
 // #include "solidChemistryModel.H"
 // #include "psiReactionThermo.H"