fluxRequired: Added setFluxRequired function to fvSchemes class

Added calls to setFluxRequired for p in all incompressible solvers which
avoids the need to add fluxRequired entries in fvSchemes dictionary.

Will add calls to setFluxRequired to the rest of the solvers.

applications/solvers/basic/potentialFoam/createFields.H

@@ -1,98 +1,99 @@
-    Info<< "Reading velocity field U\n" << endl;
+Info<< "Reading velocity field U\n" << endl;
-    volVectorField U
+volVectorField U
+(
+    IOobject
     (
-        IOobject
+        "U",
-        (
+        runTime.timeName(),
-            "U",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
-
-    U = dimensionedVector("0", U.dimensions(), vector::zero);
-
-    surfaceScalarField phi
-    (
-        IOobject
-        (
-            "phi",
-            runTime.timeName(),
-            mesh,
-            IOobject::NO_READ,
-            IOobject::AUTO_WRITE
-        ),
-        fvc::interpolate(U) & mesh.Sf()
-    );
-
-    if (args.optionFound("initialiseUBCs"))
-    {
-        U.correctBoundaryConditions();
-        phi = fvc::interpolate(U) & mesh.Sf();
-    }
-
-
-    // Default name for the pressure field
-    word pName("p");
-
-    // Update name of the pressure field from the command-line option
-    args.optionReadIfPresent("pName", pName);
-
-    // Infer the pressure BCs from the velocity BCs
-    wordList pBCTypes
-    (
-        U.boundaryField().size(),
-        fixedValueFvPatchScalarField::typeName
-    );
-
-    forAll(U.boundaryField(), patchi)
-    {
-        if (U.boundaryField()[patchi].fixesValue())
-        {
-            pBCTypes[patchi] = zeroGradientFvPatchScalarField::typeName;
-        }
-    }
-
-    Info<< "Constructing pressure field " << pName << nl << endl;
-    volScalarField p
-    (
-        IOobject
-        (
-            pName,
-            runTime.timeName(),
-            mesh,
-            IOobject::READ_IF_PRESENT,
-            IOobject::NO_WRITE
-        ),
         mesh,
-        dimensionedScalar(pName, sqr(dimVelocity), 0),
+        IOobject::MUST_READ,
-        pBCTypes
+        IOobject::AUTO_WRITE
-    );
+    ),
+    mesh
+);
 
-    Info<< "Constructing velocity potential field Phi\n" << endl;
+U = dimensionedVector("0", U.dimensions(), vector::zero);
-    volScalarField Phi
+
+surfaceScalarField phi
+(
+    IOobject
     (
-        IOobject
+        "phi",
-        (
+        runTime.timeName(),
-            "Phi",
-            runTime.timeName(),
-            mesh,
-            IOobject::READ_IF_PRESENT,
-            IOobject::NO_WRITE
-        ),
         mesh,
-        dimensionedScalar("Phi", dimLength*dimVelocity, 0),
+        IOobject::NO_READ,
-        p.boundaryField().types()
+        IOobject::AUTO_WRITE
-    );
+    ),
+    fvc::interpolate(U) & mesh.Sf()
+);
 
-    label PhiRefCell = 0;
+if (args.optionFound("initialiseUBCs"))
-    scalar PhiRefValue = 0;
+{
-    setRefCell
+    U.correctBoundaryConditions();
+    phi = fvc::interpolate(U) & mesh.Sf();
+}
+
+
+// Default name for the pressure field
+word pName("p");
+
+// Update name of the pressure field from the command-line option
+args.optionReadIfPresent("pName", pName);
+
+// Infer the pressure BCs from the velocity BCs
+wordList pBCTypes
+(
+    U.boundaryField().size(),
+    fixedValueFvPatchScalarField::typeName
+);
+
+forAll(U.boundaryField(), patchi)
+{
+    if (U.boundaryField()[patchi].fixesValue())
+    {
+        pBCTypes[patchi] = zeroGradientFvPatchScalarField::typeName;
+    }
+}
+
+Info<< "Constructing pressure field " << pName << nl << endl;
+volScalarField p
+(
+    IOobject
     (
-        Phi,
+        pName,
-        potentialFlow.dict(),
+        runTime.timeName(),
-        PhiRefCell,
+        mesh,
-        PhiRefValue
+        IOobject::READ_IF_PRESENT,
-    );
+        IOobject::NO_WRITE
+    ),
+    mesh,
+    dimensionedScalar(pName, sqr(dimVelocity), 0),
+    pBCTypes
+);
+
+Info<< "Constructing velocity potential field Phi\n" << endl;
+volScalarField Phi
+(
+    IOobject
+    (
+        "Phi",
+        runTime.timeName(),
+        mesh,
+        IOobject::READ_IF_PRESENT,
+        IOobject::NO_WRITE
+    ),
+    mesh,
+    dimensionedScalar("Phi", dimLength*dimVelocity, 0),
+    p.boundaryField().types()
+);
+
+label PhiRefCell = 0;
+scalar PhiRefValue = 0;
+setRefCell
+(
+    Phi,
+    potentialFlow.dict(),
+    PhiRefCell,
+    PhiRefValue
+);
+mesh.setFluxRequired(Phi.name());

