 |
|
|
|
|
scalarTransportFoam.CGo 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-2009 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 the 00013 Free Software Foundation; either version 2 of the License, or (at your 00014 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, write to the Free Software Foundation, 00023 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA 00024 00025 Application 00026 scalarTransportFoam 00027 00028 Description 00029 Solves a transport equation for a passive scalar 00030 00031 \*---------------------------------------------------------------------------*/ 00032 00033 #include "fvCFD.H" 00034 00035 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00036 00037 int main(int argc, char *argv[]) 00038 { 00039 00040 # include "setRootCase.H" 00041 00042 # include "createTime.H" 00043 # include "createMesh.H" 00044 # include "createFields.H" 00045 00046 00047 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00048 00049 Info<< "\nCalculating scalar transport\n" << endl; 00050 00051 # include "CourantNo.H" 00052 00053 while (runTime.loop()) 00054 { 00055 Info<< "Time = " << runTime.timeName() << nl << endl; 00056 00057 # include "readSIMPLEControls.H" 00058 00059 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) 00060 { 00061 solve 00062 ( 00063 fvm::ddt(T) 00064 + fvm::div(phi, T) 00065 - fvm::laplacian(DT, T) 00066 ); 00067 } 00068 00069 runTime.write(); 00070 } 00071 00072 Info<< "End\n" << endl; 00073 00074 return 0; 00075 } 00076 00077 00078 // ************************************************************************* // |