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

ENH: Updated film contact angle force calculation

Browse Source
This commit is contained in:
andy
2013-02-15 16:55:53 +00:00
parent 1161b53506
commit bb473d2ae3
1 changed files with 1 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -149,8 +149,6 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
volVectorField gradAlpha(fvc::grad(alpha));
scalarField nHits(owner_.regionMesh().nCells(), 0.0);
forAll(nbr, faceI)
{
const label cellO = own[faceI];
@ -173,7 +171,6 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
gradAlpha[cellI]/(mag(gradAlpha[cellI]) + ROOTVSMALL);
scalar theta = cos(degToRad(distribution_->sample()));
force[cellI] += Ccf_*n*sigma[cellI]*(1.0 - theta)/invDx;
nHits[cellI]++;
}
}
@ -200,15 +197,13 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
scalar theta = cos(degToRad(distribution_->sample()));
force[cellO] +=
Ccf_*n*sigma[cellO]*(1.0 - theta)/invDx[faceI];
nHits[cellO]++;
}
}
}
}
}
force /= (max(nHits, scalar(1.0))*magSf);
tForce().correctBoundaryConditions();
force /= magSf;
if (owner_.regionMesh().time().outputTime())
{