applications/solvers/coupled/pUCoupledFoam/pUCoupledFoam.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     pUCoupledFoam
  26
  27 Description
  28     Steady-state solver for incompressible, turbulent flow, with implicit
  29     coupling between pressure and velocity achieved by fvBlockMatrix.
  30     Turbulence is in this version solved using the existing turbulence
  31     structure.
  32
  33 Authors
  34     Klas Jareteg, Chalmers University of Technology,
  35     Vuko Vukcevic, FMENA Zagreb.
  36
  37 \*---------------------------------------------------------------------------*/
  38
  39 #include "fvCFD.H"
  40 #include "fvBlockMatrix.H"
  41 #include "singlePhaseTransportModel.H"
  42 #include "RASModel.H"
  43
  44 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  45
  46 int main(int argc, char *argv[])
  47 {
  48
  49 #   include "setRootCase.H"
  50 #   include "createTime.H"
  51 #   include "createMesh.H"
  52 #   include "createFields.H"
  53 #   include "initContinuityErrs.H"
  54 #   include "initConvergenceCheck.H"
  55
  56     Info<< "\nStarting time loop\n" << endl;
  57     while (runTime.loop())
  58     {
  59 #       include "readBlockSolverControls.H"
  60 #       include "readFieldBounds.H"
  61
  62         Info<< "Time = " << runTime.timeName() << nl << endl;
  63
  64         p.storePrevIter();
  65
  66         // Initialize the Up block system (matrix, source and reference to Up)
  67         fvBlockMatrix<vector4> UpEqn(Up);
  68
  69         // Assemble and insert momentum equation
  70 #       include "UEqn.H"
  71
  72         // Assemble and insert pressure equation
  73 #       include "pEqn.H"
  74
  75         // Assemble and insert coupling terms
  76 #       include "couplingTerms.H"
  77
  78         // Solve the block matrix
  79         maxResidual = cmptMax(UpEqn.solve().initialResidual());
  80
  81         // Retrieve solution
  82         UpEqn.retrieveSolution(0, U.internalField());
  83         UpEqn.retrieveSolution(3, p.internalField());
  84
  85         U.correctBoundaryConditions();
  86         p.correctBoundaryConditions();
  87
  88         phi = (fvc::interpolate(U) & mesh.Sf()) + pEqn.flux() + presSource;
  89
  90 #       include "continuityErrs.H"
  91
  92 #       include "boundPU.H"
  93
  94         p.relax();
  95
  96         turbulence->correct();
  97         runTime.write();
  98
  99         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
 100             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
 101             << nl << endl;
 102
 103 #       include "convergenceCheck.H"
 104     }
 105
 106     Info<< "End\n" << endl;
 107
 108     return 0;
 109 }