valueMulticomponentMixture.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Copyright (C) 2011-2022 OpenFOAM Foundation
     \\/     M anipulation  |
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.
 
    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.
 
    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
Class
    Foam::valueMulticomponentMixture
 
Description
    Thermophysical properties mixing class which applies mass-fraction weighted
    mixing to thermodynamic properties and mole-fraction weighted mixing to
    transport properties.
 
SourceFiles
    valueMulticomponentMixture.C
 
\*---------------------------------------------------------------------------*/
 
#ifndef valueMulticomponentMixture_H
#define valueMulticomponentMixture_H
 
#include "multicomponentMixture.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
/*---------------------------------------------------------------------------*\
                 Class valueMulticomponentMixture Declaration
\*---------------------------------------------------------------------------*/
 
template<class ThermoType>
class valueMulticomponentMixture
:
    public multicomponentMixture<ThermoType>
{
 
public:
 
    class thermoMixture
    {
        //- List of specie thermo
        const PtrList<ThermoType>& specieThermos_;
 
        //- List of mass fractions
        mutable List<scalar> Y_;
 
        template<class Method, class ... Args>
        scalar massWeighted(Method psiMethod, const Args& ... args) const;
 
        template<class Method, class ... Args>
        scalar harmonicMassWeighted
        (
            Method psiMethod,
            const Args& ... args
        ) const;
 
        scalar limit(const scalar T) const;
 
 
    public:
 
        friend class valueMulticomponentMixture;
 
        // Constructors
 
            thermoMixture
            (
                const PtrList<ThermoType>& specieThermos
            )
            :
                specieThermos_(specieThermos),
                Y_(specieThermos.size())
            {}
 
 
        // Fundamental properties
 
            //- Molecular weight [kg/kmol]
            scalar W() const;
 
            //- Return density [kg/m^3]
            scalar rho(scalar p, scalar T) const;
 
            //- Return compressibility [s^2/m^2]
            scalar psi(scalar p, scalar T) const;
 
            // Heat capacity at constant pressure [J/kg/K]
            scalar Cp(const scalar p, const scalar T) const;
 
            // Heat capacity at constant volume [J/kg/K]
            scalar Cv(const scalar p, const scalar T) const;
 
            // Sensible enthalpy [J/kg]
            scalar Hs(const scalar p, const scalar T) const;
 
            // Absolute enthalpy [J/kg]
            scalar Ha(const scalar p, const scalar T) const;
 
            // Enthalpy of formation [J/kg]
            scalar Hf() const;
 
 
        // Mass specific derived properties
 
            //- Heat capacity at constant pressure/volume [J/kg/K]
            scalar Cpv(const scalar p, const scalar T) const;
 
            //- Gamma = Cp/Cv []
            scalar gamma(const scalar p, const scalar T) const;
 
            //- Enthalpy/Internal energy [J/kg]
            scalar HE(const scalar p, const scalar T) const;
 
 
        // Energy->temperature  inversion functions
 
            //- Temperature from enthalpy or internal energy
            //  given an initial temperature T0
            scalar THE
            (
                const scalar he,
                const scalar p,
                const scalar T0
            ) const;
    };
 
 
    class transportMixture
    {
        //- List of specie thermo
        const PtrList<ThermoType>& specieThermos_;
 
        //- List of mole fractions
        mutable List<scalar> X_;
 
        template<class Method, class ... Args>
        scalar moleWeighted(Method psiMethod, const Args& ... args) const;
 
 
    public:
 
        friend class valueMulticomponentMixture;
 
        transportMixture
        (
            const PtrList<ThermoType>& specieThermos
        )
        :
            specieThermos_(specieThermos),
            X_(specieThermos.size())
        {}
 
 
        // Transport properties
 
            //- Dynamic viscosity [kg/m/s]
            scalar mu(const scalar p, const scalar T) const;
 
            //- Thermal conductivity [W/m/K]
            scalar kappa(const scalar p, const scalar T) const;
    };
 
 
    //- Mixing type for thermodynamic properties
    typedef thermoMixture thermoMixtureType;
 
    //- Mixing type for transport properties
    typedef transportMixture transportMixtureType;
 
 
private:
 
    // Private Data
 
        //- Mutable storage for the cell/face mixture thermo data
        mutable thermoMixtureType thermoMixture_;
 
        //- Mutable storage for the cell/face mixture transport data
        mutable transportMixtureType transportMixture_;
 
 
public:
 
    // Constructors
 
        //- Construct from dictionary, mesh and phase name
        valueMulticomponentMixture
        (
            const dictionary&,
            const fvMesh&,
            const word&
        );
 
        //- Disallow default bitwise copy construction
        valueMulticomponentMixture
        (
            const valueMulticomponentMixture<ThermoType>&
        ) = delete;
 
 
    //- Destructor
    virtual ~valueMulticomponentMixture()
    {}
 
 
    // Member Functions
 
        //- Return the instantiated type name
        static word typeName()
        {
            return
                "valueMulticomponentMixture<" + ThermoType::typeName() + '>';
        }
 
        const thermoMixtureType& cellThermoMixture(const label celli) const;
 
        const thermoMixtureType& patchFaceThermoMixture
        (
            const label patchi,
            const label facei
        ) const;
 
        const transportMixtureType& cellTransportMixture
        (
            const label celli
        ) const;
 
        const transportMixtureType& patchFaceTransportMixture
        (
            const label patchi,
            const label facei
        ) const;
 
        const transportMixtureType& cellTransportMixture
        (
            const label celli,
            const thermoMixtureType& thermoMixture
        ) const
        {
            return cellTransportMixture(celli);
        }
 
        const transportMixtureType& patchFaceTransportMixture
        (
            const label patchi,
            const label facei,
            const thermoMixtureType& thermoMixture
        ) const
        {
            return patchFaceTransportMixture(patchi, facei);
        }
};
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#ifdef NoRepository
    #include "valueMulticomponentMixture.C"
#endif
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //