|
|
|
laplacianFoam.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 laplacianFoam 00027 00028 Description 00029 Solves a simple Laplace equation, e.g. for thermal diffusion in a solid. 00030 00031 \*---------------------------------------------------------------------------*/ 00032 00033 #include "fvCFD.H" 00034 00035 00036 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00037 00038 int main(int argc, char *argv[]) 00039 { 00040 00041 # include "setRootCase.H" 00042 00043 # include "createTime.H" 00044 # include "createMesh.H" 00045 # include "createFields.H" 00046 00047 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00048 00049 Info<< "\nCalculating temperature distribution\n" << endl; 00050 00051 while (runTime.loop()) 00052 { 00053 Info<< "Time = " << runTime.timeName() << nl << endl; 00054 00055 # include "readSIMPLEControls.H" 00056 00057 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) 00058 { 00059 solve 00060 ( 00061 fvm::ddt(T) - fvm::laplacian(DT, T) 00062 ); 00063 } 00064 00065 # include "write.H" 00066 00067 Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" 00068 << " ClockTime = " << runTime.elapsedClockTime() << " s" 00069 << nl << endl; 00070 } 00071 00072 Info<< "End\n" << endl; 00073 00074 return 0; 00075 } 00076 00077 00078 // ************************************************************************* // |