FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

rhoPorousMRFPimpleFoam.C

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software: you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by
00013     the Free Software Foundation, either version 3 of the License, or
00014     (at your option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
00023 
00024 Application
00025     rhoPorousMRFPimpleFoam
00026 
00027 Description
00028     Transient solver for laminar or turbulent flow of compressible fluids
00029     with support for porous media and MRF for HVAC and similar applications.
00030 
00031     Uses the flexible PIMPLE (PISO-SIMPLE) solution for time-resolved and
00032     pseudo-transient simulations.
00033 
00034 Usage
00035     - rhoPorousMRFPimpleFoam [OPTION]
00036 
00037     @param -case <dir> \n
00038     Specify the case directory
00039 
00040     @param -parallel \n
00041     Run the case in parallel
00042 
00043     @param -help \n
00044     Display short usage message
00045 
00046     @param -doc \n
00047     Display Doxygen documentation page
00048 
00049     @param -srcDoc \n
00050     Display source code
00051 
00052 \*---------------------------------------------------------------------------*/
00053 
00054 #include <finiteVolume/fvCFD.H>
00055 #include <basicThermophysicalModels/basicPsiThermo.H>
00056 #include <compressibleTurbulenceModel/turbulenceModel.H>
00057 #include <finiteVolume/bound.H>
00058 #include <finiteVolume/MRFZones.H>
00059 #include <finiteVolume/porousZones.H>
00060 
00061 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00062 
00063 int main(int argc, char *argv[])
00064 {
00065     #include <OpenFOAM/setRootCase.H>
00066     #include <OpenFOAM/createTime.H>
00067     #include <OpenFOAM/createMesh.H>
00068     #include "createFields.H"
00069     #include <finiteVolume/initContinuityErrs.H>
00070 
00071     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00072 
00073     Info<< "\nStarting time loop\n" << endl;
00074 
00075     while (runTime.run())
00076     {
00077         #include <finiteVolume/readTimeControls.H>
00078         #include <finiteVolume/readPIMPLEControls.H>
00079         #include <finiteVolume/compressibleCourantNo.H>
00080         #include <finiteVolume/setDeltaT.H>
00081 
00082         runTime++;
00083 
00084         Info<< "Time = " << runTime.timeName() << nl << endl;
00085 
00086         if (nOuterCorr != 1)
00087         {
00088             p.storePrevIter();
00089             rho.storePrevIter();
00090         }
00091 
00092         #include <finiteVolume/rhoEqn.H>
00093 
00094         // --- Pressure-velocity PIMPLE corrector loop
00095         for (int oCorr=0; oCorr<nOuterCorr; oCorr++)
00096         {
00097             #include "UEqn.H"
00098             #include "../rhoPimpleFoam/hEqn.H"
00099 
00100             // --- PISO loop
00101             for (int corr=0; corr<nCorr; corr++)
00102             {
00103                 #include "pEqn.H"
00104             }
00105 
00106             turbulence->correct();
00107         }
00108 
00109         runTime.write();
00110 
00111         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
00112             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
00113             << nl << endl;
00114     }
00115 
00116     Info<< "End\n" << endl;
00117 
00118     return 0;
00119 }
00120 
00121 
00122 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines