applications/utilities/thermophysical/mixtureAdiabaticFlameT/mixtureAdiabaticFlameT.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 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 Application
  25     mixtureAdiabaticFlameT
  26
  27 Description
  28     Calculates the adiabatic flame temperature for a given mixture
  29     at a given temperature.
  30
  31 \*---------------------------------------------------------------------------*/
  32
  33 #include "argList.H"
  34 #include "dictionary.H"
  35 #include "IFstream.H"
  36 #include "OSspecific.H"
  37
  38 #include "specieThermo.H"
  39 #include "janafThermo.H"
  40 #include "perfectGas.H"
  41 #include "mixture.H"
  42
  43 using namespace Foam;
  44
  45 typedef specieThermo<janafThermo<perfectGas> > thermo;
  46
  47
  48 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  49
  50 int main(int argc, char *argv[])
  51 {
  52     argList::validArgs.append("controlFile");
  53     argList args(argc, argv);
  54
  55     const fileName controlFileName(args[1]);
  56
  57     // Construct control dictionary
  58     IFstream controlFile(controlFileName);
  59
  60     // Check controlFile stream is OK
  61     if (!controlFile.good())
  62     {
  63         FatalErrorIn(args.executable())
  64             << "Cannot read file " << controlFileName
  65             << abort(FatalError);
  66     }
  67
  68     dictionary control(controlFile);
  69
  70
  71     scalar T0(readScalar(control.lookup("T0")));
  72     mixture rMix(control.lookup("reactants"));
  73     mixture pMix(control.lookup("products"));
  74
  75
  76     Info<< nl << "Reading Burcat data dictionary" << endl;
  77
  78     fileName BurcatCpDataFileName(findEtcFile("thermoData/BurcatCpData"));
  79
  80     // Construct control dictionary
  81     IFstream BurcatCpDataFile(BurcatCpDataFileName);
  82
  83     // Check BurcatCpData stream is OK
  84     if (!BurcatCpDataFile.good())
  85     {
  86         FatalErrorIn(args.executable())
  87             << "Cannot read file " << BurcatCpDataFileName
  88             << abort(FatalError);
  89     }
  90
  91     dictionary CpData(BurcatCpDataFile);
  92
  93
  94     thermo reactants
  95     (
  96         rMix[0].volFrac()*thermo(CpData.lookup(rMix[0].name()))
  97     );
  98
  99     for (label i = 1; i < rMix.size(); i++)
 100     {
 101         reactants = reactants
 102             + rMix[i].volFrac()*thermo(CpData.lookup(rMix[i].name()));
 103     }
 104
 105
 106     thermo products
 107     (
 108         2*pMix[0].volFrac()*thermo(CpData.lookup(pMix[0].name()))
 109     );
 110
 111     for (label i = 1; i < pMix.size(); i++)
 112     {
 113         products = products
 114             + 2*pMix[i].volFrac()*thermo(CpData.lookup(pMix[i].name()));
 115     }
 116
 117     Info<< "Adiabatic flame temperature of mixture " << rMix.name() << " = "
 118          << products.TH(reactants.H(T0), 1000.0) << " K" << endl;
 119
 120     return 0;
 121 }
 122
 123
 124 // ************************************************************************* //