qr
codeWrite #{ const volVectorField& U(mesh().lookupObject<volVectorField>("U")); volScalarField::New ( "strainRate", sqrt(2.0)*mag(symm(fvc::grad(U))) )->write(); #};