MechanicalVentilatorModel.h

/* Distributed under the Apache License, Version 2.0.
   See accompanying NOTICE file for details.*/
 
#pragma once
 
#include "engine/common/system/Model.h"
#include "engine/common/system/Systems.h"
 
class SEEquipmentActionCollection;
class SEGasCompartment;
class SELiquidCompartment;
class SEGasSubstanceQuantity;
class SEFluidCircuitNode;
class SEFluidCircuitPath;
 
namespace pulse
{
  class PBEquipment;
  class EngineTest;
  class PULSE_DECL MechanicalVentilatorModel : public MechanicalVentilatorSystem, public Model
  {
    friend PBEquipment;//friend the serialization class
    friend EngineTest;
  public:
    MechanicalVentilatorModel(Data& pc);
    virtual ~MechanicalVentilatorModel();
 
    void Clear() override;
 
    // Set members to a stable homeostatic state
    void Initialize() override;
    // Set pointers and other member varialbes common to both homeostatic initialization and loading a state
    void SetUp() override;
 
    void StateChange() override;
 
    void AtSteadyState() override {}
    void PreProcess() override;
    void Process(bool solve_and_transport = true) override;
    void PostProcess(bool solve_and_transport = true) override;
 
  protected:
 
    void ComputeExposedModelParameters() override;
 
    void UpdateAirwayMode();
 
    void CheckInspirationTriggers();
    void CalculateInspiration();
    void CalculatePause();
    void CalculateExpiration();
    void SetVentilatorDriver();
    void CycleMode(bool patientTriggered);
    void SetLeak();
    void SetHold();
    void SetValves();
    void CalculateRespiratoryParameters();
    void CalculateInspiratoryRespiratoryParameters();
    void CalculatePauseRespiratoryParameters();
    void CalculateExpiratoryRespiratoryParameters();
    void SetResistances();
    void SetCompliance();
    void SetVolumes();
    void CheckReliefValve();
 
    // Serializable member variables (Set in Initialize and in schema)
    double                m_CurrentPeriodTime_s;
    double                m_DriverPressure_cmH2O;
    double                m_PreviousDriverFlow_L_Per_s;
    double                m_PreviousDriverPressure_cmH2O;
    double                m_DriverFlow_L_Per_s;
    double                m_CurrentVentilatorVolume_L;
    double                m_CurrentRespiratoryVolume_L;
    double                m_InspirationTime_s;
    double                m_InspiratoryFlow_L_Per_s;
    double                m_PreviousYPieceToConnectionFlow_L_Per_s;
    double                m_PreviousConnectionPressure_cmH2O;
    bool                  m_PauseOccurred;
    eBreathState          m_CurrentBreathState;
    bool                  m_Initializing;
 
    double                m_PositiveEndExpiratoryPressure_cmH2O;
    double                m_PeakExpiratoryFlow_L_Per_s;
    double                m_PeakInspiratoryFlow_L_Per_s;
    double                m_EndTidalCarbonDioxideFraction;
    double                m_EndTidalCarbonDioxidePressure_cmH2O;
    double                m_EndTidalOxygenFraction;
    double                m_EndTidalOxygenPressure_cmH2O;
 
    SERunningAverage* m_MeanAirwayPressure_cmH2O;
 
    // Stateless member variable (Set in SetUp())
    SEGasCompartment*    m_Environment;
    SEGasCompartment*    m_Ventilator;
    SEGasCompartment*    m_Connection;
    SELiquidCompartment* m_VentilatorAerosol;
    SEFluidCircuitNode*  m_VentilatorNode;
    SEFluidCircuitNode*  m_ExpiratoryLimbNode;
    SEFluidCircuitNode*  m_ExpiratoryValveNode;
    SEFluidCircuitNode*  m_InspiratoryLimbNode;
    SEFluidCircuitNode*  m_InspiratoryValveNode;
    SEFluidCircuitNode*  m_YPieceNode;
    SEFluidCircuitNode*  m_ConnectionNode;
    SEFluidCircuitNode*  m_AmbientNode;
    SEFluidCircuitPath*  m_EnvironmentToVentilator;
    SEFluidCircuitPath*  m_VentilatorToEnvironment;
    SEFluidCircuitPath*  m_YPieceToConnection;
    SEFluidCircuitPath*  m_VentilatorToExpiratoryValve;
    SEFluidCircuitPath*  m_VentilatorToInspiratoryValve;
    SEFluidCircuitPath*  m_ExpiratoryLimbToYPiece;
    SEFluidCircuitPath*  m_InspiratoryLimbToYPiece;
    SEFluidCircuitPath*  m_LeakConnectionToEnvironment;
    SEFluidCircuitPath*  m_ConnectionToReliefValve;
    SEFluidCircuitPath*  m_EnvironmentToReliefValve;
    SEFluidCircuitPath*  m_ConnectionToAirway;
    double               m_MachineClosedResistance_cmH2O_s_Per_L;
    double               m_MachineOpenResistance_cmH2O_s_Per_L;
    // Piecewise linear function points
    std::vector<std::pair<double, double>> m_leakPoints;
  };
END_NAMESPACE