solveFluid.H
Go to the documentation of this file.
1 if (pimple.frozenFlow())
2 {
3  #include "EEqn.H"
4 }
5 else
6 {
7  if (!mesh.steady() && pimples.firstPimpleIter())
8  {
9  #include "rhoEqn.H"
10  }
11 
12  #include "UEqn.H"
13  #include "YEqn.H"
14  #include "EEqn.H"
15 
16  // --- PISO loop
17  while (pimple.correct())
18  {
19  #include "../../buoyantPimpleFoam/pEqn.H"
20  }
21 
22  if (pimples.pimpleTurbCorr(i))
23  {
24  turbulence.correct();
26  }
27 
28  if (!mesh.steady() && pimples.finalPimpleIter())
29  {
30  rho = thermo.rho();
31  }
32 }
pimpleNoLoopControl & pimple
rhoReactionThermo & thermo
Definition: createFields.H:28
dynamicFvMesh & mesh
rhoReactionThermophysicalTransportModel & thermophysicalTransport
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
Definition: createFields.H:94