zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

finiteVolume: LimitedScheme: Correction when using "patchType" overrides

Browse Source
This commit is contained in:
Will Bainbridge
2021-06-08 14:59:38 +01:00
parent e6ce04f070
commit d739b7e36f
1 changed files with 5 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/LimitedScheme/LimitedScheme.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2021 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -76,6 +76,9 @@ void Foam::LimitedScheme<Type, Limiter, LimitFunc>::calcLimiter
); );
} }
const typename GeometricField<Type, fvPatchField, volMesh>::Boundary&
bPhi = phi.boundaryField();
surfaceScalarField::Boundary& bLim = surfaceScalarField::Boundary& bLim =
limiterField.boundaryFieldRef(); limiterField.boundaryFieldRef();
@ -83,7 +86,7 @@ void Foam::LimitedScheme<Type, Limiter, LimitFunc>::calcLimiter
{ {
scalarField& pLim = bLim[patchi]; scalarField& pLim = bLim[patchi];
if (bLim[patchi].coupled()) if (bPhi[patchi].coupled())
{ {
const scalarField& pCDweights = CDweights.boundaryField()[patchi]; const scalarField& pCDweights = CDweights.boundaryField()[patchi];
const scalarField& pFaceFlux = const scalarField& pFaceFlux =