applications/utilities/postProcessing/viscoelastic/PSD/PSD.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 Application
  26     PSD
  27
  28 Description
  29     Calculates and writes the Principal Stress Difference (PSD) for each time.
  30
  31 Author
  32     Jovani L. Favero.  All rights reserved
  33
  34 \*---------------------------------------------------------------------------*/
  35
  36 #include "fvCFD.H"
  37
  38 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  39
  40 int main(int argc, char *argv[])
  41 {
  42
  43 #   include "addTimeOptions.H"
  44 #   include "setRootCase.H"
  45
  46 #   include "createTime.H"
  47
  48     // Get times list
  49     instantList Times = runTime.times();
  50
  51     // set startTime and endTime depending on -time and -latestTime options
  52 #   include "checkTimeOptions.H"
  53
  54     runTime.setTime(Times[startTime], startTime);
  55
  56 #   include "createMesh.H"
  57
  58     for (label i=startTime; i<endTime; i++)
  59     {
  60         runTime.setTime(Times[i], i);
  61
  62         Info<< "Time = " << runTime.timeName() << endl;
  63
  64         mesh.readUpdate();
  65
  66         IOobject tauHeader
  67         (
  68             "tau",
  69             runTime.timeName(),
  70             mesh,
  71             IOobject::MUST_READ
  72         );
  73
  74         // Check tau exists
  75         if (tauHeader.headerOk())
  76         {
  77             mesh.readUpdate();
  78
  79             Info<< "    Reading tau" << endl;
  80             volSymmTensorField tau(tauHeader, mesh);
  81
  82             Info<< "    Calculating PSD"<< endl;
  83             volScalarField PSD
  84             (
  85                 IOobject
  86                 (
  87                     "PSD",
  88                     runTime.timeName(),
  89                     mesh,
  90                     IOobject::NO_READ
  91                 ),
  92                 Foam::sqrt
  93                 (
  94                     Foam::sqr(tau.component(symmTensor::YY)
  95                   - tau.component(symmTensor::XX))
  96                   + 4*Foam::sqr(tau.component(symmTensor::XY))
  97                 )
  98             );
  99             PSD.write();
 100         }
 101         else
 102         {
 103             Info<< "    No tau" << endl;
 104         }
 105
 106         Info<< endl;
 107
 108     } //for time
 109
 110     Info<< endl;
 111     Info<< "    End"<< nl << endl;
 112
 113     return(0);
 114 }
 115
 116
 117 // ************************************************************************* //