From 00a6e3abb43108a9c0b672818942f27b5ccce4b3 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Fri, 4 Jan 2013 13:56:47 +0000
Subject: [PATCH] ENH: Updated film contact angle force
---
.../contactAngleForce/contactAngleForce.C | 73 ++++++++++++++++---
.../contactAngleForce/contactAngleForce.H | 10 ++-
2 files changed, 71 insertions(+), 12 deletions(-)
diff --git a/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C b/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C
index 8b9f8a51dda..2909b23c4a1 100644
--- a/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C
+++ b/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -28,6 +28,7 @@ License
#include "fvcGrad.H"
#include "unitConversion.H"
#include "fvPatchField.H"
+#include "patchDist.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -43,6 +44,35 @@ namespace surfaceFilmModels
defineTypeNameAndDebug(contactAngleForce, 0);
addToRunTimeSelectionTable(force, contactAngleForce, dictionary);
+// * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * * //
+
+void contactAngleForce::initialise()
+{
+ const wordReList zeroForcePatches(coeffs_.lookup("zeroForcePatches"));
+
+ if (zeroForcePatches.size())
+ {
+ const polyBoundaryMesh& pbm = owner_.regionMesh().boundaryMesh();
+ scalar dLim = readScalar(coeffs_.lookup("zeroForceDistance"));
+
+ Info<< " Assigning zero contact force within " << dLim
+ << " of patches:" << endl;
+
+ labelHashSet patchIDs = pbm.patchSet(zeroForcePatches);
+
+ forAllConstIter(labelHashSet, patchIDs, iter)
+ {
+ label patchI = iter.key();
+ Info<< " " << pbm[patchI].name() << endl;
+ }
+
+ patchDist dist(owner_.regionMesh(), patchIDs);
+
+ mask_ = pos(dist - dimensionedScalar("dLim", dimLength, dLim));
+ }
+}
+
+
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
contactAngleForce::contactAngleForce
@@ -61,8 +91,23 @@ contactAngleForce::contactAngleForce
coeffs_.subDict("contactAngleDistribution"),
rndGen_
)
+ ),
+ mask_
+ (
+ IOobject
+ (
+ typeName + ".contactForceMask",
+ owner_.time().timeName(),
+ owner_.regionMesh(),
+ IOobject::NO_READ,
+ IOobject::NO_WRITE
+ ),
+ owner_.regionMesh(),
+ dimensionedScalar("mask", dimless, 0.0)
)
-{}
+{
+ initialise();
+}
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
@@ -121,7 +166,7 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
cellI = cellN;
}
- if (cellI != -1)
+ if (cellI != -1 && mask_[cellI] > 0)
{
const scalar dx = owner_.regionMesh().deltaCoeffs()[faceI];
const vector n =
@@ -137,20 +182,26 @@ tmp<fvVectorMatrix> contactAngleForce::correct(volVectorField& U)
if (!owner().isCoupledPatch(patchI))
{
const fvPatchField<scalar>& alphaf = alpha.boundaryField()[patchI];
+ const fvPatchField<scalar>& maskf = mask_.boundaryField()[patchI];
const scalarField& dx = alphaf.patch().deltaCoeffs();
const labelUList& faceCells = alphaf.patch().faceCells();
forAll(alphaf, faceI)
{
- label cellO = faceCells[faceI];
-
- if ((alpha[cellO] > 0.5) && (alphaf[faceI] < 0.5))
+ if (maskf[faceI] > 0)
{
- const vector n =
- gradAlpha[cellO]/(mag(gradAlpha[cellO]) + ROOTVSMALL);
- scalar theta = cos(degToRad(distribution_->sample()));
- force[cellO] += Ccf_*n*sigma[cellO]*(1.0 - theta)/dx[faceI];
- nHits[cellO]++;
+ label cellO = faceCells[faceI];
+
+ if ((alpha[cellO] > 0.5) && (alphaf[faceI] < 0.5))
+ {
+ const vector n =
+ gradAlpha[cellO]
+ /(mag(gradAlpha[cellO]) + ROOTVSMALL);
+ scalar theta = cos(degToRad(distribution_->sample()));
+ force[cellO] +=
+ Ccf_*n*sigma[cellO]*(1.0 - theta)/dx[faceI];
+ nHits[cellO]++;
+ }
}
}
}
diff --git a/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.H b/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.H
index 8a6dc7a8b3b..a299ecd1b2c 100644
--- a/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.H
+++ b/src/regionModels/surfaceFilmModels/submodels/kinematic/force/contactAngleForce/contactAngleForce.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -27,6 +27,9 @@ Class
Description
Film contact angle force
+ The effect of the contact angle force can be ignored over a specified
+ distance from patches.
+
SourceFiles
contactAngleForce.C
@@ -69,10 +72,15 @@ private:
//- Parcel size PDF model
const autoPtr<distributionModels::distributionModel> distribution_;
+ //- Mask over which force is applied
+ volScalarField mask_;
// Private member functions
+ //- Initialise
+ void initialise();
+
//- Disallow default bitwise copy construct
contactAngleForce(const contactAngleForce&);
--
GitLab