turbulentTemperatureRadFixedGradientFvPatchScalarField.C

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | Copyright (C) 2011-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/>.
 
\*---------------------------------------------------------------------------*/
 
#include "turbulentTemperatureRadFixedGradientFvPatchScalarField.H"
#include "addToRunTimeSelectionTable.H"
#include "fvPatchFieldMapper.H"
#include "volFields.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
namespace incompressible
{
// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
turbulentTemperatureRadFixedGradientFvPatchScalarField::
turbulentTemperatureRadFixedGradientFvPatchScalarField
(
    const fvPatch& p,
    const DimensionedField<scalar, volMesh>& iF
)
:
    fixedGradientFvPatchScalarField(p, iF),
    //QrName_("undefined"),
    Qs(p.size(),0.0),
    Qv(p.size(),0.0)
{
}
 
 
turbulentTemperatureRadFixedGradientFvPatchScalarField::
turbulentTemperatureRadFixedGradientFvPatchScalarField
(
    const turbulentTemperatureRadFixedGradientFvPatchScalarField& psf,
    const fvPatch& p,
    const DimensionedField<scalar, volMesh>& iF,
    const fvPatchFieldMapper& mapper
)
:
    fixedGradientFvPatchScalarField(psf, p, iF, mapper),
    //QrName_("undefined"),
    Qs(psf.Qs,mapper),
    Qv(psf.Qv,mapper)
{}
 
 
turbulentTemperatureRadFixedGradientFvPatchScalarField::
turbulentTemperatureRadFixedGradientFvPatchScalarField
(
    const fvPatch& p,
    const DimensionedField<scalar, volMesh>& iF,
    const dictionary& dict
)
:
    fixedGradientFvPatchScalarField(p, iF),
    //QrName_(dict.lookupOrDefault<word>("Qr", "Qr")),
    //Qs(dict.lookupOrDefault<scalarField>("Qs",scalarField(p.size(),0.0))),
    //Qv(dict.lookupOrDefault<scalarField>("Qv",scalarField(p.size(),0.0)))
    Qs("Qs",dict,p.size()),
    Qv("Qv",dict,p.size())
{
 
 
    if (debug)
    {
        Info << "initialeze the heat boundary conditions"<<endl;
    }
 
    if (dict.found("value"))
    {
 
        if (debug)
        {
            Info << "find value" <<endl;
        }
 
    //操作符 = 来给 value 赋值
 
        fvPatchScalarField::operator= 
        (
            scalarField("value", dict, p.size())
        );
    }
    else
    {
        if (debug)
        {
            Info << "value can not found,"<<endl;
        }
 
        fvPatchScalarField::operator=
        (
            0   
        );
    }
 
    if (dict.found("gradient"))
    {
        if (debug)
        {
            Info << "find gradient" <<endl;
        }
 
        gradient() = scalarField("gradient", dict, p.size());
        fixedGradientFvPatchScalarField::updateCoeffs();
        fixedGradientFvPatchScalarField::evaluate();
    }
    else
    {
 
        if (debug)
        {
            Info << "gardient can not found,"<< endl;
            Info << "the value is equal to the internal value near cell" << endl;
            Info << "gradient default is 0" << endl;
        }
        fvPatchField<scalar>::operator=(patchInternalField());
        gradient() = 0.0;
 
    }
 
 
    fvPatchScalarField::operator=(scalarField("value", dict, p.size()));
 
}
 
 
turbulentTemperatureRadFixedGradientFvPatchScalarField::
turbulentTemperatureRadFixedGradientFvPatchScalarField
(
    const turbulentTemperatureRadFixedGradientFvPatchScalarField& psf,
    const DimensionedField<scalar, volMesh>& iF
)
:
 
