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

Conflicts:
	src/thermophysicalModels/radiationModels/derivedFvPatchFields/radiationCoupledBase/radiationCoupledBase.C

applications/utilities/mesh/manipulation/checkMesh/checkTopology.C

@@ -7,6 +7,7 @@
 #include "pointSet.H"
 #include "IOmanip.H"
 #include "emptyPolyPatch.H"
+#include "processorPolyPatch.H"
 
 Foam::label Foam::checkTopology
 (
@@ -289,10 +290,18 @@ Foam::label Foam::checkTopology
         }
     }
 
-    if (!Pstream::parRun())
+
     {
-        Pout<< "\nChecking patch topology for multiply connected surfaces ..."
+        if (!Pstream::parRun())
-            << endl;
+        {
+            Info<< "\nChecking patch topology for multiply connected"
+                << " surfaces..." << endl;
+        }
+        else
+        {
+            Info<< "\nChecking basic patch addressing..." << endl;
+        }
+
 
         const polyBoundaryMesh& patches = mesh.boundaryMesh();
 
@@ -301,90 +310,101 @@ Foam::label Foam::checkTopology
         (
             mesh,
             "nonManifoldPoints",
-            mesh.nPoints()/100
+            mesh.nPoints()/1000
         );
 
         Pout.setf(ios_base::left);
 
-        Pout<< "    "
+        Info<< "    "
             << setw(20) << "Patch"
             << setw(9) << "Faces"
-            << setw(9) << "Points"
+            << setw(9) << "Points";
-            << setw(34) << "Surface topology";
+        if (!Pstream::parRun())
+        {
+            Info<< setw(34) << "Surface topology";
+        }
         if (allGeometry)
         {
-            Pout<< " Bounding box";
+            Info<< " Bounding box";
         }
-        Pout<< endl;
+        Info<< endl;
 
         forAll(patches, patchI)
         {
             const polyPatch& pp = patches[patchI];
 
-                Pout<< "    "
+            if (!isA<processorPolyPatch>(pp))
+            {
+                Info<< "    "
                     << setw(20) << pp.name()
-                    << setw(9) << pp.size()
+                    << setw(9) << returnReduce(pp.size(), sumOp<label>())
-                    << setw(9) << pp.nPoints();
+                    << setw(9) << returnReduce(pp.nPoints(), sumOp<label>());
 
-
+                if (!Pstream::parRun())
-            primitivePatch::surfaceTopo pTyp = pp.surfaceType();
-
-            if (pp.empty())
-            {
-                Pout<< setw(34) << "ok (empty)";
-            }
-            else if (pTyp == primitivePatch::MANIFOLD)
-            {
-                if (pp.checkPointManifold(true, &points))
                 {
-                    Pout<< setw(34) << "multiply connected (shared point)";
+                    primitivePatch::surfaceTopo pTyp = pp.surfaceType();
-                }
-                else
-                {
-                    Pout<< setw(34) << "ok (closed singly connected)";
-                }
 
-                // Add points on non-manifold edges to make set complete
+                    if (pp.empty())
-                pp.checkTopology(false, &points);
-            }
-            else
-            {
-                pp.checkTopology(false, &points);
-
-                if (pTyp == primitivePatch::OPEN)
-                {
-                    Pout<< setw(34) << "ok (non-closed singly connected)";
-                }
-                else
-                {
-                    Pout<< setw(34) << "multiply connected (shared edge)";
-                }
-            }
-
-            if (allGeometry)
-            {
-                const pointField& pts = pp.points();
-                const labelList& mp = pp.meshPoints();
-
-                if (returnReduce(mp.size(), sumOp<label>()) > 0)
-                {
-                    boundBox bb(point::max, point::min);
-                    forAll (mp, i)
                     {
-                        bb.min() = min(bb.min(), pts[mp[i]]);
+                        Info<< setw(34) << "ok (empty)";
-                        bb.max() = max(bb.max(), pts[mp[i]]);
+                    }
+                    else if (pTyp == primitivePatch::MANIFOLD)
+                    {
+                        if (pp.checkPointManifold(true, &points))
+                        {
+                            Info<< setw(34)
+                                << "multiply connected (shared point)";
+                        }
+                        else
+                        {
+                            Info<< setw(34) << "ok (closed singly connected)";
+                        }
+
+                        // Add points on non-manifold edges to make set complete
+                        pp.checkTopology(false, &points);
+                    }
+                    else
+                    {
+                        pp.checkTopology(false, &points);
+
+                        if (pTyp == primitivePatch::OPEN)
+                        {
+                            Info<< setw(34)
+                                << "ok (non-closed singly connected)";
+                        }
+                        else
+                        {
+                            Info<< setw(34)
+                                << "multiply connected (shared edge)";
+                        }
                     }
-                    reduce(bb.min(), minOp<vector>());
-                    reduce(bb.max(), maxOp<vector>());
-                    Pout<< ' ' << bb;
                 }
+
+                if (allGeometry)
+                {
+                    const pointField& pts = pp.points();
+                    const labelList& mp = pp.meshPoints();
+
+                    if (returnReduce(mp.size(), sumOp<label>()) > 0)
+                    {
+                        boundBox bb(point::max, point::min);
+                        forAll (mp, i)
+                        {
+                            bb.min() = min(bb.min(), pts[mp[i]]);
+                            bb.max() = max(bb.max(), pts[mp[i]]);
+                        }
+                        reduce(bb.min(), minOp<vector>());
+                        reduce(bb.max(), maxOp<vector>());
+                        Info<< ' ' << bb;
+                    }
+                }
+                Info<< endl;
             }
-            Pout<< endl;
         }
 
         if (points.size())
         {
-            Pout<< "  <<Writing " << points.size()
+            Info<< "  <<Writing " << returnReduce(points.size(), sumOp<label>())
                 << " conflicting points to set "
                 << points.name() << endl;
 
@@ -392,7 +412,7 @@ Foam::label Foam::checkTopology
             points.write();
         }
 
-        //Pout.setf(ios_base::right);
+        //Info.setf(ios_base::right);
     }
 
     // Force creation of all addressing if requested.

applications/utilities/mesh/manipulation/mergeMeshes/mergeMeshes.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -115,6 +115,7 @@ int main(int argc, char *argv[])
             runTimeMaster
         )
     );
+    const word oldInstance = masterMesh.pointsInstance();
 
 
     Info<< "Reading mesh to add for time = " << runTimeToAdd.timeName() << nl;
@@ -139,7 +140,13 @@ int main(int argc, char *argv[])
 
     masterMesh.addMesh(meshToAdd);
     masterMesh.merge();
-    masterMesh.polyMesh::write();
+
+    if (overwrite)
+    {
+        masterMesh.setInstance(oldInstance);
+    }
+
+    masterMesh.write();
 
     Info<< "\nEnd\n" << endl;
 

src/OSspecific/POSIX/printStack.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -167,6 +167,26 @@ void getSymbolForRaw
 }
 
 
+void error::safePrintStack(std::ostream& os)
+{
+    // Get raw stack symbols
+    void *array[100];
+    size_t size = backtrace(array, 100);
+    char **strings = backtrace_symbols(array, size);
+
+    // See if they contain function between () e.g. "(__libc_start_main+0xd0)"
+    // and see if cplus_demangle can make sense of part before +
+    for (size_t i = 0; i < size; i++)
+    {
+        string msg(strings[i]);
+        fileName programFile;
+        word address;
+
+        os  << '#' << label(i) << '\t' << msg << std::endl;
+    }
+}
+
+
 void error::printStack(Ostream& os)
 {
     // Reads the starting addresses for the dynamically linked libraries

src/OpenFOAM/db/IOstreams/Sstreams/ISstream.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -29,8 +29,7 @@ Description
 
 SourceFiles
     ISstreamI.H
-    ISread.C
+    ISstream.C
-    ISreadToken.C
 
 \*---------------------------------------------------------------------------*/
 

src/OpenFOAM/db/error/error.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -162,6 +162,10 @@ public:
         operator dictionary() const;
 
 
+        //- Helper function to print a stack (if OpenFOAM IO not yet
+        //  initialised)
+        static void safePrintStack(std::ostream&);
+
         //- Helper function to print a stack
         static void printStack(Ostream&);
 

src/OpenFOAM/fields/pointPatchFields/basic/fixedValue/fixedValuePointPatchField.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -121,6 +121,17 @@ public:
         }
 
 
+    // Member functions
+
+        // Access
+
+            //- Return true if this patch field fixes a value
+            virtual bool fixesValue() const
+            {
+                return true;
+            }
+
+
     // Member operators
 
         // Disable assignment operators

src/OpenFOAM/fields/pointPatchFields/pointPatchField/pointPatchField.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -278,6 +278,12 @@ public:
                 return internalField_;
             }
 
+            //- Return true if this patch field fixes a value
+            virtual bool fixesValue() const
+            {
+                return false;
+            }
+
             //- Return true if this patch field is coupled
             virtual bool coupled() const
             {

src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C

@@ -710,7 +710,8 @@ void Foam::autoLayerDriver::setNumLayers
     const indirectPrimitivePatch& pp,
     pointField& patchDisp,
     labelList& patchNLayers,
-    List<extrudeMode>& extrudeStatus
+    List<extrudeMode>& extrudeStatus,
+    label& nAddedCells
 ) const
 {
     const fvMesh& mesh = meshRefiner_.mesh();
@@ -797,6 +798,24 @@ void Foam::autoLayerDriver::setNumLayers
         }
     }
 
+
+    // Calculate number of cells to create
+    nAddedCells = 0;
+    forAll(pp.localFaces(), faceI)
+    {
+        const face& f = pp.localFaces()[faceI];
+
+        // Get max of extrusion per point
+        label nCells = 0;
+        forAll(f, fp)
+        {
+            nCells = max(nCells, patchNLayers[f[fp]]);
+        }
+
+        nAddedCells += nCells;
+    }
+    reduce(nAddedCells, sumOp<label>());
+
     //reduce(nConflicts, sumOp<label>());
     //
     //Info<< "Set displacement to zero for " << nConflicts
@@ -2491,9 +2510,13 @@ void Foam::autoLayerDriver::addLayers
     // extrudeStatus = EXTRUDE.
     labelList patchNLayers(pp().nPoints(), 0);
 
+    // Ideal number of cells added
+    label nIdealAddedCells = 0;
+
     // Whether to add edge for all pp.localPoints.
     List<extrudeMode> extrudeStatus(pp().nPoints(), EXTRUDE);
 
+
     {
         // Get number of layer per point from number of layers per patch
         // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -2506,7 +2529,8 @@ void Foam::autoLayerDriver::addLayers
 
             patchDisp,
             patchNLayers,
-            extrudeStatus
+            extrudeStatus,
+            nIdealAddedCells
         );
 
         // Precalculate mesh edge labels for patch edges
@@ -2776,6 +2800,7 @@ void Foam::autoLayerDriver::addLayers
         layerParams.nSmoothNormals(),
         layerParams.nSmoothSurfaceNormals(),
         layerParams.minMedianAxisAngleCos(),
+        layerParams.featureAngle(),
 
         dispVec,
         medialRatio,
@@ -3101,10 +3126,24 @@ void Foam::autoLayerDriver::addLayers
 
         label nExtruded = countExtrusion(pp, extrudeStatus);
         label nTotFaces = returnReduce(pp().size(), sumOp<label>());
+        label nAddedCells = 0;
+        {
+            forAll(flaggedCells, cellI)
+            {
+                if (flaggedCells[cellI])
+                {
+                    nAddedCells++;
+                }
+            }
+            reduce(nAddedCells, sumOp<label>());
+        }
         Info<< "Extruding " << nExtruded
             << " out of " << nTotFaces
             << " faces (" << 100.0*nExtruded/nTotFaces << "%)."
             << " Removed extrusion at " << nTotChanged << " faces."
+            << endl
+            << "Added " << nAddedCells << " out of " << nIdealAddedCells
+            << " cells (" << 100.0*nAddedCells/nIdealAddedCells << "%)."
             << endl;
 
         if (nTotChanged == 0)

src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.H

@@ -192,7 +192,8 @@ class autoLayerDriver
                     const indirectPrimitivePatch& pp,
                     pointField& patchDisp,
                     labelList& patchNLayers,
-                    List<extrudeMode>& extrudeStatus
+                    List<extrudeMode>& extrudeStatus,
+                    label& nIdealAddedCells
                 ) const;
 
                 //- Helper function to make a pointVectorField with correct
@@ -469,6 +470,7 @@ class autoLayerDriver
                     const label nSmoothNormals,
                     const label nSmoothSurfaceNormals,
                     const scalar minMedianAxisAngleCos,
+                    const scalar featureAngle,
 
                     pointVectorField& dispVec,
                     pointScalarField& medialRatio,

src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriverShrink.C

@@ -33,9 +33,10 @@ Description
 #include "motionSmoother.H"
 #include "pointData.H"
 #include "PointEdgeWave.H"
-#include "OFstream.H"
+#include "OBJstream.H"
 #include "meshTools.H"
 #include "PatchTools.H"
+#include "unitConversion.H"
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
@@ -775,6 +776,7 @@ void Foam::autoLayerDriver::medialAxisSmoothingInfo
     const label nSmoothNormals,
     const label nSmoothSurfaceNormals,
     const scalar minMedianAxisAngleCos,
+    const scalar featureAngle,
 
     pointVectorField& dispVec,
     pointScalarField& medialRatio,
@@ -904,7 +906,7 @@ void Foam::autoLayerDriver::medialAxisSmoothingInfo
                 const edge& e = edges[edgeI];
                 // Approximate medial axis location on edge.
                 //const point medialAxisPt = e.centre(points);
