BUG: heat transfer Boussinesq solvers - added fvOptions MRF support - mantis #1267

This commit is contained in:
andy
2014-04-11 12:11:36 +01:00
committed by Andrew Heather
parent a3b1a7761b
commit a840324451
2 changed files with 6 additions and 2 deletions

View File

@ -15,13 +15,15 @@
+ phig + phig
); );
fvOptions.makeRelative(phiHbyA);
// Update the fixedFluxPressure BCs to ensure flux consistency // Update the fixedFluxPressure BCs to ensure flux consistency
setSnGrad<fixedFluxPressureFvPatchScalarField> setSnGrad<fixedFluxPressureFvPatchScalarField>
( (
p_rgh.boundaryField(), p_rgh.boundaryField(),
( (
phiHbyA.boundaryField() phiHbyA.boundaryField()
- (mesh.Sf().boundaryField() & U.boundaryField()) - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField())
)/(mesh.magSf().boundaryField()*rAUf.boundaryField()) )/(mesh.magSf().boundaryField()*rAUf.boundaryField())
); );

View File

@ -14,6 +14,8 @@
(fvc::interpolate(HbyA) & mesh.Sf()) (fvc::interpolate(HbyA) & mesh.Sf())
); );
fvOptions.makeRelative(phiHbyA);
adjustPhi(phiHbyA, U, p_rgh); adjustPhi(phiHbyA, U, p_rgh);
phiHbyA += phig; phiHbyA += phig;
@ -24,7 +26,7 @@
p_rgh.boundaryField(), p_rgh.boundaryField(),
( (
phiHbyA.boundaryField() phiHbyA.boundaryField()
- (mesh.Sf().boundaryField() & U.boundaryField()) - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField())
)/(mesh.magSf().boundaryField()*rAUf.boundaryField()) )/(mesh.magSf().boundaryField()*rAUf.boundaryField())
); );