UEqn: Removed "Sp" term from the momentum equation of the steady-state solvers.

It is not clear that this "correction" term has provided any benefit and there are cases for which it has proved detrimental to stability and convergence.

UEqn: Removed "Sp" term from the momentum equation of the steady-state solvers.
It is not clear that this "correction" term has provided any benefit and there are cases for which it has proved detrimental to stability and convergence.
ed8fb0aa · henry · bb1a2476 · ed8fb0aa · ed8fb0aa · ed8fb0aa
Commit ed8fb0aa authored 15 years ago by henry
--- a/applications/solvers/compressible/rhoPorousSimpleFoam/UEqn.H
+++ b/applications/solvers/compressible/rhoPorousSimpleFoam/UEqn.H
@@ -3,7 +3,6 @@
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
-      - fvm::Sp(fvc::div(phi), U)
      + turbulence->divDevRhoReff(U)
    );


--- a/applications/solvers/compressible/rhoSimpleFoam/UEqn.H
+++ b/applications/solvers/compressible/rhoSimpleFoam/UEqn.H
@@ -3,7 +3,6 @@
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
-      - fvm::Sp(fvc::div(phi), U)
      + turbulence->divDevRhoReff(U)
    );


--- a/applications/solvers/compressible/rhoSimplecFoam/UEqn.H
+++ b/applications/solvers/compressible/rhoSimplecFoam/UEqn.H
@@ -3,7 +3,6 @@
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
-      - fvm::Sp(fvc::div(phi), U)
      + turbulence->divDevRhoReff(U)
    );


--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/UEqn.H
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/UEqn.H
@@ -3,7 +3,6 @@
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
-      - fvm::Sp(fvc::div(phi), U)
      + turbulence->divDevReff(U)
    );


--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H
@@ -3,7 +3,6 @@
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
-      - fvm::Sp(fvc::div(phi), U)
      + turbulence->divDevRhoReff(U)
    );


--- a/applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/fluid/UEqn.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/fluid/UEqn.H
@@ -2,7 +2,6 @@
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
-      - fvm::Sp(fvc::div(phi), U)
      + turb.divDevRhoReff(U)
    );


--- a/applications/solvers/incompressible/porousSimpleFoam/UEqn.H
+++ b/applications/solvers/incompressible/porousSimpleFoam/UEqn.H
@@ -3,7 +3,6 @@
    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
-      - fvm::Sp(fvc::div(phi), U)
      + turbulence->divDevReff(U)
    );


--- a/tutorials/incompressible/MRFSimpleFoam/MRFSimpleFoam/MRFSimpleFoam.C
+++ b/tutorials/incompressible/MRFSimpleFoam/MRFSimpleFoam/MRFSimpleFoam.C
@@ -64,7 +64,6 @@ int main(int argc, char *argv[])
            tmp<fvVectorMatrix> UEqn
            (
                fvm::div(phi, U)
-              - fvm::Sp(fvc::div(phi), U)
              + turbulence->divDevReff(U)
            );
            mrfZones.addCoriolis(UEqn());