From 45cfddbdcd169142e081a4ecb748ce7824c264d9 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Thu, 26 Mar 2009 18:00:35 +0000
Subject: [PATCH] added maxIter

---
 .../radiation/radiationModel/fvDOM/fvDOM/fvDOM.C             | 5 +++--
 .../radiation/radiationModel/fvDOM/fvDOM/fvDOM.H             | 3 +++
 2 files changed, 6 insertions(+), 2 deletions(-)

diff --git a/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C b/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C
index d8b00e8c34..6d07cc3bea 100644
--- a/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C
+++ b/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C
@@ -128,7 +128,8 @@ Foam::radiation::fvDOM::fvDOM(const volScalarField& T)
     aLambda_(nLambda_),
     blackBody_(nLambda_, T),
     IRay_(0),
-    convergence_(coeffs_.lookupOrDefault<scalar>("convergence", 0.0))
+    convergence_(coeffs_.lookupOrDefault<scalar>("convergence", 0.0)),
+    maxIter_(coeffs_.lookupOrDefault<label>("maxIter", 50))
 {
     if (mesh_.nSolutionD() == 3)    //3D
     {
@@ -298,7 +299,7 @@ void Foam::radiation::fvDOM::calculate()
 
         Info << "Radiation solver iter: " << radIter << endl;
 
-    } while(maxResidual > convergence_);
+    } while(maxResidual > convergence_ && radIter < maxIter_);
 
     updateG();
 }
diff --git a/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.H b/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.H
index 3565cdce86..f91163805a 100644
--- a/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.H
+++ b/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.H
@@ -117,6 +117,9 @@ class fvDOM
         //- Convergence criterion
         scalar convergence_;
 
+        //- Maximum number of iterations
+        scalar maxIter_;
+
 
     // Private member functions