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

fvMeshStitcher: Stabilise calculation of error face areas

Browse Source
This commit is contained in:
Will Bainbridge
2023-05-12 11:45:51 +01:00
parent 3479c6a37f
commit c9bf66a632
1 changed files with 6 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
src/fvMeshStitchers/moving/fvMeshStitchersMoving.C
View File

@ -987,11 +987,14 @@ void Foam::fvMeshStitchers::moving::unconformErrorFaceCorrectMeshPhi
const label errorPatchFacei0 = 2*origPatchFacei; const label errorPatchFacei0 = 2*origPatchFacei;
const label errorPatchFacei1 = 2*origPatchFacei + 1; const label errorPatchFacei1 = 2*origPatchFacei + 1;
const vector errorSf = const scalar a =
origPp.faceNormals()[origPatchFacei] phi.boundaryField()[errorPatchi][errorPatchFacei0]
*phi.boundaryField()[errorPatchi][errorPatchFacei0]
/max(meshMagUfb[origPatchi][origPatchFacei], vSmall); /max(meshMagUfb[origPatchi][origPatchFacei], vSmall);
const vector errorSf =
sign(a)*min(mag(a), origPp.magFaceAreas()[origPatchFacei])
*origPp.faceNormals()[origPatchFacei];
fvsPatchField<vector>& Sfp = fvsPatchField<vector>& Sfp =
SfSf.boundaryFieldRef()[errorPatchi]; SfSf.boundaryFieldRef()[errorPatchi];
Sfp[errorPatchFacei0] += errorSf; Sfp[errorPatchFacei0] += errorSf;