Incompressible turbulence models: Remove the correction of the laminar transport model

Explicitly correct laminar transport at the application level as is done in the multiphase solvers
2025-12-08 06:57:46 +00:00 · 2015-01-12 12:32:38 +00:00
parent 5a4a4e5e4a
commit b5277967c0
108 changed files with 133 additions and 359 deletions
--- a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/buoyantBoussinesqPimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/buoyantBoussinesqPimpleFoam.C
@ -97,6 +97,7 @@ int main(int argc, char *argv[])

            if (pimple.turbCorr())
            {
+                laminarTransport.correct();
                turbulence->correct();
            }
        }
--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
@ -81,6 +81,7 @@ int main(int argc, char *argv[])
            #include "pEqn.H"
        }

+        laminarTransport.correct();
        turbulence->correct();

        runTime.write();
--- a/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C
+++ b/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C
@ -216,6 +216,7 @@ int main(int argc, char *argv[])
            Ua.correctBoundaryConditions();
        }

+        laminarTransport.correct();
        turbulence->correct();

        runTime.write();
--- a/applications/solvers/incompressible/boundaryFoam/boundaryFoam.C
+++ b/applications/solvers/incompressible/boundaryFoam/boundaryFoam.C
@ -79,7 +79,7 @@ int main(int argc, char *argv[])
        U += (Ubar - UbarStar);
        gradP += (Ubar - UbarStar)/(1.0/UEqn.A())().weightedAverage(mesh.V());

-
+        laminarTransport.correct();
        turbulence->correct();

        Info<< "Uncorrected Ubar = " << (flowDirection & UbarStar.value())
--- a/applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/SRFPimpleFoam.C
+++ b/applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/SRFPimpleFoam.C
@ -82,6 +82,7 @@ int main(int argc, char *argv[])

            if (pimple.turbCorr())
            {
+                laminarTransport.correct();
                turbulence->correct();
            }
        }
--- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pimpleDyMFoam.C
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pimpleDyMFoam.C
@ -104,6 +104,7 @@ int main(int argc, char *argv[])

            if (pimple.turbCorr())
            {
+                laminarTransport.correct();
                turbulence->correct();
            }
        }
--- a/applications/solvers/incompressible/pimpleFoam/pimpleFoam.C
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleFoam.C
@ -83,6 +83,7 @@ int main(int argc, char *argv[])

            if (pimple.turbCorr())
            {
+                laminarTransport.correct();
                turbulence->correct();
            }
        }
--- a/applications/solvers/incompressible/pisoFoam/createFields.H
+++ b/applications/solvers/incompressible/pisoFoam/createFields.H
@ -26,7 +26,7 @@
        mesh
    );

-#   include "createPhi.H"
+    #include "createPhi.H"


    label pRefCell = 0;
--- a/applications/solvers/incompressible/pisoFoam/pisoFoam.C
+++ b/applications/solvers/incompressible/pisoFoam/pisoFoam.C
@ -130,6 +130,7 @@ int main(int argc, char *argv[])
            }
        }

+        laminarTransport.correct();
        turbulence->correct();

        runTime.write();
--- a/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C
+++ b/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C
@ -64,6 +64,7 @@ int main(int argc, char *argv[])
            #include "pEqn.H"
        }

+        laminarTransport.correct();
        turbulence->correct();

        if (runTime.outputTime())
--- a/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/porousSimpleFoam.C
+++ b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/porousSimpleFoam.C
@ -68,6 +68,7 @@ int main(int argc, char *argv[])
            #include "pEqn.H"
        }

+        laminarTransport.correct();
        turbulence->correct();

        runTime.write();
--- a/applications/solvers/incompressible/simpleFoam/simpleFoam.C
+++ b/applications/solvers/incompressible/simpleFoam/simpleFoam.C
@ -62,6 +62,7 @@ int main(int argc, char *argv[])
            #include "pEqn.H"
        }

+        laminarTransport.correct();
        turbulence->correct();

        runTime.write();
--- a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C
@ -134,6 +134,7 @@ int main(int argc, char *argv[])

            if (pimple.turbCorr())
            {
+                laminarTransport.correct();
                turbulence->correct();
            }
        }
--- a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C
@ -83,6 +83,7 @@ int main(int argc, char *argv[])

            if (pimple.turbCorr())
            {
+                laminarTransport.correct();
                turbulence->correct();
            }
        }