    fixedGradientFvPatchScalarField(psf, iF),
    //QrName_("undefined"),
    Qs(psf.Qs),
    Qv(psf.Qv)
{}
 
 
// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
void turbulentTemperatureRadFixedGradientFvPatchScalarField::updateCoeffs()
{
    // Since we're inside initEvaluate/evaluate there might be processor
    // comms underway. Change the tag we use.
    int oldTag = UPstream::msgType();
    UPstream::msgType() = oldTag+1;
 
    if (debug)
    {
        Info << "update the boundary patch field Coeffs now ... " << endl;
    }
 
    if (updated())
    {
        return;
    }
 
    const dictionary& transportProperties
        = db().lookupObject<IOdictionary>("transportProperties");
 
    dimensionedScalar Cp0(transportProperties.lookup("Cp"));
 
    dimensionedScalar rho0(transportProperties.lookup("rho"));
 
    const fvPatchField<scalar>& alphaEffW =
       patch().lookupPatchField<volScalarField, scalar>("alphaEff");
 
    const fvPatchField<scalar>& QtotalFromRadiation =
       patch().lookupPatchField<volScalarField, scalar>("QtotalRad_");
 
//    const fvPatchField<scalar>& alphaEffW =
//        patch().lookupPatchField<volScalarField, scalar>("kappaEff");
 
//    Info << "cp= "<< Cp0 << endl;
//    Info << "rho0= "<< rho0 << endl;
    //Info << "alphaEffw=" << alphaEffW << endl;
//    Info << "heatflux= " << heatFlux_ << endl;
 
            const scalar Ap = gSum(patch().magSf());
           // gradient() = q_/(Ap*rhoCp0*alphaEffp);
           // Info << "Area ="<< Ap;
           // Info << "Cp ="<< Cp0.value();
           // Info << "rho ="<< rho0.value();
           // Info << "alpha ="<< alphaEffW;
    //gradient() = heatFlux_ / (alphaEffW * Cp0.value() * rho0.value() );
 
 
    //scalarField& p = *this;
  //  scalarField Qr(p.size(), 0.0);
  //  if (QrName_ != "none" || QrName_ != "undefined")
  //  {
  //    Info << "get the " << QrName_ << " value now....." << endl;
 
  //      Qr = patch().lookupPatchField<volScalarField, scalar>(QrName_);
  //  }
    //gradient() = (QtotalFromRadiation*Ap+Qs*Ap+Qv) / (Ap*alphaEffW * Cp0.value() * rho0.value() ); //热辐射W+单位面积传热量W/m2*面积+直接发热量W
    gradient() = (Qs*Ap+Qv) / (Ap*alphaEffW * Cp0.value() * rho0.value() ); //热辐射W+单位面积传热量W/m2*面积+直接发热量W
    //if (debug)
    {
        Info << "input radiation heat flux is : " << QtotalFromRadiation << endl;
        Info << "input Qs heat flux is : " << Qs << "W/m2" << endl;
        Info << "input Qv heat flux is : " << Qv << "W" << endl;
        Info << "input Ap is : " << Ap << "m2" <<endl;
        Info << "input alphaEffw is : " << alphaEffW << "m2/s"<<endl;
//      Info << "input alpha is : " << alpha << "m2/s"<<endl;
//      Info << "input alphat is : " << alpht << "m2/s"<<endl;
        Info << "input Cp0 is : " << Cp0 << endl;
        Info << "input rho0 is : " << rho0 << endl;
        Info << "the gradient of temperature  is : " << gradient() << endl;
    }
 
    fixedGradientFvPatchScalarField::updateCoeffs();
//Info << "update coeffs" << endl;
 
 
 
 
 
 
    if (updated())
    {
        return;
    }
 
 
    // Restore tag
    UPstream::msgType() = oldTag;
}
 
 
void turbulentTemperatureRadFixedGradientFvPatchScalarField::write
(
    Ostream& os
) const
{
  //  os.writeKeyword("QrName")<< QrName_ << token::END_STATEMENT << nl;
    fvPatchScalarField::write(os);
    Qs.writeEntry("Qs", os);
    Qv.writeEntry("Qv", os);
    gradient().writeEntry("gradient", os);
    this->writeEntry("value", os);
}
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
makePatchTypeField
(
    fvPatchScalarField,
    turbulentTemperatureRadFixedGradientFvPatchScalarField
);
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} //End incompressible
} // End namespace Foam