ReactorSurface.cpp

Go to the documentation of this file.

//! @file ReactorSurface.cpp
 
// This file is part of Cantera. See License.txt in the top-level directory or
// at https://cantera.org/license.txt for license and copyright information.
 
#include "cantera/zeroD/ReactorSurface.h"
#include "cantera/zeroD/ReactorNet.h"
#include "cantera/thermo/SurfPhase.h"
#include "cantera/kinetics/Kinetics.h"
#include "cantera/kinetics/InterfaceKinetics.h"
#include "cantera/kinetics/Reaction.h"
 
namespace Cantera
{
 
ReactorSurface::ReactorSurface(shared_ptr<Solution> sol, const string& name)
    : ReactorBase(sol, name)
{
    if (!std::dynamic_pointer_cast<SurfPhase>(sol->thermo())) {
        throw CanteraError("ReactorSurface::ReactorSurface",
            "Solution object must have a SurfPhase object as the thermo manager.");
    }
 
    if (!sol->kinetics() ) {
        throw CanteraError("ReactorSurface::ReactorSurface",
            "Solution object must have kinetics manager.");
    } else if (!std::dynamic_pointer_cast<InterfaceKinetics>(sol->kinetics())) {
        throw CanteraError("ReactorSurface::ReactorSurface",
            "Kinetics manager must be an InterfaceKinetics object.");
    }
    // todo: move all member variables to use shared pointers after Cantera 3.2
    m_kinetics = m_solution->kinetics().get();
    m_thermo = m_solution->thermo().get();
    m_surf = dynamic_cast<SurfPhase*>(m_thermo);
    m_cov.resize(m_surf->nSpecies());
    m_surf->getCoverages(m_cov.data());
}
 
double ReactorSurface::area() const
{
    return m_area;
}
 
void ReactorSurface::setArea(double a)
{
    m_area = a;
}
 
void ReactorSurface::setKinetics(Kinetics* kin)
{
    warn_deprecated("ReactorSurface::setKinetics",
                    "To be removed after Cantera 3.2.");
    m_kinetics = kin;
    if (kin == nullptr) {
        m_surf = nullptr;
        return;
    }
 
    m_surf = dynamic_cast<SurfPhase*>(&kin->thermo(0));
    if (m_surf == nullptr) {
        throw CanteraError("ReactorSurface::setKinetics",
            "Specified kinetics manager does not represent a surface "
            "kinetics mechanism.");
    }
    m_cov.resize(m_surf->nSpecies());
    m_surf->getCoverages(m_cov.data());
}
 
void ReactorSurface::setKinetics(Kinetics& kin)
{
    setKinetics(&kin);
}
 
void ReactorSurface::setReactor(ReactorBase* reactor)
{
    m_reactor = reactor;
}
 
void ReactorSurface::setCoverages(const double* cov)
{
    copy(cov, cov + m_cov.size(), m_cov.begin());
}
 
void ReactorSurface::setCoverages(const Composition& cov)
{
    m_surf->setCoveragesByName(cov);
    m_surf->getCoverages(m_cov.data());
}
 
void ReactorSurface::setCoverages(const string& cov)
{
    m_surf->setCoveragesByName(cov);
    m_surf->getCoverages(m_cov.data());
}
 
void ReactorSurface::getCoverages(double* cov) const
{
    copy(m_cov.begin(), m_cov.end(), cov);
}
 
void ReactorSurface::syncState()
{
    m_surf->setTemperature(m_reactor->temperature());
    m_surf->setCoveragesNoNorm(m_cov.data());
}
 
void ReactorSurface::addSensitivityReaction(size_t rxn)
{
    if (rxn >= m_kinetics->nReactions()) {
        throw CanteraError("ReactorSurface::addSensitivityReaction",
                           "Reaction number out of range ({})", rxn);
    }
    size_t p = m_reactor->network().registerSensitivityParameter(
        m_kinetics->reaction(rxn)->equation(), 1.0, 1.0);
    m_sensParams.emplace_back(
        SensitivityParameter{rxn, p, 1.0, SensParameterType::reaction});
}
 
void ReactorSurface::setSensitivityParameters(const double* params)
{
    for (auto& p : m_sensParams) {
        p.value = m_kinetics->multiplier(p.local);
        m_kinetics->setMultiplier(p.local, p.value*params[p.global]);
    }
}
 
void ReactorSurface::resetSensitivityParameters()
{
    for (auto& p : m_sensParams) {
        m_kinetics->setMultiplier(p.local, p.value);
    }
}
 
}