zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

limitTemperature: added support for multiphase solvers

Browse Source 
Based on patch contributed by Juho Peltola, VTT

Resolves feature request https://bugs.openfoam.org/view.php?id=2572
This commit is contained in:
Henry Weller
2017-09-04 16:52:03 +01:00
parent 7fc91b4efd
commit 0ad0d34b42
5 changed files with 20 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/fvOptions/corrections/limitTemperature/limitTemperature.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2012-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2012-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -57,13 +57,16 @@ Foam::fv::limitTemperature::limitTemperature
:
cellSetOption(name, modelType, dict, mesh),
Tmin_(readScalar(coeffs_.lookup("min"))),
Tmax_(readScalar(coeffs_.lookup("max")))
Tmax_(readScalar(coeffs_.lookup("max"))),
phase_(coeffs_.lookupOrDefault<word>("phase", word::null))
{
// Set the field name to that of the energy field from which the temperature
// is obtained
const basicThermo& thermo =
mesh_.lookupObject<basicThermo>(basicThermo::dictName);
mesh_.lookupObject<basicThermo>
(
IOobject::groupName(basicThermo::dictName, phase_)
);
fieldNames_.setSize(1, thermo.he().name());
@ -92,7 +95,10 @@ bool Foam::fv::limitTemperature::read(const dictionary& dict)
void Foam::fv::limitTemperature::correct(volScalarField& he)
{
const basicThermo& thermo =
mesh_.lookupObject<basicThermo>(basicThermo::dictName);
mesh_.lookupObject<basicThermo>
(
IOobject::groupName(basicThermo::dictName, phase_)
);
scalarField Tmin(cells_.size(), Tmin_);
scalarField Tmax(cells_.size(), Tmax_);

4
src/fvOptions/corrections/limitTemperature/limitTemperature.H
View File

@ -39,6 +39,7 @@ Usage
selectionMode all;
min 200;
max 500;
phase gas; //optional
}
\endverbatim
@ -78,6 +79,9 @@ protected:
//- Maximum temperature limit [K]
scalar Tmax_;
//- Optional phase name [K]
word phase_;
private: