From e9a7c2baaef9ae64e97a21413ad42d556ba11ef0 Mon Sep 17 00:00:00 2001 From: Henry <Henry> Date: Sun, 6 Jul 2014 20:20:16 +0100 Subject: [PATCH] Revert non-orthogonality treatment for non-coupled BCs. --- .../fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.C | 9 --------- .../fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.H | 3 --- src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C | 4 +++- 3 files changed, 3 insertions(+), 13 deletions(-) diff --git a/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.C b/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.C index ddc60991217..f5ce1243ea6 100644 --- a/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.C +++ b/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.C @@ -46,15 +46,6 @@ wedgeFvPatch::wedgeFvPatch(const polyPatch& patch, const fvBoundaryMesh& bm) {} -// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // - -Foam::tmp<Foam::vectorField> Foam::wedgeFvPatch::delta() const -{ - const vectorField nHat(nf()); - return nHat*(nHat & (Cf() - Cn())); -} - - // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // } // End namespace Foam diff --git a/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.H b/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.H index 16fab3e42b8..25be556b1c1 100644 --- a/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.H +++ b/src/finiteVolume/fvMesh/fvPatches/constraint/wedge/wedgeFvPatch.H @@ -83,9 +83,6 @@ public: { return wedgePolyPatch_.cellT(); } - - //- Return cell-centre to face normal vector - virtual tmp<vectorField> delta() const; }; diff --git a/src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C b/src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C index c55ad272daa..e7dd41207e1 100644 --- a/src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C +++ b/src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C @@ -140,7 +140,9 @@ const Foam::scalarField& Foam::fvPatch::magSf() const Foam::tmp<Foam::vectorField> Foam::fvPatch::delta() const { - return Cf() - Cn(); + // Use patch-normal delta for all non-coupled BCs + const vectorField nHat(nf()); + return nHat*(nHat & (Cf() - Cn())); } -- GitLab