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("U")); + + volScalarField::New + ( + "strainRate", + sqrt(2.0)*mag(symm(fvc::grad(U))) + )->write(); +#}; + +// ************************************************************************* //