solarHeatLoad.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | Copyright (C) 2015 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::radiation::solarHeatLoad
 
Description
 
    The solar load radiation model includes Sun primary hits, their
    reflective fluxes and diffusive sky radiative fluxes.
 
    The primary hit rays are calculated using a face shading algorithm.
    The reflected fluxes are considered diffusive and use a view factors method
    to deposit the energy on "visible" walls. The sky diffusive radiation for
    horizontal and vertical walls is calculated following the Fair Weather
    Conditions Method from the ASHRAE Handbook.
 
    By default the energy is included in cells adjacent to the patches into
    the energy Equation (wallCoupled = false). On coupled patches the flux is
    by default added to the wall and considered into the solid
    (solidCoupled = true).
 
    The reflected fluxes uses a grey absoprtion/emission model wich is weighted
    by the spectral distribution. The flag useVFbeamToDiffuse should be
    switched on and the view factors should be calculated using the
    'viewFactorsGen' application.
 
    The solarHeatLoad model can be used in conjuntion with fvDOM and viewFactor
    radiation models but only using a single band spectrum. On the
    corresponding BC's for these models the flag "solarHeatLoad" must be set to
    true.
 
 
SourceFiles
    solarHeatLoad.C
 
\*---------------------------------------------------------------------------*/
 
#ifndef radiationModelsolarHeatLoad_H
#define radiationModelsolarHeatLoad_H
 
#include "radiationModel.H"
#include "singleCellFvMesh.H"
#include "scalarListIOList.H"
#include "volFields.H"
#include "faceHeatShading.H"
#include "solarCalculator.H"
#include "IOmapDistribute.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
namespace radiation
{
 
/*---------------------------------------------------------------------------*\
                           Class solarHeatLoad Declaration
\*---------------------------------------------------------------------------*/
 
class solarHeatLoad
:
    public radiationModel
{
protected:
 
    // Static data
 
        //- Static name for view factor walls
        static const word viewFactorWalls;
 
 
private:
    // Private data
 
 
        //- Agglomeration List
        labelListIOList finalAgglom_;
 
        //- Coarse mesh
        autoPtr<singleCellFvMesh> coarseMesh_;
 
        //- direct solar radiative heat flux [W/m2] 净辐射量,直射辐射
        volScalarField Qr_;
 
        //- reflect solar radiative heat flux [W/m2] 面面之间的辐射量,中间量，不用保存
        volScalarField QrefecFormOtherWallRad_;
 
        //- diffuse solar radiative heat flux [W/m2] 散射辐射量
        volScalarField QdiffRad_;
 
        //- 总辐射量total solar radiative heat flux [W/m2] 总的辐射量
        volScalarField QtotalRad_;
 
        //- Direct hit faces Ids
        autoPtr<faceHeatShading> hitFaces_;
 
        //- Constant source term
        DimensionedField<scalar, volMesh> Ru_;
 
        //- Solar calculator
        solarCalculator solarCalc_;
 
        //- Vertical direction (Default is g vector)
        vector verticalDir_;
 
        //- Include diffuse reflected heat fluxes from direct heat flux
        bool useVFbeamToDiffuse_;
 
        //- Selected patches to apply solar load
        labelList includePatches_;
 
        //- Chached coarse to fine mapping for coarse mesh
        List<labelListList> coarseToFine_;
 
        //-Number of bands
        label nBands_;
 
        //- Spectral distribution for the integrated solar heat flux
        scalarList spectralDistribution_;
 
        //- Map distribute
        autoPtr<IOmapDistribute> map_;
 
        //- Face-compact map
        labelListIOList visibleFaceFaces_;
 
        //- Couple solids through mapped boundary patch using Qr (default:true)
        bool solidCoupled_;
 
        //- Couple wall patches using Qr (default:false)
        bool wallCoupled_;
 
        //- Absorptivity list
        List<List<tmp<scalarField> > > absorptivity_;
 
        //- Update absorptivity
        bool updateAbsorptivity_;
 
        //- First iteration
        bool firstIter_;
 
        //- Update Sun position index
        label updateTimeIndex_;
 
 
    // Private Member Functions
 
 
        //- Initialise
        void initialise(const dictionary&);
 
        //- Update direct hit faces radiation
        void updateDirectHitRadiation(const labelList&, const labelHashSet&);
 
        //- Calculate diffusive heat flux
        void calculateQdiff(const labelHashSet&, const labelHashSet&);
 
        //- Update Sky diffusive radiation
        void updateSkyDiffusiveRadiation
        (
            const labelHashSet&,
            const labelHashSet&
        );
 
        //- Update hit faces
        bool updateHitFaces();
 
        //- Update absorptivity
        void updateAbsorptivity(const labelHashSet& includePatches);
 
        //- Disallow default bitwise copy construct
        solarHeatLoad(const solarHeatLoad&);
 
        //- Disallow default bitwise assignment
        void operator=(const solarHeatLoad&);
 
 
public:
 
    //- Runtime type information
    TypeName("solarHeatLoad");
 
 
    // Constructors
 
        //- Construct from volScalarField
        solarHeatLoad(const volScalarField& T);
 
        //- Construct from dictionary and volScalarField
        solarHeatLoad(const dictionary& dict, const volScalarField& T);
 
        //- Constructor from local components. Does not create a radiationModel.
        //  radWallFieldName is the solar heat field name
        solarHeatLoad
        (
            const dictionary& dict,
            const volScalarField& T,
            const word radWallFieldName
        );
 
 
    //- Destructor
    virtual ~solarHeatLoad();
 
 
    // Member functions
 
        // Edit
 
            //- Solve
            void calculate();
 
            //- Read radiation properties dictionary
            bool read();
 
            //- Source term component (for power of T^4)
            virtual tmp<volScalarField> Rp() const;
 
            //- Source term component (constant)
            virtual tmp<DimensionedField<scalar, volMesh> > Ru() const;
 
 
    // Access
 
        //- Number of bands
        label nBands() const;
};
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace radiation
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //