ReactorBase.cpp

Go to the documentation of this file.

//! @file ReactorBase.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/ReactorBase.h"
#include "cantera/zeroD/FlowDevice.h"
#include "cantera/zeroD/ReactorNet.h"
#include "cantera/zeroD/ReactorSurface.h"
#include "cantera/base/Solution.h"
#include "cantera/thermo/ThermoPhase.h"
 
namespace Cantera
{
 
ReactorBase::ReactorBase(const string& name)
{
    m_name = name;
}
 
ReactorBase::ReactorBase(shared_ptr<Solution> sol, const string& name)
{
    m_name = name;
    setSolution(sol);
}
 
void ReactorBase::setSolution(shared_ptr<Solution> sol) {
    m_solution = sol;
    setThermoMgr(*sol->thermo());
    try {
        setKineticsMgr(*sol->kinetics());
    } catch (NotImplementedError&) {
        // kinetics not used (example: Reservoir)
    }
}
 
void ReactorBase::insert(shared_ptr<Solution> sol)
{
    warn_deprecated("ReactorBase::insert",
        "To be removed after Cantera 3.1. Superseded by 'setSolution'.");
    setSolution(sol);
}
 
void ReactorBase::setThermoMgr(ThermoPhase& thermo)
{
    m_thermo = &thermo;
    m_nsp = m_thermo->nSpecies();
    m_thermo->saveState(m_state);
    m_enthalpy = m_thermo->enthalpy_mass();
    m_intEnergy = m_thermo->intEnergy_mass();
    m_pressure = m_thermo->pressure();
}
 
void ReactorBase::syncState()
{
    m_thermo->saveState(m_state);
    m_enthalpy = m_thermo->enthalpy_mass();
    m_intEnergy = m_thermo->intEnergy_mass();
    m_pressure = m_thermo->pressure();
    if (m_net) {
        m_net->setNeedsReinit();
    }
}
 
void ReactorBase::addInlet(FlowDevice& inlet)
{
    m_inlet.push_back(&inlet);
}
 
void ReactorBase::addOutlet(FlowDevice& outlet)
{
    m_outlet.push_back(&outlet);
}
 
void ReactorBase::addWall(WallBase& w, int lr)
{
    m_wall.push_back(&w);
    if (lr == 0) {
        m_lr.push_back(0);
    } else {
        m_lr.push_back(1);
    }
}
 
WallBase& ReactorBase::wall(size_t n)
{
    return *m_wall[n];
}
 
void ReactorBase::addSurface(ReactorSurface* surf)
{
    if (find(m_surfaces.begin(), m_surfaces.end(), surf) == m_surfaces.end()) {
        m_surfaces.push_back(surf);
        surf->setReactor(this);
    }
}
 
ReactorSurface* ReactorBase::surface(size_t n)
{
    return m_surfaces[n];
}
 
void ReactorBase::restoreState() {
    if (!m_thermo) {
        throw CanteraError("ReactorBase::restoreState", "No phase defined.");
    }
    m_thermo->restoreState(m_state);
}
 
ReactorNet& ReactorBase::network()
{
    if (m_net) {
        return *m_net;
    } else {
        throw CanteraError("ReactorBase::network",
                           "Reactor is not part of a ReactorNet");
    }
}
 
void ReactorBase::setNetwork(ReactorNet* net)
{
    m_net = net;
}
 
double ReactorBase::residenceTime()
{
    double mout = 0.0;
    for (size_t i = 0; i < m_outlet.size(); i++) {
        mout += m_outlet[i]->massFlowRate();
    }
    return mass()/mout;
}
 
FlowDevice& ReactorBase::inlet(size_t n)
{
    return *m_inlet[n];
}
FlowDevice& ReactorBase::outlet(size_t n)
{
    return *m_outlet[n];
}
 
}