zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: facGrad: ensure compatibility between caching and corrections (fixes #3039)

Browse Source
This commit is contained in:
Kutalmis Bercin
2023-12-04 15:13:44 +00:00
parent 4609aa38e1
commit 1482547756
1 changed files with 10 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/finiteArea/finiteArea/fac/facGrad.C
View File

@ -103,13 +103,17 @@ grad
typedef typename outerProduct<vector,Type>::type GradType;
typedef GeometricField<GradType, faPatchField, areaMesh> GradFieldType;
tmp<GradFieldType> tgGrad = fa::gradScheme<Type>::New
tmp<GradFieldType> tgGrad
(
vf.mesh(),
vf.mesh().gradScheme(name)
).ref().grad(vf);
GradFieldType& gGrad = tgGrad.ref();
fa::gradScheme<Type>::New
(
vf.mesh(),
vf.mesh().gradScheme(name)
)
.cref() // const ref to tmp contents
.grad(vf).ptr() // steal ptr or deep copy of cached gradient
);
auto& gGrad = tgGrad.ref();
const areaVectorField& n = vf.mesh().faceAreaNormals();
gGrad -= n*(n & gGrad);