BUG: LESDelta: parallel inconsistent

Internal values get set but never swapped on processor boundaries. Fixes #49.

src/TurbulenceModels/schemes/DEShybrid/DEShybrid.H

@@ -207,10 +207,10 @@ class DEShybrid
                 max(sigmaMax_*tanh(pow(A, CH1_)), sigmaMin_)
             );
 
-            // if (this->mesh().time().outputTime())
-            // {
-            //     factor.write();
-            // }
+            if (blendedSchemeBaseName::debug)
+            {
+                factor.write();
+            }
 
             return tmp<surfaceScalarField>
             (

src/TurbulenceModels/turbulenceModels/LES/LESdeltas/IDDESDelta/IDDESDelta.C

@@ -3,7 +3,7 @@
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
     \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
-     \\/     M anipulation  |
+     \\/     M anipulation  | Copyright (C) 2016 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -127,6 +127,9 @@ void Foam::LESModels::IDDESDelta::calcDelta()
             ),
             hmax
         );
+
+    // Handle coupled boundaries
+    delta_.correctBoundaryConditions();
 }
 
 

src/TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C

@@ -3,7 +3,7 @@
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
     \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
-     \\/     M anipulation  |
+     \\/     M anipulation  | Copyright (C) 2016 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -76,6 +76,9 @@ void Foam::LESModels::cubeRootVolDelta::calcDelta()
             << "Case is not 3D or 2D, LES is not applicable"
             << exit(FatalError);
     }
+
+    // Handle coupled boundaries
+    delta_.correctBoundaryConditions();
 }
 
 

src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C

@@ -3,7 +3,7 @@
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
     \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
-     \\/     M anipulation  |
+     \\/     M anipulation  | Copyright (C) 2016 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -102,6 +102,9 @@ void Foam::LESModels::maxDeltaxyz::calcDelta()
             << "Case is not 3D or 2D, LES is not applicable"
             << exit(FatalError);
     }
+
+    // Handle coupled boundaries
+    delta_.correctBoundaryConditions();
 }
 
 

src/TurbulenceModels/turbulenceModels/LES/LESdeltas/smoothDelta/smoothDelta.C

@@ -3,7 +3,7 @@
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
     \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
-     \\/     M anipulation  |
+     \\/     M anipulation  | Copyright (C) 2016 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -134,6 +134,9 @@ void Foam::LESModels::smoothDelta::calcDelta()
     {
         delta_[cellI] = cellDeltaData[cellI].delta();
     }
+
+    // Handle coupled boundaries
+    delta_.correctBoundaryConditions();
 }