faceContactSphere.C
Go to the documentation of this file.
1/*---------------------------------------------------------------------------*\
2 ========= |
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4 \\ / O peration |
5 \\ / A nd | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2011 OpenFOAM Foundation
9-------------------------------------------------------------------------------
10License
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
14 under the terms of the GNU General Public License as published by
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26Description
27 Return location of contact sphere on the face
28
29\*---------------------------------------------------------------------------*/
30
31#include "face.H"
32
33// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34
36(
37 const point& p,
38 const vector& n,
39 const UList<point>& meshPoints
40) const
41{
42 vector n1 = n/(Foam::mag(n) + SMALL);
43 vector n2 = areaNormal(meshPoints);
44
45 n2 /= Foam::mag(n2) + SMALL;
46
47 return 2*((centre(meshPoints) - p) & n2)/((n1 & n2) - 1.0);
48}
49
50
51// ************************************************************************* //
label n
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
Definition: UList.H:94
point centre(const UList< point > &points) const
Centre point of face.
Definition: face.C:514
scalar contactSphereDiameter(const point &p, const vector &n, const UList< point > &meshPoints) const
Return contact sphere diameter.
vector areaNormal(const UList< point > &p) const
The area normal - with magnitude equal to area of face.
Definition: face.C:578
volScalarField & p
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)