ENH: correcting new thermo type for reactingParcelFoam tutorials.

Correcting thermoSingleLayer.C mask field alpha to avoid heat sources where there is no film.
Tunning fvSolution for alpha for twoPhasePachuka tutorial

src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C

@@ -3,7 +3,7 @@
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
     \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
-     \\/     M anipulation  |
+     \\/     M anipulation  | Copyright (C) 2017 OpenCFD Ltd
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -269,17 +269,15 @@ void thermoSingleLayer::updateSubmodels()
 
 tmp<fvScalarMatrix> thermoSingleLayer::q(volScalarField& hs) const
 {
-    const volScalarField alpha(pos(delta_ - deltaSmall_));
-
     return
     (
         // Heat-transfer to the primary region
       - fvm::Sp(htcs_->h()/Cp_, hs)
-      + htcs_->h()*(hs/Cp_ + alpha*(TPrimary_ - T_))
+      + htcs_->h()*(hs/Cp_ + alpha_*(TPrimary_ - T_))
 
         // Heat-transfer to the wall
       - fvm::Sp(htcw_->h()/Cp_, hs)
-      + htcw_->h()*(hs/Cp_ + alpha*(Tw_- T_))
+      + htcw_->h()*(hs/Cp_ + alpha_*(Tw_- T_))
     );
 }
 

tutorials/lagrangian/reactingParcelFoam/hotBoxes/Allrun-parallel

@@ -6,8 +6,9 @@ cd ${0%/*} || exit 1    # Run from this directory
 
 ./Allrun.pre
 
-runApplication -s wallFilmRegion decomposePar -region wallFilmRegion
 runApplication -s primaryRegion  decomposePar
+runApplication -s wallFilmRegion decomposePar -region wallFilmRegion
+
 
 runParallel $(getApplication)
 

tutorials/lagrangian/reactingParcelFoam/hotBoxes/constant/chemistryProperties

@@ -18,7 +18,7 @@ FoamFile
 chemistryType
 {
     chemistrySolver   noChemistrySolver;
-    chemistryThermo   psi;
+    chemistryThermo   rho;
 }
 
 chemistry       off;

tutorials/lagrangian/reactingParcelFoam/hotBoxes/constant/combustionProperties

@@ -15,7 +15,7 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-combustionModel  laminar<psiChemistryCombustion>;
+combustionModel  laminar<rhoChemistryCombustion>;
 
 active  false;
 

tutorials/lagrangian/reactingParcelFoam/hotBoxes/constant/thermophysicalProperties

@@ -17,7 +17,7 @@ FoamFile
 
 thermoType
 {
-    type            hePsiThermo;
+    type            heRhoThermo;
     mixture         reactingMixture;
     transport       sutherland;
     thermo          janaf;

tutorials/lagrangian/reactingParcelFoam/rivuletPanel/constant/chemistryProperties

@@ -18,7 +18,7 @@ FoamFile
 chemistryType
 {
     chemistrySolver   noChemistrySolver;
-    chemistryThermo   psi;
+    chemistryThermo   rho;
 }
 
 chemistry       off;

tutorials/lagrangian/reactingParcelFoam/rivuletPanel/constant/combustionProperties

@@ -15,7 +15,7 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-combustionModel  laminar<psiChemistryCombustion>;
+combustionModel  laminar<rhoChemistryCombustion>;
 
 active  false;
 

tutorials/lagrangian/reactingParcelFoam/rivuletPanel/constant/thermophysicalProperties

@@ -17,7 +17,7 @@ FoamFile
 
 thermoType
 {
-    type            hePsiThermo;
+    type            heRhoThermo;
     mixture         reactingMixture;
     transport       sutherland;
     thermo          janaf;

tutorials/lagrangian/reactingParcelFoam/splashPanel/constant/chemistryProperties

@@ -18,7 +18,7 @@ FoamFile
 chemistryType
 {
     chemistrySolver   noChemistrySolver;
-    chemistryThermo   psi;
+    chemistryThermo   rho;
 }
 
 chemistry       off;

tutorials/lagrangian/reactingParcelFoam/splashPanel/constant/combustionProperties

@@ -15,7 +15,7 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-combustionModel  laminar<psiChemistryCombustion>;
+combustionModel  laminar<rhoChemistryCombustion>;
 
 active  false;
 

tutorials/lagrangian/reactingParcelFoam/splashPanel/constant/surfaceFilmProperties

@@ -33,7 +33,7 @@ thermoSingleLayerCoeffs
 
     filmViscosityModel liquid;
 
-    deltaWet    1e-4;
+    deltaWet    2e-4;
     hydrophilic no;
 
     turbulence  laminar;

tutorials/lagrangian/reactingParcelFoam/splashPanel/constant/thermophysicalProperties

@@ -17,7 +17,7 @@ FoamFile
 
 thermoType
 {
-    type            hePsiThermo;
+    type            heRhoThermo;
     mixture         reactingMixture;
     transport       sutherland;
     thermo          janaf;

tutorials/multiphase/MPPICInterFoam/twoPhasePachuka/system/fvSolution

@@ -24,7 +24,7 @@ solvers
         cAlpha          1;
 
         MULESCorr       yes;
-        nLimiterIter    2;
+        nLimiterIter    10;
 
         solver          smoothSolver;
         smoother        symGaussSeidel;
@@ -45,7 +45,7 @@ solvers
     {
         solver          PCG;
         preconditioner  DIC;
-        tolerance       1e-07;
+        tolerance       1e-08;
         relTol          0.05;
     }