DAPimpleFoam.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
 
    DAFoam  : Discrete Adjoint with OpenFOAM
    Version : v4
 
    Description:
        Child class for DAPimpleFoam
    
    This class is modified from OpenFOAM's source code
    applications/solvers/incompressible/pimpleFoam
 
    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/>.
 
\*---------------------------------------------------------------------------*/
 
#ifndef DAPimpleFoam_H
#define DAPimpleFoam_H
 
#include "DASolver.H"
#include "addToRunTimeSelectionTable.H"
#include "singlePhaseTransportModel.H"
#include "turbulentTransportModel.H"
#include "pimpleControlDF.H"
#include "DAFvSource.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
/*---------------------------------------------------------------------------*\
      Class DAPimpleFoam Declaration
\*---------------------------------------------------------------------------*/
 
class DAPimpleFoam
    : public DASolver
{
 
protected:
    autoPtr<pimpleControlDF> pimplePtr_;
 
    autoPtr<volScalarField> pPtr_;
 
    autoPtr<volVectorField> UPtr_;
 
    autoPtr<surfaceScalarField> phiPtr_;
 
    autoPtr<singlePhaseTransportModel> laminarTransportPtr_;
 
    autoPtr<incompressible::turbulenceModel> turbulencePtr_;
 
    autoPtr<DATurbulenceModel> daTurbulenceModelPtr_;
 
    autoPtr<DAFvSource> daFvSourcePtr_;
 
    autoPtr<volVectorField> fvSourcePtr_;
 
    label hasFvSource_ = 0;
 
    scalar cumulativeContErr_ = 0.0;
 
    label pRefCell_ = 0;
 
    scalar pRefValue_ = 0.0;
 
    label IOInterval_ = 1;
 
    label reduceIOWriteMesh_ = 0;
 
    autoPtr<dimensionedScalar> PrPtr_;
 
    autoPtr<dimensionedScalar> PrtPtr_;
 
    autoPtr<volScalarField> TPtr_;
 
    autoPtr<volScalarField> alphatPtr_;
 
    label hasTField_;
 
public:
    TypeName("DAPimpleFoam");
    // Constructors
 
    //- Construct from components
    DAPimpleFoam(
        char* argsAll,
        PyObject* pyOptions);
 
    //- Destructor
    virtual ~DAPimpleFoam()
    {
    }
 
    virtual void initSolver();
 
    virtual label solvePrimal();
 
    virtual label solveAdjointFP(
        Vec dFdW,
        Vec psi);
 
    scalar calcAdjointResiduals(
        const double* psi,
        const double* dFdW,
        double* adjRes);
};
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //