TissueModel.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 SESubstance;
class SEFluidCircuitNode;
class SEFluidCircuitPath;
class SETissueCompartment;
class SEGasCompartment;
class SELiquidCompartment;
class SELiquidSubstanceQuantity;
 
namespace pulse
{
  class PBPhysiology;
  class EngineTest;
  class PULSE_DECL TissueModel : public TissueSystem, public Model
  {
    friend PBPhysiology;//friend the serialization class
    friend EngineTest;
  public:
    TissueModel(Data& data);
    virtual ~TissueModel();
 
    void Clear() override;
 
    // Set members to a stable homeostatic state
    void Initialize() override;
    // Set pointers and other member variables common to both homeostatic initialization and loading a state
    void SetUp() 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;
 
    // Preprocess Methods
    void ProduceAlbumin(double duration_s);
    void Dehydration();
 
    /*Tissue System*/
    void CalculateMetabolicConsumptionAndProduction(double time);
    void GlucoseLipidControl(double time);
 
    /*Process Methods*/
    void CalculateDiffusion();
    void CalculatePulmonaryCapillarySubstanceTransfer();
    void CalculateVitals();
 
    /*Postprocess Methods*/
 
    /*Diffusion Utilities*/
    void DistributeMassbyVolumeWeighted(SELiquidCompartment& cmpt, const SESubstance& sub, double mass, const MassUnit& unit);
    void DistributeMassbyMassWeighted(SELiquidCompartment& cmpt, const SESubstance& sub, double mass, const MassUnit& unit);
 
    double PerfusionLimitedDiffusion(SETissueCompartment& tissue, SELiquidCompartment& vascular, const SESubstance& sub, double partitionCoeff, double timestep_s);
    void AlveolarPartialPressureGradientDiffusion(SEGasCompartment& pulmonary, SELiquidCompartment& vascular, SESubstance& sub, double DiffusingCapacityO2_mL_Per_s_mmHg, double timestep_s);
 
    double MoveMassByInstantDiffusion(SELiquidCompartment& source, SELiquidCompartment& target, const SESubstance& sub);
    double MoveMassBySimpleDiffusion(SELiquidCompartment& source, SELiquidCompartment& target, const SESubstance& sub, double permeabilityCofficient_mL_Per_s, double timestep_s);
    double MoveMassByFacilitatedDiffusion(SELiquidCompartment& source, SELiquidCompartment& target, const SESubstance& sub, double combinedCoefficient_g_Per_s, double timestep_s);
    double MoveMassByActiveTransport(SELiquidCompartment& source, SELiquidCompartment& target, const SESubstance& sub, double DiffusingCapacityO2_mL_Per_s_mmHg, double timestep_s);
 
    // Serializable member variables (Set in Initialize and in schema)
    double m_RestingTissueGlucose_g;
    double m_RestingBloodGlucose_mg_Per_mL;
    double m_RestingBloodLipid_mg_Per_mL;
    double m_RestingBloodInsulin_mg_Per_mL;
    double m_RestingFluidMass_kg;
    double m_PreviousFluidMass_kg;
    // Cache of flows to use to calculate consumption/production during cardiac arrest
    std::map<SELiquidCompartment*, double> m_CardiacArrestVascularFlows_ml_per_min;
 
    // Stateless member variable (Set in SetUp())
    double m_AlbuminProdutionRate_g_Per_s;
    SESubstance* m_Albumin;
    SESubstance* m_Glucose;
    SESubstance* m_Tristearin;
    SESubstance* m_O2;
    SESubstance* m_CO2;
    SESubstance* m_CO;
    SESubstance* m_Lactate;
    SESubstance* m_Acetoacetate;
    SESubstance* m_Creatinine;
    SESubstance* m_Sodium;
    SESubstance* m_Calcium;
    SESubstance* m_Insulin;
    SESubstanceCompound* m_Sweat;
 
    SEFluidCircuitNode* m_GutT1;
    SEFluidCircuitPath* m_GutT1ToGutT3;
 
    SELiquidSubstanceQuantity* m_LiverAcetoacetate;
    SELiquidSubstanceQuantity* m_LiverTissueAlbumin;
    SETissueCompartment* m_LeftLungTissue;
    SETissueCompartment* m_RightLungTissue;
    SETissueCompartment* m_MuscleTissue;
    SELiquidCompartment* m_MuscleIntracellular;
    SETissueCompartment* m_LiverTissue;
    SELiquidCompartment* m_LiverIntracellular;
    SETissueCompartment* m_FatTissue;
    SELiquidCompartment* m_FatIntracellular;
 
    SELiquidCompartment* m_FatVascular;
    SELiquidSubstanceQuantity* m_FatVascularLipid;
    SELiquidSubstanceQuantity* m_LiverVascularGlucose;
    SELiquidCompartment* m_MuscleVascular;
    SELiquidSubstanceQuantity* m_MuscleVascularGlucose;
 
    SEGasCompartment* m_LeftAlveoli;
    SEGasCompartment* m_RightAlveoli;
    SELiquidCompartment* m_LeftPulmonaryCapillaries;
    SELiquidCompartment* m_RightPulmonaryCapillaries;
 
    std::vector<SEFluidCircuitPath*>                     m_VascularResistances;
    std::map<SETissueCompartment*, SELiquidCompartment*> m_TissueToVascular;
    std::vector<SETissueCompartment*>                    m_ConsumptionProdutionTissues;
  };
END_NAMESPACE