applications/solvers/incompressible/adjointShapeOptimizationFoam/createFields.H

@@ -1,109 +1,111 @@
-    Info<< "Reading field p\n" << endl;
+Info<< "Reading field p\n" << endl;
-    volScalarField p
+volScalarField p
+(
+    IOobject
     (
-        IOobject
+        "p",
-        (
+        runTime.timeName(),
-            "p",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    Info<< "Reading field U\n" << endl;
+Info<< "Reading field U\n" << endl;
-    volVectorField U
+volVectorField U
+(
+    IOobject
     (
-        IOobject
+        "U",
-        (
+        runTime.timeName(),
-            "U",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    #include "createPhi.H"
+#include "createPhi.H"
 
 
-    label pRefCell = 0;
+label pRefCell = 0;
-    scalar pRefValue = 0.0;
+scalar pRefValue = 0.0;
-    setRefCell(p, simple.dict(), pRefCell, pRefValue);
+setRefCell(p, simple.dict(), pRefCell, pRefValue);
+mesh.setFluxRequired(p.name());
 
 
-    Info<< "Reading field pa\n" << endl;
+Info<< "Reading field pa\n" << endl;
-    volScalarField pa
+volScalarField pa
+(
+    IOobject
     (
-        IOobject
+        "pa",
-        (
+        runTime.timeName(),
-            "pa",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    Info<< "Reading field Ua\n" << endl;
+Info<< "Reading field Ua\n" << endl;
-    volVectorField Ua
+volVectorField Ua
+(
+    IOobject
     (
-        IOobject
+        "Ua",
-        (
+        runTime.timeName(),
-            "Ua",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    #include "createPhia.H"
+#include "createPhia.H"
 
 
-    label paRefCell = 0;
+label paRefCell = 0;
-    scalar paRefValue = 0.0;
+scalar paRefValue = 0.0;
-    setRefCell
+setRefCell
+(
+    pa,
+    simple.dict(),
+    paRefCell,
+    paRefValue
+);
+mesh.setFluxRequired(pa.name());
+
+
+singlePhaseTransportModel laminarTransport(U, phi);
+
+autoPtr<incompressible::RASModel> turbulence
+(
+    incompressible::RASModel::New(U, phi, laminarTransport)
+);
+
+
+dimensionedScalar zeroSensitivity("0", dimVelocity*dimVelocity, 0.0);
+dimensionedScalar zeroAlpha("0", dimless/dimTime, 0.0);
+
+dimensionedScalar lambda(laminarTransport.lookup("lambda"));
+dimensionedScalar alphaMax(laminarTransport.lookup("alphaMax"));
+
+const labelList& inletCells = mesh.boundary()["inlet"].faceCells();
+//const labelList& outletCells = mesh.boundary()["outlet"].faceCells();
+
+volScalarField alpha
+(
+    IOobject
     (
-        pa,
+        "alpha",
-        simple.dict(),
+        runTime.timeName(),
-        paRefCell,
+        mesh,
-        paRefValue
+        IOobject::READ_IF_PRESENT,
-    );
+        IOobject::AUTO_WRITE
-
+    ),
-
+    lambda*max(Ua & U, zeroSensitivity)
-    singlePhaseTransportModel laminarTransport(U, phi);
+);
-
+zeroCells(alpha, inletCells);
-    autoPtr<incompressible::RASModel> turbulence
+//zeroCells(alpha, outletCells);
-    (
-        incompressible::RASModel::New(U, phi, laminarTransport)
-    );
-
-
-    dimensionedScalar zeroSensitivity("0", dimVelocity*dimVelocity, 0.0);
-    dimensionedScalar zeroAlpha("0", dimless/dimTime, 0.0);
-
-    dimensionedScalar lambda(laminarTransport.lookup("lambda"));
-    dimensionedScalar alphaMax(laminarTransport.lookup("alphaMax"));
-
-    const labelList& inletCells = mesh.boundary()["inlet"].faceCells();
-    //const labelList& outletCells = mesh.boundary()["outlet"].faceCells();
-
-    volScalarField alpha
-    (
-        IOobject
-        (
-            "alpha",
-            runTime.timeName(),
-            mesh,
-            IOobject::READ_IF_PRESENT,
-            IOobject::AUTO_WRITE
-        ),
-        lambda*max(Ua & U, zeroSensitivity)
-    );
-    zeroCells(alpha, inletCells);
-    //zeroCells(alpha, outletCells);

applications/solvers/incompressible/icoFoam/createFields.H

@@ -1,55 +1,56 @@
-    Info<< "Reading transportProperties\n" << endl;
+Info<< "Reading transportProperties\n" << endl;
 
-    IOdictionary transportProperties
+IOdictionary transportProperties
+(
+    IOobject
     (
-        IOobject
+        "transportProperties",
-        (
+        runTime.constant(),
-            "transportProperties",
+        mesh,
-            runTime.constant(),
+        IOobject::MUST_READ_IF_MODIFIED,
-            mesh,
+        IOobject::NO_WRITE
-            IOobject::MUST_READ_IF_MODIFIED,
+    )
-            IOobject::NO_WRITE
+);
-        )
-    );
 
-    dimensionedScalar nu
+dimensionedScalar nu
+(
+    transportProperties.lookup("nu")
+);
+
+Info<< "Reading field p\n" << endl;
+volScalarField p
+(
+    IOobject
     (
-        transportProperties.lookup("nu")
+        "p",
-    );
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);
 
-    Info<< "Reading field p\n" << endl;
+
-    volScalarField p
+Info<< "Reading field U\n" << endl;
+volVectorField U
+(
+    IOobject
     (
-        IOobject
+        "U",
-        (
+        runTime.timeName(),
-            "p",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
 
-    Info<< "Reading field U\n" << endl;
+#include "createPhi.H"
-    volVectorField U
-    (
-        IOobject
-        (
-            "U",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
 
 
-    #include "createPhi.H"
+label pRefCell = 0;
-
+scalar pRefValue = 0.0;
-
+setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
-    label pRefCell = 0;
+mesh.setFluxRequired(p.name());
-    scalar pRefValue = 0.0;
-    setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);

applications/solvers/incompressible/nonNewtonianIcoFoam/createFields.H

@@ -1,39 +1,40 @@
-    Info<< "Reading field p\n" << endl;
+Info<< "Reading field p\n" << endl;
-    volScalarField p
+volScalarField p
+(
+    IOobject
     (
-        IOobject
+        "p",
-        (
+        runTime.timeName(),
-            "p",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
 
-    Info<< "Reading field U\n" << endl;
+Info<< "Reading field U\n" << endl;
-    volVectorField U
+volVectorField U
+(
+    IOobject
     (
-        IOobject
+        "U",
-        (
+        runTime.timeName(),
-            "U",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
 
-    #include "createPhi.H"
+#include "createPhi.H"
 
 
-    singlePhaseTransportModel fluid(U, phi);
+singlePhaseTransportModel fluid(U, phi);
 
 
-    label pRefCell = 0;
+label pRefCell = 0;
-    scalar pRefValue = 0.0;
+scalar pRefValue = 0.0;
-    setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
+setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
+mesh.setFluxRequired(p.name());

applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/createFields.H

@@ -43,6 +43,7 @@ surfaceScalarField phi
 label pRefCell = 0;
 scalar pRefValue = 0.0;
 setRefCell(p, pimple.dict(), pRefCell, pRefValue);
+mesh.setFluxRequired(p.name());
 
 singlePhaseTransportModel laminarTransport(Urel, phi);
 

applications/solvers/incompressible/pimpleFoam/createFields.H

@@ -32,6 +32,7 @@ volVectorField U
 label pRefCell = 0;
 scalar pRefValue = 0.0;
 setRefCell(p, pimple.dict(), pRefCell, pRefValue);
+mesh.setFluxRequired(p.name());
 
 
 singlePhaseTransportModel laminarTransport(U, phi);

applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/createFields.H

@@ -0,0 +1,2 @@
+#include "../createFields.H"
+mesh.setFluxRequired("pcorr");

applications/solvers/incompressible/pisoFoam/createFields.H

@@ -1,42 +1,43 @@
-    Info<< "Reading field p\n" << endl;
+Info<< "Reading field p\n" << endl;
-    volScalarField p
+volScalarField p
+(
+    IOobject
     (
-        IOobject
+        "p",
-        (
+        runTime.timeName(),
-            "p",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    Info<< "Reading field U\n" << endl;
+Info<< "Reading field U\n" << endl;
-    volVectorField U
+volVectorField U
+(
+    IOobject
     (
-        IOobject
+        "U",
-        (
+        runTime.timeName(),
-            "U",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    #include "createPhi.H"
+#include "createPhi.H"
 
 
-    label pRefCell = 0;
+label pRefCell = 0;
-    scalar pRefValue = 0.0;
+scalar pRefValue = 0.0;
-    setRefCell(p, piso.dict(), pRefCell, pRefValue);
+setRefCell(p, piso.dict(), pRefCell, pRefValue);
+mesh.setFluxRequired(p.name());
 
 
-    singlePhaseTransportModel laminarTransport(U, phi);
+singlePhaseTransportModel laminarTransport(U, phi);
 
-    autoPtr<incompressible::turbulenceModel> turbulence
+autoPtr<incompressible::turbulenceModel> turbulence
-    (
+(
-        incompressible::turbulenceModel::New(U, phi, laminarTransport)
+    incompressible::turbulenceModel::New(U, phi, laminarTransport)
-    );
+);

applications/solvers/incompressible/shallowWaterFoam/createFields.H

@@ -73,3 +73,5 @@ hTotal.write();
 
 Info<< "Creating Coriolis Force" << endl;
 const dimensionedVector F("F", ((2.0*Omega) & gHat)*gHat);
+
+mesh.setFluxRequired(h.name());

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

@@ -1,55 +1,56 @@
-    Info<< "Reading field p\n" << endl;
+Info<< "Reading field p\n" << endl;
-    volScalarField p
+volScalarField p
+(
+    IOobject
     (
-        IOobject
+        "p",
-        (
+        runTime.timeName(),
-            "p",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    Info<< "Reading field Urel\n" << endl;
+Info<< "Reading field Urel\n" << endl;
-    volVectorField Urel
+volVectorField Urel
+(
+    IOobject
     (
-        IOobject
+        "Urel",
-        (
+        runTime.timeName(),
-            "Urel",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    Info<< "Reading/calculating face flux field phi\n" << endl;
+Info<< "Reading/calculating face flux field phi\n" << endl;
-    surfaceScalarField phi
+surfaceScalarField phi
+(
+    IOobject
     (
-        IOobject
+        "phi",
-        (
+        runTime.timeName(),
-            "phi",
+        mesh,
-            runTime.timeName(),
+        IOobject::READ_IF_PRESENT,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::READ_IF_PRESENT,
+    ),
-            IOobject::AUTO_WRITE
+    linearInterpolate(Urel) & mesh.Sf()
-        ),
+);
-        linearInterpolate(Urel) & mesh.Sf()
-    );
 
-    label pRefCell = 0;
+label pRefCell = 0;
-    scalar pRefValue = 0.0;
+scalar pRefValue = 0.0;
-    setRefCell(p, simple.dict(), pRefCell, pRefValue);
+setRefCell(p, simple.dict(), pRefCell, pRefValue);
+mesh.setFluxRequired(p.name());
 
-    singlePhaseTransportModel laminarTransport(Urel, phi);
+singlePhaseTransportModel laminarTransport(Urel, phi);
 
-    autoPtr<incompressible::RASModel> turbulence
+autoPtr<incompressible::RASModel> turbulence
-    (
+(
-        incompressible::RASModel::New(Urel, phi, laminarTransport)
+    incompressible::RASModel::New(Urel, phi, laminarTransport)
-    );
+);
 
-    Info<< "Creating SRF model\n" << endl;
+Info<< "Creating SRF model\n" << endl;
-    autoPtr<SRF::SRFModel> SRF(SRF::SRFModel::New(Urel));
+autoPtr<SRF::SRFModel> SRF(SRF::SRFModel::New(Urel));

applications/solvers/incompressible/simpleFoam/createFields.H

@@ -1,42 +1,43 @@
-    Info<< "Reading field p\n" << endl;
+Info<< "Reading field p\n" << endl;
-    volScalarField p
+volScalarField p
+(
+    IOobject
     (
-        IOobject
+        "p",
-        (
+        runTime.timeName(),
-            "p",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    Info<< "Reading field U\n" << endl;
+Info<< "Reading field U\n" << endl;
-    volVectorField U
+volVectorField U
+(
+    IOobject
     (
-        IOobject
+        "U",
-        (
+        runTime.timeName(),
-            "U",
+        mesh,
-            runTime.timeName(),
+        IOobject::MUST_READ,
-            mesh,
+        IOobject::AUTO_WRITE
-            IOobject::MUST_READ,
+    ),
-            IOobject::AUTO_WRITE
+    mesh
-        ),
+);
-        mesh
-    );
 
-    #include "createPhi.H"
+#include "createPhi.H"
 
 
-    label pRefCell = 0;
+label pRefCell = 0;
-    scalar pRefValue = 0.0;
+scalar pRefValue = 0.0;
-    setRefCell(p, simple.dict(), pRefCell, pRefValue);
+setRefCell(p, simple.dict(), pRefCell, pRefValue);
+mesh.setFluxRequired(p.name());
 
 
-    singlePhaseTransportModel laminarTransport(U, phi);
+singlePhaseTransportModel laminarTransport(U, phi);
 
-    autoPtr<incompressible::turbulenceModel> turbulence
+autoPtr<incompressible::turbulenceModel> turbulence
-    (
+(
-        incompressible::turbulenceModel::New(U, phi, laminarTransport)
+    incompressible::turbulenceModel::New(U, phi, laminarTransport)
-    );
+);

src/finiteVolume/finiteVolume/fvSchemes/fvSchemes.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -42,13 +42,13 @@ void Foam::fvSchemes::clear()
     defaultD2dt2Scheme_.clear();
     interpolationSchemes_.clear();
     defaultInterpolationScheme_.clear();
-    divSchemes_.clear(); // optional
+    divSchemes_.clear();
     defaultDivScheme_.clear();
-    gradSchemes_.clear(); // optional
+    gradSchemes_.clear();
     defaultGradScheme_.clear();
     snGradSchemes_.clear();
     defaultSnGradScheme_.clear();
-    laplacianSchemes_.clear(); // optional
+    laplacianSchemes_.clear();
     defaultLaplacianScheme_.clear();
     fluxRequired_.clear();
     defaultFluxRequired_ = false;
@@ -558,6 +558,17 @@ Foam::ITstream& Foam::fvSchemes::laplacianScheme(const word& name) const
 }
 
 
+void Foam::fvSchemes::setFluxRequired(const word& name)
+{
+    if (debug)
+    {
+        Info<< "Setting fluxRequired for " << name << endl;
+    }
+
+    fluxRequired_.add(name, true);
+}
+
+
 bool Foam::fvSchemes::fluxRequired(const word& name) const
 {
     if (debug)

src/finiteVolume/finiteVolume/fvSchemes/fvSchemes.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -130,6 +130,8 @@ public:
 
             ITstream& laplacianScheme(const word& name) const;
 
+            void setFluxRequired(const word& name);
+
             bool fluxRequired(const word& name) const;
 
             //- Return true if the default ddtScheme is steadyState

tutorials/basic/potentialFoam/cylinder/system/fvSchemes

@@ -47,11 +47,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    Phi             ;
-}
-
 
 // ************************************************************************* //

tutorials/basic/potentialFoam/pitzDaily/system/fvSchemes

@@ -47,11 +47,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    Phi             ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/SRFPimpleFoam/rotor2D/system/fvSchemes

@@ -49,11 +49,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/SRFSimpleFoam/mixer/system/fvSchemes

@@ -54,12 +54,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/adjointShapeOptimizationFoam/pitzDaily/system/fvSchemes

@@ -53,11 +53,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p;
-    pa;
-}
 
 // ************************************************************************* //

tutorials/incompressible/icoFoam/cavity/system/fvSchemes

@@ -47,11 +47,5 @@ snGradSchemes
     default         orthogonal;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/icoFoam/cavityClipped/system/fvSchemes

@@ -46,11 +46,5 @@ snGradSchemes
     default         orthogonal;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/icoFoam/cavityGrade/system/fvSchemes

@@ -46,11 +46,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/icoFoam/elbow/system/fvSchemes

@@ -46,11 +46,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/nonNewtonianIcoFoam/offsetCylinder/system/fvSchemes

@@ -47,11 +47,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleDyMFoam/mixerVesselAMI2D/system/fvSchemes

@@ -50,12 +50,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    pcorr           ;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleDyMFoam/movingCone/system/fvSchemes

@@ -50,12 +50,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    pcorr           ;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleDyMFoam/oscillatingInletACMI2D/system/fvSchemes

@@ -52,12 +52,5 @@ snGradSchemes
     default         limited corrected 0.33;
 }
 
-fluxRequired
-{
-    default         no;
-    pcorr           ;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleDyMFoam/propeller/system/fvSchemes

@@ -52,12 +52,5 @@ snGradSchemes
     default         limited corrected 0.33;
 }
 
-fluxRequired
-{
-    default         no;
-    pcorr           ;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam/system/fvSchemes

@@ -50,13 +50,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    pcorr           ;
-    p;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_simpleFoam/system/fvSchemes

@@ -50,12 +50,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/pimpleDyMFoam/wingMotion/wingMotion_snappyHexMesh/system/fvSchemes

@@ -54,11 +54,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleFoam/TJunction/system/fvSchemes

@@ -52,11 +52,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleFoam/TJunctionFan/system/fvSchemes

@@ -52,11 +52,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pimpleFoam/channel395/system/fvSchemes

@@ -51,12 +51,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/pimpleFoam/elipsekkLOmega/system/fvSchemes

@@ -50,12 +50,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/pimpleFoam/pitzDaily/system/fvSchemes

@@ -61,11 +61,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pisoFoam/les/motorBike/lesFiles/fvSchemes

@@ -57,12 +57,6 @@ snGradSchemes
     default         limited corrected 0.33;
 }
 
-fluxRequired
-{
-    default         no;
-    p;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/pisoFoam/les/motorBike/motorBike/system/fvSchemes

@@ -52,13 +52,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p;
-    Phi;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/pisoFoam/les/pitzDaily/system/fvSchemes

@@ -48,11 +48,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pisoFoam/les/pitzDailyMapped/system/fvSchemes

@@ -51,11 +51,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pisoFoam/ras/cavity/system/fvSchemes

@@ -53,11 +53,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/pisoFoam/ras/cavityCoupledU/system/fvSchemes

@@ -52,11 +52,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/porousSimpleFoam/angledDuctExplicit/system/fvSchemes

@@ -50,11 +50,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/porousSimpleFoam/angledDuctImplicit/system/fvSchemes

@@ -50,11 +50,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/porousSimpleFoam/straightDuctImplicit/system/fvSchemes

@@ -50,11 +50,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/shallowWaterFoam/squareBump/system/fvSchemes

@@ -45,10 +45,5 @@ snGradSchemes
     default         uncorrected;
 }
 
-fluxRequired
-{
-    h;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/simpleFoam/airFoil2D/system/fvSchemes

@@ -48,12 +48,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/simpleFoam/mixerVessel2D/system/fvSchemes

@@ -49,11 +49,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/simpleFoam/motorBike/system/fvSchemes

@@ -49,13 +49,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p;
-    Phi;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/simpleFoam/pipeCyclic/system/fvSchemes

@@ -52,11 +52,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/simpleFoam/pitzDaily/system/fvSchemes

@@ -52,12 +52,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/simpleFoam/pitzDailyExptInlet/system/fvSchemes

@@ -52,11 +52,5 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/simpleFoam/rotorDisk/system/fvSchemes

@@ -58,12 +58,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 wallDist
 {
     method meshWave;

tutorials/incompressible/simpleFoam/turbineSiting/system/fvSchemes

@@ -50,11 +50,5 @@ snGradSchemes
     default             limited corrected 0.333;
 }
 
-fluxRequired
-{
-    default             no;
-    p;
-}
-
 
 // ************************************************************************* //

tutorials/incompressible/simpleFoam/windAroundBuildings/system/fvSchemes

@@ -57,12 +57,6 @@ snGradSchemes
     default         corrected;
 }
 
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
 wallDist
 {
     method meshWave;