applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/createFields.H

   1     Info<< "\nReading transportProperties\n" << endl;
   2
   3     IOdictionary transportProperties
   4     (
   5         IOobject
   6         (
   7             "transportProperties",
   8             runTime.constant(),
   9             mesh,
  10             IOobject::MUST_READ,
  11             IOobject::NO_WRITE
  12         )
  13     );
  14
  15
  16     dimensionedScalar nu
  17     (
  18         transportProperties.lookup("nu")
  19     );
  20
  21
  22     dimensionedScalar rhoFluid
  23     (
  24         transportProperties.lookup("rho")
  25     );
  26
  27
  28     Info<< "Reading field p\n" << endl;
  29     volScalarField p
  30     (
  31         IOobject
  32         (
  33             "p",
  34             runTime.timeName(),
  35             mesh,
  36             IOobject::MUST_READ,
  37             IOobject::AUTO_WRITE
  38         ),
  39         mesh
  40     );
  41
  42
  43     Info<< "Reading field U\n" << endl;
  44     volVectorField U
  45     (
  46         IOobject
  47         (
  48             "U",
  49             runTime.timeName(),
  50             mesh,
  51             IOobject::MUST_READ,
  52             IOobject::AUTO_WRITE
  53         ),
  54         mesh
  55     );
  56
  57 #   include "createPhi.H"
  58
  59     label pRefCell = 0;
  60     scalar pRefValue = 0.0;
  61     setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
  62
  63 //     pRefCell = mesh.findNearestCell(vector(0,0.205,0));
  64 //     Info << "Ref. cell: " << pRefCell << endl;
  65 //     Info << "Ref. cell position: " << mesh.C()[pRefCell] << endl;