applications/test/nearWallDist-wave/testYPlus.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 Description
  26     Calculate distance to wall.
  27
  28 \*---------------------------------------------------------------------------*/
  29
  30 #include <finiteVolume/fvCFD.H>
  31 #include <finiteVolume/wallDistData.H>
  32 #include <finiteVolume/wallPointYPlus.H>
  33
  34 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  35 // Main program:
  36
  37 int main(int argc, char *argv[])
  38 {
  39 #   include <OpenFOAM/setRootCase.H>
  40 #   include <OpenFOAM/createTime.H>
  41 #   include <OpenFOAM/createMesh.H>
  42
  43     Info<< "Mesh read in = "
  44         << runTime.cpuTimeIncrement()
  45         << " s\n" << endl << endl;
  46
  47
  48     Info<< "Time now = " << runTime.timeName() << endl;
  49
  50     // wall distance and yStar
  51
  52     volScalarField yStar
  53     (
  54         IOobject
  55         (
  56             "yStar",
  57             mesh.time().timeName(),
  58             mesh
  59         ),
  60         mesh,
  61         dimensionedScalar("yStar", dimless, 1.0)
  62     );
  63
  64     // Fill wall patches of yStar with some value.
  65     forAll(mesh.boundary(), patchI)
  66     {
  67         const fvPatch& patch = mesh.boundary()[patchI];
  68
  69         if (typeid(patch) == typeid(wallFvPatch))
  70         {
  71             fvPatchScalarField& wallData = yStar.boundaryField()[patchI];
  72
  73             forAll(patch, patchFaceI)
  74             {
  75 // Hack. Just some value.
  76                 wallData[patchFaceI] = 1/2500.0;
  77             }
  78         }
  79     }
  80
  81
  82     // Do distance calculation, transporting values of yStar
  83     wallPointYPlus::yPlusCutOff = 200;
  84     wallDistData<wallPointYPlus> y(mesh, yStar, true);
  85
  86     if (y.nUnset() != 0)
  87     {
  88         WarningIn(args.executable())
  89             << "There are " << y.nUnset()
  90             << " remaining unset cells and/or boundary values" << endl;
  91     }
  92
  93
  94     y.write();
  95
  96     y.data().write();
  97
  98     volScalarField yPlus
  99     (
 100         IOobject
 101         (
 102             "yPlus",
 103             mesh.time().timeName(),
 104             mesh
 105         ),
 106         y/y.data()
 107     );
 108
 109     yPlus.write();
 110
 111     Info<< "End\n" << endl;
 112
 113     return 0;
 114
 115
 116
 117 }
 118
 119
 120 // ************************ vim: set sw=4 sts=4 et: ************************ //