 |
|
|
|
|
rhoEqn.HGo 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 Global 00026 rhoEqn 00027 00028 Description 00029 Solve the continuity for density. 00030 00031 \*---------------------------------------------------------------------------*/ 00032 00033 volScalarField Sevap 00034 ( 00035 IOobject 00036 ( 00037 "Sevap", 00038 runTime.timeName(), 00039 mesh, 00040 IOobject::NO_READ, 00041 IOobject::NO_WRITE 00042 ), 00043 mesh, 00044 dimensionedScalar("zero", dimensionSet(1, -3, -1, 0, 0), 0.0) 00045 ); 00046 00047 for (label i=0; i<Y.size(); i++) 00048 { 00049 if (dieselSpray.isLiquidFuel()[i]) 00050 { 00051 Sevap += dieselSpray.evaporationSource(i); 00052 } 00053 } 00054 00055 { 00056 solve 00057 ( 00058 fvm::ddt(rho) 00059 + fvc::div(phi) 00060 == 00061 Sevap 00062 ); 00063 } 00064 00065 // ************************************************************************* // |