PoissonPatchDistMethod.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 | Copyright (C) 2015-2016 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "PoissonPatchDistMethod.H"
27 #include "fvcGrad.H"
28 #include "fvmLaplacian.H"
30 
31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
32 
33 namespace Foam
34 {
35 namespace patchDistMethods
36 {
37  defineTypeNameAndDebug(Poisson, 0);
38  addToRunTimeSelectionTable(patchDistMethod, Poisson, dictionary);
39 }
40 }
41 
42 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
43 
44 Foam::patchDistMethods::Poisson::Poisson
45 (
46  const dictionary& dict,
47  const fvMesh& mesh,
48  const labelHashSet& patchIDs
49 )
50 :
51  patchDistMethod(mesh, patchIDs)
52 {}
53 
54 
55 Foam::patchDistMethods::Poisson::Poisson
56 (
57  const fvMesh& mesh,
58  const labelHashSet& patchIDs
59 )
60 :
61  patchDistMethod(mesh, patchIDs)
62 {}
63 
64 
65 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
66 
68 {
69  return correct(y, const_cast<volVectorField&>(volVectorField::null()));
70 }
71 
72 
74 (
75  volScalarField& y,
77 )
78 {
79  if (!tyPsi_.valid())
80  {
81  tyPsi_ = tmp<volScalarField>
82  (
83  new volScalarField
84  (
85  IOobject
86  (
87  "yPsi",
88  mesh_.time().timeName(),
89  mesh_
90  ),
91  mesh_,
92  dimensionedScalar("yPsi", sqr(dimLength), 0.0),
93  y.boundaryFieldRef().types()
94  )
95  );
96  }
97  volScalarField& yPsi = tyPsi_.ref();
98 
99  solve(fvm::laplacian(yPsi) == dimensionedScalar("1", dimless, -1.0));
100 
101  volVectorField gradyPsi(fvc::grad(yPsi));
102  volScalarField magGradyPsi(mag(gradyPsi));
103 
104  y = sqrt(magSqr(gradyPsi) + 2*yPsi) - magGradyPsi;
105 
106  // Cache yPsi if the mesh is moving otherwise delete
107  if (!mesh_.changing())
108  {
109  tyPsi_.clear();
110  }
111 
112  // Only calculate n if the field is defined
113  if (notNull(n))
114  {
115  n =
116  -gradyPsi
117  /max
118  (
119  magGradyPsi,
120  dimensionedScalar("smallMagGradyPsi", dimLength, SMALL)
121  );
122  }
123 
124  return true;
125 }
126 
127 
128 // ************************************************************************* //
virtual bool correct(volScalarField &y)
Correct the given distance-to-patch field.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
addToRunTimeSelectionTable(patchDistMethod, advectionDiffusion, dictionary)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Calculate the matrix for the laplacian of the field.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
Calculate the gradient of the given field.
void clear()
Clear the list, i.e. set size to zero.
Definition: ListI.H:124
SolverPerformance< Type > solve(fvMatrix< Type > &, const dictionary &)
Solve returning the solution statistics given convergence tolerance.
dimensioned< scalar > magSqr(const dimensioned< Type > &)
Info<< "Predicted p max-min : "<< max(p).value()<< " "<< min(p).value()<< endl;rho==max(psi *p+alphal *rhol0+((alphav *psiv+alphal *psil) - psi) *pSat, rhoMin);# 1 "/home/ubuntu/OpenFOAM-5.0/applications/solvers/multiphase/cavitatingFoam/alphavPsi.H" 1{ alphav=max(min((rho - rholSat)/(rhovSat - rholSat), scalar(1)), scalar(0));alphal=1.0 - alphav;Info<< "max-min alphav: "<< max(alphav).value()<< " "<< min(alphav).value()<< endl;psiModel-> correct()
Definition: pEqn.H:72
Internal & ref()
Return a reference to the dimensioned internal field.
bool notNull(const T &t)
Return true if t is not a reference to the nullObject of type T.
Definition: nullObjectI.H:46
static const GeometricField< vector, fvPatchField, volMesh > & null()
Return a null geometric field.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
Definition: dimensionSets.H:47
Boundary & boundaryFieldRef()
Return a reference to the boundary field.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
Definition: dimensionSets.H:50
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvmLaplacian.C:46
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
defineTypeNameAndDebug(advectionDiffusion, 0)
dimensioned< scalar > mag(const dimensioned< Type > &)
A class for managing temporary objects.
Definition: PtrList.H:53
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:92
Specialisation of patchDist for wall distance calculation.
Namespace for OpenFOAM.