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-2021 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 typename CloudType::parcelType& p,
86  const typename CloudType::parcelType::trackingData& td,
87  const scalar dt,
88  const scalar mass0,
89  const scalar mass,
90  const scalar T,
91  const scalarField& YGasEff,
92  const scalarField& YLiquidEff,
93  const scalarField& YSolidEff,
94  label& canCombust,
95  scalarField& dMassDV
96 ) const
97 {
98  bool done = true;
99  forAll(volatileData_, i)
100  {
101  const label id = volatileToGasMap_[i];
102  const scalar massVolatile0 = mass0*YVolatile0_[i];
103  const scalar massVolatile = mass*YGasEff[id];
104 
105  // Combustion allowed once all volatile components evolved
106  done = done && (massVolatile <= residualCoeff_*massVolatile0);
107 
108  // Model coefficients
109  const scalar A1 = volatileData_[i].A1();
110  const scalar E = volatileData_[i].E();
111 
112  // Kinetic rate
113  const scalar kappa = A1*exp(-E/(RR*T));
114 
115  // Mass transferred from particle to carrier gas phase
116  dMassDV[id] = min(dt*kappa*massVolatile, massVolatile);
117  }
118 
119  if (done && canCombust != -1)
120  {
121  canCombust = 1;
122  }
123 }
124 
125 
126 // ************************************************************************* //
dictionary dict
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
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
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Definition: DSMCCloud.H:221
static const char nl
Definition: Ostream.H:260
virtual void calculate(const typename CloudType::parcelType &p, const typename CloudType::parcelType::trackingData &td, const scalar dt, 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.
#define WarningInFunction
Report a warning using Foam::Warning.
messageStream Info
Templated devolatilisation model class.
Templated base class for dsmc cloud.
Definition: DSMCCloud.H:75
SingleKineticRateDevolatilisation(const dictionary &dict, CloudType &owner)
Construct from dictionary.