1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
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 the
13 Free Software Foundation; either version 2 of the License, or (at your
14 option) any later version.
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
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM; if not, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
25 \*---------------------------------------------------------------------------*/
27 #include "gradientRhoFvPatchScalarField.H"
28 #include "fvPatchFieldMapper.H"
29 #include "volFields.H"
30 #include "surfaceFields.H"
31 #include "addToRunTimeSelectionTable.H"
33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
38 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
40 gradientRhoFvPatchScalarField::gradientRhoFvPatchScalarField
43 const DimensionedField<scalar, volMesh>& iF
46 fixedGradientFvPatchScalarField(p, iF)
50 gradientRhoFvPatchScalarField::gradientRhoFvPatchScalarField
52 const gradientRhoFvPatchScalarField& ptf,
54 const DimensionedField<scalar, volMesh>& iF,
55 const fvPatchFieldMapper& mapper
58 fixedGradientFvPatchScalarField(ptf, p, iF, mapper)
62 gradientRhoFvPatchScalarField::gradientRhoFvPatchScalarField
65 const DimensionedField<scalar, volMesh>& iF,
66 const dictionary& dict
69 fixedGradientFvPatchScalarField(p, iF)
71 if (dict.found("gradient"))
73 gradient() = scalarField("gradient", dict, p.size());
74 fixedGradientFvPatchScalarField::updateCoeffs();
75 fixedGradientFvPatchScalarField::evaluate();
79 fvPatchField<scalar>::operator=(patchInternalField());
85 gradientRhoFvPatchScalarField::gradientRhoFvPatchScalarField
87 const gradientRhoFvPatchScalarField& wbppsf
90 fixedGradientFvPatchScalarField(wbppsf)
94 gradientRhoFvPatchScalarField::gradientRhoFvPatchScalarField
96 const gradientRhoFvPatchScalarField& wbppsf,
97 const DimensionedField<scalar, volMesh>& iF
100 fixedGradientFvPatchScalarField(wbppsf, iF)
104 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
106 void gradientRhoFvPatchScalarField::updateCoeffs()
113 const fvPatchField<scalar>& psip =
114 patch().lookupPatchField<volScalarField, scalar>("psi");
116 const fvPatchField<scalar>& pp =
117 patch().lookupPatchField<volScalarField, scalar>("p");
119 gradient() = psip*pp.snGrad() + psip.snGrad()*pp;
121 fixedGradientFvPatchScalarField::updateCoeffs();
125 void gradientRhoFvPatchScalarField::write(Ostream& os) const
127 fixedGradientFvPatchScalarField::write(os);
128 writeEntry("value", os);
132 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
134 makePatchTypeField(fvPatchScalarField, gradientRhoFvPatchScalarField);
136 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
138 } // End namespace Foam
140 // ************************************************************************* //