-                vector eVec = e.vec(mesh.points());
+                vector eVec = e.vec(points);
                 scalar eMag = mag(eVec);
                 if (eMag > VSMALL)
                 {
@@ -961,42 +963,111 @@ void Foam::autoLayerDriver::medialAxisSmoothingInfo
 
         labelHashSet adaptPatches(meshMover.adaptPatchIDs());
 
+
         forAll(patches, patchI)
         {
             const polyPatch& pp = patches[patchI];
-            //const pointPatchVectorField& pvf =
+            const pointPatchVectorField& pvf =
-            //    meshMover.displacement().boundaryField()[patchI];
+                meshMover.displacement().boundaryField()[patchI];
 
             if
             (
                 !pp.coupled()
              && !isA<emptyPolyPatch>(pp)
-             //&&  pvf.constraintType() != word::null //exclude fixedValue
              && !adaptPatches.found(patchI)
             )
             {
-                Info<< "Inserting points on patch " << pp.name() << endl;
-
                 const labelList& meshPoints = pp.meshPoints();
 
-                forAll(meshPoints, i)
+                // Check the type of the patchField. The types are
+                //  - fixedValue (0 or more layers) but the >0 layers have
+                //    already been handled in the adaptPatches loop
+                //  - constraint (but not coupled) types, e.g. symmetryPlane,
+                //    slip.
+                if (pvf.fixesValue())
                 {
-                    label pointI = meshPoints[i];
+                    // Disable all movement on fixedValue patchFields
+                    Info<< "Inserting all points on patch " << pp.name()
+                        << endl;
 
-                    if (!pointMedialDist[pointI].valid(dummyTrackData))
+                    forAll(meshPoints, i)
                     {
-                        maxPoints.append(pointI);
+                        label pointI = meshPoints[i];
-                        maxInfo.append
+                        if (!pointMedialDist[pointI].valid(dummyTrackData))
-                        (
+                        {
-                            pointData
+                            maxPoints.append(pointI);
+                            maxInfo.append
                             (
-                                points[pointI],
+                                pointData
-                                0.0,
+                                (
-                                pointI,         // passive data
+                                    points[pointI],
-                                vector::zero    // passive data
+                                    0.0,
-                            )
+                                    pointI,         // passive data
-                        );
+                                    vector::zero    // passive data
-                        pointMedialDist[pointI] = maxInfo.last();
+                                )
+                            );
+                            pointMedialDist[pointI] = maxInfo.last();
+                        }
+                    }
+                }
+                else
+                {
+                    // Based on geometry: analyse angle w.r.t. nearest moving
+                    // point. In the pointWallDist we transported the
+                    // normal as the passive vector. Note that this points
+                    // out of the originating wall so inside of the domain
+                    // on this patch.
+                    Info<< "Inserting points on patch " << pp.name()
+                        << " if angle to nearest layer patch > "
+                        << featureAngle << " degrees." << endl;
+
+                    scalar featureAngleCos = Foam::cos(degToRad(featureAngle));
+                    pointField pointNormals
+                    (
+                        PatchTools::pointNormals
+                        (
+                            mesh,
+                            pp,
+                            identity(pp.size())+pp.start()
+                        )
+                    );
+
+                    forAll(meshPoints, i)
+                    {
+                        label pointI = meshPoints[i];
+
+                        if (!pointMedialDist[pointI].valid(dummyTrackData))
+                        {
+                            // Check if angle not too large.
+                            scalar cosAngle =
+                            (
+                               -pointWallDist[pointI].v()
+                              & pointNormals[i]
+                            );
+                            if (cosAngle > featureAngleCos)
+                            {
+                                // Extrusion direction practically perpendicular
+                                // to the patch. Disable movement at the patch.
+
+                                maxPoints.append(pointI);
+                                maxInfo.append
+                                (
+                                    pointData
+                                    (
+                                        points[pointI],
+                                        0.0,
+                                        pointI,         // passive data
+                                        vector::zero    // passive data
+                                    )
+                                );
+                                pointMedialDist[pointI] = maxInfo.last();
+                            }
+                            else
+                            {
+                                // Extrusion direction makes angle with patch
+                                // so allow sliding - don't insert zero points
+                            }
+                        }
                     }
                 }
             }
@@ -1141,13 +1212,12 @@ void Foam::autoLayerDriver::shrinkMeshMedialDistance
     // reduce thickness where thickness/medial axis distance large
     // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-    autoPtr<OFstream> str;
+    autoPtr<OBJstream> str;
-    label vertI = 0;
     if (debug)
     {
         str.reset
         (
-            new OFstream
+            new OBJstream
             (
                 mesh.time().path()
               / "thicknessRatioExcludePoints_"
@@ -1159,13 +1229,12 @@ void Foam::autoLayerDriver::shrinkMeshMedialDistance
             << str().name() << endl;
     }
 
-    autoPtr<OFstream> medialVecStr;
+    autoPtr<OBJstream> medialVecStr;
-    label medialVertI = 0;
     if (debug)
     {
         medialVecStr.reset
         (
-            new OFstream
+            new OBJstream
             (
                 mesh.time().path()
               / "thicknessRatioExcludeMedialVec_"
@@ -1227,21 +1296,19 @@ void Foam::autoLayerDriver::shrinkMeshMedialDistance
                 if (str.valid())
                 {
                     const point& pt = mesh.points()[pointI];
-                    meshTools::writeOBJ(str(), pt);
+                    str().write(linePointRef(pt, pt+patchDisp[patchPointI]));
-                    vertI++;
-                    meshTools::writeOBJ(str(), pt+patchDisp[patchPointI]);
-                    vertI++;
-                    str()<< "l " << vertI-1 << ' ' << vertI << nl;
                 }
                 if (medialVecStr.valid())
                 {
                     const point& pt = mesh.points()[pointI];
-                    meshTools::writeOBJ(medialVecStr(), pt);
+                    medialVecStr().write
-                    medialVertI++;
+                    (
-                    meshTools::writeOBJ(medialVecStr(), medialVec[pointI]);
+                        linePointRef
-                    medialVertI++;
+                        (
-                    medialVecStr()<< "l " << medialVertI-1
+                            pt,
-                        << ' ' << medialVertI << nl;
+                            medialVec[pointI]
+                        )
+                    );
                 }
             }
         }

src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoSnapDriver.H

@@ -178,6 +178,24 @@ class autoSnapDriver
                     vector& edgeOffset  // offset from pt to point on edge
                 ) const;
 
+
+                //- For any reverse (so from feature back to mesh) attraction:
+                //  add attraction if diagonal points on face attracted
+                void stringFeatureEdges
+                (
+                    const label iter,
+                    const scalar featureCos,
+
+                    const indirectPrimitivePatch& pp,
+                    const scalarField& snapDist,
+
+                    const vectorField& rawPatchAttraction,
+                    const List<pointConstraint>& rawPatchConstraints,
+
+                    vectorField& patchAttraction,
+                    List<pointConstraint>& patchConstraints
+                ) const;
+
                 //- Return hit if on multiple points
                 pointIndexHit findMultiPatchPoint
                 (

src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoSnapDriverFeature.C

@@ -848,17 +848,17 @@ void Foam::autoSnapDriver::featureAttractionUsingReconstruction
         vector d = r.refPoint()-pt;
         d -= (d&n)*n;
 
-        // Correct for attraction to non-dominant face
+        //// Correct for attraction to non-dominant face
-        correctAttraction
+        //correctAttraction
-        (
+        //(
-            surfacePoints,
+        //    surfacePoints,
-            surfaceCount,
+        //    surfaceCount,
-            r.refPoint(),
+        //    r.refPoint(),
-            n,                  // normalised normal
+        //    n,                  // normalised normal
-            pt,
+        //    pt,
-
+        //
-            d                   // perpendicular offset vector
+        //    d                   // perpendicular offset vector
-        );
+        //);
 
         // Trim to snap distance
         if (magSqr(d) > sqr(snapDist[pointI]))
@@ -893,6 +893,15 @@ void Foam::autoSnapDriver::featureAttractionUsingReconstruction
         patchConstraint.applyConstraint(surfaceNormals[0]);
         patchConstraint.applyConstraint(surfaceNormals[1]);
         patchConstraint.applyConstraint(surfaceNormals[2]);
+
+        //Pout<< "# Feature point " << pt << nl;
+        //meshTools::writeOBJ(Pout, pt);
+        //meshTools::writeOBJ(Pout, surfacePoints[0]);
+        //meshTools::writeOBJ(Pout, surfacePoints[1]);
+        //meshTools::writeOBJ(Pout, surfacePoints[2]);
+        //Pout<< "l 1 2" << nl
+        //    << "l 1 3" << nl
+        //    << "l 1 4" << nl;
     }
 }
 
@@ -1001,6 +1010,195 @@ void Foam::autoSnapDriver::featureAttractionUsingReconstruction
 }
 
 
+void Foam::autoSnapDriver::stringFeatureEdges
+(
+    const label iter,
+    const scalar featureCos,
+
+    const indirectPrimitivePatch& pp,
+    const scalarField& snapDist,
+
+    const vectorField& rawPatchAttraction,
+    const List<pointConstraint>& rawPatchConstraints,
+
+    vectorField& patchAttraction,
+    List<pointConstraint>& patchConstraints
+) const
+{
+    // Snap edges to feature edges
+    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~
+    // Walk existing edges and snap remaining ones (that are marked as
+    // feature edges in rawPatchConstraints)
+
+    // What this does is fill in any faces where not all points
+    // on the face are being attracted:
+    /*
+           +
+          / \
+         /   \
+      ---+    +---
+         \   /
+          \ /
+           +
+    */
+    // so the top and bottom will never get attracted since the nearest
+    // back from the feature edge will always be one of the left or right
+    // points since the face is diamond like. So here we walk the feature edges
+    // and add any non-attracted points.
+
+
+    while (true)
+    {
+        label nChanged = 0;
+
+        const labelListList& pointEdges = pp.pointEdges();
+        forAll(pointEdges, pointI)
+        {
+            if (patchConstraints[pointI].first() == 2)
+            {
+                const point& pt = pp.localPoints()[pointI];
+                const labelList& pEdges = pointEdges[pointI];
+                const vector& featVec = patchConstraints[pointI].second();
+
+                // Detect whether there are edges in both directions.
+                // (direction along the feature edge that is)
+                bool hasPos = false;
+                bool hasNeg = false;
+
+                forAll(pEdges, pEdgeI)
+                {
+                    const edge& e = pp.edges()[pEdges[pEdgeI]];
+                    label nbrPointI = e.otherVertex(pointI);
+
+                    if (patchConstraints[nbrPointI].first() > 1)
+                    {
+                        const point& nbrPt = pp.localPoints()[nbrPointI];
+                        const point featPt =
+                            nbrPt + patchAttraction[nbrPointI];
+                        const scalar cosAngle = (featVec & (featPt-pt));
+
+                        if (cosAngle > 0)
+                        {
+                            hasPos = true;
+                        }
+                        else
+                        {
+                            hasNeg = true;
+                        }
+                    }
+                }
+
+                if (!hasPos || !hasNeg)
+                {
+                    //Pout<< "**Detected feature string end at  "
+                    //    << pp.localPoints()[pointI] << endl;
+
+                    // No string. Assign best choice on either side
+                    label bestPosPointI = -1;
+                    scalar minPosDistSqr = GREAT;
+                    label bestNegPointI = -1;
+                    scalar minNegDistSqr = GREAT;
+
+                    forAll(pEdges, pEdgeI)
+                    {
+                        const edge& e = pp.edges()[pEdges[pEdgeI]];
+                        label nbrPointI = e.otherVertex(pointI);
+
+                        if
+                        (
+                            patchConstraints[nbrPointI].first() <= 1
+                         && rawPatchConstraints[nbrPointI].first() > 1
+                        )
+                        {
+                            const vector& nbrFeatVec =
+                                rawPatchConstraints[pointI].second();
+
+                            if (mag(featVec&nbrFeatVec) > featureCos)
+                            {
+                                // nbrPointI attracted to sameish feature
+                                // Note: also check on position.
+
+                                scalar d2 = magSqr
+                                (
+                                    rawPatchAttraction[nbrPointI]
+                                );
+
+                                const point featPt =
+                                    pp.localPoints()[nbrPointI]
+                                  + rawPatchAttraction[nbrPointI];
+                                const scalar cosAngle =
+                                    (featVec & (featPt-pt));
+
+                                if (cosAngle > 0)
+                                {
+                                    if (!hasPos && d2 < minPosDistSqr)
+                                    {
+                                        minPosDistSqr = d2;
+                                        bestPosPointI = nbrPointI;
+                                    }
+                                }
+                                else
+                                {
+                                    if (!hasNeg && d2 < minNegDistSqr)
+                                    {
+                                        minNegDistSqr = d2;
+                                        bestNegPointI = nbrPointI;
+                                    }
+                                }
+                            }
+                        }
+                    }
+
+                    if (bestPosPointI != -1)
+                    {
+                        // Use reconstructed-feature attraction. Use only
+                        // part of it since not sure...
+                        //const point& bestPt =
+                        //    pp.localPoints()[bestPosPointI];
+                        //Pout<< "**Overriding point " << bestPt
+                        //    << " on reconstructed feature edge at "
+                        //    << rawPatchAttraction[bestPosPointI]+bestPt
+                        //    << " to attracted-to-feature-edge." << endl;
+                        patchAttraction[bestPosPointI] =
+                            0.5*rawPatchAttraction[bestPosPointI];
+                        patchConstraints[bestPosPointI] =
+                            rawPatchConstraints[bestPosPointI];
+
+                        nChanged++;
+                    }
+                    if (bestNegPointI != -1)
+                    {
+                        // Use reconstructed-feature attraction. Use only
+                        // part of it since not sure...
+                        //const point& bestPt =
+                        //    pp.localPoints()[bestNegPointI];
+                        //Pout<< "**Overriding point " << bestPt
+                        //    << " on reconstructed feature edge at "
+                        //    << rawPatchAttraction[bestNegPointI]+bestPt
+                        //    << " to attracted-to-feature-edge." << endl;
+                        patchAttraction[bestNegPointI] =
+                            0.5*rawPatchAttraction[bestNegPointI];
+                        patchConstraints[bestNegPointI] =
+                            rawPatchConstraints[bestNegPointI];
+
+                        nChanged++;
+                    }
+                }
+            }
+        }
+
+
+        reduce(nChanged, sumOp<label>());
+        Info<< "Stringing feature edges : changed " << nChanged << " points"
+            << endl;
+        if (nChanged == 0)
+        {
+            break;
+        }
+    }
+}
+
+
 Foam::pointIndexHit Foam::autoSnapDriver::findNearFeatureEdge
 (
     const indirectPrimitivePatch& pp,
@@ -2113,156 +2311,20 @@ void Foam::autoSnapDriver::featureAttractionUsingFeatureEdges
     // Walk existing edges and snap remaining ones (that are marked as
     // feature edges in allPatchConstraints)
 
-    while (true)
+    stringFeatureEdges
-    {
+    (
-        label nChanged = 0;
+        iter,
+        featureCos,
 
-        const labelListList& pointEdges = pp.pointEdges();
+        pp,
-        forAll(pointEdges, pointI)
+        snapDist,
-        {
-            if (patchConstraints[pointI].first() == 2)
-            {
-                const point& pt = pp.localPoints()[pointI];
-                const labelList& pEdges = pointEdges[pointI];
-                const vector& featVec = patchConstraints[pointI].second();
 
-                // Detect whether there are edges in both directions.
+        allPatchAttraction,
-                // (direction along the feature edge that is)
+        allPatchConstraints,
-                bool hasPos = false;
-                bool hasNeg = false;
-
-                forAll(pEdges, pEdgeI)
-                {
-                    const edge& e = pp.edges()[pEdges[pEdgeI]];
-                    label nbrPointI = e.otherVertex(pointI);
-
-                    if (patchConstraints[nbrPointI].first() > 1)
-                    {
-                        const point& nbrPt = pp.localPoints()[nbrPointI];
-                        const point featPt =
-                            nbrPt + patchAttraction[nbrPointI];
-                        const scalar cosAngle = (featVec & (featPt-pt));
-
-                        if (cosAngle > 0)
-                        {
-                            hasPos = true;
-                        }
-                        else
-                        {
-                            hasNeg = true;
-                        }
-                    }
-                }
-
-                if (!hasPos || !hasNeg)
-                {
-                    //Pout<< "**Detected feature string end at  "
-                    //    << pp.localPoints()[pointI] << endl;
-
-                    // No string. Assign best choice on either side
-                    label bestPosPointI = -1;
-                    scalar minPosDistSqr = GREAT;
-                    label bestNegPointI = -1;
-                    scalar minNegDistSqr = GREAT;
-
-                    forAll(pEdges, pEdgeI)
-                    {
-                        const edge& e = pp.edges()[pEdges[pEdgeI]];
-                        label nbrPointI = e.otherVertex(pointI);
-
-                        if
-                        (
-                            patchConstraints[nbrPointI].first() <= 1
-                         && allPatchConstraints[nbrPointI].first() > 1
-                        )
-                        {
-                            const vector& nbrFeatVec =
-                                allPatchConstraints[pointI].second();
-
-                            if (mag(featVec&nbrFeatVec) > featureCos)
-                            {
-                                // nbrPointI attracted to sameish feature
-                                // Note: also check on position.
-
-                                scalar d2 = magSqr
-                                (
-                                    allPatchAttraction[nbrPointI]
-                                );
-
-                                const point featPt =
-                                    pp.localPoints()[nbrPointI]
-                                  + allPatchAttraction[nbrPointI];
-                                const scalar cosAngle =
-                                    (featVec & (featPt-pt));
-
-                                if (cosAngle > 0)
-                                {
-                                    if (!hasPos && d2 < minPosDistSqr)
-                                    {
-                                        minPosDistSqr = d2;
-                                        bestPosPointI = nbrPointI;
-                                    }
-                                }
-                                else
-                                {
-                                    if (!hasNeg && d2 < minNegDistSqr)
-                                    {
-                                        minNegDistSqr = d2;
-                                        bestNegPointI = nbrPointI;
-                                    }
-                                }
-                            }
-                        }
-                    }
-
-                    if (bestPosPointI != -1)
-                    {
-                        // Use reconstructed-feature attraction. Use only
-                        // part of it since not sure...
-                        //const point& bestPt =
-                        //    pp.localPoints()[bestPosPointI];
-                        //Pout<< "**Overriding point " << bestPt
-                        //    << " on reconstructed feature edge at "
-                        //    << allPatchAttraction[bestPosPointI]+bestPt
-                        //    << " to attracted-to-feature-edge." << endl;
-                        patchAttraction[bestPosPointI] =
-                            0.5*allPatchAttraction[bestPosPointI];
-                        patchConstraints[bestPosPointI] =
-                            allPatchConstraints[bestPosPointI];
-
-                        nChanged++;
-                    }
-                    if (bestNegPointI != -1)
-                    {
-                        // Use reconstructed-feature attraction. Use only
-                        // part of it since not sure...
-                        //const point& bestPt =
-                        //    pp.localPoints()[bestNegPointI];
-                        //Pout<< "**Overriding point " << bestPt
-                        //    << " on reconstructed feature edge at "
-                        //    << allPatchAttraction[bestNegPointI]+bestPt
-                        //    << " to attracted-to-feature-edge." << endl;
-                        patchAttraction[bestNegPointI] =
-                            0.5*allPatchAttraction[bestNegPointI];
-                        patchConstraints[bestNegPointI] =
-                            allPatchConstraints[bestNegPointI];
-
-                        nChanged++;
-                    }
-                }
-            }
-        }
-
-
-        reduce(nChanged, sumOp<label>());
-        Info<< "Stringing feature edges : changed " << nChanged << " points"
-            << endl;
-        if (nChanged == 0)
-        {
-            break;
-        }
-    }
 
+        patchAttraction,
+        patchConstraints
+    );
 
 
     // Avoid diagonal attraction

src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementRefine.C

@@ -312,7 +312,8 @@ Foam::label Foam::meshRefinement::markFeatureRefinement
                     {
                         Pout<< "Adding particle from point:" << pointI
                             << " coord:" << featureMesh.points()[pointI]
-                            << " pEdges:" << pointEdges[pointI]
+                            << " since number of emanating edges:"
+                            << pointEdges[pointI].size()
                             << endl;
                     }
 

src/regionModels/thermoBaffleModels/derivedFvPatchFields/temperatureThermoBaffle/temperatureThermoBaffleFvPatchScalarField.C

@@ -46,7 +46,7 @@ temperatureThermoBaffleFvPatchScalarField
     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField(p, iF),
     owner_(false),
     baffle_(),
-    solidThermoType_("undefined"),
+    solidThermoType_(new primitiveEntry("thermoType", "undefined")),
     dict_(dictionary::null)
 {}
 
@@ -69,7 +69,7 @@ temperatureThermoBaffleFvPatchScalarField
     ),
     owner_(ptf.owner_),
     baffle_(ptf.baffle_),
-    solidThermoType_(ptf.solidThermoType_),
+    solidThermoType_(ptf.solidThermoType_().clone()),
     dict_(ptf.dict_)
 {}
 
@@ -85,7 +85,7 @@ temperatureThermoBaffleFvPatchScalarField
     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField(p, iF, dict),
     owner_(false),
     baffle_(),
-    solidThermoType_(),
+    solidThermoType_(new primitiveEntry("thermoType", "undefined")),
     dict_(dict)
 {
     if (!isA<mappedPatchBase>(patch().patch()))
@@ -126,7 +126,7 @@ temperatureThermoBaffleFvPatchScalarField
         Info << "Creating thermal baffle" <<  nbrMesh << endl;
         baffle_.reset(baffle::New(thisMesh, dict).ptr());
         owner_ = true;
-        dict.lookup("thermoType") >> solidThermoType_;
+        solidThermoType_ = dict.lookupEntry("thermoType", false, false).clone();
         baffle_->rename(nbrMesh);
     }
 }
