Buoyant solvers: Add special handling for ghRef in the case g = (0 0 0)

This commit is contained in:
Henry
2015-03-19 15:26:06 +00:00
parent 29b9eaaf4d
commit 28abb8cda5
15 changed files with 24 additions and 71 deletions

View File

@ -55,11 +55,7 @@
#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);
#include "gh.H"
volScalarField p