zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

TDACChemistryModel: Removed unnecessary warning messages

Browse Source 
Patch contributed by Francesco Contino
This commit is contained in:
Henry Weller
2017-12-04 12:28:59 +00:00
parent 026804c0e0
commit 745b20acb6
3 changed files with 8 additions and 20 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/thermophysicalModels/chemistryModel/chemistryModel/TDACChemistryModel/reduction/DAC/DAC.C
View File

@ -642,17 +642,15 @@ void Foam::chemistryReductionMethods::DAC<CompType, ThermoType>::reduceMechanism
{
label u = Q.pop();
scalar Den = max(PA[u],CA[u]);
if (Den!=0.0)
if (Den != 0)
{
for (label v=0; v<NbrABInit[u]; v++)
{
label otherSpec = rABOtherSpec(u, v);
scalar rAB = mag(rABNum(u, v))/Den;
if (rAB>1)
if (rAB > 1)
{
Info<< "Badly Conditioned rAB : " << rAB
<< "species involved : "<<u << "," << otherSpec << endl;
rAB=1.0;
rAB = 1;
}
// The direct link is weaker than the user-defined tolerance
if (rAB >= this->tolerance())

5
src/thermophysicalModels/chemistryModel/chemistryModel/TDACChemistryModel/reduction/DRG/DRG.C
View File

@ -253,11 +253,8 @@ void Foam::chemistryReductionMethods::DRG<CompType, ThermoType>::reduceMechanism
label otherSpec = rABOtherSpec(u, v);
scalar rAB = rABNum(u, v)/Den;
if (rAB>1)
if (rAB > 1)
{
Info<< "Badly Conditioned rAB : " << rAB
<< "species involved : " << u << "," << otherSpec
<< endl;
rAB = 1;
}

15
src/thermophysicalModels/chemistryModel/chemistryModel/TDACChemistryModel/reduction/DRGEP/DRGEP.C
View File

@ -442,11 +442,9 @@ reduceMechanism
{
label otherSpec = rABOtherSpec(u, v);
scalar rAB = mag(rABNum(u, v))/Den;
if (rAB>1)
if (rAB > 1)
{
Info<< "Badly Conditioned rAB : " << rAB
<< "species involved : "<<u << "," << otherSpec << endl;
rAB=1.0;
rAB = 1;
}
scalar Rtemp = Rvalue[u]*rAB;
@ -647,14 +645,9 @@ reduceMechanism
if (!disabledSpecies[otherSpec])
{
scalar rAB = mag(rABNum(u, v))/Den;
if (rAB>1.0)
if (rAB > 1)
{
Info<< "Badly Conditioned rAB : " << rAB
<< "species involved : "
<<this->chemistry_.Y()[u].member() << ","
<< this->chemistry_.Y()[otherSpec].member()
<< endl;
rAB=1.0;
rAB = 1;
}
scalar Rtemp = Rvalue[u]*rAB;