@@ -142,7 +142,7 @@ temperatureThermoBaffleFvPatchScalarField
     turbulentTemperatureCoupledBaffleMixedFvPatchScalarField(ptf, iF),
     owner_(ptf.owner_),
     baffle_(ptf.baffle_),
-    solidThermoType_(ptf.solidThermoType_),
+    solidThermoType_(ptf.solidThermoType_().clone()),
     dict_(ptf.dict_)
 {}
 
@@ -219,8 +219,7 @@ void temperatureThermoBaffleFvPatchScalarField::write(Ostream& os) const
         os.writeKeyword(word(thermoModel + "Coeffs"));
         os << dict_.subDict(thermoModel + "Coeffs") << nl;
 
-        os.writeKeyword("thermoType") << solidThermoType_
+        os << solidThermoType_() << nl;
-            << token::END_STATEMENT << nl;
 
         os.writeKeyword("mixture");
         os << dict_.subDict("mixture") << nl;

src/regionModels/thermoBaffleModels/derivedFvPatchFields/temperatureThermoBaffle/temperatureThermoBaffleFvPatchScalarField.H

@@ -140,7 +140,7 @@ class temperatureThermoBaffleFvPatchScalarField
         autoPtr<regionModels::thermoBaffleModels::thermoBaffleModel> baffle_;
 
         //- Solid thermo type
-        word solidThermoType_;
+        autoPtr<entry> solidThermoType_;
 
         //- Dictionary
         dictionary dict_;

src/regionModels/thermoBaffleModels/thermoBaffle2D/thermoBaffle2D.C

@@ -318,7 +318,13 @@ const tmp<volScalarField> thermoBaffle2D::Cp() const
 
 const volScalarField& thermoBaffle2D::kappaRad() const
 {
-    return thermo_->kappaRad();
+    // ***HGW return thermo_->kappaRad();
+    FatalErrorIn
+    (
+        "thermoBaffle2D::kappaRad()"
+    )   << "not currently supported"
+        << exit(FatalError);
+    return volScalarField::null();
 }
 
 

src/thermophysicalModels/basic/Make/files

@@ -1,11 +1,11 @@
 basicThermo/basicThermo.C
 fluidThermo/fluidThermo.C
 
-psiThermo/psiThermo/psiThermo.C
+psiThermo/psiThermo.C
-psiThermo/hePsiThermo/hePsiThermos.C
+psiThermo/psiThermos.C
 
-rhoThermo/rhoThermo/rhoThermo.C
+rhoThermo/rhoThermo.C
-rhoThermo/heRhoThermo/heRhoThermos.C
+rhoThermo/rhoThermos.C
 
 derivedFvPatchFields/fixedEnergy/fixedEnergyFvPatchScalarField.C
 derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.C

src/thermophysicalModels/basic/basicThermo/basicThermo.H

@@ -105,10 +105,22 @@ public:
 
     // Selectors
 
+        //- Generic lookup for each of the related thermodynamics packages
+        template<class Thermo, class Table>
+        static typename Table::iterator lookupThermo
+        (
+            const dictionary& thermoDict,
+            Table* tablePtr
+        );
+
         //- Generic New for each of the related thermodynamics packages
         template<class Thermo>
         static autoPtr<Thermo> New(const fvMesh&);
 
+        //- Generic New for each of the related thermodynamics packages
+        template<class Thermo>
+        static autoPtr<Thermo> New(const fvMesh&, const dictionary&);
+
         //- Specialisation of the Generic New for basicThermo
         static autoPtr<basicThermo> New(const fvMesh&);
 

src/thermophysicalModels/basic/basicThermo/basicThermoTemplates.C

@@ -27,25 +27,13 @@ License
 
 // * * * * * * * * * * * * * * * * Selectors * * * * * * * * * * * * * * * * //
 
-template<class Thermo>
+template<class Thermo, class Table>
-Foam::autoPtr<Thermo> Foam::basicThermo::New
+typename Table::iterator Foam::basicThermo::lookupThermo
 (
-    const fvMesh& mesh
+    const dictionary& thermoDict,
+    Table* tablePtr
 )
 {
-    IOdictionary thermoDict
-    (
-        IOobject
-        (
-            "thermophysicalProperties",
-            mesh.time().constant(),
-            mesh,
-            IOobject::MUST_READ_IF_MODIFIED,
-            IOobject::NO_WRITE,
-            false
-        )
-    );
-
     word thermoTypeName;
 
     if (thermoDict.isDict("thermoType"))
@@ -76,16 +64,13 @@ Foam::autoPtr<Thermo> Foam::basicThermo::New
           + word(thermoTypeDict.lookup("specie")) + ">>,"
           + word(thermoTypeDict.lookup("energy")) + ">>>";
 
-        Info<< thermoTypeName << endl;
-
         // Lookup the thermo package
-        typename Thermo::fvMeshConstructorTable::iterator cstrIter =
+        typename Table::iterator cstrIter = tablePtr->find(thermoTypeName);
-            Thermo::fvMeshConstructorTablePtr_->find(thermoTypeName);
 
         // Print error message if package not found in the table
-        if (cstrIter == Thermo::fvMeshConstructorTablePtr_->end())
+        if (cstrIter == tablePtr->end())
         {
-            FatalErrorIn(Thermo::typeName + "::New(const fvMesh&)")
+            FatalErrorIn(Thermo::typeName + "::New")
                 << "Unknown " << Thermo::typeName << " type " << nl
                 << "thermoType" << thermoTypeDict << nl << nl
                 << "Valid " << Thermo::typeName << " types are:" << nl << nl;
@@ -93,7 +78,7 @@ Foam::autoPtr<Thermo> Foam::basicThermo::New
             // Get the list of all the suitable thermo packages available
             wordList validThermoTypeNames
             (
-                Thermo::fvMeshConstructorTablePtr_->sortedToc()
+                tablePtr->sortedToc()
             );
 
             // Build a table of the thermo packages constituent parts
@@ -123,7 +108,7 @@ Foam::autoPtr<Thermo> Foam::basicThermo::New
             FatalError<< exit(FatalError);
         }
 
-        return autoPtr<Thermo>(cstrIter()(mesh));
+        return cstrIter;
     }
     else
     {
@@ -131,22 +116,69 @@ Foam::autoPtr<Thermo> Foam::basicThermo::New
 
         Info<< "Selecting thermodynamics package " << thermoTypeName << endl;
 
-        typename Thermo::fvMeshConstructorTable::iterator cstrIter =
+        typename Table::iterator cstrIter = tablePtr->find(thermoTypeName);
-            Thermo::fvMeshConstructorTablePtr_->find(thermoTypeName);
 
-        if (cstrIter == Thermo::fvMeshConstructorTablePtr_->end())
+        if (cstrIter == tablePtr->end())
         {
-            FatalErrorIn(Thermo::typeName + "::New(const fvMesh&)")
+            FatalErrorIn(Thermo::typeName + "::New")
                 << "Unknown " << Thermo::typeName << " type "
                 << thermoTypeName << nl << nl
                 << "Valid " << Thermo::typeName << " types are:" << nl
-                << Thermo::fvMeshConstructorTablePtr_->sortedToc() << nl
+                << tablePtr->sortedToc() << nl
                 << exit(FatalError);
         }
 
-        return autoPtr<Thermo>(cstrIter()(mesh));
+        return cstrIter;
     }
 }
 
 
+template<class Thermo>
+Foam::autoPtr<Thermo> Foam::basicThermo::New
+(
+    const fvMesh& mesh
+)
+{
+    IOdictionary thermoDict
+    (
+        IOobject
+        (
+            "thermophysicalProperties",
+            mesh.time().constant(),
+            mesh,
+            IOobject::MUST_READ_IF_MODIFIED,
+            IOobject::NO_WRITE,
+            false
+        )
+    );
+
+    typename Thermo::fvMeshConstructorTable::iterator cstrIter =
+        lookupThermo<Thermo, typename Thermo::fvMeshConstructorTable>
+        (
+            thermoDict,
+            Thermo::fvMeshConstructorTablePtr_
+        );
+
+    return autoPtr<Thermo>(cstrIter()(mesh));
+}
+
+
+template<class Thermo>
+Foam::autoPtr<Thermo> Foam::basicThermo::New
+(
+    const fvMesh& mesh,
+    const dictionary& dict
+)
+{
+    typename Thermo::dictionaryConstructorTable::iterator cstrIter =
+        lookupThermo<Thermo, typename Thermo::dictionaryConstructorTable>
+        (
+            dict,
+            Thermo::dictionaryConstructorTablePtr_
+        );
+
+    return autoPtr<Thermo>(cstrIter()(mesh, dict));
+}
+
+
 // ************************************************************************* //

src/thermophysicalModels/basic/fluidThermo/makeThermo.H

@@ -37,7 +37,7 @@ Description
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-#define makeThermo(BaseThermo,Cthermo,Mixture,Transport,Type,Thermo,EqnOfState,Specie)\
+#define makeThermoTypedefs(BaseThermo,Cthermo,Mixture,Transport,Type,Thermo,EqnOfState,Specie)\
                                                                               \
 typedef                                                                       \
     Transport                                                                 \
@@ -58,10 +58,8 @@ typedef                                                                       \
 typedef                                                                       \
     Cthermo                                                                   \
     <                                                                         \
-        Mixture                                                               \
+        BaseThermo,                                                           \
-        <                                                                     \
+        Mixture<Transport##Type##Thermo##EqnOfState##Specie>                  \
-            Transport##Type##Thermo##EqnOfState##Specie                       \
-        >                                                                     \
     > Cthermo##Mixture##Transport##Type##Thermo##EqnOfState##Specie;          \
                                                                               \
 defineTemplateTypeNameAndDebugWithName                                        \
@@ -73,7 +71,22 @@ defineTemplateTypeNameAndDebugWithName                                        \
       + ">>"                                                                  \
     ).c_str(),                                                                \
     0                                                                         \
