CircuitCoupling.h

Go to the documentation of this file.

 
#ifndef RIZER_CIRCUIT_COUPLING_H
#define RIZER_CIRCUIT_COUPLING_H
 
#include "PlasmaChannel1D.h"
 
#include <cstddef>
#include <functional>
#include <vector>
 
namespace rizer
{

    struct CircuitSample
    {
        double Rp = 0.0;
        double Rp_state = 0.0;
        double Vp = 0.0;
        double current = 0.0;
        double electric_field = 0.0;
    };

    class ChannelCircuit
    {
    public:
        ChannelCircuit(const PlasmaChannel1D &channel, double gap,
                       const std::function<double(double, double)> &plasma_voltage,
                       std::size_t channel_offset, std::size_t n_vars);

        bool active() const { return static_cast<bool>(*m_plasma_voltage); }

        void seedResistance(const std::vector<double> &x);

        void closeOver(double t0, double t1, const std::vector<double> &x);

        double electricField(double t) const;

        CircuitSample sample(const std::vector<double> &x) const;
 
    private:
        double conductance(const std::vector<double> &x) const;

        double resistanceFrom(double conductance_value) const;
 
        const PlasmaChannel1D *m_channel; // Non-owning
        double m_gap;                     // Inter-electrode gap [m]
        const std::function<double(double, double)> *m_plasma_voltage; // Non-owning
        std::size_t m_channel_offset;    // Offset into the global state [-]
        std::size_t m_n_vars;            // Components per node, 2+K [-]
        std::size_t m_n_points;          // Radial grid points [-]
        std::vector<double> m_cell_area; // Per-node cross-section [m^2]
 
        // The circuit's state over the step currently being advanced.
        double m_t0 = 0.0;  // Start of that step [s]
        double m_dt = 1.0;  // Its length [s]
        double m_Vp0 = 0.0; // Terminal voltage at m_t0 [V]
        double m_Vp1 = 0.0; // Terminal voltage at m_t0 + m_dt [V]
        double m_Rp;        // Resistance the step is driven with [Ohm]
    };
 
} // namespace rizer
 
#endif // RIZER_CIRCUIT_COUPLING_H