Boundary1D.h

Go to the documentation of this file.

/**
 * @file Boundary1D.h
 *
 * Boundary objects for one-dimensional simulations.
 */
 
// 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.
 
#ifndef CT_BOUNDARY1D_H
#define CT_BOUNDARY1D_H
 
#include "Domain1D.h"
#include "cantera/thermo/SurfPhase.h"
#include "cantera/kinetics/InterfaceKinetics.h"
#include "StFlow.h"
 
namespace Cantera
{
 
const int LeftInlet = 1;
const int RightInlet = -1;
 
/**
 * The base class for boundaries between one-dimensional spatial domains. The
 * boundary may have its own internal variables, such as surface species
 * coverages.
 *
 * The boundary types are an inlet, an outlet, a symmetry plane, and a surface.
 *
 * The public methods are all virtual, and the base class implementations throw
 * exceptions.
 * @ingroup onedim
 */
class Boundary1D : public Domain1D
{
public:
    Boundary1D();
 
    virtual void init() {
        _init(1);
    }
 
    /// Set the temperature.
    virtual void setTemperature(double t) {
        m_temp = t;
    }
 
    /// Temperature [K].
    virtual double temperature() {
        return m_temp;
    }
 
    virtual size_t nSpecies() {
        return 0;
    }
 
    /// Set the mole fractions by specifying a std::string.
    virtual void setMoleFractions(const std::string& xin) {
        throw NotImplementedError("Boundary1D::setMoleFractions");
    }
 
    /// Set the mole fractions by specifying an array.
    virtual void setMoleFractions(const double* xin) {
        throw NotImplementedError("Boundary1D::setMoleFractions");
    }
 
    /// Mass fraction of species k.
    virtual double massFraction(size_t k) {
        throw NotImplementedError("Boundary1D::massFraction");
    }
 
    /// Set the total mass flow rate.
    virtual void setMdot(double mdot) {
        m_mdot = mdot;
    }
 
    /// The total mass flow rate [kg/m2/s].
    virtual double mdot() {
        return m_mdot;
    }
 
    virtual void setupGrid(size_t n, const double* z) {}
 
protected:
    void _init(size_t n);
 
    StFlow* m_flow_left, *m_flow_right;
    size_t m_ilr, m_left_nv, m_right_nv;
    size_t m_left_loc, m_right_loc;
    size_t m_left_points;
    size_t m_left_nsp, m_right_nsp;
    size_t m_sp_left, m_sp_right;
    size_t m_start_left, m_start_right;
    ThermoPhase* m_phase_left, *m_phase_right;
    double m_temp, m_mdot;
};
 
 
/**
 * An inlet.
 * @ingroup onedim
 */
class Inlet1D : public Boundary1D
{
public:
    Inlet1D();
 
    /// set spreading rate
    virtual void setSpreadRate(double V0) {
        m_V0 = V0;
        needJacUpdate();
    }
 
    /// spreading rate
    virtual double spreadRate() {
        return m_V0;
    }
 
    virtual void showSolution(const double* x);
 
    virtual size_t nSpecies() {
        return m_nsp;
    }
 
    virtual void setMoleFractions(const std::string& xin);
    virtual void setMoleFractions(const double* xin);
    virtual double massFraction(size_t k) {
        return m_yin[k];
    }
    virtual void init();
    virtual void eval(size_t jg, double* xg, double* rg,
                      integer* diagg, double rdt);
    virtual XML_Node& save(XML_Node& o, const double* const soln);
    virtual void restore(const XML_Node& dom, double* soln, int loglevel);
    virtual AnyMap serialize(const double* soln) const;
    virtual void restore(const AnyMap& state, double* soln, int loglevel);
 
protected:
    int m_ilr;
    double m_V0;
    size_t m_nsp;
    vector_fp m_yin;
    std::string m_xstr;
    StFlow* m_flow;
};
 
/**
 * A terminator that does nothing.
 * @ingroup onedim
 */
class Empty1D : public Boundary1D
{
public:
    Empty1D() : Boundary1D() {
        m_type = cEmptyType;
    }
 
    virtual void showSolution(const double* x) {}
 
    virtual void init();
 
    virtual void eval(size_t jg, double* xg, double* rg,
                      integer* diagg, double rdt);
 