-);                                                                            \
+);
+
+
+#define makeThermo(BaseThermo,Cthermo,Mixture,Transport,Type,Thermo,EqnOfState,Specie)\
+                                                                              \
+makeThermoTypedefs                                                            \
+(                                                                             \
+    BaseThermo,                                                               \
+    Cthermo,                                                                  \
+    Mixture,                                                                  \
+    Transport,                                                                \
+    Type,                                                                     \
+    Thermo,                                                                   \
+    EqnOfState,                                                               \
+    Specie                                                                    \
+)                                                                             \
                                                                               \
 addToRunTimeSelectionTable                                                    \
 (                                                                             \

src/thermophysicalModels/basic/psiThermo/hePsiThermo.C

@@ -27,8 +27,8 @@ License
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-void Foam::hePsiThermo<MixtureType>::calculate()
+void Foam::hePsiThermo<BasicPsiThermo, MixtureType>::calculate()
 {
     const scalarField& hCells = this->he_.internalField();
     const scalarField& pCells = this->p_.internalField();
@@ -101,10 +101,10 @@ void Foam::hePsiThermo<MixtureType>::calculate()
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-Foam::hePsiThermo<MixtureType>::hePsiThermo(const fvMesh& mesh)
+Foam::hePsiThermo<BasicPsiThermo, MixtureType>::hePsiThermo(const fvMesh& mesh)
 :
-    heThermo<psiThermo, MixtureType>(mesh)
+    heThermo<BasicPsiThermo, MixtureType>(mesh)
 {
     calculate();
 
@@ -115,19 +115,20 @@ Foam::hePsiThermo<MixtureType>::hePsiThermo(const fvMesh& mesh)
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-Foam::hePsiThermo<MixtureType>::~hePsiThermo()
+Foam::hePsiThermo<BasicPsiThermo, MixtureType>::~hePsiThermo()
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-void Foam::hePsiThermo<MixtureType>::correct()
+void Foam::hePsiThermo<BasicPsiThermo, MixtureType>::correct()
 {
     if (debug)
     {
-        Info<< "entering hePsiThermo<MixtureType>::correct()" << endl;
+        Info<< "entering hePsiThermo<BasicPsiThermo, MixtureType>::correct()"
+            << endl;
     }
 
     // force the saving of the old-time values
@@ -137,7 +138,8 @@ void Foam::hePsiThermo<MixtureType>::correct()
 
     if (debug)
     {
-        Info<< "exiting hePsiThermo<MixtureType>::correct()" << endl;
+        Info<< "exiting hePsiThermo<BasicPsiThermo, MixtureType>::correct()"
+            << endl;
     }
 }
 

src/thermophysicalModels/basic/psiThermo/hePsiThermo.H

@@ -25,7 +25,7 @@ Class
     Foam::hePsiThermo
 
 Description
-    Enthalpy for a mixture based on compressibility
+    Energy for a mixture based on compressibility
 
 SourceFiles
     hePsiThermo.C
@@ -47,10 +47,10 @@ namespace Foam
                          Class hePsiThermo Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 class hePsiThermo
 :
-    public heThermo<psiThermo, MixtureType>
+    public heThermo<BasicPsiThermo, MixtureType>
 {
     // Private Member Functions
 
@@ -58,7 +58,7 @@ class hePsiThermo
         void calculate();
 
         //- Construct as copy (not implemented)
-        hePsiThermo(const hePsiThermo<MixtureType>&);
+        hePsiThermo(const hePsiThermo<BasicPsiThermo, MixtureType>&);
 
 public:
 

src/thermophysicalModels/basic/rhoThermo/heRhoThermo.C

@@ -27,8 +27,8 @@ License
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-void Foam::heRhoThermo<MixtureType>::calculate()
+void Foam::heRhoThermo<BasicPsiThermo, MixtureType>::calculate()
 {
     const scalarField& hCells = this->he().internalField();
     const scalarField& pCells = this->p_.internalField();
@@ -106,10 +106,10 @@ void Foam::heRhoThermo<MixtureType>::calculate()
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-Foam::heRhoThermo<MixtureType>::heRhoThermo(const fvMesh& mesh)
+Foam::heRhoThermo<BasicPsiThermo, MixtureType>::heRhoThermo(const fvMesh& mesh)
 :
-    heThermo<rhoThermo, MixtureType>(mesh)
+    heThermo<BasicPsiThermo, MixtureType>(mesh)
 {
     calculate();
 }
@@ -117,15 +117,15 @@ Foam::heRhoThermo<MixtureType>::heRhoThermo(const fvMesh& mesh)
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-Foam::heRhoThermo<MixtureType>::~heRhoThermo()
+Foam::heRhoThermo<BasicPsiThermo, MixtureType>::~heRhoThermo()
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-void Foam::heRhoThermo<MixtureType>::correct()
+void Foam::heRhoThermo<BasicPsiThermo, MixtureType>::correct()
 {
     if (debug)
     {

src/thermophysicalModels/basic/rhoThermo/heRhoThermo.H

@@ -25,7 +25,7 @@ Class
     Foam::heRhoThermo
 
 Description
-    Enthalpy for a mixture based on density
+    Energy for a mixture based on density
 
 SourceFiles
     heRhoThermo.C
@@ -47,10 +47,10 @@ namespace Foam
                          Class heRhoThermo Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 class heRhoThermo
 :
-    public heThermo<rhoThermo, MixtureType>
+    public heThermo<BasicPsiThermo, MixtureType>
 {
     // Private Member Functions
 
@@ -58,7 +58,7 @@ class heRhoThermo
         void calculate();
 
         //- Construct as copy (not implemented)
-        heRhoThermo(const heRhoThermo<MixtureType>&);
+        heRhoThermo(const heRhoThermo<BasicPsiThermo, MixtureType>&);
 
 
 public:

src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel/basicChemistryModelTemplates.C

@@ -58,12 +58,11 @@ Foam::autoPtr<ChemistryModel> Foam::basicChemistryModel::New
 
         Info<< "Selecting chemistry type " << chemistryTypeDict << endl;
 
-        const int nCmpt = 8;
+        const int nCmpt = 7;
         const char* cmptNames[nCmpt] =
         {
             "chemistrySolver",
-            "chemistryModel",
+            "chemistryThermo",
-            "???ChemistryModel",
             "transport",
             "thermo",
             "equationOfState",
@@ -109,9 +108,8 @@ Foam::autoPtr<ChemistryModel> Foam::basicChemistryModel::New
         // Construct the name of the chemistry type from the components
         chemistryTypeName =
             word(chemistryTypeDict.lookup("chemistrySolver")) + '<'
-          + word(chemistryTypeDict.lookup("chemistryModel")) + '<'
+          + word(chemistryTypeDict.lookup("chemistryThermo")) + ','
-          + ChemistryModel::typeName + ','
+          + thermoTypeName + ">";
-          + thermoTypeName + ">>";
 
         typename ChemistryModel::fvMeshConstructorTable::iterator cstrIter =
             ChemistryModel::fvMeshConstructorTablePtr_->find(chemistryTypeName);
@@ -160,8 +158,7 @@ Foam::autoPtr<ChemistryModel> Foam::basicChemistryModel::New
             FatalError<< exit(FatalError);
         }
 
-        return autoPtr<ChemistryModel>
+        return autoPtr<ChemistryModel>(cstrIter()(mesh));
-            (cstrIter()(mesh, typeName, chemistryTypeName));
     }
     else
     {
@@ -183,8 +180,7 @@ Foam::autoPtr<ChemistryModel> Foam::basicChemistryModel::New
                 << exit(FatalError);
         }
 
-        return autoPtr<ChemistryModel>
+        return autoPtr<ChemistryModel>(cstrIter()(mesh));
-            (cstrIter()(mesh, typeName, chemistryTypeName));
     }
 }
 

src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C

@@ -23,21 +23,18 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-#include "ODEChemistryModel.H"
+#include "chemistryModel.H"
-#include "chemistrySolver.H"
 #include "reactingMixture.H"
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
 template<class CompType, class ThermoType>
-Foam::ODEChemistryModel<CompType, ThermoType>::ODEChemistryModel
+Foam::chemistryModel<CompType, ThermoType>::chemistryModel
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& ODEModelName,
-    const word& thermoTypeName
 )
 :
-    CompType(mesh, thermoTypeName),
+    CompType(mesh),
 
     ODE(),
 
@@ -80,7 +77,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::ODEChemistryModel
         );
     }
 
-    Info<< "ODEChemistryModel: Number of species = " << nSpecie_
+    Info<< "chemistryModel: Number of species = " << nSpecie_
         << " and reactions = " << nReaction_ << endl;
 }
 
@@ -88,7 +85,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::ODEChemistryModel
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
 template<class CompType, class ThermoType>
-Foam::ODEChemistryModel<CompType, ThermoType>::~ODEChemistryModel()
+Foam::chemistryModel<CompType, ThermoType>::~chemistryModel()
 {}
 
 
@@ -96,7 +93,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::~ODEChemistryModel()
 
 template<class CompType, class ThermoType>
 Foam::tmp<Foam::scalarField>
-Foam::ODEChemistryModel<CompType, ThermoType>::omega
+Foam::chemistryModel<CompType, ThermoType>::omega
 (
     const scalarField& c,
     const scalar T,
@@ -138,7 +135,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::omega
 
 
 template<class CompType, class ThermoType>
-Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::omegaI
+Foam::scalar Foam::chemistryModel<CompType, ThermoType>::omegaI
 (
     const label index,
     const scalarField& c,
@@ -160,7 +157,7 @@ Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::omegaI
 
 
 template<class CompType, class ThermoType>
-void Foam::ODEChemistryModel<CompType, ThermoType>::updateConcsInReactionI
+void Foam::chemistryModel<CompType, ThermoType>::updateConcsInReactionI
 (
     const label index,
     const scalar dt,
@@ -189,7 +186,7 @@ void Foam::ODEChemistryModel<CompType, ThermoType>::updateConcsInReactionI
 
 
 template<class CompType, class ThermoType>
-void Foam::ODEChemistryModel<CompType, ThermoType>::updateRRInReactionI
+void Foam::chemistryModel<CompType, ThermoType>::updateRRInReactionI
 (
     const label index,
     const scalar pr,
@@ -220,7 +217,7 @@ void Foam::ODEChemistryModel<CompType, ThermoType>::updateRRInReactionI
 
 
 template<class CompType, class ThermoType>
-Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::omega
+Foam::scalar Foam::chemistryModel<CompType, ThermoType>::omega
 (
     const Reaction<ThermoType>& R,
     const scalarField& c,
@@ -338,7 +335,7 @@ Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::omega
 
 
 template<class CompType, class ThermoType>
-void Foam::ODEChemistryModel<CompType, ThermoType>::derivatives
+void Foam::chemistryModel<CompType, ThermoType>::derivatives
 (
     const scalar time,
     const scalarField &c,
@@ -389,7 +386,7 @@ void Foam::ODEChemistryModel<CompType, ThermoType>::derivatives
 
 
 template<class CompType, class ThermoType>
-void Foam::ODEChemistryModel<CompType, ThermoType>::jacobian
+void Foam::chemistryModel<CompType, ThermoType>::jacobian
 (
     const scalar t,
     const scalarField& c,
@@ -532,7 +529,7 @@ void Foam::ODEChemistryModel<CompType, ThermoType>::jacobian
 
 template<class CompType, class ThermoType>
 Foam::tmp<Foam::volScalarField>
-Foam::ODEChemistryModel<CompType, ThermoType>::tc() const
+Foam::chemistryModel<CompType, ThermoType>::tc() const
 {
     scalar pf, cf, pr, cr;
     label lRef, rRef;
@@ -570,6 +567,8 @@ Foam::ODEChemistryModel<CompType, ThermoType>::tc() const
     );
 
     scalarField& tc = ttc();
+    const scalarField& T = this->thermo().T();
+    const scalarField& p = this->thermo().p();
 
     const label nReaction = reactions_.size();
 
@@ -578,8 +577,8 @@ Foam::ODEChemistryModel<CompType, ThermoType>::tc() const
         forAll(rho, celli)
         {
             scalar rhoi = rho[celli];
-            scalar Ti = this->thermo().T()[celli];
+            scalar Ti = T[celli];
-            scalar pi = this->thermo().p()[celli];
+            scalar pi = p[celli];
             scalarField c(nSpecie_);
             scalar cSum = 0.0;
 
@@ -615,7 +614,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::tc() const
 
 template<class CompType, class ThermoType>
 Foam::tmp<Foam::volScalarField>
-Foam::ODEChemistryModel<CompType, ThermoType>::Sh() const
+Foam::chemistryModel<CompType, ThermoType>::Sh() const
 {
     tmp<volScalarField> tSh
     (
@@ -657,7 +656,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::Sh() const
 
 template<class CompType, class ThermoType>
 Foam::tmp<Foam::volScalarField>
-Foam::ODEChemistryModel<CompType, ThermoType>::dQ() const
+Foam::chemistryModel<CompType, ThermoType>::dQ() const
 {
     tmp<volScalarField> tdQ
     (
@@ -689,7 +688,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::dQ() const
 
 
 template<class CompType, class ThermoType>
-Foam::label Foam::ODEChemistryModel<CompType, ThermoType>::nEqns() const
+Foam::label Foam::chemistryModel<CompType, ThermoType>::nEqns() const
 {
     // nEqns = number of species + temperature + pressure
     return nSpecie_ + 2;
@@ -697,7 +696,7 @@ Foam::label Foam::ODEChemistryModel<CompType, ThermoType>::nEqns() const
 
 
 template<class CompType, class ThermoType>
-void Foam::ODEChemistryModel<CompType, ThermoType>::calculate()
+void Foam::chemistryModel<CompType, ThermoType>::calculate()
 {
     if (!this->chemistry_)
     {
@@ -718,11 +717,14 @@ void Foam::ODEChemistryModel<CompType, ThermoType>::calculate()
         this->thermo().rho()
     );
 
+    const scalarField& T = this->thermo().T();
+    const scalarField& p = this->thermo().p();
+
     forAll(rho, celli)
     {
         const scalar rhoi = rho[celli];
-        const scalar Ti = this->thermo().T()[celli];
+        const scalar Ti = T[celli];
-        const scalar pi = this->thermo().p()[celli];
+        const scalar pi = p[celli];
 
         scalarField c(nSpecie_, 0.0);
         for (label i=0; i<nSpecie_; i++)
@@ -742,7 +744,7 @@ void Foam::ODEChemistryModel<CompType, ThermoType>::calculate()
 
 
 template<class CompType, class ThermoType>
-Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::solve
+Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve
 (
     const scalar t0,
     const scalar deltaT
@@ -773,13 +775,16 @@ Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::solve
 
     tmp<volScalarField> thc = this->thermo().hc();
     const scalarField& hc = thc();
+    const scalarField& he = this->thermo().he();
+    const scalarField& T = this->thermo().T();
+    const scalarField& p = this->thermo().p();
 
     forAll(rho, celli)
     {
         const scalar rhoi = rho[celli];
-        const scalar hi = this->thermo().he()[celli] + hc[celli];
+        const scalar hi = he[celli] + hc[celli];
-        const scalar pi = this->thermo().p()[celli];
+        const scalar pi = p[celli];
-        scalar Ti = this->thermo().T()[celli];
+        scalar Ti = T[celli];
 
         scalarField c(nSpecie_, 0.0);
         scalarField c0(nSpecie_, 0.0);
@@ -833,7 +838,7 @@ Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::solve
 
 
 template<class CompType, class ThermoType>
-Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::solve
+Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve
 (
     scalarField &c,
     const scalar T,
@@ -844,7 +849,7 @@ Foam::scalar Foam::ODEChemistryModel<CompType, ThermoType>::solve
 {
     notImplemented
     (
-        "ODEChemistryModel::solve"
+        "chemistryModel::solve"
         "("
             "scalarField&, "
             "const scalar, "

src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.H

@@ -22,7 +22,7 @@ License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
 Class
-    Foam::ODEChemistryModel
+    Foam::chemistryModel
 
 Description
     Extends base chemistry model by adding a thermo package, and ODE functions.
@@ -30,13 +30,13 @@ Description
     terms.
 
 SourceFiles
-    ODEChemistryModelI.H
+    chemistryModelI.H
-    ODEChemistryModel.C
+    chemistryModel.C
 
 \*---------------------------------------------------------------------------*/
 
-#ifndef ODEChemistryModel_H
+#ifndef chemistryModel_H
-#define ODEChemistryModel_H
+#define chemistryModel_H
 
 #include "Reaction.H"
 #include "ODE.H"
@@ -53,11 +53,11 @@ namespace Foam
 class fvMesh;
 
 /*---------------------------------------------------------------------------*\
-                      Class ODEChemistryModel Declaration
+                      Class chemistryModel Declaration
 \*---------------------------------------------------------------------------*/
 
 template<class CompType, class ThermoType>
-class ODEChemistryModel
+class chemistryModel
 :
     public CompType,
     public ODE
@@ -65,10 +65,10 @@ class ODEChemistryModel
     // Private Member Functions
 
         //- Disallow copy constructor
-        ODEChemistryModel(const ODEChemistryModel&);
+        chemistryModel(const chemistryModel&);
 
         //- Disallow default bitwise assignment
-        void operator=(const ODEChemistryModel&);
+        void operator=(const chemistryModel&);
 
 
 protected:
@@ -104,22 +104,17 @@ protected:
 public:
 
     //- Runtime type information
-    TypeName("ODEChemistryModel");
+    TypeName("chemistryModel");
 
 
     // Constructors
 
-        //- Construct from components
+        //- Construct from mesh
-        ODEChemistryModel
+        chemistryModel(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& ODEModelName,
-            const word& thermoTypeName
-        );
 
 
     //- Destructor
-    virtual ~ODEChemistryModel();
+    virtual ~chemistryModel();
 
 
     // Member Functions
@@ -263,12 +258,12 @@ public:
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-#include "ODEChemistryModelI.H"
+#include "chemistryModelI.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 #ifdef NoRepository
-#   include "ODEChemistryModel.C"
+#   include "chemistryModel.C"
 #endif
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModelI.H

@@ -30,7 +30,7 @@ License
 
 template<class CompType, class ThermoType>
 inline Foam::PtrList<Foam::DimensionedField<Foam::scalar, Foam::volMesh> >&
-Foam::ODEChemistryModel<CompType, ThermoType>::RR()
+Foam::chemistryModel<CompType, ThermoType>::RR()
 {
     return RR_;
 }
@@ -38,7 +38,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::RR()
 
 template<class CompType, class ThermoType>
 inline const Foam::PtrList<Foam::Reaction<ThermoType> >&
-Foam::ODEChemistryModel<CompType, ThermoType>::reactions() const
+Foam::chemistryModel<CompType, ThermoType>::reactions() const
 {
     return reactions_;
 }
@@ -46,7 +46,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::reactions() const
 
 template<class CompType, class ThermoType>
 inline const Foam::PtrList<ThermoType>&
-Foam::ODEChemistryModel<CompType, ThermoType>::specieThermo() const
+Foam::chemistryModel<CompType, ThermoType>::specieThermo() const
 {
     return specieThermo_;
 }
@@ -54,7 +54,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::specieThermo() const
 
 template<class CompType, class ThermoType>
 inline Foam::label
-Foam::ODEChemistryModel<CompType, ThermoType>::nSpecie() const
+Foam::chemistryModel<CompType, ThermoType>::nSpecie() const
 {
     return nSpecie_;
 }
@@ -62,7 +62,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::nSpecie() const
 
 template<class CompType, class ThermoType>
 inline Foam::label
-Foam::ODEChemistryModel<CompType, ThermoType>::nReaction() const
+Foam::chemistryModel<CompType, ThermoType>::nReaction() const
 {
     return nReaction_;
 }
@@ -70,7 +70,7 @@ Foam::ODEChemistryModel<CompType, ThermoType>::nReaction() const
 
 template<class CompType, class ThermoType>
 inline const Foam::DimensionedField<Foam::scalar, Foam::volMesh>&
-Foam::ODEChemistryModel<CompType, ThermoType>::RR
+Foam::chemistryModel<CompType, ThermoType>::RR
 (
     const label i
 ) const

src/thermophysicalModels/chemistryModel/chemistryModel/psiChemistryModel/psiChemistryModel.C

@@ -25,7 +25,6 @@ License
 
 #include "psiChemistryModel.H"
 #include "fvMesh.H"
-#include "Time.H"
 
 /* * * * * * * * * * * * * * * private static data * * * * * * * * * * * * * */
 
@@ -35,12 +34,12 @@ namespace Foam
     defineRunTimeSelectionTable(psiChemistryModel, fvMesh);
 }
 
+
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
 Foam::psiChemistryModel::psiChemistryModel
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& thermoTypeName
 )
 :
     basicChemistryModel(mesh),

src/thermophysicalModels/chemistryModel/chemistryModel/psiChemistryModel/psiChemistryModel.H

@@ -78,7 +78,7 @@ protected:
 public:
 
     //- Runtime type information
-    TypeName("psiChemistryModel");
+    TypeName("psi");
 
 
     //- Declare run-time constructor selection tables
@@ -87,19 +87,15 @@ public:
         autoPtr,
         psiChemistryModel,
         fvMesh,
-        (
+        (const fvMesh& mesh),
-            const fvMesh& mesh,
+        (mesh)
-            const word& compTypeName,
-            const word& thermoTypeName
-        ),
-        (mesh, compTypeName, thermoTypeName)
     );
 
 
     // Constructors
 
         //- Construct from mesh
-        psiChemistryModel(const fvMesh& mesh, const word& thermoTypeName);
+        psiChemistryModel(const fvMesh& mesh);
 
 
     //- Selector

src/thermophysicalModels/chemistryModel/chemistryModel/psiChemistryModel/psiChemistryModels.C

@@ -32,7 +32,7 @@ Description
 #include "makeChemistryModel.H"
 
 #include "psiChemistryModel.H"
-#include "ODEChemistryModel.H"
+#include "chemistryModel.H"
 #include "thermoPhysicsTypes.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -41,35 +41,35 @@ namespace Foam
 {
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         psiChemistryModel,
         constGasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         psiChemistryModel,
         gasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         psiChemistryModel,
         constIncompressibleGasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         psiChemistryModel,
         incompressibleGasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         psiChemistryModel,
         icoPoly8ThermoPhysics
     );

src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/rhoChemistryModel.C

@@ -25,7 +25,6 @@ License
 
 #include "rhoChemistryModel.H"
 #include "fvMesh.H"
-#include "Time.H"
 
 /* * * * * * * * * * * * * * * private static data * * * * * * * * * * * * * */
 
@@ -35,12 +34,12 @@ namespace Foam
     defineRunTimeSelectionTable(rhoChemistryModel, fvMesh);
 }
 
+
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
 Foam::rhoChemistryModel::rhoChemistryModel
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& thermoTypeName
 )
 :
     basicChemistryModel(mesh),

src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/rhoChemistryModel.H

@@ -78,7 +78,7 @@ protected:
 public:
 
     //- Runtime type information
-    TypeName("rhoChemistryModel");
+    TypeName("rho");
 
 
     //- Declare run-time constructor selection tables
@@ -88,18 +88,16 @@ public:
         rhoChemistryModel,
         fvMesh,
         (
-            const fvMesh& mesh,
+            const fvMesh& mesh
-            const word& compTypeName,
-            const word& thermoTypeName
         ),
-        (mesh, compTypeName, thermoTypeName)
+        (mesh)
     );
 
 
     // Constructors
 
         //- Construct from mesh
-        rhoChemistryModel(const fvMesh& mesh, const word& thermoTypeName);
+        rhoChemistryModel(const fvMesh& mesh);
 
 
     //- Selector

src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/rhoChemistryModels.C

@@ -32,7 +32,7 @@ Description
 #include "makeChemistryModel.H"
 
 #include "rhoChemistryModel.H"
-#include "ODEChemistryModel.H"
+#include "chemistryModel.H"
 #include "thermoPhysicsTypes.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -41,35 +41,35 @@ namespace Foam
 {
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         rhoChemistryModel,
         constGasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         rhoChemistryModel,
         gasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         rhoChemistryModel,
         constIncompressibleGasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         rhoChemistryModel,
         incompressibleGasThermoPhysics
     );
 
     makeChemistryModel
     (
-        ODEChemistryModel,
+        chemistryModel,
         rhoChemistryModel,
         icoPoly8ThermoPhysics
     );

src/thermophysicalModels/chemistryModel/chemistrySolver/EulerImplicit/EulerImplicit.C

@@ -29,15 +29,13 @@ License
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::EulerImplicit<ODEChemistryType>::EulerImplicit
+Foam::EulerImplicit<ChemistryModel>::EulerImplicit
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& ODEModelName,
-    const word& thermoType
 )
 :
-    chemistrySolver<ODEChemistryType>(mesh, ODEModelName, thermoType),
+    chemistrySolver<ChemistryModel>(mesh),
     coeffsDict_(this->subDict("EulerImplicitCoeffs")),
     cTauChem_(readScalar(coeffsDict_.lookup("cTauChem"))),
     eqRateLimiter_(coeffsDict_.lookup("equilibriumRateLimiter"))
@@ -46,15 +44,15 @@ Foam::EulerImplicit<ODEChemistryType>::EulerImplicit
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::EulerImplicit<ODEChemistryType>::~EulerImplicit()
+Foam::EulerImplicit<ChemistryModel>::~EulerImplicit()
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::scalar Foam::EulerImplicit<ODEChemistryType>::solve
+Foam::scalar Foam::EulerImplicit<ChemistryModel>::solve
 (
     scalarField &c,
     const scalar T,

src/thermophysicalModels/chemistryModel/chemistrySolver/EulerImplicit/EulerImplicit.H

@@ -46,10 +46,10 @@ namespace Foam
                        Class EulerImplicit Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
 class EulerImplicit
 :
-    public chemistrySolver<ODEChemistryType>
+    public chemistrySolver<ChemistryModel>
 {
     // Private data
 
@@ -74,13 +74,8 @@ public:
 
     // Constructors
 
-        //- Construct from components
+        //- Construct from mesh
-        EulerImplicit
+        EulerImplicit(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& ODEmodelName,
-            const word& thermoType
-        );
 
 
     //- Destructor

src/thermophysicalModels/chemistryModel/chemistrySolver/chemistrySolver/chemistrySolver.C

@@ -27,23 +27,20 @@ License
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::chemistrySolver<ODEChemistryType>::chemistrySolver
+Foam::chemistrySolver<ChemistryModel>::chemistrySolver
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& ODEModelName,
-    const word& thermoType
 )
 :
-    ODEChemistryType(mesh, ODEModelName, thermoType),
+    ChemistryModel(mesh)
-    name_(ODEModelName)
 {}
 
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::chemistrySolver<ODEChemistryType>::~chemistrySolver()
+Foam::chemistrySolver<ChemistryModel>::~chemistrySolver()
 {}
 
 

src/thermophysicalModels/chemistryModel/chemistrySolver/chemistrySolver/chemistrySolver.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -35,7 +35,7 @@ SourceFiles
 #ifndef chemistrySolver_H
 #define chemistrySolver_H
 
-#include "ODEChemistryModel.H"
+#include "chemistryModel.H"
 #include "IOdictionary.H"
 #include "scalarField.H"
 
@@ -48,30 +48,18 @@ namespace Foam
                       Class chemistrySolver Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
 class chemistrySolver
 :
-    public ODEChemistryType
+    public ChemistryModel
 {
-protected:
-
-    // Protected data
-
-        //- Name of the chemistry solver
-        const word name_;
-
 
 public:
 
     // Constructors
 
         //- Construct from components
-        chemistrySolver
+        chemistrySolver(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& ODEModelName,
-            const word& thermoType
-        );
 
 
     //- Destructor

src/thermophysicalModels/chemistryModel/chemistrySolver/chemistrySolver/makeChemistrySolverTypes.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -28,7 +28,7 @@ License
 
 #include "chemistrySolver.H"
 
-#include "ODEChemistryModel.H"
+#include "chemistryModel.H"
 
 #include "noChemistrySolver.H"
 #include "EulerImplicit.H"
@@ -37,21 +37,22 @@ License
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-#define makeChemistrySolverType(SS, ODEChem, Comp, Thermo)                    \
+#define makeChemistrySolverType(SS, Comp, Thermo)                             \
                                                                               \
-    typedef SS<ODEChem<Comp, Thermo> > SS##ODEChem##Comp##Thermo;             \
+    typedef SS<chemistryModel<Comp, Thermo> > SS##Comp##Thermo;               \
                                                                               \
     defineTemplateTypeNameAndDebugWithName                                    \
     (                                                                         \
-        SS##ODEChem##Comp##Thermo,                                            \
+        SS##Comp##Thermo,                                                     \
-        (#SS"<"#ODEChem"<"#Comp"," + Thermo::typeName() + ">>").c_str(),      \
+        (#SS"<" + word(Comp::typeName_())                                     \
+      + "," + Thermo::typeName() + ">").c_str(),                              \
         0                                                                     \
     );                                                                        \
                                                                               \
     addToRunTimeSelectionTable                                                \
     (                                                                         \
         Comp,                                                                 \
-        SS##ODEChem##Comp##Thermo,                                            \
+        SS##Comp##Thermo,                                                     \
         fvMesh                                                                \
     );
 
@@ -61,7 +62,6 @@ License
     makeChemistrySolverType                                                   \
     (                                                                         \
         noChemistrySolver,                                                    \
-        ODEChemistryModel,                                                    \
         CompChemModel,                                                        \
         Thermo                                                                \
     );                                                                        \
@@ -69,7 +69,6 @@ License
     makeChemistrySolverType                                                   \
     (                                                                         \
         EulerImplicit,                                                        \
-        ODEChemistryModel,                                                    \
         CompChemModel,                                                        \
         Thermo                                                                \
     );                                                                        \
@@ -77,7 +76,6 @@ License
     makeChemistrySolverType                                                   \
     (                                                                         \
         ode,                                                                  \
-        ODEChemistryModel,                                                    \
         CompChemModel,                                                        \
         Thermo                                                                \
     );                                                                        \
@@ -85,7 +83,6 @@ License
     makeChemistrySolverType                                                   \
     (                                                                         \
         sequential,                                                           \
-        ODEChemistryModel,                                                    \
         CompChemModel,                                                        \
         Thermo                                                                \
     );

src/thermophysicalModels/chemistryModel/chemistrySolver/noChemistrySolver/noChemistrySolver.C

@@ -28,29 +28,27 @@ License
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::noChemistrySolver<ODEChemistryType>::noChemistrySolver
+Foam::noChemistrySolver<ChemistryModel>::noChemistrySolver
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& ODEModelName,
-    const word& thermoType
 )
 :
-    chemistrySolver<ODEChemistryType>(mesh, ODEModelName, thermoType)
+    chemistrySolver<ChemistryModel>(mesh)
 {}
 
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::noChemistrySolver<ODEChemistryType>::~noChemistrySolver()
+Foam::noChemistrySolver<ChemistryModel>::~noChemistrySolver()
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::scalar Foam::noChemistrySolver<ODEChemistryType>::solve
+Foam::scalar Foam::noChemistrySolver<ChemistryModel>::solve
 (
     scalarField&,
     const scalar,

src/thermophysicalModels/chemistryModel/chemistrySolver/noChemistrySolver/noChemistrySolver.H

@@ -47,10 +47,10 @@ namespace Foam
                      Class noChemistrySolver Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
 class noChemistrySolver
 :
-    public chemistrySolver<ODEChemistryType>
+    public chemistrySolver<ChemistryModel>
 {
 
 public:
@@ -61,14 +61,8 @@ public:
 
     // Constructors
 
-
         //- Construct from components
-        noChemistrySolver
+        noChemistrySolver(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& ODEModelName,
-            const word& thermoType
-        );
 
 
     //- Destructor

src/thermophysicalModels/chemistryModel/chemistrySolver/ode/ode.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -24,19 +24,17 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "ode.H"
-#include "ODEChemistryModel.H"
+#include "chemistryModel.H"
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::ode<ODEChemistryType>::ode
+Foam::ode<ChemistryModel>::ode
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& ODEModelName,
-    const word& thermoType
 )
 :
-    chemistrySolver<ODEChemistryType>(mesh, ODEModelName, thermoType),
+    chemistrySolver<ChemistryModel>(mesh),
     coeffsDict_(this->subDict("odeCoeffs")),
     solverName_(coeffsDict_.lookup("solver")),
     odeSolver_(ODESolver::New(solverName_, *this)),
@@ -46,15 +44,15 @@ Foam::ode<ODEChemistryType>::ode
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::ode<ODEChemistryType>::~ode()
+Foam::ode<ChemistryModel>::~ode()
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::scalar Foam::ode<ODEChemistryType>::solve
+Foam::scalar Foam::ode<ChemistryModel>::solve
 (
     scalarField& c,
     const scalar T,

src/thermophysicalModels/chemistryModel/chemistrySolver/ode/ode.H

@@ -47,10 +47,10 @@ namespace Foam
                             Class ode Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
 class ode
 :
-    public chemistrySolver<ODEChemistryType>
+    public chemistrySolver<ChemistryModel>
 {
     // Private data
 
@@ -66,18 +66,13 @@ class ode
 public:
 
     //- Runtime type information
-    TypeName("ODE");
+    TypeName("ode");
 
 
     // Constructors
 
-        //- Construct from components
+        //- Construct from mesh
-        ode
+        ode(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& ODEModeNewlName,
-            const word& thermoType
-        );
 
 
     //- Destructor

src/thermophysicalModels/chemistryModel/chemistrySolver/sequential/sequential.C

@@ -28,15 +28,13 @@ License
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::sequential<ODEChemistryType>::sequential
+Foam::sequential<ChemistryModel>::sequential
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& ODEModelName,
-    const word& thermoType
 )
 :
-    chemistrySolver<ODEChemistryType>(mesh, ODEModelName, thermoType),
+    chemistrySolver<ChemistryModel>(mesh),
     coeffsDict_(this->subDict("sequentialCoeffs")),
     cTauChem_(readScalar(coeffsDict_.lookup("cTauChem"))),
     eqRateLimiter_(coeffsDict_.lookup("equilibriumRateLimiter"))
@@ -45,15 +43,15 @@ Foam::sequential<ODEChemistryType>::sequential
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::sequential<ODEChemistryType>::~sequential()
+Foam::sequential<ChemistryModel>::~sequential()
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
-Foam::scalar Foam::sequential<ODEChemistryType>::solve
+Foam::scalar Foam::sequential<ChemistryModel>::solve
 (
     scalarField &c,
     const scalar T,

src/thermophysicalModels/chemistryModel/chemistrySolver/sequential/sequential.H

@@ -48,17 +48,16 @@ namespace Foam
                          Class sequential Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class ODEChemistryType>
+template<class ChemistryModel>
 class sequential
 :
-    public chemistrySolver<ODEChemistryType>
+    public chemistrySolver<ChemistryModel>
 {
     // Private data
 
         //- Coefficients dictionary
         dictionary coeffsDict_;
 
-
         // Model constants
 
             //- Chemistry time scale
@@ -76,13 +75,8 @@ public:
 
     // Constructors
 
-        //- Construct from components
+        //- Construct from mesh
-        sequential
+        sequential(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& ODEModelName,
-            const word& thermoType
-        );
 
 
     //- Destructor

src/thermophysicalModels/radiationModels/derivedFvPatchFields/radiationCoupledBase/radiationCoupledBase.C

@@ -149,6 +149,7 @@ Foam::scalarField Foam::radiationCoupledBase::emissivity() const
                     "radiationProperties"
                 );
 
+
             // Force recalculation of mapping and schedule
             const mapDistribute& distMap = mpp.map();
 
@@ -157,6 +158,7 @@ Foam::scalarField Foam::radiationCoupledBase::emissivity() const
                 nbrMesh
             ).boundary()[mpp.samplePolyPatch().index()];
 
+
             scalarField emissivity
             (
                 radiation.absorptionEmission().e()().boundaryField()
@@ -182,8 +184,7 @@ Foam::scalarField Foam::radiationCoupledBase::emissivity() const
             FatalErrorIn
             (
                 "radiationCoupledBase::emissivity(const scalarField&)"
-            )
+            )   << "Unimplemented method " << method_ << endl
-                << "Unimplemented method " << method_ << endl
                 << "Please set 'emissivity' to one of "
                 << emissivityMethodTypeNames_.toc()
                 << " and 'emissivityName' to the name of the volScalar"
@@ -191,6 +192,7 @@ Foam::scalarField Foam::radiationCoupledBase::emissivity() const
         }
         break;
     }
+
     return scalarField(0);
 }
 

src/thermophysicalModels/reactionThermo/Make/files

@@ -4,14 +4,14 @@ chemistryReaders/chemistryReader/makeChemistryReaders.C
 
 mixtures/basicMultiComponentMixture/basicMultiComponentMixture.C
 
-psiReactionThermo/psiReactionThermo/psiReactionThermo.C
+psiReactionThermo/psiReactionThermo.C
-psiReactionThermo/hePsiReactionThermo/hePsiReactionThermos.C
+psiReactionThermo/psiReactionThermos.C
 
-psiReactionThermo/psiuReactionThermo/psiuReactionThermo.C
+psiuReactionThermo/psiuReactionThermo.C
-psiReactionThermo/heheuReactionThermo/heheuReactionThermos.C
+psiuReactionThermo/psiuReactionThermos.C
 
-rhoReactionThermo/rhoReactionThermo/rhoReactionThermo.C
+rhoReactionThermo/rhoReactionThermo.C
-rhoReactionThermo/heRhoReactionThermo/heRhoReactionThermos.C
+rhoReactionThermo/rhoReactionThermos.C
 
 derivedFvPatchFields/fixedUnburntEnthalpy/fixedUnburntEnthalpyFvPatchScalarField.C
 derivedFvPatchFields/gradientUnburntEnthalpy/gradientUnburntEnthalpyFvPatchScalarField.C

src/thermophysicalModels/reactionThermo/makeReactionThermo.H

@@ -35,6 +35,7 @@ License
                                                                                \
 typedef MixtureThermo                                                          \
 <                                                                              \
+    CThermo,                                                                   \
     SpecieMixture                                                              \
     <                                                                          \
         Mixture                                                                \

src/thermophysicalModels/reactionThermo/psiReactionThermo/hePsiReactionThermo/hePsiReactionThermo.C

@@ -1,147 +0,0 @@
-/*---------------------------------------------------------------------------*\
-  =========                 |
-  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
-   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
-     \\/     M anipulation  |
--------------------------------------------------------------------------------
-License
-    This file is part of OpenFOAM.
-
-    OpenFOAM is free software: you can redistribute it and/or modify it
-    under the terms of the GNU General Public License as published by
-    the Free Software Foundation, either version 3 of the License, or
-    (at your option) any later version.
-
-    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
-    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
-    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
-    for more details.
-
-    You should have received a copy of the GNU General Public License
-    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
-
-\*---------------------------------------------------------------------------*/
-
-#include "hePsiReactionThermo.H"
-#include "fvMesh.H"
-#include "fixedValueFvPatchFields.H"
-
-// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
-
-template<class MixtureType>
-void Foam::hePsiReactionThermo<MixtureType>::calculate()
-{
-    const scalarField& hCells = this->he_.internalField();
-    const scalarField& pCells = this->p_.internalField();
-
-    scalarField& TCells = this->T_.internalField();
-    scalarField& psiCells = this->psi_.internalField();
-    scalarField& muCells = this->mu_.internalField();
-    scalarField& alphaCells = this->alpha_.internalField();
-
-    forAll(TCells, celli)
-    {
-        const typename MixtureType::thermoType& mixture =
-            this->cellMixture(celli);
-
-        TCells[celli] = mixture.THE
-        (
-            hCells[celli],
-            pCells[celli],
-            TCells[celli]
-        );
-
-        psiCells[celli] = mixture.psi(pCells[celli], TCells[celli]);
-
-        muCells[celli] = mixture.mu(pCells[celli], TCells[celli]);
-        alphaCells[celli] = mixture.alphah(pCells[celli], TCells[celli]);
-    }
-
-    forAll(this->T_.boundaryField(), patchi)
-    {
-        fvPatchScalarField& pp = this->p_.boundaryField()[patchi];
-        fvPatchScalarField& pT = this->T_.boundaryField()[patchi];
-        fvPatchScalarField& ppsi = this->psi_.boundaryField()[patchi];
-
-        fvPatchScalarField& ph = this->he_.boundaryField()[patchi];
-
-        fvPatchScalarField& pmu_ = this->mu_.boundaryField()[patchi];
-        fvPatchScalarField& palpha_ = this->alpha_.boundaryField()[patchi];
-
-        if (pT.fixesValue())
-        {
-            forAll(pT, facei)
-            {
-                const typename MixtureType::thermoType& mixture =
-                    this->patchFaceMixture(patchi, facei);
-
-                ph[facei] = mixture.HE(pp[facei], pT[facei]);
-
-                ppsi[facei] = mixture.psi(pp[facei], pT[facei]);
-                pmu_[facei] = mixture.mu(pp[facei], pT[facei]);
-                palpha_[facei] = mixture.alphah(pp[facei], pT[facei]);
-            }
-        }
-        else
-        {
-            forAll(pT, facei)
-            {
-                const typename MixtureType::thermoType& mixture =
-                    this->patchFaceMixture(patchi, facei);
-
-                pT[facei] = mixture.THE(ph[facei], pp[facei], pT[facei]);
-
-                ppsi[facei] = mixture.psi(pp[facei], pT[facei]);
-                pmu_[facei] = mixture.mu(pp[facei], pT[facei]);
-                palpha_[facei] = mixture.alphah(pp[facei], pT[facei]);
-            }
-        }
-    }
-}
-
-
-// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
-
-template<class MixtureType>
-Foam::hePsiReactionThermo<MixtureType>::hePsiReactionThermo(const fvMesh& mesh)
-:
-    heThermo<psiReactionThermo, MixtureType>(mesh)
-{
-    calculate();
-
-    // Switch on saving old time
-    this->psi_.oldTime();
-}
-
-
-// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
-
-template<class MixtureType>
-Foam::hePsiReactionThermo<MixtureType>::~hePsiReactionThermo()
-{}
-
-
-// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
-
-template<class MixtureType>
-void Foam::hePsiReactionThermo<MixtureType>::correct()
-{
-    if (debug)
-    {
-        Info<< "entering hePsiReactionThermo<MixtureType>::correct()" << endl;
-    }
-
-    // force the saving of the old-time values
-    this->psi_.oldTime();
-
-    calculate();
-
-    if (debug)
-    {
-        Info<< "exiting hePsiReactionThermo<MixtureType>::correct()" << endl;
-    }
-}
-
-
-// ************************************************************************* //

src/thermophysicalModels/reactionThermo/psiReactionThermo/psiReactionThermos.C

@@ -26,7 +26,7 @@ License
 #include "makeReactionThermo.H"
 
 #include "psiReactionThermo.H"
-#include "hePsiReactionThermo.H"
+#include "hePsiThermo.H"
 
 #include "specie.H"
 #include "perfectGas.H"
@@ -59,7 +59,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     homogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -72,7 +72,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     inhomogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -85,7 +85,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     veryInhomogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -101,7 +101,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     homogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -114,7 +114,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     inhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -127,7 +127,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     veryInhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -143,7 +143,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     homogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -156,7 +156,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     inhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -169,7 +169,7 @@ makeReactionThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     veryInhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -185,7 +185,7 @@ makeReactionMixtureThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     multiComponentMixture,
     constGasThermoPhysics
 );
@@ -194,7 +194,7 @@ makeReactionMixtureThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     multiComponentMixture,
     gasThermoPhysics
 );
@@ -206,7 +206,7 @@ makeReactionMixtureThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     reactingMixture,
     constGasThermoPhysics
 );
@@ -215,7 +215,7 @@ makeReactionMixtureThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     reactingMixture,
     gasThermoPhysics
 );
@@ -224,7 +224,7 @@ makeReactionMixtureThermo
 (
     psiThermo,
     psiReactionThermo,
-    hePsiReactionThermo,
+    hePsiThermo,
     singleStepReactingMixture,
     gasThermoPhysics
 );

src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C

@@ -23,14 +23,14 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-#include "heheuReactionThermo.H"
+#include "heheuPsiThermo.H"
 #include "fvMesh.H"
 #include "fixedValueFvPatchFields.H"
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-void Foam::heheuReactionThermo<MixtureType>::calculate()
+void Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::calculate()
 {
     const scalarField& hCells = this->he_.internalField();
     const scalarField& heuCells = this->heu_.internalField();
@@ -118,8 +118,11 @@ void Foam::heheuReactionThermo<MixtureType>::calculate()
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-Foam::heheuReactionThermo<MixtureType>::heheuReactionThermo(const fvMesh& mesh)
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::heheuPsiThermo
+(
+    const fvMesh& mesh
+)
 :
     heThermo<psiuReactionThermo, MixtureType>(mesh),
     Tu_
@@ -188,19 +191,20 @@ Foam::heheuReactionThermo<MixtureType>::heheuReactionThermo(const fvMesh& mesh)
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-Foam::heheuReactionThermo<MixtureType>::~heheuReactionThermo()
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::~heheuPsiThermo()
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-void Foam::heheuReactionThermo<MixtureType>::correct()
+void Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::correct()
 {
     if (debug)
     {
-        Info<< "entering heheuReactionThermo<MixtureType>::correct()" << endl;
+        Info<< "entering heheuPsiThermo<BasicPsiThermo, MixtureType>::correct()"
+            << endl;
     }
 
     // force the saving of the old-time values
@@ -210,14 +214,15 @@ void Foam::heheuReactionThermo<MixtureType>::correct()
 
     if (debug)
     {
-        Info<< "exiting heheuReactionThermo<MixtureType>::correct()" << endl;
+        Info<< "exiting heheuPsiThermo<BasicPsiThermo, MixtureType>::correct()"
+            << endl;
     }
 }
 
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 Foam::tmp<Foam::scalarField>
-Foam::heheuReactionThermo<MixtureType>::heu
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::heu
 (
     const scalarField& p,
     const scalarField& Tu,
@@ -236,9 +241,9 @@ Foam::heheuReactionThermo<MixtureType>::heu
 }
 
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 Foam::tmp<Foam::scalarField>
-Foam::heheuReactionThermo<MixtureType>::heu
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::heu
 (
     const scalarField& p,
     const scalarField& Tu,
@@ -258,9 +263,9 @@ Foam::heheuReactionThermo<MixtureType>::heu
 }
 
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 Foam::tmp<Foam::volScalarField>
-Foam::heheuReactionThermo<MixtureType>::Tb() const
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::Tb() const
 {
     tmp<volScalarField> tTb
     (
@@ -314,9 +319,9 @@ Foam::heheuReactionThermo<MixtureType>::Tb() const
 }
 
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 Foam::tmp<Foam::volScalarField>
-Foam::heheuReactionThermo<MixtureType>::psiu() const
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::psiu() const
 {
     tmp<volScalarField> tpsiu
     (
@@ -365,9 +370,9 @@ Foam::heheuReactionThermo<MixtureType>::psiu() const
 }
 
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 Foam::tmp<Foam::volScalarField>
-Foam::heheuReactionThermo<MixtureType>::psib() const
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::psib() const
 {
     tmp<volScalarField> tpsib
     (
@@ -417,9 +422,9 @@ Foam::heheuReactionThermo<MixtureType>::psib() const
 }
 
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 Foam::tmp<Foam::volScalarField>
-Foam::heheuReactionThermo<MixtureType>::muu() const
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::muu() const
 {
     tmp<volScalarField> tmuu
     (
@@ -472,9 +477,9 @@ Foam::heheuReactionThermo<MixtureType>::muu() const
 }
 
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
 Foam::tmp<Foam::volScalarField>
-Foam::heheuReactionThermo<MixtureType>::mub() const
+Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::mub() const
 {
     tmp<volScalarField> tmub
     (

src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.H

@@ -36,7 +36,6 @@ SourceFiles
 #define heheuReactionThermo_H
 
 #include "heThermo.H"
-#include "psiuReactionThermo.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -44,13 +43,13 @@ namespace Foam
 {
 
 /*---------------------------------------------------------------------------*\
-                      Class heheuReactionThermo Declaration
+                      Class heheuPsiThermo Declaration
 \*---------------------------------------------------------------------------*/
 
-template<class MixtureType>
+template<class BasicPsiThermo, class MixtureType>
-class heheuReactionThermo
+class heheuPsiThermo
 :
-    public heThermo<psiuReactionThermo, MixtureType>
+    public heThermo<BasicPsiThermo, MixtureType>
 {
     // Private Member Functions
 
@@ -63,23 +62,26 @@ class heheuReactionThermo
         void calculate();
 
         //- Construct as copy (not implemented)
-        heheuReactionThermo(const heheuReactionThermo<MixtureType>&);
+        heheuPsiThermo
+        (
+            const heheuPsiThermo<BasicPsiThermo, MixtureType>&
+        );
 
 
 public:
 
     //- Runtime type information
-    TypeName("heheuReactionThermo");
+    TypeName("heheuPsiThermo");
 
 
     // Constructors
 
         //- Construct from mesh
-        heheuReactionThermo(const fvMesh&);
+        heheuPsiThermo(const fvMesh&);
 
 
     //- Destructor
-    virtual ~heheuReactionThermo();
+    virtual ~heheuPsiThermo();
 
 
     // Member functions
@@ -155,7 +157,7 @@ public:
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
 
 #ifdef NoRepository
-#   include "heheuReactionThermo.C"
+#   include "heheuPsiThermo.C"
 #endif
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

src/thermophysicalModels/reactionThermo/psiuReactionThermo/psiuReactionThermos.C

@@ -24,7 +24,7 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "psiuReactionThermo.H"
-#include "heheuReactionThermo.H"
+#include "heheuPsiThermo.H"
 
 #include "makeReactionThermo.H"
 #include "addToRunTimeSelectionTable.H"
@@ -56,7 +56,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     homogeneousMixture,
     constTransport,
     absoluteEnthalpy,
@@ -69,7 +69,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     inhomogeneousMixture,
     constTransport,
     absoluteEnthalpy,
@@ -82,7 +82,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     veryInhomogeneousMixture,
     constTransport,
     absoluteEnthalpy,
@@ -95,7 +95,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     homogeneousMixture,
     sutherlandTransport,
     absoluteEnthalpy,
@@ -108,7 +108,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     inhomogeneousMixture,
     sutherlandTransport,
     absoluteEnthalpy,
@@ -121,7 +121,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     veryInhomogeneousMixture,
     sutherlandTransport,
     absoluteEnthalpy,
@@ -134,7 +134,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     egrMixture,
     constTransport,
     absoluteEnthalpy,
@@ -147,7 +147,7 @@ makeReactionThermo
 (
     psiThermo,
     psiuReactionThermo,
-    heheuReactionThermo,
+    heheuPsiThermo,
     egrMixture,
     sutherlandTransport,
     absoluteEnthalpy,

src/thermophysicalModels/reactionThermo/rhoReactionThermo/heRhoReactionThermo/heRhoReactionThermo.C

@@ -1,145 +0,0 @@
-/*---------------------------------------------------------------------------*\
-  =========                 |
-  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
-   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
-     \\/     M anipulation  |
--------------------------------------------------------------------------------
-License
-    This file is part of OpenFOAM.
-
-    OpenFOAM is free software: you can redistribute it and/or modify it
-    under the terms of the GNU General Public License as published by
-    the Free Software Foundation, either version 3 of the License, or
-    (at your option) any later version.
-
-    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
-    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
-    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
-    for more details.
-
-    You should have received a copy of the GNU General Public License
-    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
-
-\*---------------------------------------------------------------------------*/
-
-#include "heRhoReactionThermo.H"
-#include "fvMesh.H"
-#include "fixedValueFvPatchFields.H"
-
-// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
-
-template<class MixtureType>
-void Foam::heRhoReactionThermo<MixtureType>::calculate()
-{
-    const scalarField& hCells = this->he_.internalField();
-    const scalarField& pCells = this->p_.internalField();
-
-    scalarField& TCells = this->T_.internalField();
-    scalarField& psiCells = this->psi_.internalField();
-    scalarField& rhoCells = this->rho_.internalField();
-    scalarField& muCells = this->mu_.internalField();
-    scalarField& alphaCells = this->alpha_.internalField();
-
-    forAll(TCells, celli)
-    {
-        const typename MixtureType::thermoType& mixture =
-            this->cellMixture(celli);
-
-        TCells[celli] = mixture.THE
-        (
-            hCells[celli],
-            pCells[celli],
-            TCells[celli]
-        );
-        psiCells[celli] = mixture.psi(pCells[celli], TCells[celli]);
-        rhoCells[celli] = mixture.rho(pCells[celli], TCells[celli]);
-
-        muCells[celli] = mixture.mu(pCells[celli], TCells[celli]);
-        alphaCells[celli] = mixture.alphah(pCells[celli], TCells[celli]);
-    }
-
-    forAll(this->T_.boundaryField(), patchi)
-    {
-        fvPatchScalarField& pp = this->p_.boundaryField()[patchi];
-        fvPatchScalarField& pT = this->T_.boundaryField()[patchi];
-        fvPatchScalarField& ppsi = this->psi_.boundaryField()[patchi];
-        fvPatchScalarField& prho = this->rho_.boundaryField()[patchi];
-
-        fvPatchScalarField& ph = this->he_.boundaryField()[patchi];
-
-        fvPatchScalarField& pmu_ = this->mu_.boundaryField()[patchi];
-        fvPatchScalarField& palpha_ = this->alpha_.boundaryField()[patchi];
-
-        if (pT.fixesValue())
-        {
-            forAll(pT, facei)
-            {
-                const typename MixtureType::thermoType& mixture =
-                    this->patchFaceMixture(patchi, facei);
-
-                ph[facei] = mixture.HE(pp[facei], pT[facei]);
-
-                ppsi[facei] = mixture.psi(pp[facei], pT[facei]);
-                prho[facei] = mixture.rho(pp[facei], pT[facei]);
-                pmu_[facei] = mixture.mu(pp[facei], pT[facei]);
-                palpha_[facei] = mixture.alphah(pp[facei], pT[facei]);
-            }
-        }
-        else
-        {
-            forAll(pT, facei)
-            {
-                const typename MixtureType::thermoType& mixture =
-                    this->patchFaceMixture(patchi, facei);
-
-                pT[facei] = mixture.THE(ph[facei], pp[facei], pT[facei]);
-
-                ppsi[facei] = mixture.psi(pp[facei], pT[facei]);
-                prho[facei] = mixture.rho(pp[facei], pT[facei]);
-                pmu_[facei] = mixture.mu(pp[facei], pT[facei]);
-                palpha_[facei] = mixture.alphah(pp[facei], pT[facei]);
-            }
-        }
-    }
-}
-
-
-// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
-
-template<class MixtureType>
-Foam::heRhoReactionThermo<MixtureType>::heRhoReactionThermo(const fvMesh& mesh)
-:
-    heThermo<rhoReactionThermo, MixtureType>(mesh)
-{
-    calculate();
-}
-
-
-// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
-
-template<class MixtureType>
-Foam::heRhoReactionThermo<MixtureType>::~heRhoReactionThermo()
-{}
-
-
-// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
-
-template<class MixtureType>
-void Foam::heRhoReactionThermo<MixtureType>::correct()
-{
-    if (debug)
-    {
-        Info<< "entering heRhoReactionThermo<MixtureType>::correct()" << endl;
-    }
-
-    calculate();
-
-    if (debug)
-    {
-        Info<< "exiting heRhoReactionThermo<MixtureType>::correct()" << endl;
-    }
-}
-
-
-// ************************************************************************* //

src/thermophysicalModels/reactionThermo/rhoReactionThermo/heRhoReactionThermo/heRhoReactionThermo.H

@@ -1,100 +0,0 @@
-/*---------------------------------------------------------------------------*\
-  =========                 |
-  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
-   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
-     \\/     M anipulation  |
--------------------------------------------------------------------------------
-License
-    This file is part of OpenFOAM.
-
-    OpenFOAM is free software: you can redistribute it and/or modify it
-    under the terms of the GNU General Public License as published by
-    the Free Software Foundation, either version 3 of the License, or
-    (at your option) any later version.
-
-    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
-    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
-    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
-    for more details.
-
-    You should have received a copy of the GNU General Public License
-    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
-
-Class
-    Foam::heRhoReactionThermo
-
-Description
-    Foam::heRhoReactionThermo
-
-SourceFiles
-    heRhoReactionThermo.C
-
-\*---------------------------------------------------------------------------*/
-
-#ifndef heRhoReactionThermo_H
-#define heRhoReactionThermo_H
-
-#include "heThermo.H"
-#include "rhoReactionThermo.H"
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-namespace Foam
-{
-
-/*---------------------------------------------------------------------------*\
-                      Class heRhoReactionThermo Declaration
-\*---------------------------------------------------------------------------*/
-
-template<class MixtureType>
-class heRhoReactionThermo
-:
-    public heThermo<rhoReactionThermo, MixtureType>
-{
-    // Private Member Functions
-
-        void calculate();
-
-        //- Construct as copy (not implemented)
-        heRhoReactionThermo(const heRhoReactionThermo<MixtureType>&);
-
-
-public:
-
-    //- Runtime type information
-    TypeName("heRhoReactionThermo");
-
-
-    // Constructors
-
-        //- Construct from mesh
-        heRhoReactionThermo(const fvMesh&);
-
-
-    //- Destructor
-    virtual ~heRhoReactionThermo();
-
-
-    // Member functions
-
-        //- Update properties
-        virtual void correct();
-};
-
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-} // End namespace Foam
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-#ifdef NoRepository
-#   include "heRhoReactionThermo.C"
-#endif
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-#endif
-
-// ************************************************************************* //

src/thermophysicalModels/reactionThermo/rhoReactionThermo/rhoReactionThermos.C

@@ -26,7 +26,7 @@ License
 #include "makeReactionThermo.H"
 
 #include "rhoReactionThermo.H"
-#include "heRhoReactionThermo.H"
+#include "heRhoThermo.H"
 
 #include "specie.H"
 #include "perfectGas.H"
@@ -59,7 +59,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     homogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -72,7 +72,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     inhomogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -85,7 +85,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     veryInhomogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -98,7 +98,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     homogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -111,7 +111,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     inhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -124,7 +124,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     veryInhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -138,7 +138,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     homogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -151,7 +151,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     inhomogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -164,7 +164,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     veryInhomogeneousMixture,
     constTransport,
     sensibleEnthalpy,
@@ -177,7 +177,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     homogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -190,7 +190,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     inhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -203,7 +203,7 @@ makeReactionThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     veryInhomogeneousMixture,
     sutherlandTransport,
     sensibleEnthalpy,
@@ -219,7 +219,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     multiComponentMixture,
     constGasThermoPhysics
 );
@@ -228,7 +228,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     multiComponentMixture,
     gasThermoPhysics
 );
@@ -237,7 +237,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     multiComponentMixture,
     constIncompressibleGasThermoPhysics
 );
@@ -246,7 +246,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     multiComponentMixture,
     incompressibleGasThermoPhysics
 );
@@ -255,7 +255,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     multiComponentMixture,
     icoPoly8ThermoPhysics
 );
@@ -267,7 +267,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     reactingMixture,
     constGasThermoPhysics
 );
@@ -276,7 +276,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     reactingMixture,
     gasThermoPhysics
 );
@@ -285,7 +285,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     reactingMixture,
     constIncompressibleGasThermoPhysics
 );
@@ -294,7 +294,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     reactingMixture,
     incompressibleGasThermoPhysics
 );
@@ -303,7 +303,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     reactingMixture,
     icoPoly8ThermoPhysics
 );
@@ -312,7 +312,7 @@ makeReactionMixtureThermo
 (
     rhoThermo,
     rhoReactionThermo,
-    heRhoReactionThermo,
+    heRhoThermo,
     singleStepReactingMixture,
     gasThermoPhysics
 );

src/thermophysicalModels/solidChemistryModel/ODESolidChemistryModel/ODESolidChemistryModel.C

@@ -32,12 +32,10 @@ template<class CompType, class SolidThermo, class GasThermo>
 Foam::ODESolidChemistryModel<CompType, SolidThermo, GasThermo>::
 ODESolidChemistryModel
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& compTypeName,
-    const word& solidThermoName
 )
 :
-    CompType(mesh, solidThermoName),
+    CompType(mesh),
     ODE(),
     Ys_(this->solid().composition().Y()),
     pyrolisisGases_

src/thermophysicalModels/solidChemistryModel/ODESolidChemistryModel/ODESolidChemistryModel.H

@@ -138,13 +138,8 @@ public:
 
     // Constructors
 
-        //- Construct from components
+        //- Construct from mesh
-        ODESolidChemistryModel
+        ODESolidChemistryModel(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& compTypeName,
-            const word& SolidThermoName
-        );
 
 
     //- Destructor

src/thermophysicalModels/solidChemistryModel/solidChemistryModel/solidChemistryModel.C

@@ -39,12 +39,11 @@ namespace Foam
 
 Foam::solidChemistryModel::solidChemistryModel
 (
-    const fvMesh& mesh,
+    const fvMesh& mesh
-    const word& solidThermoTypeName
 )
 :
     basicChemistryModel(mesh),
-    solidThermo_(solidReactionThermo::New(mesh)) //, thermoTypeName))
+    solidThermo_(solidReactionThermo::New(mesh))
 {}
 
 

src/thermophysicalModels/solidChemistryModel/solidChemistryModel/solidChemistryModel.H

@@ -88,23 +88,15 @@ public:
         autoPtr,
         solidChemistryModel,
         fvMesh,
-        (
+        (const fvMesh& mesh),
-            const fvMesh& mesh,
+        (mesh)
-            const word& compTypeName,
-            const word& solidThermoTypeName
-        ),
-        (mesh, compTypeName, solidThermoTypeName)
     );
 
 
     // Constructors
 
-        //- Construct from mesh and thermo type name
+        //- Construct from mesh
-        solidChemistryModel
+        solidChemistryModel(const fvMesh& mesh);
-        (
-            const fvMesh& mesh,
-            const word& solidThermoTypeName
-        );
 
 
     //- Selector

src/thermophysicalModels/solidChemistryModel/solidChemistryModel/solidChemistryModelNew.C

@@ -104,8 +104,7 @@ Foam::autoPtr<Foam::solidChemistryModel> Foam::solidChemistryModel::New
         }
     }
 
-    return autoPtr<solidChemistryModel>
+    return autoPtr<solidChemistryModel>(cstrIter()(mesh));
-        (cstrIter()(mesh, ODEModelName, thermoTypeName));
 }
 
 

src/thermophysicalModels/solidSpecie/Make/files

@@ -1,4 +1,4 @@
 reaction/Reactions/solidReaction/solidReaction.C
-reaction/reactions/makeSolidReactionThermoReactions.C
+reaction/reactions/makeSolidReactions.C
 
 LIB = $(FOAM_LIBBIN)/libsolidSpecie

src/thermophysicalModels/solidSpecie/reaction/reactions/makeSolidReaction.H

@@ -22,10 +22,10 @@ License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
 InClass
-    Foam::makeSolidReactionThermo
+    Foam::makeSolidReaction
 
 Description
-    Macros for instantiating reactions on given solid thermo packages
+    Macros for instantiating reactions for solid
 
 \*---------------------------------------------------------------------------*/
 
@@ -79,5 +79,4 @@ namespace Foam
 
 #endif
 
-
 // ************************************************************************* //

src/thermophysicalModels/solidSpecie/reaction/reactions/makeSolidReactions.C

@@ -23,7 +23,7 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-#include "makeSolidReactionThermo.H"
+#include "makeSolidReaction.H"
 #include "solidArrheniusReactionRate.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

src/thermophysicalModels/solidThermo/Make/files

@@ -1,10 +1,11 @@
 solidThermo/solidThermo.C
 solidThermo/solidThermoNew.C
+solidThermo/solidThermos.C
 
 mixtures/basicSolidMixture/basicSolidMixture.C
 
 solidReactionThermo/solidReactionThermo.C
 solidReactionThermo/solidReactionThermoNew.C
-heSolidThermo/heSolidThermos.C
+solidReactionThermo/solidReactionThermos.C
 
 LIB = $(FOAM_LIBBIN)/libsolidThermo

src/thermophysicalModels/solidThermo/mixtures/basicSolidMixture/basicSolidMixture.H

@@ -191,21 +191,9 @@ public:
             //- Density
             virtual scalar rho(scalar p, scalar T, label celli) const = 0;
 
-            //- Absorption coefficient
-            virtual scalar kappaRad(scalar p, scalar T, label celli) const = 0;
-
-            //- Scatter coefficient
-            virtual scalar sigmaS(scalar p, scalar T, label celli) const = 0;
-
             //- Thermal conductivity
             virtual scalar kappa(scalar p, scalar T, label celli) const = 0;
 
-            //- Emissivity coefficient
-            virtual scalar emissivity
-            (
-                scalar p, scalar T, label celli
-            ) const = 0;
-
             //- Specific heat capacity
             virtual scalar Cp(scalar p, scalar T, label celli) const = 0;
 };
@@ -217,7 +205,7 @@ public:
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-#   include "basicSolidMixtureI.H"
+#include "basicSolidMixtureI.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 

src/thermophysicalModels/solidThermo/mixtures/basicSolidMixture/basicSolidMixtures.C

@@ -29,7 +29,6 @@ Description
 #include "basicMixture.H"
 #include "makeBasicMixture.H"
 
-
 #include "rhoConst.H"
 
 #include "hConstThermo.H"
@@ -39,8 +38,6 @@ Description
 #include "constAnIsoSolidTransport.H"
 #include "exponentialSolidTransport.H"
 
-#include "constSolidRad.H"
-
 #include "sensibleInternalEnergy.H"
 #include "sensibleEnthalpy.H"
 #include "thermo.H"

src/thermophysicalModels/solidThermo/mixtures/multiComponentSolidMixture/multiComponentSolidMixture.C

@@ -316,36 +316,6 @@ Foam::scalar Foam::multiComponentSolidMixture<ThermoType>::rho
 }
 
 
-template<class ThermoType>
-Foam::scalar Foam::multiComponentSolidMixture<ThermoType>::kappaRad
-(
-    scalar p, scalar T, label celli
-) const
-{
-    scalar tmp = 0.0;
-    forAll(solidData_, i)
-    {
-        tmp += solidData_[i].kappaRad(T)*X(i, celli, p, T);
-    }
-    return tmp;
-}
-
-
-template<class ThermoType>
-Foam::scalar Foam::multiComponentSolidMixture<ThermoType>::sigmaS
-(
-    scalar p, scalar T, label celli
-) const
-{
-    scalar tmp = 0.0;
-    forAll(solidData_, i)
-    {
-        tmp += solidData_[i].sigmaS(T)*X(i, celli, p, T);
-    }
-    return tmp;
-}
-
-
 template<class ThermoType>
 Foam::scalar Foam::multiComponentSolidMixture<ThermoType>::kappa
 (
@@ -361,21 +331,6 @@ Foam::scalar Foam::multiComponentSolidMixture<ThermoType>::kappa
 }
 
 
-template<class ThermoType>
-Foam::scalar Foam::multiComponentSolidMixture<ThermoType>::emissivity
-(
-    scalar p, scalar T, label celli
-) const
-{
-    scalar tmp = 0.0;
-    forAll(solidData_, i)
-    {
-        tmp += solidData_[i].emissivity(T)*X(i, celli, p, T);
-    }
-    return tmp;
-}
-
-
 template<class ThermoType>
 Foam::scalar Foam::multiComponentSolidMixture<ThermoType>::Cp
 (

src/thermophysicalModels/solidThermo/mixtures/multiComponentSolidMixture/multiComponentSolidMixture.H

@@ -208,21 +208,11 @@ public:
             //- Density
             virtual scalar rho(scalar p, scalar T, label celli) const;
 
-            //- Absorption coefficient
-            virtual scalar kappaRad(scalar p, scalar T, label celli) const;
-
-            //- Scatter coefficient
-            virtual scalar sigmaS(scalar p, scalar T, label celli) const;
-
             //- Thermal conductivity
             virtual scalar kappa(scalar p, scalar T, label celli) const;
 
-            //- Emissivity coefficient
-            virtual scalar emissivity(scalar p, scalar T, label celli) const;
-
             //- Specific heat capacity
             virtual scalar Cp(scalar p, scalar T, label celli) const;
-
 };
 
 

src/thermophysicalModels/solidThermo/solidReactionThermo/solidReactionThermo.C

@@ -31,7 +31,7 @@ License
 namespace Foam
 {
     defineTypeNameAndDebug(solidReactionThermo, 0);
-    defineRunTimeSelectionTable(solidReactionThermo, mesh);
+    defineRunTimeSelectionTable(solidReactionThermo, fvMesh);
     defineRunTimeSelectionTable(solidReactionThermo, dictionary);
 }
 

src/thermophysicalModels/solidThermo/solidReactionThermo/solidReactionThermo.H

@@ -66,7 +66,7 @@ public:
     (
         autoPtr,
         solidReactionThermo,
-        mesh,
+        fvMesh,
         (const fvMesh& mesh),
         (mesh)
     );

src/thermophysicalModels/solidThermo/solidReactionThermo/solidReactionThermoNew.C

@@ -53,16 +53,16 @@ Foam::autoPtr<Foam::solidReactionThermo> Foam::solidReactionThermo::New
 
     Info<< "Selecting thermodynamics package " << modelType << endl;
 
-    meshConstructorTable::iterator cstrIter =
+    fvMeshConstructorTable::iterator cstrIter =
-        meshConstructorTablePtr_->find(modelType);
+        fvMeshConstructorTablePtr_->find(modelType);
 
-    if (cstrIter == meshConstructorTablePtr_->end())
+    if (cstrIter == fvMeshConstructorTablePtr_->end())
     {
         FatalErrorIn("solidReactionThermo::New(const fvMesh&)")
             << "Unknown solidReactionThermo type "
             << modelType << nl << nl
             << "Valid solidReactionThermo types:" << nl
-            << meshConstructorTablePtr_->sortedToc() << nl
+            << fvMeshConstructorTablePtr_->sortedToc() << nl
             << exit(FatalError);
     }
 

src/thermophysicalModels/solidThermo/solidReactionThermo/solidReactionThermos.C

@@ -21,80 +21,46 @@ License
     You should have received a copy of the GNU General Public License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
-Class
-    Foam::hePsiReactionThermo
-
-Description
-    Foam::hePsiReactionThermo
-
-SourceFiles
-    hePsiReactionThermo.C
-
 \*---------------------------------------------------------------------------*/
 
-#ifndef hePsiReactionThermo_H
+#include "makeSolidThermo.H"
-#define hePsiReactionThermo_H
+#include "solidReactionThermo.H"
+#include "heSolidThermo.H"
+
+#include "specie.H"
+#include "rhoConst.H"
+#include "hConstThermo.H"
+#include "hExponentialThermo.H"
+#include "constIsoSolidTransport.H"
+#include "constAnIsoSolidTransport.H"
+#include "exponentialSolidTransport.H"
+#include "reactingSolidMixture.H"
+#include "sensibleEnthalpy.H"
+#include "thermo.H"
 
-#include "heThermo.H"
-#include "psiReactionThermo.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
 
-/*---------------------------------------------------------------------------*\
+/* * * * * * * * * * * * * * * * * Enthalpy-based * * * * * * * * * * * * * */
-                      Class hePsiReactionThermo Declaration
-\*---------------------------------------------------------------------------*/
 
-template<class MixtureType>
+makeSolidThermo
-class hePsiReactionThermo
+(
-:
+    solidReactionThermo,
-    public heThermo<psiReactionThermo, MixtureType>
+    heSolidThermo,
-{
+    reactingSolidMixture,
-    // Private Member Functions
+    constIsoSolidTransport,
-
+    sensibleEnthalpy,
-        void calculate();
+    hConstThermo,
-
+    rhoConst,
-        //- Construct as copy (not implemented)
+    specie
-        hePsiReactionThermo(const hePsiReactionThermo<MixtureType>&);
+);
-
-
-public:
-
-    //- Runtime type information
-    TypeName("hePsiReactionThermo");
-
-
-    // Constructors
-
-        //- Construct from mesh
-        hePsiReactionThermo(const fvMesh&);
-
-
-    //- Destructor
-    virtual ~hePsiReactionThermo();
-
-
-    // Member functions
-
-        //- Update properties
-        virtual void correct();
-};
 
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 } // End namespace Foam
 
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-#ifdef NoRepository
-#   include "hePsiReactionThermo.C"
-#endif
-
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-#endif
-
 // ************************************************************************* //

src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C

@@ -230,192 +230,6 @@ Foam::heSolidThermo<BasicSolidThermo, MixtureType>::Kappa() const
 }
 
 
-template<class BasicSolidThermo, class MixtureType>
-Foam::tmp<Foam::volScalarField>
-Foam::heSolidThermo<BasicSolidThermo, MixtureType>::kappaRad() const
-{
-    const fvMesh& mesh = this->T_.mesh();
-
-    tmp<volScalarField> tkappaRad
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "kappaRad",
-                mesh.time().timeName(),
-                mesh,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE
-            ),
-            mesh,
-            inv(dimLength)
-        )
-    );
-
-    volScalarField& kappaRad = tkappaRad();
-    scalarField& kappaRadCells = kappaRad.internalField();
-    const scalarField& TCells = this->T_.internalField();
-    const scalarField& pCells = this->p_.internalField();
-
-    forAll(kappaRadCells, celli)
-    {
-        kappaRadCells[celli] =
-            this->cellVolMixture
-            (
-                pCells[celli],
-                TCells[celli],
-                celli
-            ).kappaRad(TCells[celli]);
-    }
-
-    forAll(kappaRad.boundaryField(), patchi)
-    {
-        scalarField& kappaRadp = kappaRad.boundaryField()[patchi];
-        const scalarField& pT = this->T_.boundaryField()[patchi];
-        const scalarField& pp = this->p_.boundaryField()[patchi];
-
-        forAll(kappaRadp, facei)
-        {
-            kappaRadp[facei] =
-                this->patchFaceVolMixture
-                (
-                    pp[facei],
-                    pT[facei],
-                    patchi,
-                    facei
-                ).kappaRad(pT[facei]);
-        }
-    }
-
-    return tkappaRad;
-}
-
-
-template<class BasicSolidThermo, class MixtureType>
-Foam::tmp<Foam::volScalarField>
-Foam::heSolidThermo<BasicSolidThermo, MixtureType>::sigmaS() const
-{
-    const fvMesh& mesh = this->T_.mesh();
-
-    tmp<volScalarField> tsigmaS
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "sigmaS",
-                mesh.time().timeName(),
-                mesh,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE
-            ),
-            mesh,
-            inv(dimLength)
-        )
-    );
-
-    volScalarField& sigmaS = tsigmaS();
-    scalarField& sigmaSCells = sigmaS.internalField();
-    const scalarField& TCells = this->T_.internalField();
-    const scalarField& pCells = this->p_.internalField();
-
-    forAll(sigmaSCells, celli)
-    {
-        sigmaSCells[celli] =
-            this->cellVolMixture
-            (
-                pCells[celli],
-                TCells[celli],
-                celli
-            ).sigmaS(TCells[celli]);
-    }
-
-    forAll(sigmaS.boundaryField(), patchi)
-    {
-        scalarField& sigmaSp = sigmaS.boundaryField()[patchi];
-        const scalarField& pT = this->T_.boundaryField()[patchi];
-        const scalarField& pp = this->p_.boundaryField()[patchi];
-
-        forAll(sigmaSp, facei)
-        {
-            sigmaSp[facei] =
-                this->patchFaceVolMixture
-                (
-                    pp[facei],
-                    pT[facei],
-                    patchi,
-                    facei
-                ).sigmaS(pT[facei]);
-        }
-    }
-
-    return tsigmaS;
-}
-
-
-template<class BasicSolidThermo, class MixtureType>
-Foam::tmp<Foam::volScalarField>
-Foam::heSolidThermo<BasicSolidThermo, MixtureType>::emissivity() const
-{
-    const fvMesh& mesh = this->T_.mesh();
-
-    tmp<volScalarField> temissivity
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "emissivity",
-                mesh.time().timeName(),
-                mesh,
-                IOobject::NO_READ,
-                IOobject::NO_WRITE
-            ),
-            mesh,
-            inv(dimLength)
-        )
-    );
-
-    volScalarField& emissivity = temissivity();
-    scalarField& emissivityCells = emissivity.internalField();
-    const scalarField& TCells = this->T_.internalField();
-    const scalarField& pCells = this->p_.internalField();
-
-    forAll(emissivityCells, celli)
-    {
-        emissivityCells[celli] =
-            this->cellVolMixture
-            (
-                pCells[celli],
-                TCells[celli],
-                celli
-            ).emissivity(TCells[celli]);
-    }
-
-    forAll(emissivity.boundaryField(), patchi)
-    {
-        scalarField& emissivityp = emissivity.boundaryField()[patchi];
-        const scalarField& pT = this->T_.boundaryField()[patchi];
-        const scalarField& pp = this->p_.boundaryField()[patchi];
-
-        forAll(emissivityp, facei)
-        {
-            emissivityp[facei] =
-                this->patchFaceVolMixture
-                (
-                    pp[facei],
-                    pT[facei],
-                    patchi,
-                    facei
-                ).emissivity(pT[facei]);
-        }
-    }
-
-    return temissivity;
-}
-
-
 template<class BasicSolidThermo, class MixtureType>
 Foam::tmp<Foam::vectorField>
 Foam::heSolidThermo<BasicSolidThermo, MixtureType>::Kappa
