continuityErrsPython.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
 
    DAFoam  : Discrete Adjoint with OpenFOAM
    Version : v3
 
    This file is modified from OpenFOAM's source code
    src/finiteVolume/cfdTools/incompressible/continuityErrs.H
 
    OpenFOAM: The Open Source CFD Toolbox
 
    Copyright (C): 2011-2016 OpenFOAM Foundation
 
    OpenFOAM License:
 
        OpenFOAM is free software: you can redistribute it and/or modify it
        under the terms of the GNU General Public License as published by
        the Free Software Foundation, either version 3 of the License, or
        (at your option) any later version.
    
        OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
        ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
        FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
        for more details.
    
        You should have received a copy of the GNU General Public License
        along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
\*---------------------------------------------------------------------------*/
 
// Calculates and prints the continuity errors.
 
{
    volScalarField contErr(fvc::div(phi));
 
    scalar sumLocalContErr = runTime.deltaTValue() * mag(contErr)().weightedAverage(mesh.V()).value();
 
    scalar globalContErr = runTime.deltaTValue() * contErr.weightedAverage(mesh.V()).value();
    cumulativeContErr += globalContErr;
 
    Info << "Time step continuity errors : sum local = " << sumLocalContErr << endl
         << "                                 global = " << globalContErr << endl
         << "                             cumulative = " << cumulativeContErr << endl;
}