Bugfix: consistent notation of 1/aP

This commit is contained in:
Hrvoje Jasak 2016-08-01 14:38:56 +01:00
parent 5e7de9eb8b
commit 2e1ab51dbe

View file

@ -2,7 +2,7 @@
// Prepare clean 1/Ap without contribution from under-relaxation
// HJ, 26/Oct/2015
volScalarField rUA
volScalarField rAU
(
"(1|A(U))",
1/HUEqn().A()
@ -12,7 +12,7 @@
// the flux precursor
U.storePrevIter();
U = rUA*HUEqn().H();
U = rAU*HUEqn().H();
HUEqn.clear();
phi = fvc::interpolate(U) & mesh.Sf();
adjustPhi(phi, U, p);