@@ -445,89 +259,4 @@ Foam::heSolidThermo<BasicSolidThermo, MixtureType>::Kappa
 }
 
 
-template<class BasicSolidThermo, class MixtureType>
-Foam::tmp<Foam::scalarField>
-Foam::heSolidThermo<BasicSolidThermo, MixtureType>::kappaRad
-(
-    const label patchi
-) const
-{
-    const scalarField& Tp = this->T_.boundaryField()[patchi];
-    tmp<scalarField> tKappaRad(new scalarField(Tp.size()));
-    scalarField& KappapRadp = tKappaRad();
-    const scalarField& pp = this->p_.boundaryField()[patchi];
-
-    forAll(Tp, facei)
-    {
-        KappapRadp[facei] =
-            this->patchFaceVolMixture
-            (
-                pp[facei],
-                Tp[facei],
-                patchi,
-                facei
-            ).kappaRad(Tp[facei]);
-    }
-
-    return tKappaRad;
-}
-
-
-template<class BasicSolidThermo, class MixtureType>
-Foam::tmp<Foam::scalarField>
-Foam::heSolidThermo<BasicSolidThermo, MixtureType>::sigmaS
-(
-    const label patchi
-) const
-{
-    const scalarField& Tp = this->T_.boundaryField()[patchi];
-    tmp<scalarField> tsigmaS(new scalarField(Tp.size()));
-    scalarField& sigmaSp = tsigmaS();
-    const scalarField& pp = this->p_.boundaryField()[patchi];
-
-
-    forAll(Tp, facei)
-    {
-        sigmaSp[facei] =
-            this->patchFaceVolMixture
-            (
-                pp[facei],
-                Tp[facei],
-                patchi,
-                facei
-            ).sigmaS(Tp[facei]);
-    }
-
-    return tsigmaS;
-}
-
-
-template<class BasicSolidThermo, class MixtureType>
-Foam::tmp<Foam::scalarField>
-Foam::heSolidThermo<BasicSolidThermo, MixtureType>::emissivity
-(
-    const label patchi
-) const
-{
-    const scalarField& Tp = this->T_.boundaryField()[patchi];
-    tmp<scalarField> temissivity(new scalarField(Tp.size()));
-    scalarField& emissivity = temissivity();
-    const scalarField& pp = this->p_.boundaryField()[patchi];
-
-    forAll(Tp, facei)
-    {
-        emissivity[facei] =
-            this->patchFaceVolMixture
-            (
-                pp[facei],
-                Tp[facei],
-                patchi,
-                facei
-            ).emissivity(Tp[facei]);
-    }
-
-    return temissivity;
-}
-
-
 // ************************************************************************* //

