00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 2010-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 Global 00025 setMultiRegionDeltaT 00026 00027 Description 00028 Reset the timestep to maintain a constant maximum courant and diffusion 00029 Numbers. Reduction of time-step is immediate, but increase is damped to avoid 00030 unstable oscillations. 00031 00032 \*---------------------------------------------------------------------------*/ 00033 00034 if (adjustTimeStep) 00035 { 00036 if (CoNum == -GREAT) 00037 { 00038 CoNum = SMALL; 00039 } 00040 00041 if (DiNum == -GREAT) 00042 { 00043 DiNum = SMALL; 00044 } 00045 00046 scalar maxDeltaTFluid = maxCo/(CoNum + SMALL); 00047 scalar maxDeltaTSolid = maxDi/(DiNum + SMALL); 00048 00049 scalar deltaTFluid = 00050 min 00051 ( 00052 min(maxDeltaTFluid, 1.0 + 0.1*maxDeltaTFluid), 00053 1.2 00054 ); 00055 00056 runTime.setDeltaT 00057 ( 00058 min 00059 ( 00060 min(deltaTFluid, maxDeltaTSolid)*runTime.deltaT().value(), 00061 maxDeltaT 00062 ) 00063 ); 00064 00065 Info<< "deltaT = " << runTime.deltaT().value() << endl; 00066 } 00067 00068 // ************************ vim: set sw=4 sts=4 et: ************************ //
