SingleKineticRateDevolatilisation.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2020 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 
27 
28 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 
30 template<class CloudType>
33 (
34  const dictionary& dict,
35  CloudType& owner
36 )
37 :
38  DevolatilisationModel<CloudType>(dict, owner, typeName),
39  volatileData_(this->coeffDict().lookup("volatileData")),
40  YVolatile0_(volatileData_.size()),
41  volatileToGasMap_(volatileData_.size()),
42  residualCoeff_(this->coeffDict().template lookup<scalar>("residualCoeff"))
43 {
44  if (volatileData_.empty())
45  {
47  << "Devolatilisation model selected, but no volatiles defined"
48  << nl << endl;
49  }
50  else
51  {
52  Info<< "Participating volatile species:" << endl;
53 
54  // Determine mapping between active volatiles and cloud gas components
55  const label idGas = owner.composition().idGas();
56  const scalar YGasTot = owner.composition().YMixture0()[idGas];
57  const scalarField& YGas = owner.composition().Y0(idGas);
58  forAll(volatileData_, i)
59  {
60  const word& specieName = volatileData_[i].name();
61  const label id = owner.composition().localId(idGas, specieName);
62  volatileToGasMap_[i] = id;
63  YVolatile0_[i] = YGasTot*YGas[id];
64 
65  Info<< " " << specieName << ": particle mass fraction = "
66  << YVolatile0_[i] << endl;
67  }
68  }
69 }
70 
71 
72 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
73 
74 template<class CloudType>
77 {}
78 
79 
80 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
81 
82 template<class CloudType>
84 (
85  const scalar dt,
86  const scalar age,
87  const scalar mass0,
88  const scalar mass,
89  const scalar T,
90  const scalarField& YGasEff,
91  const scalarField& YLiquidEff,
92  const scalarField& YSolidEff,
93  label& canCombust,
94  scalarField& dMassDV
95 ) const
96 {
97  bool done = true;
98  forAll(volatileData_, i)
99  {
100  const label id = volatileToGasMap_[i];
101  const scalar massVolatile0 = mass0*YVolatile0_[i];
102  const scalar massVolatile = mass*YGasEff[id];
103 
104  // Combustion allowed once all volatile components evolved
105  done = done && (massVolatile <= residualCoeff_*massVolatile0);
106 
107  // Model coefficients
108  const scalar A1 = volatileData_[i].A1();
109  const scalar E = volatileData_[i].E();
110 
111  // Kinetic rate
112  const scalar kappa = A1*exp(-E/(RR*T));
113 
114  // Mass transferred from particle to carrier gas phase
115  dMassDV[id] = min(dt*kappa*massVolatile, massVolatile);
116  }
117 
118  if (done && canCombust != -1)
119  {
120  canCombust = 1;
121  }
122 }
123 
124 
125 // ************************************************************************* //
dictionary dict
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:156
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
const dimensionedScalar RR
Universal gas constant: default SI units: [J/kmol/K].
const dimensionedScalar kappa
Coulomb constant: default SI units: [N.m2/C2].
dimensionedScalar exp(const dimensionedScalar &ds)
A class for handling words, derived from string.
Definition: word.H:59
static const char nl
Definition: Ostream.H:260
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
#define WarningInFunction
Report a warning using Foam::Warning.
messageStream Info
Templated devolatilisation model class.
virtual void calculate(const scalar dt, const scalar age, const scalar mass0, const scalar mass, const scalar T, const scalarField &YGasEff, const scalarField &YLiquidEff, const scalarField &YSolidEff, label &canCombust, scalarField &dMassDV) const
Update model.
Templated base class for dsmc cloud.
Definition: DSMCCloud.H:75
SingleKineticRateDevolatilisation(const dictionary &dict, CloudType &owner)
Construct from dictionary.