diff --git a/applications/utilities/mesh/generation/extrudeMesh/extrudeMesh.C b/applications/utilities/mesh/generation/extrudeMesh/extrudeMesh.C
index 38bce28409..b351f01e37 100644
--- a/applications/utilities/mesh/generation/extrudeMesh/extrudeMesh.C
+++ b/applications/utilities/mesh/generation/extrudeMesh/extrudeMesh.C
@@ -25,10 +25,7 @@ Application
     extrudeMesh
 
 Description
-    Extrude mesh from existing patch (by default outwards facing normals;
-    optional flips faces) or from patch read from file.
-
-    Note: Merges close points so be careful.
+    Extrude mesh from existing patch or from patch read from file.
 
     Type of extrusion prescribed by run-time selectable model.
 
@@ -96,25 +93,70 @@ label findIndex(const polyBoundaryMesh& patches, const word& name)
 }
 
 
-labelList patchFaces(const polyBoundaryMesh& patches, const wordList& names)
+labelList findPatchZones(polyMesh& mesh, const wordList& zoneNames)
 {
+    labelList patchZones(zoneNames.size());
+
+    // Add optional cellZone generation for the extrusion layer
+    forAll(zoneNames, i)
+    {
+        const label cellZonei = mesh.cellZones().findIndex(zoneNames[i]);
+        if (cellZonei >= 0)
+        {
+            patchZones[i] = cellZonei;
+        }
+        else
+        {
+            label nZones = mesh.cellZones().size();
+            mesh.cellZones().setSize(nZones + 1);
+            mesh.cellZones().set
+            (
+                nZones,
+                new cellZone
+                (
+                    zoneNames[i],
+                    labelList(),
+                    mesh.cellZones()
+                )
+            );
+            patchZones[i] = nZones;
+        }
+    }
+
+    return patchZones;
+}
+
+
+labelList findPatchFaces
+(
+    polyMesh& mesh,
+    const polyBoundaryMesh& patches,
+    const wordList& patchNames,
+    const labelList& patchZones,
+    labelList& faceCellZones
+)
+{
+    // Count of the number of faces in the named patches
     label n = 0;
 
-    forAll(names, i)
+    forAll(patchNames, i)
     {
-        const polyPatch& pp = patches[findIndex(patches, names[i])];
-
-        n += pp.size();
+        n += patches[findIndex(patches, patchNames[i])].size();
     }
+
     labelList faceLabels(n);
+    faceCellZones.setSize(n);
     n = 0;
-    forAll(names, i)
+
+    forAll(patchNames, i)
     {
-        const polyPatch& pp = patches[findIndex(patches, names[i])];
+        const polyPatch& pp = patches[findIndex(patches, patchNames[i])];
 
         forAll(pp, j)
         {
-            faceLabels[n++] = pp.start()+j;
+            faceLabels[n] = pp.start()+j;
+            faceCellZones[n] = patchZones[i];
+            n++;
         }
     }
 
@@ -122,6 +164,24 @@ labelList patchFaces(const polyBoundaryMesh& patches, const wordList& names)
 }
 
 
+labelList findPatchFaces
+(
+    const polyBoundaryMesh& patches,
+    const word& name
+)
+{
+    const polyPatch& pp = patches[findIndex(patches, name)];
+    labelList faceLabels(pp.size());
+
+    forAll(pp, i)
+    {
+        faceLabels[i] = pp.start() + i;
+    }
+
+    return faceLabels;
+}
+
+
 void updateFaceLabels(const polyTopoChangeMap& map, labelList& faceLabels)
 {
     const labelList& reverseMap = map.reverseFaceMap();
@@ -282,14 +342,33 @@ int main(int argc, char *argv[])
                 sourceCaseDir
                /"processor" + Foam::name(Pstream::myProcNo());
         }
-        wordList sourcePatches;
-        dict.lookup("sourcePatches") >> sourcePatches;
+        wordList sourcePatches(dict.lookup("sourcePatches"));
 
         if (sourcePatches.size() == 1)
         {
             frontPatchName = sourcePatches[0];
         }
 
+        wordList zoneNames(dict.lookupOrDefault("zoneNames", wordList::null()));
+        if (zoneNames.size() == 1)
+        {
+            if (zoneNames[0] == "patchNames")
+            {
+                zoneNames = sourcePatches;
+            }
+            else
+            {
+                zoneNames = wordList(sourcePatches.size(), zoneNames[0]);
+            }
+        }
+        else if (zoneNames.size() && zoneNames.size() != sourcePatches.size())
+        {
+            FatalErrorInFunction
+                << "Number of zoneNames "
+                   "does not equal the number of sourcePatches"
+                << exit(FatalError);
+        }
+
         Info<< "Extruding patches " << sourcePatches
             << " on mesh " << sourceCasePath << nl
             << endl;