    virtual XML_Node& save(XML_Node& o, const double* const soln);
    virtual void restore(const XML_Node& dom, double* soln, int loglevel);
    virtual AnyMap serialize(const double* soln) const;
};
 
/**
 * A symmetry plane. The axial velocity u = 0, and all other components have
 * zero axial gradients.
 * @ingroup onedim
 */
class Symm1D : public Boundary1D
{
public:
    Symm1D() : Boundary1D() {
        m_type = cSymmType;
    }
 
    virtual void init();
 
    virtual void eval(size_t jg, double* xg, double* rg,
                      integer* diagg, double rdt);
 
    virtual XML_Node& save(XML_Node& o, const double* const soln);
    virtual void restore(const XML_Node& dom, double* soln, int loglevel);
    virtual AnyMap serialize(const double* soln) const;
};
 
 
/**
 * An outlet.
 * @ingroup onedim
 */
class Outlet1D : public Boundary1D
{
public:
    Outlet1D() : Boundary1D() {
        m_type = cOutletType;
    }
 
    virtual void init();
 
    virtual void eval(size_t jg, double* xg, double* rg,
                      integer* diagg, double rdt);
 
    virtual XML_Node& save(XML_Node& o, const double* const soln);
    virtual void restore(const XML_Node& dom, double* soln, int loglevel);
    virtual AnyMap serialize(const double* soln) const;
};
 
 
/**
 * An outlet with specified composition.
 * @ingroup onedim
 */
class OutletRes1D : public Boundary1D
{
public:
    OutletRes1D();
 
    virtual void showSolution(const double* x) {}
 
    virtual size_t nSpecies() {
        return m_nsp;
    }
 
    virtual void setMoleFractions(const std::string& xin);
    virtual void setMoleFractions(const double* xin);
    virtual double massFraction(size_t k) {
        return m_yres[k];
    }
    virtual void init();
    virtual void eval(size_t jg, double* xg, double* rg,
                      integer* diagg, double rdt);
    virtual XML_Node& save(XML_Node& o, const double* const soln);
    virtual void restore(const XML_Node& dom, double* soln, int loglevel);
    virtual AnyMap serialize(const double* soln) const;
    virtual void restore(const AnyMap& state, double* soln, int loglevel);
 
protected:
    size_t m_nsp;
    vector_fp m_yres;
    std::string m_xstr;
    StFlow* m_flow;
};
 
/**
 * A non-reacting surface. The axial velocity is zero (impermeable), as is the
 * transverse velocity (no slip). The temperature is specified, and a zero flux
 * condition is imposed for the species.
 * @ingroup onedim
 */
class Surf1D : public Boundary1D
{
public:
    Surf1D() : Boundary1D() {
        m_type = cSurfType;
    }
 
    virtual void init();
 
    virtual void eval(size_t jg, double* xg, double* rg,
                      integer* diagg, double rdt);
 
    virtual XML_Node& save(XML_Node& o, const double* const soln);
    virtual void restore(const XML_Node& dom, double* soln, int loglevel);
    virtual AnyMap serialize(const double* soln) const;
    virtual void restore(const AnyMap& state, double* soln, int loglevel);
 
    virtual void showSolution_s(std::ostream& s, const double* x);
 
    virtual void showSolution(const double* x);
};
 
/**
 * A reacting surface.
 * @ingroup onedim
 */
class ReactingSurf1D : public Boundary1D
{
public:
    ReactingSurf1D();
 
    void setKineticsMgr(InterfaceKinetics* kin);
 
    void enableCoverageEquations(bool docov) {
        m_enabled = docov;
    }
 
    bool coverageEnabled() {
        return m_enabled;
    }
 
    virtual std::string componentName(size_t n) const;
 
    virtual void init();
    virtual void resetBadValues(double* xg);
 
    virtual void eval(size_t jg, double* xg, double* rg,
                      integer* diagg, double rdt);
 
    virtual XML_Node& save(XML_Node& o, const double* const soln);
    virtual void restore(const XML_Node& dom, double* soln, int loglevel);
    virtual AnyMap serialize(const double* soln) const;
    virtual void restore(const AnyMap& state, double* soln, int loglevel);
 
    virtual void _getInitialSoln(double* x) {
        m_sphase->getCoverages(x);
    }
 
    virtual void _finalize(const double* x) {
        std::copy(x, x+m_nsp, m_fixed_cov.begin());
    }
 
    virtual void showSolution(const double* x);
 
protected:
    InterfaceKinetics* m_kin;
    SurfPhase* m_sphase;
    size_t m_surfindex, m_nsp;
    bool m_enabled;
    vector_fp m_work;
    vector_fp m_fixed_cov;
};
 
}
 
#endif