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: ************************ //