--- a/applications/test/PisoFoam/Make/options
+++ b/applications/test/PisoFoam/Make/options
@ -1,4 +1,5 @@
 EXE_INC = \
+    -I$(FOAM_SOLVERS)/incompressible/pisoFoam \
    -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/incompressible/lnInclude \
    -I$(LIB_SRC)/transportModels \
--- a/applications/test/PisoFoam/PisoFoam.C
+++ b/applications/test/PisoFoam/PisoFoam.C
@ -33,7 +33,7 @@ Description

 #include "fvCFD.H"
 #include "singlePhaseTransportModel.H"
-#include "IncompressibleTurbulenceModel.H"
+#include "turbulenceModel.H"

 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

--- a/applications/test/PisoFoam/createFields.H
+++ b/applications/test/PisoFoam/createFields.H
@ -1,48 +0,0 @@
-    Info<< "Reading field p\n" << endl;
-    volScalarField p
-    (
-        IOobject
-        (
-            "p",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
-
-    Info<< "Reading field U\n" << endl;
-    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);
-
-
-    singlePhaseTransportModel laminarTransport(U, phi);
-
-    autoPtr<IncompressibleTurbulenceModel<transportModel> >
-    turbulence
-    (
-        IncompressibleTurbulenceModel<transportModel>::New
-        (
-            U,
-            phi,
-            laminarTransport
-        )
-    );
--- a/applications/test/PisoFoam/turbulenceModel.H
+++ b/applications/test/PisoFoam/turbulenceModel.H
@ -0,0 +1,9 @@
+#include "IncompressibleTurbulenceModel.H"
+
+namespace Foam
+{
+    namespace incompressible
+    {
+        typedef IncompressibleTurbulenceModel<transportModel> turbulenceModel;
+    }
+}
--- a/applications/test/RhoPimpleFoam/Make/options
+++ b/applications/test/RhoPimpleFoam/Make/options
@ -1,4 +1,5 @@
 EXE_INC = \
+    -I$(FOAM_SOLVERS)/compressible/rhoPimpleFoam \
    -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
--- a/applications/test/RhoPimpleFoam/UEqn.H
+++ b/applications/test/RhoPimpleFoam/UEqn.H
@ -1,22 +0,0 @@
-// Solve the Momentum equation
-
-tmp<fvVectorMatrix> UEqn
-(
-    fvm::ddt(rho, U)
-  + fvm::div(phi, U)
-  + turbulence->divDevRhoReff(U)
- ==
-    fvOptions(rho, U)
-);
-
-UEqn().relax();
-
-fvOptions.constrain(UEqn());
-
-if (pimple.momentumPredictor())
-{
-    solve(UEqn() == -fvc::grad(p));
-
-    fvOptions.correct(U);
-    K = 0.5*magSqr(U);
-}
--- a/applications/test/RhoPimpleFoam/createFields.H
+++ b/applications/test/RhoPimpleFoam/createFields.H
@ -1,73 +0,0 @@
-    Info<< "Reading thermophysical properties\n" << endl;
-
-    autoPtr<psiThermo> pThermo
-    (
-        psiThermo::New(mesh)
-    );
-    psiThermo& thermo = pThermo();
-    thermo.validate(args.executable(), "h", "e");
-
-    volScalarField& p = thermo.p();
-    const volScalarField& psi = thermo.psi();
-
-    volScalarField rho
-    (
-        IOobject
-        (
-            "rho",
-            runTime.timeName(),
-            mesh,
-            IOobject::READ_IF_PRESENT,
-            IOobject::AUTO_WRITE
-        ),
-        thermo.rho()
-    );
-
-    Info<< "Reading field U\n" << endl;
-    volVectorField U
-    (
-        IOobject
-        (
-            "U",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
-
-    #include "compressibleCreatePhi.H"
-
-    dimensionedScalar rhoMax(pimple.dict().lookup("rhoMax"));
-    dimensionedScalar rhoMin(pimple.dict().lookup("rhoMin"));
-
-    Info<< "Creating turbulence model\n" << endl;
-
-    autoPtr<CompressibleTurbulenceModel<fluidThermo> >
-    turbulence
-    (
-        CompressibleTurbulenceModel<fluidThermo>::New
-        (
-            rho,
-            U,
-            phi,
-            thermo
-        )
-    );
-
-    Info<< "Creating field dpdt\n" << endl;
-    volScalarField dpdt
-    (
-        IOobject
-        (
-            "dpdt",
-            runTime.timeName(),
-            mesh
-        ),
-        mesh,
-        dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
-    );
-
-    Info<< "Creating field kinetic energy K\n" << endl;
-    volScalarField K("K", 0.5*magSqr(U));
--- a/applications/test/RhoPimpleFoam/pEqn.H
+++ b/applications/test/RhoPimpleFoam/pEqn.H
@ -1,113 +0,0 @@
-rho = thermo.rho();
-rho = max(rho, rhoMin);
-rho = min(rho, rhoMax);
-rho.relax();
-
-volScalarField rAU(1.0/UEqn().A());
-volVectorField HbyA("HbyA", U);
-HbyA = rAU*UEqn().H();
-
-if (pimple.nCorrPISO() <= 1)
-{
-    UEqn.clear();
-}
-
-if (pimple.transonic())
-{
-    surfaceScalarField phid
-    (
-        "phid",
-        fvc::interpolate(psi)
-       *(
-            (fvc::interpolate(HbyA) & mesh.Sf())
-          + fvc::ddtCorr(rAU, rho, U, phi)
-        )
-    );
-
-    fvOptions.makeRelative(fvc::interpolate(psi), phid);
-
-    volScalarField Dp("Dp", rho*rAU);
-
-    while (pimple.correctNonOrthogonal())
-    {
-        fvScalarMatrix pEqn
-        (
-            fvm::ddt(psi, p)
-          + fvm::div(phid, p)
-          - fvm::laplacian(Dp, p)
-          ==
-            fvOptions(psi, p, rho.name())
-        );
-
-        fvOptions.constrain(pEqn);
-
-        pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
-
-        if (pimple.finalNonOrthogonalIter())
-        {
-            phi == pEqn.flux();
-        }
-    }
-}
-else
-{
-    surfaceScalarField phiHbyA
-    (
-        "phiHbyA",
-        fvc::interpolate(rho)
-       *(
-            (fvc::interpolate(HbyA) & mesh.Sf())
-          + fvc::ddtCorr(rAU, rho, U, phi)
-        )
-    );
-
-    fvOptions.makeRelative(fvc::interpolate(rho), phiHbyA);
-
-    volScalarField Dp("Dp", rho*rAU);
-
-    while (pimple.correctNonOrthogonal())
-    {
-        // Pressure corrector
-        fvScalarMatrix pEqn
-        (
-            fvm::ddt(psi, p)
-          + fvc::div(phiHbyA)
-          - fvm::laplacian(Dp, p)
-          ==
-            fvOptions(psi, p, rho.name())
-        );
-
-        fvOptions.constrain(pEqn);
-
-        pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
-
-        if (pimple.finalNonOrthogonalIter())
-        {
-            phi = phiHbyA + pEqn.flux();
-        }
-    }
-}
-
-#include "rhoEqn.H"
-#include "compressibleContinuityErrs.H"
-
-// Explicitly relax pressure for momentum corrector
-p.relax();
-
-// Recalculate density from the relaxed pressure
-rho = thermo.rho();
-rho = max(rho, rhoMin);
-rho = min(rho, rhoMax);
-rho.relax();
-Info<< "rho max/min : " << max(rho).value()
-    << " " << min(rho).value() << endl;
-
-U = HbyA - rAU*fvc::grad(p);
-U.correctBoundaryConditions();
-fvOptions.correct(U);
-K = 0.5*magSqr(U);
-
-if (thermo.dpdt())
-{
-    dpdt = fvc::ddt(p);
-}
--- a/applications/test/RhoPimpleFoam/rhoPimpleFoam.C
+++ b/applications/test/RhoPimpleFoam/rhoPimpleFoam.C
@ -35,7 +35,7 @@ Description

 #include "fvCFD.H"
 #include "psiThermo.H"
-#include "CompressibleTurbulenceModel.H"
+#include "turbulenceModel.H"
 #include "bound.H"
 #include "pimpleControl.H"
 #include "fvIOoptionList.H"
--- a/applications/test/RhoPimpleFoam/turbulenceModel.H
+++ b/applications/test/RhoPimpleFoam/turbulenceModel.H
@ -0,0 +1,9 @@
+#include "CompressibleTurbulenceModel.H"
+
+namespace Foam
+{
+    namespace compressible
+    {
+        typedef CompressibleTurbulenceModel<fluidThermo> turbulenceModel;
+    }
+}
--- a/src/turbulenceModels/incompressible/LES/DESModel/DESModel.C
+++ b/src/turbulenceModels/incompressible/LES/DESModel/DESModel.C
@ -39,7 +39,7 @@ DESModel::DESModel
    const word& type,
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName
 )
 :
--- a/src/turbulenceModels/incompressible/LES/DESModel/DESModel.H
+++ b/src/turbulenceModels/incompressible/LES/DESModel/DESModel.H
@ -84,7 +84,7 @@ public:
            const word& type,
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName
        );

--- a/src/turbulenceModels/incompressible/LES/DeardorffDiffStress/DeardorffDiffStress.C
+++ b/src/turbulenceModels/incompressible/LES/DeardorffDiffStress/DeardorffDiffStress.C
@ -55,7 +55,7 @@ DeardorffDiffStress::DeardorffDiffStress
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/DeardorffDiffStress/DeardorffDiffStress.H
+++ b/src/turbulenceModels/incompressible/LES/DeardorffDiffStress/DeardorffDiffStress.H
@ -105,7 +105,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/GenEddyVisc/GenEddyVisc.C
+++ b/src/turbulenceModels/incompressible/LES/GenEddyVisc/GenEddyVisc.C
@ -45,7 +45,7 @@ GenEddyVisc::GenEddyVisc
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/GenEddyVisc/GenEddyVisc.H
+++ b/src/turbulenceModels/incompressible/LES/GenEddyVisc/GenEddyVisc.H
@ -91,7 +91,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/GenSGSStress/GenSGSStress.C
+++ b/src/turbulenceModels/incompressible/LES/GenSGSStress/GenSGSStress.C
@ -45,7 +45,7 @@ GenSGSStress::GenSGSStress
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/GenSGSStress/GenSGSStress.H
+++ b/src/turbulenceModels/incompressible/LES/GenSGSStress/GenSGSStress.H
@ -91,7 +91,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/LESModel/LESModel.C
+++ b/src/turbulenceModels/incompressible/LES/LESModel/LESModel.C
@ -57,7 +57,7 @@ LESModel::LESModel
    const word& type,
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName
 )
 :