@@ -305,42 +384,59 @@ int main(int argc, char *argv[])
 
         const polyBoundaryMesh& patches = mesh.boundaryMesh();
 
+        const labelList patchZones
+        (
+            zoneNames.size()
+          ? findPatchZones(mesh, zoneNames)
+          : labelList(sourcePatches.size(), -1)
+        );
+
+        labelList faceCellZones;
+        const labelList patchFaces
+        (
+            findPatchFaces
+            (
+                mesh,
+                patches,
+                sourcePatches,
+                patchZones,
+                faceCellZones
+            )
+        );
 
         // Extrusion engine. Either adding to existing mesh or
         // creating separate mesh.
         addPatchCellLayer layerExtrude(mesh, (mode == MESH));
 
-        const labelList meshFaces(patchFaces(patches, sourcePatches));
-
         if (mode == PATCH && flipNormals)
         {
             // Cheat. Flip patch faces in mesh. This invalidates the
             // mesh (open cells) but does produce the correct extrusion.
             polyTopoChange meshMod(mesh);
-            forAll(meshFaces, i)
+            forAll(patchFaces, i)
             {
-                label meshFacei = meshFaces[i];
+                label patchFacei = patchFaces[i];
 
-                label patchi = patches.whichPatch(meshFacei);
-                label own = mesh.faceOwner()[meshFacei];
+                label patchi = patches.whichPatch(patchFacei);
+                label own = mesh.faceOwner()[patchFacei];
                 label nei = -1;
                 if (patchi == -1)
                 {
-                    nei = mesh.faceNeighbour()[meshFacei];
+                    nei = mesh.faceNeighbour()[patchFacei];
                 }
 
-                label zoneI = mesh.faceZones().whichZone(meshFacei);
+                label zoneI = mesh.faceZones().whichZone(patchFacei);
                 bool zoneFlip = false;
                 if (zoneI != -1)
                 {
-                    label index = mesh.faceZones()[zoneI].whichFace(meshFacei);
+                    label index = mesh.faceZones()[zoneI].whichFace(patchFacei);
                     zoneFlip = mesh.faceZones()[zoneI].flipMap()[index];
                 }
 
                 meshMod.modifyFace
                 (
-                    mesh.faces()[meshFacei].reverseFace(),  // modified face
-                    meshFacei,                      // label of face
+                    mesh.faces()[patchFacei].reverseFace(), // modified face
+                    patchFacei,                     // label of face
                     own,                            // owner
                     nei,                            // neighbour
                     true,                           // face flip
@@ -363,7 +459,7 @@ int main(int argc, char *argv[])
             IndirectList<face>
             (
                 mesh.faces(),
-                meshFaces
+                patchFaces
             ),
             mesh.points()
         );
@@ -428,7 +524,7 @@ int main(int argc, char *argv[])
         );
 
 
-        // Add any patches.
+        // Add any patches
 
         label nAdded = nPatches - mesh.boundaryMesh().size();
         reduce(nAdded, sumOp<label>());
@@ -572,6 +668,7 @@ int main(int argc, char *argv[])
             nFaceLayers,
             nPointLayers,
             firstLayerDisp,
+            faceCellZones,
             meshMod()
         );
 
@@ -718,16 +815,16 @@ int main(int argc, char *argv[])
 
         // Get the faces on front and back
         frontPatchName = "originalPatch";
-        frontPatchFaces = patchFaces
+        frontPatchFaces = findPatchFaces
         (
             meshFromSurface().boundaryMesh(),
-            wordList(1, frontPatchName)
+            frontPatchName
         );
         backPatchName = "otherSide";
-        backPatchFaces = patchFaces
+        backPatchFaces = findPatchFaces
         (
             meshFromSurface().boundaryMesh(),
-            wordList(1, backPatchName)
+            backPatchName
         );
     }
 
diff --git a/etc/caseDicts/annotated/extrudeMeshDict b/etc/caseDicts/annotated/extrudeMeshDict
index 50bd77c28f..0c146bd943 100644
--- a/etc/caseDicts/annotated/extrudeMeshDict
+++ b/etc/caseDicts/annotated/extrudeMeshDict
@@ -36,6 +36,14 @@ surface "movingWall.stl";
 // patch.
 flipNormals false;
 
