Files
openfoam/src/transportModels/interfaceProperties/interfaceProperties.H
2008-06-25 15:01:46 +02:00

161 lines
4.3 KiB
C++

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 1991-2008 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation; either version 2 of the License, or (at your
option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM; if not, write to the Free Software Foundation,
Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
Class
Foam::interfaceProperties
Description
Contains the interface properties.
Properties to aid interFoam:
-# Correct the gamma boundary condition for dynamic contact angle.
-# Calculate interface curvature.
SourceFiles
interfaceProperties.C
\*---------------------------------------------------------------------------*/
#ifndef interfaceProperties_H
#define interfaceProperties_H
#include "IOdictionary.H"
#include "volFields.H"
#include "surfaceFields.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
/*---------------------------------------------------------------------------*\
Class interfaceProperties Declaration
\*---------------------------------------------------------------------------*/
class interfaceProperties
{
// Private data
//- Keep a reference to the transportProperties dictionary
const dictionary& transportPropertiesDict_;
//- Compression coefficient
scalar cGamma_;
//- Surface tension
dimensionedScalar sigma_;
//- Stabilisation for normalisation of the interface normal
const dimensionedScalar deltaN_;
const volScalarField& gamma_;
const volVectorField& U_;
surfaceScalarField nHatf_;
volScalarField K_;
// Private Member Functions
//- Disallow default bitwise copy construct and assignment
interfaceProperties(const interfaceProperties&);
void operator=(const interfaceProperties&);
//- Correction for the boundary condition on the unit normal nHat on
// walls to produce the correct contact dynamic angle
// calculated from the component of U parallel to the wall
void correctContactAngle
(
surfaceVectorField::GeometricBoundaryField& nHat
) const;
//- Re-calculate the interface curvature
void calculateK();
public:
//- Conversion factor for degrees into radians
static const scalar convertToRad;
// Constructors
//- Construct from volume fraction field gamma and IOdictionary
interfaceProperties
(
const volScalarField& gamma,
const volVectorField& U,
const IOdictionary&
);
// Member Functions
scalar cGamma() const
{
return cGamma_;
}
const dimensionedScalar& deltaN() const
{
return deltaN_;
}
const surfaceScalarField& nHatf() const
{
return nHatf_;
}
const volScalarField& K() const
{
return K_;
}
const dimensionedScalar& sigma() const
{
return sigma_;
}
tmp<volScalarField> sigmaK() const
{
return sigma_*K_;
}
void correct()
{
calculateK();
}
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************************************************************* //