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

movingMappedWallVelocity: Fix calculation of mapped shear velocity

Browse Source
This commit is contained in:
Will Bainbridge
2024-09-11 16:52:25 +01:00
parent f1068151ee
commit 2853f9ddfd
1 changed files with 7 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/finiteVolume/fields/fvPatchFields/derived/movingMappedWallVelocity/movingMappedWallVelocityFvPatchVectorField.C
View File

@ -107,18 +107,19 @@ void Foam::movingMappedWallVelocityFvPatchVectorField::updateCoeffs()
: tmp<scalarField>(new scalarField(fvp.size(), Zero)); : tmp<scalarField>(new scalarField(fvp.size(), Zero));
// Calculate and map the mesh velocity from the neighbour // Calculate and map the mesh velocity from the neighbour
vectorField nbrCf0(nbrFvp.size()); vectorField nbrCf(nbrFvp.size()), nbrCf0(nbrFvp.size());
forAll(nbrCf0, nbrPatchFacei) forAll(nbrCf0, nbrPatchFacei)
{ {
const label nbrPolyFacei =
nbrMesh.polyFacesBf()[nbrFvp.index()][nbrPatchFacei];
nbrCf[nbrPatchFacei] =
nbrMesh.faceCentres()[nbrPolyFacei];
nbrCf0[nbrPatchFacei] = nbrCf0[nbrPatchFacei] =
nbrMesh.faces() nbrMesh.faces()[nbrPolyFacei].centre(nbrMesh.oldPoints());
[
nbrMesh.polyFacesBf()[nbrFvp.index()][nbrPatchFacei]
].centre(nbrMesh.oldPoints());
} }
const vectorField nbrUp const vectorField nbrUp
( (
(nbrFvp.Cf() - nbrCf0)/db().time().deltaTValue() (nbrCf - nbrCf0)/db().time().deltaTValue()
); );
const vectorField Up(mapper.fromNeighbour(nbrUp)); const vectorField Up(mapper.fromNeighbour(nbrUp));