BUG: Correct treatment of contactAngleForce near boundaries

This commit is contained in:
sergio
2022-04-14 11:08:17 -07:00
parent c78101a9ae
commit 8b50f8db25
2 changed files with 3 additions and 16 deletions

View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2020 OpenCFD Ltd.
Copyright (C) 2020-2022 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -153,19 +153,6 @@ tmp<faVectorMatrix> contactAngleForce::correct(areaVectorField& U)
}
}
forAll(sigma.boundaryField(), patchi)
{
const faPatchField<scalar>& alphaPf = sigma.boundaryField()[patchi];
const faPatchField<scalar>& sigmaPf = sigma.boundaryField()[patchi];
const labelUList& faces = alphaPf.patch().edgeFaces();
forAll(sigmaPf, edgei)
{
label face0 = faces[edgei];
force[face0] = vector::zero;
}
}
force /= magSf.field();

View File

@ -5,9 +5,9 @@ cd "${0%/*}" || exit # Run from this directory
runApplication blockMesh
runApplication decomposePar
runApplication makeFaMesh
runParallel makeFaMesh
runApplication decomposePar
runParallel $(getApplication)