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

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