@ -95,7 +95,7 @@ autoPtr<LESModel> LESModel::New
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName
 )
 {
--- a/src/turbulenceModels/incompressible/LES/LESModel/LESModel.H
+++ b/src/turbulenceModels/incompressible/LES/LESModel/LESModel.H
@ -122,7 +122,7 @@ public:
            (
                const volVectorField& U,
                const surfaceScalarField& phi,
-                transportModel& transport,
+                const transportModel& transport,
                const word& turbulenceModelName
            ),
            (U, phi, transport, turbulenceModelName)
@ -137,7 +137,7 @@ public:
            const word& type,
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName
        );

@ -149,7 +149,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName
        );

--- a/src/turbulenceModels/incompressible/LES/LRRDiffStress/LRRDiffStress.C
+++ b/src/turbulenceModels/incompressible/LES/LRRDiffStress/LRRDiffStress.C
@ -55,7 +55,7 @@ LRRDiffStress::LRRDiffStress
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/LRRDiffStress/LRRDiffStress.H
+++ b/src/turbulenceModels/incompressible/LES/LRRDiffStress/LRRDiffStress.H
@ -102,7 +102,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/Smagorinsky/Smagorinsky.C
+++ b/src/turbulenceModels/incompressible/LES/Smagorinsky/Smagorinsky.C
@ -55,7 +55,7 @@ Smagorinsky::Smagorinsky
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/Smagorinsky/Smagorinsky.H
+++ b/src/turbulenceModels/incompressible/LES/Smagorinsky/Smagorinsky.H
@ -99,7 +99,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.C
+++ b/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.C
@ -46,7 +46,7 @@ Smagorinsky2::Smagorinsky2
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.H
+++ b/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.H
@ -96,7 +96,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.C
@ -180,7 +180,7 @@ SpalartAllmaras::SpalartAllmaras
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.H
+++ b/src/turbulenceModels/incompressible/LES/SpalartAllmaras/SpalartAllmaras.H
@ -157,7 +157,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/SpalartAllmarasDDES/SpalartAllmarasDDES.C
+++ b/src/turbulenceModels/incompressible/LES/SpalartAllmarasDDES/SpalartAllmarasDDES.C
@ -99,7 +99,7 @@ SpalartAllmarasDDES::SpalartAllmarasDDES
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/SpalartAllmarasDDES/SpalartAllmarasDDES.H
+++ b/src/turbulenceModels/incompressible/LES/SpalartAllmarasDDES/SpalartAllmarasDDES.H
@ -100,7 +100,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/SpalartAllmarasIDDES.C
+++ b/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/SpalartAllmarasIDDES.C
@ -146,7 +146,7 @@ SpalartAllmarasIDDES::SpalartAllmarasIDDES
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/SpalartAllmarasIDDES.H
+++ b/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/SpalartAllmarasIDDES.H
@ -109,7 +109,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.C
+++ b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.C
@ -58,7 +58,7 @@ dynLagrangian::dynLagrangian
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
+++ b/src/turbulenceModels/incompressible/LES/dynLagrangian/dynLagrangian.H
@ -137,7 +137,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
+++ b/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.C
@ -106,7 +106,7 @@ dynOneEqEddy::dynOneEqEddy
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.H
+++ b/src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.H
@ -135,7 +135,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/homogeneousDynOneEqEddy/homogeneousDynOneEqEddy.C
+++ b/src/turbulenceModels/incompressible/LES/homogeneousDynOneEqEddy/homogeneousDynOneEqEddy.C
@ -122,7 +122,7 @@ homogeneousDynOneEqEddy::homogeneousDynOneEqEddy
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/homogeneousDynOneEqEddy/homogeneousDynOneEqEddy.H
+++ b/src/turbulenceModels/incompressible/LES/homogeneousDynOneEqEddy/homogeneousDynOneEqEddy.H
@ -113,7 +113,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.C
+++ b/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.C
@ -110,7 +110,7 @@ homogeneousDynSmagorinsky::homogeneousDynSmagorinsky
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.H
+++ b/src/turbulenceModels/incompressible/LES/homogeneousDynSmagorinsky/homogeneousDynSmagorinsky.H
@ -126,7 +126,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/kOmegaSSTSAS/kOmegaSSTSAS.C
+++ b/src/turbulenceModels/incompressible/LES/kOmegaSSTSAS/kOmegaSSTSAS.C
@ -120,7 +120,7 @@ kOmegaSSTSAS::kOmegaSSTSAS
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/kOmegaSSTSAS/kOmegaSSTSAS.H
+++ b/src/turbulenceModels/incompressible/LES/kOmegaSSTSAS/kOmegaSSTSAS.H
@ -199,7 +199,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/laminar/laminar.C
+++ b/src/turbulenceModels/incompressible/LES/laminar/laminar.C
@ -47,7 +47,7 @@ laminar::laminar
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/laminar/laminar.H
+++ b/src/turbulenceModels/incompressible/LES/laminar/laminar.H
@ -76,7 +76,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/mixedSmagorinsky/mixedSmagorinsky.C
+++ b/src/turbulenceModels/incompressible/LES/mixedSmagorinsky/mixedSmagorinsky.C
@ -46,7 +46,7 @@ mixedSmagorinsky::mixedSmagorinsky
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/mixedSmagorinsky/mixedSmagorinsky.H
+++ b/src/turbulenceModels/incompressible/LES/mixedSmagorinsky/mixedSmagorinsky.H
@ -102,7 +102,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.C
+++ b/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.C
@ -56,7 +56,7 @@ oneEqEddy::oneEqEddy
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.H
+++ b/src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.H
@ -105,7 +105,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/scaleSimilarity/scaleSimilarity.C
+++ b/src/turbulenceModels/incompressible/LES/scaleSimilarity/scaleSimilarity.C
@ -44,7 +44,7 @@ scaleSimilarity::scaleSimilarity
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/scaleSimilarity/scaleSimilarity.H
+++ b/src/turbulenceModels/incompressible/LES/scaleSimilarity/scaleSimilarity.H
@ -87,7 +87,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/LES/spectEddyVisc/spectEddyVisc.C
+++ b/src/turbulenceModels/incompressible/LES/spectEddyVisc/spectEddyVisc.C
@ -65,7 +65,7 @@ spectEddyVisc::spectEddyVisc
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/LES/spectEddyVisc/spectEddyVisc.H
+++ b/src/turbulenceModels/incompressible/LES/spectEddyVisc/spectEddyVisc.H
@ -107,7 +107,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/LRR/LRR.C
+++ b/src/turbulenceModels/incompressible/RAS/LRR/LRR.C
@ -49,7 +49,7 @@ LRR::LRR
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/LRR/LRR.H
+++ b/src/turbulenceModels/incompressible/RAS/LRR/LRR.H
@ -115,7 +115,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.C
+++ b/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.C
@ -48,7 +48,7 @@ LamBremhorstKE::LamBremhorstKE
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.H
+++ b/src/turbulenceModels/incompressible/RAS/LamBremhorstKE/LamBremhorstKE.H
@ -96,7 +96,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
+++ b/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
@ -49,7 +49,7 @@ LaunderGibsonRSTM::LaunderGibsonRSTM
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.H
+++ b/src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.H
@ -127,7 +127,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
+++ b/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C
@ -64,7 +64,7 @@ LaunderSharmaKE::LaunderSharmaKE
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.H
+++ b/src/turbulenceModels/incompressible/RAS/LaunderSharmaKE/LaunderSharmaKE.H
@ -108,7 +108,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.C
+++ b/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.C
@ -48,7 +48,7 @@ LienCubicKE::LienCubicKE
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.H
+++ b/src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.H
@ -104,7 +104,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.C
+++ b/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.C
@ -49,7 +49,7 @@ LienCubicKELowRe::LienCubicKELowRe
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.H
+++ b/src/turbulenceModels/incompressible/RAS/LienCubicKELowRe/LienCubicKELowRe.H
@ -134,7 +134,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.C
+++ b/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.C
@ -49,7 +49,7 @@ LienLeschzinerLowRe::LienLeschzinerLowRe
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.H
+++ b/src/turbulenceModels/incompressible/RAS/LienLeschzinerLowRe/LienLeschzinerLowRe.H
@ -104,7 +104,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.C
+++ b/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.C
@ -48,7 +48,7 @@ NonlinearKEShih::NonlinearKEShih
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.H
+++ b/src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.H
@ -107,7 +107,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C
+++ b/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C
@ -57,7 +57,7 @@ RASModel::RASModel
    const word& type,
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName
 )
 :