src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.H

@@ -25,7 +25,7 @@ Class
     Foam::heSolidThermo
 
 Description
-    Enthalpy, internal energy based for a solid mixture.
+    Energy for a solid mixture
 
 SourceFiles
     heSolidThermo.C
@@ -51,7 +51,6 @@ class heSolidThermo
 :
     public heThermo<BasicSolidThermo, MixtureType>
 {
-
     // Private Member Functions
 
         //- Calculate the thermo variables
@@ -91,29 +90,11 @@ public:
             //- Anisotropic thermal conductivity [W/m/K]
             virtual tmp<volVectorField> Kappa() const;
 
-            //- Absorption coefficient [1/m]
-            virtual tmp<volScalarField> kappaRad() const;
-
-            //- Scatter coefficient
-            virtual tmp<volScalarField> sigmaS() const;
-
-            //- Emissivity coefficient [1/m]
-            virtual tmp<volScalarField> emissivity() const;
-
 
         // Per patch calculation
 
             //- Anisotropic thermal conductivity [W/m/K]
             virtual tmp<vectorField> Kappa(const label patchI) const;
-
-             //- Absorption coefficient [1/m]
-            virtual tmp<scalarField> kappaRad(const label patchI) const;
-
-            //- Scatter coefficient
-            virtual tmp<scalarField> sigmaS(const label patchI) const;
-
-            //- Emissivity coefficient [1/m]
-            virtual tmp<scalarField> emissivity(const label patchI) const;
 };
 
 

