Commit c51fc017 authored by sergio's avatar sergio Committed by Andrew Heather
Browse files

ENH: exponential: new nucleate flux model

parent 1f067169
...@@ -231,6 +231,7 @@ $(TDNBModels)/Schroeder/Schroeder.C ...@@ -231,6 +231,7 @@ $(TDNBModels)/Schroeder/Schroeder.C
nucleateFluxModels = $(wallBoilingSubModels)/nucleateFluxModels nucleateFluxModels = $(wallBoilingSubModels)/nucleateFluxModels
$(nucleateFluxModels)/nucleateFluxModel/nucleateFluxModel.C $(nucleateFluxModels)/nucleateFluxModel/nucleateFluxModel.C
$(nucleateFluxModels)/Kutadeladze/Kutadeladze.C $(nucleateFluxModels)/Kutadeladze/Kutadeladze.C
$(nucleateFluxModels)/exponential/exponential.C
/* Turbulence */ /* Turbulence */
turbulence/multiphaseCompressibleTurbulenceModels.C turbulence/multiphaseCompressibleTurbulenceModels.C
......
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2021 OpenCFD Ltd
-------------------------------------------------------------------------------
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 "exponential.H"
#include "addToRunTimeSelectionTable.H"
#include "uniformDimensionedFields.H"
#include "phasePairKey.H"
#include "phaseSystem.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
namespace Foam
{
namespace wallBoilingModels
{
namespace nucleateFluxModels
{
defineTypeNameAndDebug(exponential, 0);
addToRunTimeSelectionTable
(
nucleateFluxModel,
exponential,
dictionary
);
}
}
}
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
Foam::wallBoilingModels::nucleateFluxModels::exponential::exponential
(
const dictionary& dict
)
:
nucleateFluxModel(),
a_(dict.getOrDefault<scalar>("a", 6309)),
b_(dict.getOrDefault<scalar>("b", 2.52))
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
Foam::tmp<Foam::scalarField>
Foam::wallBoilingModels::nucleateFluxModels::exponential::qNucleate
(
const phaseModel& liquid,
const phaseModel& vapor,
const label patchi,
const scalarField& Tl,
const scalarField& Tsatw,
const scalarField& L
) const
{
const fvPatchScalarField& Tw =
liquid.thermo().T().boundaryField()[patchi];
return a_*pow(max((Tw-Tsatw), scalar(0)), b_);
}
void Foam::wallBoilingModels::nucleateFluxModels::exponential::write
(
Ostream& os
) const
{
nucleateFluxModel::write(os);
os.writeEntry("a", a_);
os.writeEntry("b", b_);
}
// ************************************************************************* //
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2021 OpenCFD Ltd
-------------------------------------------------------------------------------
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::wallBoilingModels::nucleateFluxModels::exponential
Description
Nucleate flux sub-cooling correlation
References:
\verbatim
Frost W. & Dzakowic G. S. (1967)
An extension of the methods of predicting incipient
boiling on commercially finished surfaces.
ASME AIChE Heat Transfer Conf 67-HT-61, Seattle.
Shirai, Y., Tatsumoto, H., Shiotsu, M., Hata, K.,
Kobayashi, H., Naruo, Y., & Inatani, Y. (2010).
Boiling heat transfer from a horizontal flat
plate in a pool of liquid hydrogen.
Cryogenics, 50(6-7), 410-416.
DOI:10.1016/j.cryogenics.2010.04.001
\endverbatim
Usage
Example of the model specification:
\verbatim
nucleateFluxModel
{
// Mandatory entries
type exponential;
// Optional entries
a <scalar>;
b <scalar>;
}
\endverbatim
where the entries mean:
\table
Property | Description | Type | Reqd | Deflt
type | Type name: exponential | word | yes | -
a | Pre-factor coefficient | scalar | no | 6309
b | Exponent coefficient | scalar | no | 2.52
\endtable
SourceFiles
exponential.C
\*---------------------------------------------------------------------------*/
#ifndef exponential_H
#define exponential_H
#include "nucleateFluxModel.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
namespace wallBoilingModels
{
namespace nucleateFluxModels
{
/*---------------------------------------------------------------------------*\
Class exponential Declaration
\*---------------------------------------------------------------------------*/
class exponential
:
public nucleateFluxModel
{
// Private Data
//- Pre-factor coefficient
scalar a_;
//- Exponent coefficient
scalar b_;
// Private Member Functions
//- No copy construct
exponential(const exponential&) = delete;
//- No copy assignment
void operator=(const exponential&) = delete;
public:
//- Runtime type information
TypeName("exponential");
// Constructors
//- Construct from a dictionary
exponential(const dictionary& dict);
//- Destructor
virtual ~exponential() = default;
// Member Functions
//- Calculate and return the nucleation-site density
virtual tmp<scalarField> qNucleate
(
const phaseModel& liquid,
const phaseModel& vapor,
const label patchi,
const scalarField& Tl,
const scalarField& Tsatw,
const scalarField& L
) const;
//- Write
virtual void write(Ostream& os) const;
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace nucleateFluxModels
} // End namespace wallBoilingModels
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************************************************************* //
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment