flowControllers.cpp

Go to the documentation of this file.

//! @file flowControllers.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/flowControllers.h"
#include "cantera/zeroD/ReactorBase.h"
#include "cantera/numerics/Func1.h"
 
namespace Cantera
{
 
void MassFlowController::setMassFlowRate(double mdot)
{
    m_tfunc.reset();
    m_coeff = mdot;
}
 
void MassFlowController::updateMassFlowRate(double time)
{
    if (!ready()) {
        throw CanteraError("MassFlowController::updateMassFlowRate",
                           "Device is not ready; some parameters have not been set.");
    }
    double mdot = m_coeff;
    if (m_tfunc) {
        mdot *= m_tfunc->eval(time);
    }
    m_mdot = std::max(mdot, 0.0);
}
 
void PressureController::updateMassFlowRate(double time)
{
    if (!ready()) {
        throw CanteraError("PressureController::updateMassFlowRate",
                           "Device is not ready; some parameters have not been set.");
    }
    double mdot = m_coeff;
    double delta_P = in().pressure() - out().pressure();
    if (m_pfunc) {
        mdot *= m_pfunc->eval(delta_P);
    } else {
        mdot *= delta_P;
    }
    m_primary->updateMassFlowRate(time);
    mdot += m_primary->massFlowRate();
    m_mdot = std::max(mdot, 0.0);
}
 
double PressureController::massFlowRate_ddP() const
{
    if (m_mdot <= 0.0) {
        return 0.0;
    }
    double delta_P = in().pressure() - out().pressure();
    return m_coeff * pressureFunction_ddP(delta_P);
}
 
void PressureController::addMassFlowRateJacobian(
    SparseTriplets& trips, size_t row, double coeff, bool includePressureSpecies)
{
    FlowDevice::addMassFlowRateJacobian(trips, row, coeff, includePressureSpecies);
    if (m_primary) {
        m_primary->addMassFlowRateJacobian(trips, row, coeff, includePressureSpecies);
    }
}
 
void PressureController::setPrimary(shared_ptr<ConnectorNode> primary)
{
    auto dev = std::dynamic_pointer_cast<FlowDevice>(primary);
    if (!dev) {
        throw CanteraError("PressureController::setPrimary",
                           "Invalid primary mass flow controller with type {}.",
                           primary->type());
    }
    m_primary = dev.get();
}
 
void Valve::updateMassFlowRate(double time)
{
    if (!ready()) {
        throw CanteraError("Valve::updateMassFlowRate",
                           "Device is not ready; some parameters have not been set.");
    }
    double mdot = m_coeff;
    if (m_tfunc) {
        mdot *= m_tfunc->eval(time);
    }
    double delta_P = in().pressure() - out().pressure();
    if (m_pfunc) {
        mdot *= m_pfunc->eval(delta_P);
    } else {
        mdot *= delta_P;
    }
    m_mdot = std::max(mdot, 0.0);
}
 
double Valve::massFlowRate_ddP() const
{
    if (m_mdot <= 0.0) {
        return 0.0;
    }
    double alpha = m_coeff;
    if (m_tfunc) {
        alpha *= m_tfunc->eval(m_time);
    }
    double delta_P = in().pressure() - out().pressure();
    alpha *= pressureFunction_ddP(delta_P);
    return alpha;
}
 
}