src/thermophysicalModels/solidThermo/solidThermo/makeSolidThermo.H

@@ -32,79 +32,42 @@ Description
 #ifndef makeSolidThermo_H
 #define makesolidThermo_H
 
-#include "addToRunTimeSelectionTable.H"
+#include "makeThermo.H"
-#include "basicThermo.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-#define makeSolidThermo(BaseThermo,Cthermo,Mixture,Transport,Radiation,Type,Thermo,EqnOfState,Specie)\
+#define makeSolidThermo(BaseThermo,Cthermo,Mixture,Transport,Type,Thermo,EqnOfState,Specie)\
                                                                               \
-typedef                                                                       \
+makeThermoTypedefs                                                            \
-    Transport                                                                 \
-    <                                                                         \
-        Radiation                                                             \
-        <                                                                     \
-            species::thermo                                                   \
-            <                                                                 \
-                Thermo                                                        \
-                <                                                             \
-                    EqnOfState                                                \
-                    <                                                         \
-                        Specie                                                \
-                    >                                                         \
-                >,                                                            \
-                Type                                                          \
-            >                                                                 \
-        >                                                                     \
-    > Transport##Radiation##Type##Thermo##EqnOfState##Specie;                 \
-                                                                              \
-typedef                                                                       \
-    Mixture                                                                   \
-    <                                                                         \
-        Transport##Radiation##Type##Thermo##EqnOfState##Specie                \
-    > Mixture##Transport##Radiation##Type##Thermo##EqnOfState##Specie;        \
-                                                                              \
-typedef                                                                       \
-    Cthermo                                                                   \
-    <                                                                         \
-        BaseThermo,                                                           \
-        Mixture##Transport##Radiation##Type##Thermo##EqnOfState##Specie       \
-    > Cthermo##Mixture##Transport##Radiation##Type##Thermo                    \
-          ##EqnOfState##Specie;                                               \
-                                                                              \
-defineTemplateTypeNameAndDebugWithName                                        \
 (                                                                             \
-    Cthermo##Mixture##Transport##Radiation##Type##Thermo                      \
+    BaseThermo,                                                               \
-         ##EqnOfState##Specie,                                                \
+    Cthermo,                                                                  \
-    (                                                                         \
+    Mixture,                                                                  \
-        #Cthermo"<"#Mixture"<"                                                \
+    Transport,                                                                \
-      + Transport##Radiation##Type##Thermo##EqnOfState##Specie::typeName()    \
+    Type,                                                                     \
-      + ">>"                                                                  \
+    Thermo,                                                                   \
-    ).c_str(),                                                                \
+    EqnOfState,                                                               \
-    0                                                                         \
+    Specie                                                                    \
-);                                                                            \
+)                                                                             \
                                                                               \
 addToRunTimeSelectionTable                                                    \
 (                                                                             \
     basicThermo,                                                              \
-    Cthermo##Mixture##Transport##Radiation##Type##Thermo                      \
+    Cthermo##Mixture##Transport##Type##Thermo##EqnOfState##Specie,            \
-        ##EqnOfState##Specie,                                                 \
     fvMesh                                                                    \
 );                                                                            \
                                                                               \
 addToRunTimeSelectionTable                                                    \
 (                                                                             \
     BaseThermo,                                                               \
-    Cthermo##Mixture##Transport##Radiation##Type##Thermo                      \
+    Cthermo##Mixture##Transport##Type##Thermo##EqnOfState##Specie,            \
-        ##EqnOfState##Specie,                                                 \
     fvMesh                                                                    \
 );                                                                            \
                                                                               \
 addToRunTimeSelectionTable                                                    \
 (                                                                             \
     BaseThermo,                                                               \
-    Cthermo##Mixture##Transport##Radiation##Type##Thermo                      \
+    Cthermo##Mixture##Transport##Type##Thermo##EqnOfState##Specie,            \
-        ##EqnOfState##Specie,                                                 \
     dictionary                                                                \
 );
 

