ENH: Introduced filed sources to remaining ico SIMPLE solvers

applications/solvers/incompressible/simpleFoam/MRFSimpleFoam/MRFSimpleFoam.C

@@ -35,6 +35,7 @@ Description
 #include "RASModel.H"
 #include "MRFZones.H"
 #include "simpleControl.H"
+#include "IObasicSourceList.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 

applications/solvers/incompressible/simpleFoam/MRFSimpleFoam/Make/options

@@ -4,9 +4,11 @@ EXE_INC = \
     -I$(LIB_SRC)/turbulenceModels/incompressible/RAS/RASModel \
     -I$(LIB_SRC)/transportModels \
     -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \
-    -I$(LIB_SRC)/finiteVolume/lnInclude
+    -I$(LIB_SRC)/finiteVolume/lnInclude \
+    -I$(LIB_SRC)/meshTools/lnInclude
 
 EXE_LIBS = \
     -lincompressibleRASModels \
     -lincompressibleTransportModels \
-    -lfiniteVolume
+    -lfiniteVolume \
+    -lmeshTools

applications/solvers/incompressible/simpleFoam/MRFSimpleFoam/UEqn.H

@@ -3,10 +3,13 @@
     tmp<fvVectorMatrix> UEqn
     (
         fvm::div(phi, U)
-        + turbulence->divDevReff(U)
+      + turbulence->divDevReff(U)
     );
+
     mrfZones.addCoriolis(UEqn());
 
+    sources.apply(UEqn());
+
     UEqn().relax();
 
     solve(UEqn() == -fvc::grad(p));

applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/Make/options

@@ -4,7 +4,8 @@ EXE_INC = \
     -I$(LIB_SRC)/turbulenceModels \
     -I$(LIB_SRC)/turbulenceModels/incompressible/RAS/RASModel \
     -I$(LIB_SRC)/transportModels \
-    -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel
+    -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \
+    -I$(LIB_SRC)/meshTools/lnInclude
 
 EXE_LIBS = \
     -lincompressibleRASModels \

applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C

@@ -35,6 +35,7 @@ Description
 #include "RASModel.H"
 #include "SRFModel.H"
 #include "simpleControl.H"
+#include "IObasicSourceList.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 

applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/UrelEqn.H

@@ -7,6 +7,8 @@
       + SRF->Su()
     );
 
+    sources.apply(UrelEqn());
+
     UrelEqn().relax();
 
     solve(UrelEqn() == -fvc::grad(p));

applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/createFields.H

@@ -56,3 +56,5 @@
     (
         SRF::SRFModel::New(Urel)
     );
+
+    IObasicSourceList sources(mesh);

applications/solvers/incompressible/simpleFoam/porousSimpleFoam/UEqn.H

@@ -6,6 +6,8 @@
       + turbulence->divDevReff(U)
     );
 
+    sources.apply(UEqn());
+
     UEqn().relax();
 
     // Include the porous media resistance and solve the momentum equation

applications/solvers/incompressible/simpleFoam/porousSimpleFoam/porousSimpleFoam.C

@@ -35,6 +35,7 @@ Description
 #include "RASModel.H"
 #include "porousZones.H"
 #include "simpleControl.H"
+#include "IObasicSourceList.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //