ENH: added fields sources to simpleFoam

applications/solvers/incompressible/simpleFoam/Make/options

@@ -3,10 +3,12 @@ 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 = \
     -lincompressibleTurbulenceModel \
     -lincompressibleRASModels \
     -lincompressibleTransportModels \
-    -lfiniteVolume
+    -lfiniteVolume \
+    -lmeshTools

applications/solvers/incompressible/simpleFoam/UEqn.H

@@ -6,6 +6,8 @@
       + turbulence->divDevReff(U)
     );
 
+    sources.apply(UEqn());
+
     UEqn().relax();
 
     solve(UEqn() == -fvc::grad(p));

applications/solvers/incompressible/simpleFoam/createFields.H

@@ -39,3 +39,5 @@
     (
         incompressible::RASModel::New(U, phi, laminarTransport)
     );
+
+    IObasicSourceList sources(mesh);

applications/solvers/incompressible/simpleFoam/simpleFoam.C

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