diff --git a/applications/solvers/incompressible/potentialFreeSurfaceFoam/pEqn.H b/applications/solvers/incompressible/potentialFreeSurfaceFoam/pEqn.H index 6f6b9834e4..888d06a889 100644 --- a/applications/solvers/incompressible/potentialFreeSurfaceFoam/pEqn.H +++ b/applications/solvers/incompressible/potentialFreeSurfaceFoam/pEqn.H @@ -1,3 +1,5 @@ +#include "resetPhiPatches.H" + volScalarField rAU(1.0/UEqn().A()); surfaceScalarField rAUf("Dp", fvc::interpolate(rAU)); diff --git a/applications/solvers/incompressible/potentialFreeSurfaceFoam/resetPhiPatches.H b/applications/solvers/incompressible/potentialFreeSurfaceFoam/resetPhiPatches.H new file mode 100644 index 0000000000..fe428b6e3b --- /dev/null +++ b/applications/solvers/incompressible/potentialFreeSurfaceFoam/resetPhiPatches.H @@ -0,0 +1,8 @@ +FieldField& phibf = phi.boundaryField(); +const FieldField& Ubf = U.boundaryField(); +const FieldField& Sfbf = mesh.Sf().boundaryField(); + +forAll(phibf, patchI) +{ + phibf[patchI] = (Ubf[patchI] & Sfbf[patchI]); +} diff --git a/applications/test/GAMGAgglomeration/Test-GAMGAgglomeration.C b/applications/test/GAMGAgglomeration/Test-GAMGAgglomeration.C index 82ccd8ce33..2dd7905294 100644 --- a/applications/test/GAMGAgglomeration/Test-GAMGAgglomeration.C +++ b/applications/test/GAMGAgglomeration/Test-GAMGAgglomeration.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -117,6 +117,30 @@ int main(int argc, char *argv[]) << " agglomerated size : " << returnReduce(coarseSize, sumOp