48 void Foam::interfaceProperties::correctContactAngle
50 surfaceVectorField::GeometricBoundaryField& nHatb,
51 surfaceVectorField::GeometricBoundaryField& gradAlphaf
55 const volScalarField::GeometricBoundaryField& abf = alpha1_.
boundaryField();
57 const fvBoundaryMesh&
boundary = mesh.boundary();
61 if (isA<alphaContactAngleFvPatchScalarField>(abf[
patchi]))
63 alphaContactAngleFvPatchScalarField& acap =
64 const_cast<alphaContactAngleFvPatchScalarField&
>
66 refCast<const alphaContactAngleFvPatchScalarField>
88 b2[facei] =
cos(
acos(a12[facei]) - theta[facei]);
96 nHatp = a*nf + b*nHatp;
98 nHatp /= (
mag(nHatp) + deltaN_.
value());
100 acap.gradient() = (nf & nHatp)*
mag(gradAlphaf[patchi]);
107 void Foam::interfaceProperties::calculateK()
109 const fvMesh& mesh = alpha1_.mesh();
123 correctContactAngle(nHatfv.boundaryField(), gradAlphaf.boundaryField());
126 nHatf_ = nHatfv & Sf;
147 Foam::interfaceProperties::interfaceProperties
151 const IOdictionary& dict
154 transportPropertiesDict_(dict),
159 alpha1.mesh().
solutionDict().subDict(
"PISO").lookup(
"cAlpha")
162 sigma_(dict.lookup(
"sigma")),