Model Mathematics

Component: environment

Component: heart_parameters

Component: driver_function

tau=timeiftimeperiodtime-periodiftimeperiod2time-period2iftimeperiod3time-period3iftimeperiod4time-period4iftimeperiod5time-period5iftimeperiod6time-period6iftimeperiod7time-period7iftimeperiod8time-period8iftimeperiod9time-period9iftimeperiod10time-period10iftimeperiod11time-period11iftimeperiod12time-period12iftimeperiod13 e_t=A-Btau-C2

Component: pericardium

V_pcd=V_lv+V_rv P_pcd=P_0_pcdlambda_pcdV_pcd-V_0_pcd-1 P_peri=P_pcd+P_th

Component: left_ventricle

V_lvf=V_lv-V_spt P_lvf=e_tP_es_lvf+1-e_tP_ed_lvf P_lv=P_lvf+P_peri ddtimeV_lv=0ifQ_mt<0Q_av<0-Q_avifQ_mt<0Q_mtifQ_av<0Q_mt-Q_avotherwise

Component: lvf_calculator

P_es_lvf=E_es_lvfV_lvf-V_d_lvf P_ed_lvf=P_0_lvflambda_lvfV_lvf-V_0_lvf-1

Component: right_ventricle

V_rvf=V_rv+V_spt P_rvf=e_tP_es_rvf+1-e_tP_ed_rvf P_rv=P_rvf+P_peri ddtimeV_rv=0ifQ_tc<0Q_pv<0-Q_pvifQ_tc<0Q_tcifQ_pv<0Q_tc-Q_pvotherwise

Component: rvf_calculator

P_es_rvf=E_es_rvfV_rvf-V_d_rvf P_ed_rvf=P_0_rvflambda_rvfV_rvf-V_0_rvf-1

Component: septum

0=e_tE_es_sptV_spt-V_d_spt+one-e_tP_0_sptlambda_sptV_spt-V_0_spt-one-e_tE_es_lvfV_lv-V_spt-1-e_tP_0_lvflambda_lvfV_lv-V_spt-1+e_tE_es_rvfV_rv+V_spt+1-e_tP_0_rvflambda_rvfV_rv+V_spt-1

Component: pulmonary_artery

P_pa=E_es_paV_pa-V_d_pa+P_th ddtimeV_pa=-Q_pulifQ_pv<0Q_pv-Q_pulotherwise

Component: pulmonary_vein

P_pu=E_es_puV_pu-V_d_pu+P_th ddtimeV_pu=Q_pulifQ_mt<0Q_pul-Q_mtotherwise

Component: aorta

P_ao=E_es_aoV_ao-V_d_ao ddtimeV_ao=-Q_sysifQ_av<0Q_av-Q_sysotherwise

Component: vena_cava

P_vc=E_es_vcV_vc-V_d_vc ddtimeV_vc=Q_sysifQ_tc<0Q_sys-Q_tcotherwise

Component: flow

Q_sys=P_ao-P_vcR_sys Q_pul=P_pa-P_puR_pul ddtimeQ_mt=0ifP_pu-P_lv<0Q_mt<0P_pu-P_lv-Q_mtR_mtL_mtotherwise ddtimeQ_av=0ifP_lv-P_ao<0Q_av<0P_lv-P_ao-Q_avR_avL_avotherwise ddtimeQ_tc=0ifP_vc-P_rv<0Q_tc<0P_vc-P_rv-Q_tcR_tcL_tcotherwise ddtimeQ_pv=0ifP_rv-P_pa<0Q_pv<0P_rv-P_pa-Q_pvR_pvL_pvotherwise