Merge branch 'master' of /home/dm4/OpenFOAM/OpenFOAM-dev

src/regionModels/pyrolysisModels/noPyrolysis/noPyrolysis.C

@@ -42,9 +42,18 @@ defineTypeNameAndDebug(noPyrolysis, 0);
 addToRunTimeSelectionTable(pyrolysisModel, noPyrolysis, mesh);
 addToRunTimeSelectionTable(pyrolysisModel, noPyrolysis, dictionary);
 
-
 // * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
 
+void noPyrolysis::constructThermoChemistry()
+{
+    solidChemistry_.reset
+    (
+        solidChemistryModel::New(regionMesh()).ptr()
+    );
+
+    solidThermo_.reset(&solidChemistry_->solidThermo());
+}
+
 bool noPyrolysis::read()
 {
     if (pyrolysisModel::read())
@@ -77,10 +86,15 @@ bool noPyrolysis::read(const dictionary& dict)
 
 noPyrolysis::noPyrolysis(const word& modelType, const fvMesh& mesh)
 :
-    pyrolysisModel(modelType, mesh),
+    pyrolysisModel(mesh),
-    solidChemistry_(solidChemistryModel::New(regionMesh())),
+    solidChemistry_(NULL),
-    solidThermo_(solidChemistry_->solidThermo())
+    solidThermo_(NULL)
-{}
+{
+    if (active())
+    {
+        constructThermoChemistry();
+    }
+}
 
 
 noPyrolysis::noPyrolysis
@@ -89,11 +103,15 @@ noPyrolysis::noPyrolysis
     const fvMesh& mesh,
     const dictionary& dict
 ):
-    pyrolysisModel(modelType, mesh, dict),
+    pyrolysisModel(mesh),
-    solidChemistry_(solidChemistryModel::New(regionMesh())),
+    solidChemistry_(NULL),
-    solidThermo_(solidChemistry_->solidThermo())
+    solidThermo_(NULL)
-{}
+{
-
+    if (active())
+    {
+        constructThermoChemistry();
+    }
+}
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
@@ -114,33 +132,34 @@ void noPyrolysis::evolveRegion()
     //Do nothing
 }
 
+
 const volScalarField& noPyrolysis::rho() const
 {
-    return (solidThermo_.rho());
+    return (solidThermo_->rho());
 }
 
 
 const volScalarField& noPyrolysis::T() const
 {
-    return (solidThermo_.T());
+    return (solidThermo_->T());
 }
 
 
 const tmp<volScalarField> noPyrolysis::Cp() const
 {
-    return (solidThermo_.Cp());
+    return (solidThermo_->Cp());
 }
 
 
 const volScalarField& noPyrolysis::kappa() const
 {
-    return (solidThermo_.kappa());
+    return (solidThermo_->kappa());
 }
 
 
 const volScalarField& noPyrolysis::K() const
 {
-    return (solidThermo_.K());
+     return (solidThermo_->K());
 }
 
 
@@ -154,7 +173,7 @@ const surfaceScalarField& noPyrolysis::phiGas() const
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-} // End namespace pyrolysisFilmModels
+} // End namespace surfaceFilmModels
 } // End namespace regionModels
 } // End namespace Foam
 

src/regionModels/pyrolysisModels/noPyrolysis/noPyrolysis.H

@@ -76,11 +76,14 @@ protected:
         //- Read control parameters from dictionary
         virtual bool read(const dictionary& dict);
 
+        //- Reset solidChemistryModel and basicSolidThermo pointers
+        void constructThermoChemistry();
+
         //- Reference to the solid chemistry model
         autoPtr<solidChemistryModel> solidChemistry_;
 
         //- Reference to solid thermo
-        basicSolidThermo& solidThermo_;
+        autoPtr<basicSolidThermo> solidThermo_;
 
 
 public:
@@ -103,7 +106,6 @@ public:
         );
 
 
-
     //- Destructor
     virtual ~noPyrolysis();
 

src/regionModels/pyrolysisModels/pyrolysisModel/pyrolysisModelCollection.C

@@ -142,30 +142,16 @@ scalar pyrolysisModelCollection::maxDiff() const
 
 scalar pyrolysisModelCollection::solidRegionDiffNo() const
 {
-    scalar regionDiNum = 0.0;
+    scalar totalDiNum = GREAT;
-    scalar totalDiNum = 0.0;
 
     forAll(*this, i)
     {
-        const pyrolysisModel& pyrolysis = this->operator[](i);
+        if
-
+        (
-        if (pyrolysis.regionMesh().nInternalFaces() > 0)
+            totalDiNum > this->operator[](i).solidRegionDiffNo()
+        )
         {
-            surfaceScalarField KrhoCpbyDelta
+            totalDiNum = this->operator[](i).solidRegionDiffNo();
-            (
-                pyrolysis.regionMesh().surfaceInterpolation::deltaCoeffs()
-            * fvc::interpolate(pyrolysis.K())
-            / fvc::interpolate(pyrolysis.Cp()*pyrolysis.rho())
-            );
-
-            regionDiNum =
-                 max(KrhoCpbyDelta.internalField())
-                *pyrolysis.time().deltaTValue();
-        }
-
-        if (regionDiNum > totalDiNum)
-        {
-            totalDiNum = regionDiNum;
         }
     }
 

tutorials/combustion/fireFoam/les/smallPoolFire2D/constant/pyrolysisZones

@@ -11,27 +11,34 @@ FoamFile
     format      binary;
     class       dictionary;
     location    "constant";
-    object      pyrolysisProperties;
+    object      pyrolysisZones;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-active          false;
+1
+(
+    pyrolysis
+    {
 
-pyrolysisModel  none;
+        active          false;
 
-regionName      panelRegion;
+        pyrolysisModel  none;
 
-reactingOneDimCoeffs
+        regionName      panelRegion;
-{
-    filmCoupled     false;
 
-    radFluxName     Qr;
+        reactingOneDimCoeffs
+        {
+            filmCoupled     false;
 
-    minimumDelta    1e-8;
+            radFluxName     Qr;
 
-    moveMesh        false;
+            minimumDelta    1e-8;
-}
 
-infoOutput      false;
+            moveMesh        false;
+        }
+
+        infoOutput      false;
+    }
+)
 
 // ************************************************************************* //

tutorials/combustion/fireFoam/les/smallPoolFire3D/constant/pyrolysisZones

@@ -11,27 +11,34 @@ FoamFile
     format      binary;
     class       dictionary;
     location    "constant";
-    object      pyrolysisProperties;
+    object      pyrolysisZones;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-active          false;
+1
+(
+    pyrolysis
+    {
 
-pyrolysisModel  none;
+        active          false;
 
-regionName      panelRegion;
+        pyrolysisModel  none;
 
-reactingOneDimCoeffs
+        regionName      panelRegion;
-{
-    filmCoupled     false;
 
-    radFluxName     Qr;
+        reactingOneDimCoeffs
+        {
+            filmCoupled     false;
 
-    minimumDelta    1e-8;
+            radFluxName     Qr;
 
-    moveMesh        false;
+            minimumDelta    1e-8;
-}
 
-infoOutput      false;
+            moveMesh        false;
+        }
+
+        infoOutput      false;
+    }
+)
 
 // ************************************************************************* //