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: Corrected surface film contact angle force model

Browse Source
This commit is contained in:
andy
2013-02-15 11:35:01 +00:00
parent 124ad071a7
commit 54c55d9f41
1 changed files with 7 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C
View File

@ -103,7 +103,7 @@ contactAngleForce::contactAngleForce
IOobject::NO_WRITE
),
owner_.regionMesh(),
dimensionedScalar("mask", dimless, 0.0)
dimensionedScalar("mask", dimless, 1.0)
)
{
initialise();
@ -166,13 +166,13 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
cellI = cellN;
}
if (cellI != -1 && mask_[cellI] > 0)
if (cellI != -1 && mask_[cellI] > 0.5)
{
const scalar dx = owner_.regionMesh().deltaCoeffs()[faceI];
const scalar invDx = owner_.regionMesh().deltaCoeffs()[faceI];
const vector n =
gradAlpha[cellI]/(mag(gradAlpha[cellI]) + ROOTVSMALL);
scalar theta = cos(degToRad(distribution_->sample()));
force[cellI] += Ccf_*n*sigma[cellI]*(1.0 - theta)/dx;
force[cellI] += Ccf_*n*sigma[cellI]*(1.0 - theta)/invDx;
nHits[cellI]++;
}
}
@ -183,12 +183,12 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
{
const fvPatchField<scalar>& alphaf = alpha.boundaryField()[patchI];
const fvPatchField<scalar>& maskf = mask_.boundaryField()[patchI];
const scalarField& dx = alphaf.patch().deltaCoeffs();
const scalarField& invDx = alphaf.patch().deltaCoeffs();
const labelUList& faceCells = alphaf.patch().faceCells();
forAll(alphaf, faceI)
{
if (maskf[faceI] > 0)
if (maskf[faceI] > 0.5)
{
label cellO = faceCells[faceI];
@ -199,7 +199,7 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
/(mag(gradAlpha[cellO]) + ROOTVSMALL);
scalar theta = cos(degToRad(distribution_->sample()));
force[cellO] +=
Ccf_*n*sigma[cellO]*(1.0 - theta)/dx[faceI];
Ccf_*n*sigma[cellO]*(1.0 - theta)/invDx[faceI];
nHits[cellO]++;
}
}