src/thermophysicalModels/solidThermo/solidThermo/solidThermo.H

@@ -127,49 +127,20 @@ public:
             //- Thermal conductivity [W/m/K]
             virtual tmp<volVectorField> Kappa() const = 0;
 
-            //- Absorption coefficient [1/m]
-            virtual tmp<volScalarField> kappaRad() const = 0;
-
-            //- Scatter coefficient
-            virtual tmp<volScalarField> sigmaS() const = 0;
-
-            //- Emissivity coefficient [1/m]
-            virtual tmp<volScalarField> emissivity() const = 0;
-
 
         // Per patch calculation
 
-
             //- Anisotropic thermal conductivity [W/m/K]
             virtual tmp<vectorField> Kappa
             (
                 const label patchI
             ) const = 0;
 
-             //- Absorption coefficient [1/m]
-            virtual tmp<scalarField> kappaRad
-            (
-                const label patchI
-            ) const = 0;
-
-            //- Scatter coefficient
-            virtual tmp<scalarField> sigmaS
-            (
-                const label patchI
-            ) const = 0;
-
-            //- Emissivity coefficient [1/m]
-            virtual tmp<scalarField> emissivity
-            (
-                const label patchI
-            ) const = 0;
-
 
     // I-O
 
         //- Read thermophysicalProperties dictionary
         virtual bool read();
-
 };
 
 

src/thermophysicalModels/solidThermo/solidThermo/solidThermoNew.C

@@ -32,48 +32,20 @@ Foam::autoPtr<Foam::solidThermo> Foam::solidThermo::New
     const fvMesh& mesh
 )
 {
-    if (debug)
+    return basicThermo::New<solidThermo>(mesh);
-    {
-        Info<< "solidThermo::New(const fvMesh&): "
-            << "constructing solidThermo"
-            << endl;
-    }
-
-    const word thermoType
-    (
-        IOdictionary
-        (
-            IOobject
-            (
-                "thermophysicalProperties",
-                mesh.time().constant(),
-                mesh,
-                IOobject::MUST_READ_IF_MODIFIED,
-                IOobject::NO_WRITE,
-                false
-            )
-        ).lookup("thermoType")
-    );
-
-    fvMeshConstructorTable::iterator cstrIter =
-        fvMeshConstructorTablePtr_->find(thermoType);
-
-    if (cstrIter == fvMeshConstructorTablePtr_->end())
-    {
-        FatalErrorIn
-        (
-            "solidThermo::New(const fvMesh&)"
-        )   << "Unknown solidThermo type " << thermoType
-            << endl << endl
-            << "Valid solidThermo types are :" << endl
-            << fvMeshConstructorTablePtr_->toc()
-            << exit(FatalError);
-    }
-
-    return autoPtr<solidThermo>(cstrIter()(mesh));
 }
 
 
+Foam::autoPtr<Foam::solidThermo> Foam::solidThermo::New
+(
+    const fvMesh& mesh,
+    const dictionary& dict
+)
+{
+    return basicThermo::New<solidThermo>(mesh, dict);
+}
+
+/*
 Foam::autoPtr<Foam::solidThermo> Foam::solidThermo::New
 (
     const fvMesh& mesh,
@@ -106,6 +78,7 @@ Foam::autoPtr<Foam::solidThermo> Foam::solidThermo::New
 
     return autoPtr<solidThermo>(cstrIter()(mesh, dict));
 }
+*/
 
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

src/thermophysicalModels/solidThermo/solidThermo/solidThermos.C

@@ -23,8 +23,9 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-
 #include "makeSolidThermo.H"
+#include "solidThermo.H"
+#include "heSolidThermo.H"
 
 #include "specie.H"
 #include "rhoConst.H"
@@ -33,7 +34,6 @@ License
 #include "constIsoSolidTransport.H"
 #include "constAnIsoSolidTransport.H"
 #include "exponentialSolidTransport.H"
-#include "constSolidRad.H"
 #include "pureSolidMixture.H"
 #include "multiComponentSolidMixture.H"
 #include "reactingSolidMixture.H"
@@ -41,11 +41,6 @@ License
 #include "sensibleInternalEnergy.H"
 #include "thermo.H"
 
-#include "heThermo.H"
-
-#include "solidThermo.H"
-#include "solidReactionThermo.H"
-#include "heSolidThermo.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -60,7 +55,6 @@ makeSolidThermo
     heSolidThermo,
     pureSolidMixture,
     constIsoSolidTransport,
-    constSolidRad,
     sensibleEnthalpy,
     hConstThermo,
     rhoConst,
@@ -73,7 +67,6 @@ makeSolidThermo
     heSolidThermo,
     pureSolidMixture,
     constAnIsoSolidTransport,
-    constSolidRad,
     sensibleEnthalpy,
     hConstThermo,
     rhoConst,
@@ -86,7 +79,6 @@ makeSolidThermo
     heSolidThermo,
     pureSolidMixture,
     exponentialSolidTransport,
-    constSolidRad,
     sensibleEnthalpy,
     hExponentialThermo,
     rhoConst,
@@ -99,20 +91,6 @@ makeSolidThermo
     heSolidThermo,
     multiComponentSolidMixture,
     constIsoSolidTransport,
-    constSolidRad,
-    sensibleEnthalpy,
-    hConstThermo,
-    rhoConst,
-    specie
-);
-
-makeSolidThermo
-(
-    solidReactionThermo,
-    heSolidThermo,
-    reactingSolidMixture,
-    constIsoSolidTransport,
-    constSolidRad,
     sensibleEnthalpy,
     hConstThermo,
     rhoConst,

src/thermophysicalModels/specie/Make/files

@@ -1,11 +1,6 @@
-atomicWeights = atomicWeights
+atomicWeights/atomicWeights.C
-specie = specie
+specie/specie.C
-equationOfState = equationOfState
+reaction/reactions/makeReactions.C
-reactions = reaction/reactions
+reaction/reactions/makeLangmuirHinshelwoodReactions.C
-
-$(atomicWeights)/atomicWeights.C
-$(specie)/specie.C
-$(reactions)/makeReactionThermoReactions.C
-$(reactions)/makeLangmuirHinshelwoodReactions.C
 
 LIB = $(FOAM_LIBBIN)/libspecie

tutorials/combustion/PDRFoam/flamePropagationWithObstacles/constant/thermophysicalProperties

@@ -17,7 +17,7 @@ FoamFile
 
 thermoType
 {
-    type            heheuReactionThermo;
+    type            heheuPsiThermo;
     mixture         inhomogeneousMixture;
     transport       sutherland;
     thermo          janaf;

tutorials/combustion/XiFoam/ras/moriyoshiHomogeneous/constant/thermophysicalProperties

@@ -17,7 +17,7 @@ FoamFile
 
 thermoType
 {
-    type            heheuReactionThermo;
+    type            heheuPsiThermo;
     mixture         homogeneousMixture;
     transport       sutherland;
     thermo          janaf;

tutorials/combustion/chemFoam/gri/constant/chemistryProperties

@@ -17,8 +17,8 @@ FoamFile
 
 chemistryType
 {
-    chemistryModel    ODEChemistryModel;
     chemistrySolver   ode;
+    chemistryThermo   psi;
 }
 
 chemistry       on;

tutorials/combustion/chemFoam/gri/constant/thermophysicalProperties

@@ -17,7 +17,7 @@ FoamFile
 
 thermoType
 {
-    type            hePsiReactionThermo;
+    type            hePsiThermo;
     mixture         reactingMixture;
     transport       sutherland;
     thermo          janaf;

tutorials/combustion/chemFoam/h2/constant/chemistryProperties

@@ -17,8 +17,8 @@ FoamFile
 
 chemistryType
 {
-    chemistryModel    ODEChemistryModel;
     chemistrySolver   ode;
+    chemistryThermo   psi;
 }
 
 chemistry       on;

tutorials/combustion/chemFoam/h2/constant/thermophysicalProperties

@@ -17,7 +17,7 @@ FoamFile
 
 thermoType
 {
-    type            hePsiReactionThermo;
+    type            hePsiThermo;
     mixture         reactingMixture;
     transport       sutherland;
     thermo          janaf;

tutorials/combustion/chemFoam/ic8h18/constant/chemistryProperties

@@ -17,8 +17,8 @@ FoamFile
 
 chemistryType
 {
-    chemistryModel    ODEChemistryModel;
     chemistrySolver   ode;
+    chemistryThermo   psi;
 }
 
 chemistry       on;