Skip to content
Snippets Groups Projects
Commit afcd290d authored by Henry's avatar Henry
Browse files

Solvers based on p_rgh: Added support for optional hRef 

Allows the specification of a reference height, for example the height
of the free-surface in a VoF simulation, which reduces the range of p_rgh.

hRef is a uniformDimensionedScalarField specified via the constant/hRef
file, equivalent to the way in which g is specified, so that it can be
looked-up from the database.  For example see the constant/hRef file in
the DTCHull LTSInterFoam and interDyMFoam cases.
parent 930d3d43
Branches
Tags
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 79 additions and 35 deletions
applications/solvers/combustion/fireFoam/createFields.H
+ 6
− 2
View file @ afcd290d
......@@ -101,10 +101,14 @@
volScalarField K("K", 0.5*magSqr(U));
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
surfaceScalarField ghf("ghf", g & mesh.Cf());
volScalarField p_rgh
(
......
applications/solvers/combustion/fireFoam/fireFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -54,7 +54,6 @@ int main(int argc, char *argv[])
pimpleControl pimple(mesh);
#include "readGravitationalAcceleration.H"
#include "createFields.H"
#include "createFvOptions.H"
#include "createClouds.H"
......
applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/createFields.H
+ 7
− 2
View file @ afcd290d
......@@ -63,9 +63,14 @@ autoPtr<compressible::turbulenceModel> turbulence
reaction->setTurbulence(turbulence());
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
Info<< "Reading field p_rgh\n" << endl;
volScalarField p_rgh
......
applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/rhoReactingBuoyantFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -48,7 +48,6 @@ int main(int argc, char *argv[])
pimpleControl pimple(mesh);
#include "readGravitationalAcceleration.H"
#include "createFields.H"
#include "createFvOptions.H"
#include "initContinuityErrs.H"
......
applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/buoyantBoussinesqPimpleFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -63,7 +63,6 @@ int main(int argc, char *argv[])
pimpleControl pimple(mesh);
#include "readGravitationalAcceleration.H"
#include "createFields.H"
#include "createIncompressibleRadiationModel.H"
#include "createFvOptions.H"
......
applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H
+ 8
− 2
View file @ afcd290d
......@@ -79,9 +79,15 @@
mesh
);
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
volScalarField p
(
......
applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -62,7 +62,6 @@ int main(int argc, char *argv[])
simpleControl simple(mesh);
#include "readGravitationalAcceleration.H"
#include "createFields.H"
#include "createFvOptions.H"
#include "initContinuityErrs.H"
......
applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H
+ 8
− 2
View file @ afcd290d
......@@ -79,9 +79,15 @@
mesh
);
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
volScalarField p
(
......
applications/solvers/lagrangian/reactingParcelFilmFoam/createFields.H
+ 7
− 2
View file @ afcd290d
......@@ -80,10 +80,15 @@
Info<< "Creating field kinetic energy K\n" << endl;
volScalarField K("K", 0.5*magSqr(U));
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
surfaceScalarField ghf("ghf", g & mesh.Cf());
volScalarField p_rgh
(
......
applications/solvers/lagrangian/reactingParcelFilmFoam/reactingParcelFilmFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -52,7 +52,6 @@ int main(int argc, char *argv[])
pimpleControl pimple(mesh);
#include "readGravitationalAcceleration.H"
#include "createFields.H"
#include "createFvOptions.H"
#include "createClouds.H"
......
applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
+ 2
− 3
View file @ afcd290d
......@@ -56,7 +56,6 @@ int main(int argc, char *argv[])
#include "setRootCase.H"
#include "createTime.H"
#include "createDynamicFvMesh.H"
#include "readGravitationalAcceleration.H"
#include "initContinuityErrs.H"
pimpleControl pimple(mesh);
......@@ -98,8 +97,8 @@ int main(int argc, char *argv[])
<< runTime.elapsedCpuTime() - timeBeforeMeshUpdate
<< " s" << endl;
gh = g & mesh.C();
ghf = g & mesh.Cf();
gh = (g & mesh.C()) - ghRef;
ghf = (g & mesh.Cf()) - ghRef;
}
if (mesh.changing() && correctPhi)
......
applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -53,7 +53,6 @@ int main(int argc, char *argv[])
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"
#include "readGravitationalAcceleration.H"
pimpleControl pimple(mesh);
......
applications/solvers/multiphase/compressibleInterFoam/createFields.H
+ 14
− 3
View file @ afcd290d
......@@ -57,11 +57,22 @@
);
dimensionedScalar pMin(mixture.lookup("pMin"));
dimensionedScalar pMin
(
"pMin",
dimPressure,
mixture.lookup("pMin")
);
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
// Mass flux
// Initialisation does not matter because rhoPhi is reset after the
......
applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseInterFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -48,7 +48,6 @@ int main(int argc, char *argv[])
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"
#include "readGravitationalAcceleration.H"
pimpleControl pimple(mesh);
......
applications/solvers/multiphase/compressibleMultiphaseInterFoam/createFields.H
+ 7
− 2
View file @ afcd290d
......@@ -49,9 +49,14 @@
dimensionedScalar pMin(mixture.lookup("pMin"));
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
// Construct compressible turbulence model
autoPtr<compressible::turbulenceModel> turbulence
......
applications/solvers/multiphase/driftFluxFoam/createFields.H
+ 9
− 3
View file @ afcd290d
......@@ -90,9 +90,15 @@
::New(rho, U, rhoPhi, mixture)
);
Info<< "Calculating field (g.h)f\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
volScalarField p
(
......
applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
+ 0
− 1
View file @ afcd290d
......@@ -57,7 +57,6 @@ int main(int argc, char *argv[])
pimpleControl pimple(mesh);
#include "readGravitationalAcceleration.H"
#include "createFields.H"
#include "createFvOptions.H"
#include "initContinuityErrs.H"
......
applications/solvers/multiphase/interFoam/createFields.H
+ 4
− 2
View file @ afcd290d
......@@ -77,10 +77,12 @@
);
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
volScalarField p
(
......
applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
+ 2
− 2
View file @ afcd290d
......@@ -110,8 +110,8 @@ int main(int argc, char *argv[])
<< runTime.elapsedCpuTime() - timeBeforeMeshUpdate
<< " s" << endl;
gh = g & mesh.C();
ghf = g & mesh.Cf();
gh = (g & mesh.C()) - ghRef;
ghf = (g & mesh.Cf()) - ghRef;
}
if (mesh.changing() && correctPhi)
......
applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H
+ 5
− 2
View file @ afcd290d
......@@ -79,10 +79,13 @@
incompressible::turbulenceModel::New(U, phi, mixture)
);
#include "readGravitationalAcceleration.H"
#include "readhRef.H"
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
dimensionedScalar ghRef(g & (cmptMag(g.value())/mag(g.value()))*hRef);
volScalarField gh("gh", (g & mesh.C()) - ghRef);
surfaceScalarField ghf("ghf", (g & mesh.Cf()) - ghRef);
volScalarField p
(
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment