NusseltNumber.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) 2021 OpenCFD Ltd.
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
26\*---------------------------------------------------------------------------*/
27
28#include "NusseltNumber.H"
29#include "ThermoCloud.H"
30
31// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
32
33template<class CloudType>
35(
36 const dictionary& dict,
37 CloudType& owner,
38 const word& modelName
39)
40:
41 CloudFunctionObject<CloudType>(dict, owner, modelName, typeName)
42{}
43
44
45template<class CloudType>
47(
49)
50:
52{}
53
54
55// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
56
57template<class CloudType>
59(
60 const typename parcelType::trackingData& td
61)
62{
63 auto& c = this->owner();
64 const auto& tc =
65 static_cast<const ThermoCloud<KinematicCloud<Cloud<parcelType>>>&>(c);
66
67 if (!c.template foundObject<IOField<scalar>>("Nu"))
68 {
69 auto* NuPtr =
71 (
73 (
74 "Nu",
75 c.time().timeName(),
76 c,
78 )
79 );
80
81 NuPtr->store();
82 }
83
84 auto& Nu = c.template lookupObjectRef<IOField<scalar>>("Nu");
85 Nu.setSize(c.size());
86
87 const auto& heatTransfer = tc.heatTransfer();
88 typename parcelType::trackingData& nctd =
89 const_cast<typename parcelType::trackingData&>(td);
90
91 label parceli = 0;
92 forAllConstIters(c, parcelIter)
93 {
94 const parcelType& p = parcelIter();
95
96 scalar Ts, rhos, mus, Pr, kappas;
97 p.template calcSurfaceValues<CloudType>
98 (
99 c, nctd, p.T(), Ts, rhos, mus, Pr, kappas
100 );
101 const scalar Re = p.Re(rhos, p.U(), td.Uc(), p.d(), mus);
102
103 Nu[parceli++] = heatTransfer.Nu(Re, Pr);
104 }
105
106
107 if (c.size() && c.time().writeTime())
108 {
109 Nu.write();
110 }
111}
112
113
114// ************************************************************************* //
Templated cloud function object base class.
Templated base class for dsmc cloud.
Definition: DSMCCloud.H:75
tmp< GeometricField< Type, PatchField, GeoMesh > > T() const
Return transpose (only if it is a tensor field)
A primitive field of type <T> with automated input and output.
Definition: IOField.H:58
Defines the attributes of an object for which implicit objectRegistry management is supported,...
Definition: IOobject.H:170
virtual void postEvolve()
Post-evolve hook.
Calculates and writes particle Nusselt number field on the cloud.
Templated base class for thermodynamic cloud.
Definition: ThermoCloud.H:70
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:126
Class used to pass data into container.
virtual bool write(const bool valid=true) const
Write using setting from DB.
A class for handling words, derived from Foam::string.
Definition: word.H:68
volScalarField & p
scalarField Re(const UList< complex > &cf)
Extract real component.
Definition: complexField.C:159
volScalarField & nu
dictionary dict
dimensionedScalar Pr("Pr", dimless, laminarTransport)
#define forAllConstIters(container, iter)
Iterate across all elements of the container object with const access.
Definition: stdFoam.H:278