@ -99,7 +99,7 @@ autoPtr<RASModel> RASModel::New
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName
 )
 {
--- a/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.H
+++ b/src/turbulenceModels/incompressible/RAS/RASModel/RASModel.H
@ -130,7 +130,7 @@ public:
            (
                const volVectorField& U,
                const surfaceScalarField& phi,
-                transportModel& transport,
+                const transportModel& transport,
                const word& turbulenceModelName
            ),
            (U, phi, transport, turbulenceModelName)
@ -145,7 +145,7 @@ public:
            const word& type,
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName
        );

@ -157,7 +157,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName
        );

--- a/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C
+++ b/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C
@ -48,7 +48,7 @@ RNGkEpsilon::RNGkEpsilon
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.H
+++ b/src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.H
@ -114,7 +114,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
+++ b/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.C
@ -143,7 +143,7 @@ SpalartAllmaras::SpalartAllmaras
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.H
+++ b/src/turbulenceModels/incompressible/RAS/SpalartAllmaras/SpalartAllmaras.H
@ -159,7 +159,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C
+++ b/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C
@ -48,7 +48,7 @@ kEpsilon::kEpsilon
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.H
+++ b/src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.H
@ -111,7 +111,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.C
+++ b/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.C
@ -48,7 +48,7 @@ kOmega::kOmega
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.H
+++ b/src/turbulenceModels/incompressible/RAS/kOmega/kOmega.H
@ -114,7 +114,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
+++ b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
@ -117,7 +117,7 @@ kOmegaSST::kOmegaSST
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.H
+++ b/src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.H
@ -206,7 +206,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C
+++ b/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C
@ -206,7 +206,7 @@ kkLOmega::kkLOmega
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.H
+++ b/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.H
@ -200,7 +200,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/laminar/laminar.C
+++ b/src/turbulenceModels/incompressible/RAS/laminar/laminar.C
@ -46,7 +46,7 @@ laminar::laminar
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/laminar/laminar.H
+++ b/src/turbulenceModels/incompressible/RAS/laminar/laminar.H
@ -70,7 +70,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.C
+++ b/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.C
@ -74,7 +74,7 @@ qZeta::qZeta
 (
    const volVectorField& U,
    const surfaceScalarField& phi,
-    transportModel& transport,
+    const transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
 )
--- a/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.H
+++ b/src/turbulenceModels/incompressible/RAS/qZeta/qZeta.H
@ -119,7 +119,7 @@ public:
        (
            const volVectorField& U,
            const surfaceScalarField& phi,
-            transportModel& transport,
+            const transportModel& transport,
            const word& turbulenceModelName = turbulenceModel::typeName,
            const word& modelName = typeName
        );
--- a/Show More
+++ b/Show More