76 lines
1.8 KiB
C++
76 lines
1.8 KiB
C++
//! @file FlowDevice.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/FlowDevice.h"
|
|
#include "cantera/zeroD/ReactorBase.h"
|
|
#include "cantera/numerics/Func1.h"
|
|
|
|
namespace Cantera
|
|
{
|
|
|
|
FlowDevice::FlowDevice() : m_mdot(0.0), m_pfunc(0), m_tfunc(0),
|
|
m_coeff(1.0), m_type(0),
|
|
m_nspin(0), m_nspout(0),
|
|
m_in(0), m_out(0) {}
|
|
|
|
bool FlowDevice::install(ReactorBase& in, ReactorBase& out)
|
|
{
|
|
if (m_in || m_out) {
|
|
throw CanteraError("FlowDevice::install", "Already installed");
|
|
}
|
|
m_in = ∈
|
|
m_out = &out;
|
|
m_in->addOutlet(*this);
|
|
m_out->addInlet(*this);
|
|
|
|
// construct adapters between inlet and outlet species
|
|
const ThermoPhase& mixin = m_in->contents();
|
|
const ThermoPhase& mixout = m_out->contents();
|
|
|
|
m_nspin = mixin.nSpecies();
|
|
m_nspout = mixout.nSpecies();
|
|
std::string nm;
|
|
size_t ki, ko;
|
|
for (ki = 0; ki < m_nspin; ki++) {
|
|
nm = mixin.speciesName(ki);
|
|
ko = mixout.speciesIndex(nm);
|
|
m_in2out.push_back(ko);
|
|
}
|
|
for (ko = 0; ko < m_nspout; ko++) {
|
|
nm = mixout.speciesName(ko);
|
|
ki = mixin.speciesIndex(nm);
|
|
m_out2in.push_back(ki);
|
|
}
|
|
return true;
|
|
}
|
|
|
|
void FlowDevice::setPressureFunction(Func1* f)
|
|
{
|
|
m_pfunc = f;
|
|
}
|
|
|
|
void FlowDevice::setTimeFunction(Func1* g)
|
|
{
|
|
m_tfunc = g;
|
|
}
|
|
|
|
double FlowDevice::outletSpeciesMassFlowRate(size_t k)
|
|
{
|
|
if (k >= m_nspout) {
|
|
return 0.0;
|
|
}
|
|
size_t ki = m_out2in[k];
|
|
if (ki == npos) {
|
|
return 0.0;
|
|
}
|
|
return m_mdot * m_in->massFraction(ki);
|
|
}
|
|
|
|
double FlowDevice::enthalpy_mass()
|
|
{
|
|
return m_in->enthalpy_mass();
|
|
}
|
|
|
|
}
|