src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "forceCoeffs.H"
  27 #include "dictionary.H"
  28 #include "Time.H"
  29 #include "Pstream.H"
  30
  31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  32
  33 defineTypeNameAndDebug(Foam::forceCoeffs, 0);
  34
  35
  36 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  37
  38 Foam::forceCoeffs::forceCoeffs
  39 (
  40     const word& name,
  41     const objectRegistry& obr,
  42     const dictionary& dict,
  43     const bool loadFromFiles
  44 )
  45 :
  46     forces(name, obr, dict, loadFromFiles),
  47     liftDir_(vector::zero),
  48     dragDir_(vector::zero),
  49     pitchAxis_(vector::zero),
  50     magUInf_(0.0),
  51     lRef_(0.0),
  52     Aref_(0.0)
  53 {
  54     read(dict);
  55 }
  56
  57
  58 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  59
  60 Foam::forceCoeffs::~forceCoeffs()
  61 {}
  62
  63
  64 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  65
  66 void Foam::forceCoeffs::read(const dictionary& dict)
  67 {
  68     if (active_)
  69     {
  70         forces::read(dict);
  71
  72         // Directions for lift and drag forces, and pitch moment
  73         dict.lookup("liftDir") >> liftDir_;
  74         dict.lookup("dragDir") >> dragDir_;
  75         dict.lookup("pitchAxis") >> pitchAxis_;
  76
  77         // Free stream velocity magnitude
  78         dict.lookup("magUInf") >> magUInf_;
  79
  80         // Reference length and area scales
  81         dict.lookup("lRef") >> lRef_;
  82         dict.lookup("Aref") >> Aref_;
  83     }
  84 }
  85
  86
  87 void Foam::forceCoeffs::writeFileHeader()
  88 {
  89     if (forcesFilePtr_.valid())
  90     {
  91         forcesFilePtr_()
  92             << "# Time" << tab << "Cd" << tab << "Cl" << tab << "Cm" << endl;
  93     }
  94 }
  95
  96
  97 void Foam::forceCoeffs::execute()
  98 {
  99     // Do nothing - only valid on write
 100 }
 101
 102
 103 void Foam::forceCoeffs::end()
 104 {
 105     // Do nothing - only valid on write
 106 }
 107
 108
 109 void Foam::forceCoeffs::write()
 110 {
 111     if (active_)
 112     {
 113         // Create the forces file if not already created
 114         makeFile();
 115
 116         forcesMoments fm = forces::calcForcesMoment();
 117
 118         scalar pDyn = 0.5*rhoRef_*magUInf_*magUInf_;
 119
 120         vector totForce = fm.first().first() + fm.first().second();
 121         vector totMoment = fm.second().first() + fm.second().second();
 122
 123         scalar liftForce = totForce & liftDir_;
 124         scalar dragForce = totForce & dragDir_;
 125         scalar pitchMoment = totMoment & pitchAxis_;
 126
 127         scalar Cl = liftForce/(Aref_*pDyn);
 128         scalar Cd = dragForce/(Aref_*pDyn);
 129         scalar Cm = pitchMoment/(Aref_*lRef_*pDyn);
 130
 131         if (Pstream::master())
 132         {
 133             forcesFilePtr_()
 134                 << obr_.time().value() << tab
 135                 << Cd << tab << Cl << tab << Cm << endl;
 136
 137             if (log_)
 138             {
 139                 Info<< "forceCoeffs output:" << nl
 140                     << "    Cd = " << Cd << nl
 141                     << "    Cl = " << Cl << nl
 142                     << "    Cm = " << Cm << nl
 143                     << endl;
 144             }
 145         }
 146     }
 147 }
 148
 149
 150 // ************************************************************************* //