From 58cb08d8cd9f01bc40ff3282aaf4e17c0361165d Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Mon, 18 Sep 2023 11:38:54 +0100
Subject: [PATCH] tutorials/incompressibleFluid/pitzDailySteady: Added
 strainRate coded functionObject

to demonstrate how easy it is to add field post-processing to cases.
---
 .../pitzDailySteady/system/controlDict        |  2 ++
 .../pitzDailySteady/system/strainRate         | 32 +++++++++++++++++++
 2 files changed, 34 insertions(+)
 create mode 100644 tutorials/incompressibleFluid/pitzDailySteady/system/strainRate

diff --git a/tutorials/incompressibleFluid/pitzDailySteady/system/controlDict b/tutorials/incompressibleFluid/pitzDailySteady/system/controlDict
index ca48a991eb..d51067c3bb 100644
--- a/tutorials/incompressibleFluid/pitzDailySteady/system/controlDict
+++ b/tutorials/incompressibleFluid/pitzDailySteady/system/controlDict
@@ -63,6 +63,8 @@ functions
     )
 
     #includeFunc writeObjects(kEpsilon:G)
+
+    #includeFunc strainRate(writeControl=writeTime)
 }
 
 // ************************************************************************* //
diff --git a/tutorials/incompressibleFluid/pitzDailySteady/system/strainRate b/tutorials/incompressibleFluid/pitzDailySteady/system/strainRate
new file mode 100644
index 0000000000..900cbddae1
--- /dev/null
+++ b/tutorials/incompressibleFluid/pitzDailySteady/system/strainRate
@@ -0,0 +1,32 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+
+type            coded;
+
+// Load the library containing the 'coded' functionObject
+libs            ("libutilityFunctionObjects.so");
+
+name            strainRate;
+
+codeInclude
+#{
+    #include "fvcGrad.H"
+#};
+
+codeWrite
+#{
+    const volVectorField& U(mesh().lookupObject<volVectorField>("U"));
+
+    volScalarField::New
+    (
+        "strainRate",
+        sqrt(2.0)*mag(symm(fvc::grad(U)))
+    )->write();
+#};
+
+// ************************************************************************* //