From 258091faee6bdc7d8037d0b552658db24b2ec541 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Fri, 13 Jan 2017 14:11:28 +0000
Subject: [PATCH] species::thermo<Thermo, Type>::T: Check for negative initial
temperature
Starting from an negative initial temperature causes non-convergence and a
misleading error, now a specific message is generated.
---
.../specie/thermo/thermo/thermoI.H | 11 +++++++++--
1 file changed, 9 insertions(+), 2 deletions(-)
diff --git a/src/thermophysicalModels/specie/thermo/thermo/thermoI.H b/src/thermophysicalModels/specie/thermo/thermo/thermoI.H
index 44197fd18f..669d7e7a21 100644
--- a/src/thermophysicalModels/specie/thermo/thermo/thermoI.H
+++ b/src/thermophysicalModels/specie/thermo/thermo/thermoI.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -49,6 +49,13 @@ inline Foam::scalar Foam::species::thermo<Thermo, Type>::T
scalar (thermo<Thermo, Type>::*limit)(const scalar) const
) const
{
+ if (T0 < 0)
+ {
+ FatalErrorInFunction
+ << "Negative initial temperature T0: " << T0
+ << abort(FatalError);
+ }
+
scalar Test = T0;
scalar Tnew = T0;
scalar Ttol = T0*tol_;
@@ -64,7 +71,7 @@ inline Foam::scalar Foam::species::thermo<Thermo, Type>::T
if (iter++ > maxIter_)
{
FatalErrorInFunction
- << "Maximum number of iterations exceeded"
+ << "Maximum number of iterations exceeded: " << maxIter_
<< abort(FatalError);
}
--
GitLab