diff --git a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C index c79ac392d..6248707db 100644 --- a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C +++ b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C @@ -189,7 +189,10 @@ void Foam::movingImmersedBoundaryVelocityFvPatchVectorField::updateCoeffs() scalarField phip = p.patchField(fvc::meshPhi(U)); - vectorField n = p.nf(); + // Warning: cannot use patch normal but the real face normal + // THEY MAY NOT BE THE SAME! HJ, 28/Mar/2019 + vectorField n = p.Sf()/(p.magSf()); + const scalarField& magSf = p.magSf(); scalarField Un = phip/(magSf + VSMALL);