SIMPLEControlOpt.C

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
    Copyright (C) 2007-2023 PCOpt/NTUA
    Copyright (C) 2013-2023 FOSS GP
    Copyright (C) 2019-2020 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.
 
    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/>.
 
\*---------------------------------------------------------------------------*/
 
#include "SIMPLEControlOpt.H"
#include "addToRunTimeSelectionTable.H"
#include "Time.H"
 
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
namespace Foam
{
    defineTypeNameAndDebug(SIMPLEControlOpt, 1);
    addToRunTimeSelectionTable( SIMPLEControl, SIMPLEControlOpt, dictionary);
}
 
// * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
 
bool Foam::SIMPLEControlOpt::read()
{
    nInitialIters_ = dict().getOrDefault<label>("nInitialIters", nIters_);
    return SIMPLEControl::read();
}
 
 
bool Foam::SIMPLEControlOpt::criteriaSatisfied()
{
    bool satisfied(false);
 
    // Do not check criteria in the first iteration of the algorithm.
    // Used to avoid stopping the solution of the flow equations
    // due to a converged solution in the previous optimisation cycle
    if (subCycledTimePtr_().index() == 1)
    {
        satisfied = false;
    }
    else
    {
        satisfied = simpleControl::criteriaSatisfied();
    }
 
    return satisfied;
}
 
 
const Foam::label& Foam::SIMPLEControlOpt::nIters() const
{
    if (mesh_.time().timeIndex() == mesh_.time().startTimeIndex() + 1)
    {
        return nInitialIters_;
    }
    else
    {
        return nIters_;
    }
}
 
 
void Foam::SIMPLEControlOpt::resetDeltaT()
{
    Time& runTime = const_cast<Time&>(mesh_.time());
    if (runTime.deltaTValue() != deltaTSubSycle_)
    {
        runTime.setDeltaT(deltaTSubSycle_, false);
    }
}
 
 
// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
Foam::SIMPLEControlOpt::SIMPLEControlOpt
(
    fvMesh& mesh,
    const word& managerType,
    const solver& solver
)
:
    SIMPLEControl(mesh, managerType, solver),
    nInitialIters_(0),
    subCycledTimePtr_(nullptr),
    deltaTSubSycle_(Zero)
{
    this->read();
}
 
 
// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
bool Foam::SIMPLEControlOpt::write(const bool valid) const
{
    // Intentionally does nothing
    // Write is called only at the last iteration of the cycle
    return false;
}
 
 
bool Foam::SIMPLEControlOpt::loop()
{
    this->read();
 
    Time& runTime = const_cast<Time&>(mesh_.time());
 
    bool doNextIter(true);
 
    if (nIters() > 0)
    {
        // Sub-cycle time if this is the first iter
        if (!subCycledTimePtr_)
        {
            subCycledTimePtr_.reset(new subCycleTime(runTime, nIters()));
            Info<< "Solving equations for solver "
                << solver_.solverName() << "\n" << endl;
            deltaTSubSycle_ = runTime.deltaTValue();
 
            // Reset iteration count to zero
            iter_ = 0;
 
            // Reset previous time index of fvMesh::data, to avoid the rare
            // occurance of a solver satisfying the convergence criteria at the
            // first iteration, which then causes all subsequent optimisation
            // cycles to be seen as converged, irrespective of the residual
            // level, since the data::prevTimeIndex_ is not updated
            //mesh_.data::setPreviousTimeIndex(0);
        }
 
        // Increase index
        subCycledTimePtr_()++;
        iter_ = subCycledTimePtr_().index();
 
 
        if (criteriaSatisfied())
        {
            Info<< nl
                << solver_.solverName()
                << " solution converged in "
                << subCycledTimePtr_->index() << " iterations" << nl << endl;
 
            subCycledTimePtr_->endSubCycle();
            subCycledTimePtr_.clear();
 
            // Write solution before continuing to next solver
            runTime.write();
            solver_.write();
 
            // Check whether mean fields have not been computed due to an
            // unexpectedly early convergence
            checkMeanSolution();
 
            doNextIter = false;
        }
        else if (subCycledTimePtr_->end())
        {
            Info<< nl
                << solver_.solverName()
                << " solution reached max. number of iterations "
                << subCycledTimePtr_().nSubCycles() << nl << endl;
 
            subCycledTimePtr_->endSubCycle();
            subCycledTimePtr_.clear();
 
            // Write solution before continuing to next solver
            runTime.write();
            solver_.write();
 
            doNextIter = false;
        }
        else
        {
            // Since dicts are not updated when Time is sub-cycled,
            // do it manually here
            runTime.readModifiedObjects();
            resetDeltaT();
 
            DebugInfo
                << "Iteration " << subCycledTimePtr_().index()
                << "|" << subCycledTimePtr_().nSubCycles() << endl;
 
            simpleControl::storePrevIterFields();
 
            doNextIter = true;
        }
    }
    else
    {
        WarningInFunction
            << "Number of iterations is non-positive (" << nIters() << ")."
            << nl
            << "Skipping the solution of the equations corresponding to solver "
            << solver_.solverName()
            << nl << endl;
 
        doNextIter = false;
    }
 
    return doNextIter;
}
 
 
// ************************************************************************* //