+// Optional zone name/names for the cells in the extruded regions
+// The cells can be in a single zone or a separate zone for each extruded patch
+//zoneNames (movingWallCells)
+
+// ... or in a separate zone for each extruded patch
+// with the same name as the corresponding patch
+//zoneNames (patchNames)
+
 //- Linear extrusion in point-normal direction
 //extrudeModel        linearNormal;
 
diff --git a/src/OpenFOAM/meshes/polyMesh/zones/cellZone/cellZone.C b/src/OpenFOAM/meshes/polyMesh/zones/cellZone/cellZone.C
index 21d16fd9b9..54bbf2d23b 100644
--- a/src/OpenFOAM/meshes/polyMesh/zones/cellZone/cellZone.C
+++ b/src/OpenFOAM/meshes/polyMesh/zones/cellZone/cellZone.C
@@ -134,6 +134,7 @@ void Foam::cellZone::topoChange(const polyTopoChangeMap& map)
 
     labelHashSet newIndices;
     const labelList& cellMap = map.cellMap();
+    const labelList& reverseCellMap = map.reverseCellMap();
 
     forAll(cellMap, celli)
     {
@@ -143,6 +144,14 @@ void Foam::cellZone::topoChange(const polyTopoChangeMap& map)
         }
     }
 
+    forAll(reverseCellMap, celli)
+    {
+        if (reverseCellMap[celli] >= 0 && localIndex(celli) != -1)
+        {
+            newIndices.insert(reverseCellMap[celli]);
+        }
+    }
+
     labelList::operator=(newIndices.sortedToc());
 }
 
diff --git a/src/polyTopoChange/polyTopoChange/addPatchCellLayer.C b/src/polyTopoChange/polyTopoChange/addPatchCellLayer.C
index bd2acf7461..d48235e446 100644
--- a/src/polyTopoChange/polyTopoChange/addPatchCellLayer.C
+++ b/src/polyTopoChange/polyTopoChange/addPatchCellLayer.C
@@ -698,6 +698,7 @@ void Foam::addPatchCellLayer::setRefinement
     const labelList& nFaceLayers,
     const labelList& nPointLayers,
     const vectorField& firstLayerDisp,
+    const labelList& faceCellZones,
     polyTopoChange& meshMod
 )
 {
@@ -1077,7 +1078,7 @@ void Foam::addPatchCellLayer::setRefinement
 
             const labelList ownZones
             (
-                mesh_.cellZones().whichZones( mesh_.faceOwner()[meshFacei])
+                mesh_.cellZones().whichZones(mesh_.faceOwner()[meshFacei])
             );
 
             for (label i = 0; i < nFaceLayers[patchFacei]; i++)
@@ -1096,6 +1097,14 @@ void Foam::addPatchCellLayer::setRefinement
                         addedCells[patchFacei][i]
                     );
                 }
+
+                if (faceCellZones[patchFacei] != -1)
+                {
+                    cellZonesAddedCells_[faceCellZones[patchFacei]].insert
+                    (
+                        addedCells[patchFacei][i]
+                    );
+                }
             }
         }
     }
diff --git a/src/polyTopoChange/polyTopoChange/addPatchCellLayer.H b/src/polyTopoChange/polyTopoChange/addPatchCellLayer.H
index e95af7822a..109a5fb558 100644
--- a/src/polyTopoChange/polyTopoChange/addPatchCellLayer.H
+++ b/src/polyTopoChange/polyTopoChange/addPatchCellLayer.H
@@ -350,6 +350,7 @@ public:
                 const labelList& nFaceLayers,
                 const labelList& nPointLayers,
                 const vectorField& firstLayerDisp,
+                const labelList& faceCellZones,
                 polyTopoChange& meshMod
             );
 
@@ -363,6 +364,7 @@ public:
                 const indirectPrimitivePatch& pp,
                 const labelList& sidePatchID,
                 const vectorField& overallDisplacement,
+                const labelList& faceCellZones,
                 polyTopoChange& meshMod
             )
             {
@@ -377,6 +379,7 @@ public:
                     labelList(pp.size(), nLayers),      // nFaceLayers
                     labelList(pp.nPoints(), nLayers),   // nPointLayers
                     overallDisplacement / nLayers,      // firstLayerDisp
+                    faceCellZones,
                     meshMod
                 );
             }