From 1a56cdfd14c912d213f5be09d0c96e258eaecda2 Mon Sep 17 00:00:00 2001 From: Henry Date: Mon, 25 Nov 2013 10:10:00 +0000 Subject: [PATCH] interfaceProperties: Name grad(alpha1) "nHat" so that a specific scheme can be selected for it --- .../interfaceProperties/interfaceProperties.C | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/src/transportModels/interfaceProperties/interfaceProperties.C b/src/transportModels/interfaceProperties/interfaceProperties.C index cfa20630ae..9387769cae 100644 --- a/src/transportModels/interfaceProperties/interfaceProperties.C +++ b/src/transportModels/interfaceProperties/interfaceProperties.C @@ -112,7 +112,7 @@ void Foam::interfaceProperties::calculateK() const surfaceVectorField& Sf = mesh.Sf(); // Cell gradient of alpha - const volVectorField gradAlpha(fvc::grad(alpha1_)); + const volVectorField gradAlpha(fvc::grad(alpha1_, "nHat")); // Interpolated face-gradient of alpha surfaceVectorField gradAlphaf(fvc::interpolate(gradAlpha)); @@ -123,6 +123,11 @@ void Foam::interfaceProperties::calculateK() // Face unit interface normal surfaceVectorField nHatfv(gradAlphaf/(mag(gradAlphaf) + deltaN_)); + // surfaceVectorField nHatfv + // ( + // (gradAlphaf + deltaN_*vector(0, 0, 1) + // *sign(gradAlphaf.component(vector::Z)))/(mag(gradAlphaf) + deltaN_) + // ); correctContactAngle(nHatfv.boundaryField(), gradAlphaf.boundaryField()); // Face unit interface normal flux