FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

homogeneousMixture.C

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software: you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by
00013     the Free Software Foundation, either version 3 of the License, or
00014     (at your option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
00023 
00024 \*---------------------------------------------------------------------------*/
00025 
00026 #include "homogeneousMixture.H"
00027 #include <finiteVolume/fvMesh.H>
00028 
00029 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
00030 
00031 template<class ThermoType>
00032 const char* Foam::homogeneousMixture<ThermoType>::specieNames_[1] = {"bprog"};
00033 
00034 
00035 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
00036 
00037 template<class ThermoType>
00038 Foam::homogeneousMixture<ThermoType>::homogeneousMixture
00039 (
00040     const dictionary& thermoDict,
00041     const fvMesh& mesh
00042 )
00043 :
00044     basicMultiComponentMixture
00045     (
00046         thermoDict,
00047         speciesTable(nSpecies_, specieNames_),
00048         mesh
00049     ),
00050 
00051     reactants_(thermoDict.lookup("reactants")),
00052     products_(thermoDict.lookup("products")),
00053     mixture_("mixture", reactants_),
00054     b_(Y("bprog"))
00055 {}
00056 
00057 
00058 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
00059 
00060 template<class ThermoType>
00061 const ThermoType& Foam::homogeneousMixture<ThermoType>::mixture
00062 (
00063     const scalar b
00064 ) const
00065 {
00066     if (b > 0.999)
00067     {
00068         return reactants_;
00069     }
00070     else if (b < 0.001)
00071     {
00072         return products_;
00073     }
00074     else
00075     {
00076         mixture_ = b/reactants_.W()*reactants_;
00077         mixture_ += (1 - b)/products_.W()*products_;
00078 
00079         return mixture_;
00080     }
00081 }
00082 
00083 
00084 template<class ThermoType>
00085 void Foam::homogeneousMixture<ThermoType>::read(const dictionary& thermoDict)
00086 {
00087     reactants_ = ThermoType(thermoDict.lookup("reactants"));
00088     products_ = ThermoType(thermoDict.lookup("products"));
00089 }
00090 
00091 
00092 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines