radialActuationDiskSource.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-2016 OpenFOAM Foundation
9 Copyright (C) 2020 OpenCFD Ltd.
10-------------------------------------------------------------------------------
11License
12 This file is part of OpenFOAM.
13
14 OpenFOAM is free software: you can redistribute it and/or modify it
15 under the terms of the GNU General Public License as published by
16 the Free Software Foundation, either version 3 of the License, or
17 (at your option) any later version.
18
19 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
20 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
21 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
22 for more details.
23
24 You should have received a copy of the GNU General Public License
25 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
26
27\*---------------------------------------------------------------------------*/
28
30#include "geometricOneField.H"
32
33// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
34
35namespace Foam
36{
37namespace fv
38{
41}
42}
43
44// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
45
47(
48 const word& name,
49 const word& modelType,
50 const dictionary& dict,
51 const fvMesh& mesh
52)
53:
54 actuationDiskSource(name, modelType, dict, mesh),
55 radialCoeffs_(coeffs_.lookup("coeffs"))
56{
57 Info<< " - creating radial actuation disk zone: " << name_ << endl;
58}
59
60
61// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
62
64(
66 const label fieldi
67)
68{
69 const scalarField& cellsV = mesh_.V();
70 vectorField& Usource = eqn.source();
71 const vectorField& U = eqn.psi();
72
73 if (V_ > VSMALL)
74 {
75 addRadialActuationDiskAxialInertialResistance
76 (
77 Usource,
78 cells_,
79 cellsV,
81 U
82 );
83 }
84}
85
86
88(
89 const volScalarField& rho,
91 const label fieldi
92)
93{
94 const scalarField& cellsV = mesh_.V();
95 vectorField& Usource = eqn.source();
96 const vectorField& U = eqn.psi();
97
98 if (V_ > VSMALL)
99 {
100 addRadialActuationDiskAxialInertialResistance
101 (
102 Usource,
103 cells_,
104 cellsV,
105 rho,
106 U
107 );
108 }
109}
110
111
113{
115 {
116 coeffs_.readEntry("coeffs", radialCoeffs_);
117
118 return true;
119 }
120
121 return false;
122}
123
124
125// ************************************************************************* //
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
virtual bool read()
Re-read model coefficients if they have changed.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:126
A special matrix type and solver, designed for finite volume solutions of scalar equations....
Definition: fvMatrix.H:121
const GeometricField< Type, fvPatchField, volMesh > & psi(const label i=0) const
Return psi.
Definition: fvMatrix.H:412
Field< Type > & source() noexcept
Definition: fvMatrix.H:458
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:91
Applies sources on velocity (i.e. U) within a specified region to enable actuator disk models for aer...
Base abstract class for handling finite volume options (i.e. fvOption).
Definition: fvOption.H:127
const word name_
Source name.
Definition: fvOption.H:133
Applies sources on velocity field (i.e. U) within a specified region to enable actuator disk models w...
virtual void addSup(fvMatrix< vector > &eqn, const label fieldi)
Add explicit contribution to momentum equation.
A class representing the concept of a GeometricField of 1 used to avoid unnecessary manipulations for...
Lookup type of boundary radiation properties.
Definition: lookup.H:66
A class for handling words, derived from Foam::string.
Definition: word.H:68
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
Definition: className.H:121
U
Definition: pEqn.H:72
dynamicFvMesh & mesh
Namespace for OpenFOAM.
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
Definition: exprTraits.C:59
labelList fv(nPoints)
dictionary dict