-
Notifications
You must be signed in to change notification settings - Fork 1
/
main_ADAN-86-Heart_cellml.txt
5456 lines (4627 loc) · 167 KB
/
main_ADAN-86-Heart_cellml.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
def model CardiovascularSystem as
def import using "Units.cellml" for
unit m3 using unit m3;
unit m3_per_s using unit m3_per_s;
unit m3_per_s2 using unit m3_per_s2;
unit m6_per_J using unit m6_per_J;
unit per_m using unit per_m;
unit per_m3 using unit per_m3;
unit J_per_m3 using unit J_per_m3;
unit J_s_per_m6 using unit J_s_per_m6;
unit J_s2_per_m6 using unit J_s2_per_m6;
unit J_s2_per_m5 using unit J_s2_per_m5;
unit J_s2_per_m9 using unit J_s2_per_m9;
unit J_per_s using unit J_per_s;
unit J_per_m6 using unit J_per_m6;
unit J_s_per_m3 using unit J_s_per_m3;
unit UnitValve using unit UnitValve;
unit s2 using unit s2;
unit s_per_m3 using unit s_per_m3;
enddef;
def import using "Parameters.cellml" for
comp Parameters_Systemic using comp Parameters_Systemic;
comp Parameters_Pulmonary using comp Parameters_Pulmonary;
// comp Parameters_Coronary using comp Parameters_Coronary;
comp Parameters_Heart using comp Parameters_Heart;
enddef;
def comp environment as
var time: second {pub: out};
enddef;
//------------ Systemic Circulation ------------//
def comp Systemic as
var t: second {pub: in, priv: out};
// Constitutive Parameters
var C_svl: m6_per_J {priv: in};
var C_svn: m6_per_J {priv: in};
var C_svc: m6_per_J {priv: in};
var C_ivl: m6_per_J {priv: in};
var C_ivn: m6_per_J {priv: in};
var C_ivc: m6_per_J {priv: in};
var R_svl: J_s_per_m6 {priv: in};
var R_svn: J_s_per_m6 {priv: in};
var R_svc: J_s_per_m6 {priv: in};
var R_ivl: J_s_per_m6 {priv: in};
var R_ivn: J_s_per_m6 {priv: in};
var R_ivc: J_s_per_m6 {priv: in};
var I_svl: J_s2_per_m6 {priv: in};
var I_svn: J_s2_per_m6 {priv: in};
var I_svc: J_s2_per_m6 {priv: in};
var I_ivl: J_s2_per_m6 {priv: in};
var I_ivn: J_s2_per_m6 {priv: in};
var I_ivc: J_s2_per_m6 {priv: in};
// State Variables
var u_sas: J_per_m3 {pub: out};
var u_root: J_per_m3 {pub: out, priv: in};
var v_sas: m3_per_s {priv: out};
// var v_coronary: m3_per_s {pub: in};
//
var v_aov: m3_per_s {pub: in};
var u_ra: J_per_m3 {pub: in};
var u_svl: J_per_m3 {init: 100, priv: out};
var u_svn: J_per_m3 {init: 0.0, priv: out};
var u_svc: J_per_m3 {init: 0.0, priv: out};
var u_ivl: J_per_m3 {init: 100, priv: out};
var u_ivn: J_per_m3 {init: 0.0, priv: out};
var u_ivc: J_per_m3 {init: 0.0, priv: out};
var v_svl: m3_per_s {init: 0.0, priv: out};
var v_svn: m3_per_s {init: 0.0, priv: out};
var v_svc: m3_per_s {init: 0.0, priv: out};
var v_ivl: m3_per_s {init: 0.0, priv: out};
var v_ivn: m3_per_s {init: 0.0, priv: out};
var v_ivc: m3_per_s {init: 0.0, priv: out};
var v_sup_venacava: m3_per_s {pub: out};
var v_inf_venacava: m3_per_s {pub: out};
var v_chest: m3_per_s {pub: out};
var v_trunk: m3_per_s {pub: out};
var v_kidneys: m3_per_s {pub: out};
var v_intestines: m3_per_s {pub: out};
var v_legs: m3_per_s {pub: out};
var v_arms: m3_per_s {pub: out};
var v_head: m3_per_s {pub: out};
var v_posterior_intercostal_T1_R98: m3_per_s {priv: in};
var v_posterior_intercostal_T1_L102: m3_per_s {priv: in};
var v_posterior_intercostal_T2_R106: m3_per_s {priv: in};
var v_posterior_intercostal_T2_L110: m3_per_s {priv: in};
var v_left_gastric_T3_C120: m3_per_s {priv: in};
var v_dorsal_pancreatic_T1_C124: m3_per_s {priv: in};
var v_splenic_T2_C126: m3_per_s {priv: in};
var v_hepatic_artery_proper_left_branch_C132: m3_per_s {priv: in};
var v_hepatic_artery_proper_right_branch_C134: m3_per_s {priv: in};
var v_middle_colic_T8_C140: m3_per_s {priv: in};
var v_jejunal_3_T10_C144: m3_per_s {priv: in};
var v_jejunal_6_T11_C148: m3_per_s {priv: in};
var v_ileocolic_T9_C152: m3_per_s {priv: in};
var v_ileal_4_T12_C156: m3_per_s {priv: in};
var v_ileal_6_T13_C160: m3_per_s {priv: in};
var v_superior_mesenteric_T4_C162: m3_per_s {priv: in};
var v_inferior_segmental_T5_L170: m3_per_s {priv: in};
var v_superior_segmental_T4_L172: m3_per_s {priv: in};
var v_renal_posterior_branch_T3_L174: m3_per_s {priv: in};
var v_superior_segmental_T4_R182: m3_per_s {priv: in};
var v_inferior_segmental_T5_R184: m3_per_s {priv: in};
var v_renal_posterior_branch_T3_R186: m3_per_s {priv: in};
var v_inferior_mesenteric_T5_C190: m3_per_s {priv: in};
var v_internal_iliac_T1_R218: m3_per_s {priv: in};
var v_profundus_T2_R224: m3_per_s {priv: in};
var v_anterior_tibial_T3_R230: m3_per_s {priv: in};
var v_posterior_tibial_T4_R236: m3_per_s {priv: in};
var v_internal_iliac_T1_L196: m3_per_s {priv: in};
var v_profundus_T2_L202: m3_per_s {priv: in};
var v_anterior_tibial_T3_L208: m3_per_s {priv: in};
var v_posterior_tibial_T4_L214: m3_per_s {priv: in};
var v_posterior_interosseous_T3_R40: m3_per_s {priv: in};
var v_ulnar_T2_R42: m3_per_s {priv: in};
var v_radial_T1_R44: m3_per_s {priv: in};
var v_posterior_interosseous_T3_L88: m3_per_s {priv: in};
var v_ulnar_T2_L90: m3_per_s {priv: in};
var v_radial_T1_L92: m3_per_s {priv: in};
var v_external_carotid_T2_R26: m3_per_s {priv: in};
var v_external_carotid_T2_L62: m3_per_s {priv: in};
var v_internal_carotid_R8_C: m3_per_s {priv: in};
var v_internal_carotid_L50_C: m3_per_s {priv: in};
var v_vertebral_R272: m3_per_s {priv: in};
var v_vertebral_L2: m3_per_s {priv: in};
// Constitutive Relations
u_sas = u_root;
v_sas = v_aov;
// v_sas = v_aov+v_coronary;
// ode(u_svl, t) = (v_sup_venacava-v_svl)/C_svl;
ode(u_svn, t) = (v_svl-v_svn)/C_svn;
ode(u_svc, t) = (v_svn-v_svc)/C_svc;
// ode(u_ivl, t) = (v_inf_venacava-v_ivl)/C_ivl;
ode(u_ivn, t) = (v_ivl-v_ivn)/C_ivn;
ode(u_ivc, t) = (v_ivn-v_ivc)/C_ivc;
ode(v_svl, t) = (u_svl-u_svn-v_svl*R_svl)/I_svl;
ode(v_svn, t) = (u_svn-u_svc-v_svn*R_svn)/I_svn;
ode(v_svc, t) = (u_svc-u_ra-v_svc*R_svc)/I_svc;
ode(v_ivl, t) = (u_ivl-u_ivn-v_ivl*R_ivl)/I_ivl;
ode(v_ivn, t) = (u_ivn-u_ivc-v_ivn*R_ivn)/I_ivn;
ode(v_ivc, t) = (u_ivc-u_ra-v_ivc*R_ivc)/I_ivc;
// Conservation Laws
v_chest = v_posterior_intercostal_T1_R98+v_posterior_intercostal_T1_L102+v_posterior_intercostal_T2_R106+v_posterior_intercostal_T2_L110;
v_trunk = v_left_gastric_T3_C120+v_dorsal_pancreatic_T1_C124+v_splenic_T2_C126+v_hepatic_artery_proper_left_branch_C132+v_hepatic_artery_proper_right_branch_C134;
v_intestines = v_middle_colic_T8_C140+v_jejunal_3_T10_C144+v_jejunal_6_T11_C148+v_ileocolic_T9_C152+v_ileal_4_T12_C156+v_ileal_6_T13_C160+v_superior_mesenteric_T4_C162+v_inferior_mesenteric_T5_C190;
v_kidneys = v_inferior_segmental_T5_L170+v_superior_segmental_T4_L172+v_renal_posterior_branch_T3_L174+v_superior_segmental_T4_R182+v_inferior_segmental_T5_R184+v_renal_posterior_branch_T3_R186;
v_legs = v_internal_iliac_T1_R218+v_profundus_T2_R224+v_anterior_tibial_T3_R230+v_posterior_tibial_T4_R236+v_internal_iliac_T1_L196+v_profundus_T2_L202+v_anterior_tibial_T3_L208+v_posterior_tibial_T4_L214;
v_arms = v_posterior_interosseous_T3_R40+v_ulnar_T2_R42+v_radial_T1_R44+v_posterior_interosseous_T3_L88+v_ulnar_T2_L90+v_radial_T1_L92;
v_head = v_external_carotid_T2_R26+v_external_carotid_T2_L62+v_internal_carotid_R8_C+v_internal_carotid_L50_C+v_vertebral_R272+v_vertebral_L2;
v_sup_venacava = v_arms+v_head;
v_inf_venacava = v_chest+v_trunk+v_intestines+v_kidneys+v_legs;
enddef;
//------------ Pulmonary Circulation ------------//
def comp Pulmonary as
var t: second {pub: in};
// Constitutive Parameters
var C_pas: m6_per_J {priv: in};
var C_pat: m6_per_J {priv: in};
var C_par: m6_per_J {priv: in};
var C_pcp: m6_per_J {priv: in};
var C_pvn: m6_per_J {priv: in};
var C_pvc: m6_per_J {priv: in};
var R_pas: J_s_per_m6 {priv: in};
var R_pat: J_s_per_m6 {priv: in};
var R_par: J_s_per_m6 {priv: in};
var R_pcp: J_s_per_m6 {priv: in};
var R_psh: J_s_per_m6 {priv: in};
var R_pvn: J_s_per_m6 {priv: in};
var I_pas: J_s2_per_m6 {priv: in};
var I_pat: J_s2_per_m6 {priv: in};
var I_par: J_s2_per_m6 {priv: in};
var I_pcp: J_s2_per_m6 {priv: in};
var I_pvn: J_s2_per_m6 {priv: in};
var I_psh: J_s2_per_m6 {priv: in};
// State Variables
var u_pas: J_per_m3 {init: 4000.0, pub: out};
var u_pat: J_per_m3 {init: 0.0, pub: out};
var u_par: J_per_m3 {init: 0.0, pub: out};
var u_pcp: J_per_m3 {init: 0.0, pub: out};
var u_pvn: J_per_m3 {init: 0.0, pub: out};
var u_la: J_per_m3 {pub: in};
var v_pas: m3_per_s {init: 0.0, pub: out};
var v_pat: m3_per_s {init: 0.0, pub: out};
var v_par: m3_per_s {init: 0.0, pub: out};
var v_pcp: m3_per_s {init: 0.0, pub: out};
var v_psh: m3_per_s {init: 0.0, pub: out};
var v_pvn: m3_per_s {init: 0.0, pub: out};
var v_puv: m3_per_s {pub: in};
// Constitutive Relations
ode(u_pas, t) = (v_puv-v_pas)/C_pas;
ode(u_pat, t) = (v_pas-v_pat)/C_pat;
ode(u_par, t) = (v_pat-v_psh-v_par)/C_par;
ode(u_pcp, t) = (v_par-v_pcp)/C_pcp;
ode(u_pvn, t) = (v_pcp+v_psh-v_pvn)/C_pvn;
ode(v_pas, t) = (u_pas-u_pat-v_pas*R_pas)/I_pas;
ode(v_pat, t) = (u_pat-u_par-v_pat*R_pat)/I_pat;
ode(v_par, t) = (u_par-u_pcp-v_par*R_par)/I_par;
ode(v_pcp, t) = (u_pcp-u_pvn-v_pcp*R_pcp)/I_pcp;
ode(v_pvn, t) = (u_pvn-u_la-v_pvn*R_pvn)/I_pvn;
ode(v_psh, t) = (u_par-u_pvn-v_psh*R_psh)/I_psh;
enddef;
//------------ Coronary ------------//
// def comp Coronary as
// var t: second {pub: in};
// // Constitutive Parameters
// var C_corepi: m6_per_J {priv: in};
// var C_corint: m6_per_J {priv: in};
// var C_corcap: m6_per_J {priv: in};
// var C_corvns: m6_per_J {priv: in};
// var R_corepi: J_s_per_m6 {priv: in};
// var R_corint: J_s_per_m6 {priv: in};
// var R_corcap: J_s_per_m6 {priv: in};
// var R_corvns: J_s_per_m6 {priv: in};
// var I_corepi: J_s2_per_m6 {priv: in};
// var I_corint: J_s2_per_m6 {priv: in};
// var I_corcap: J_s2_per_m6 {priv: in};
// var I_corvns: J_s2_per_m6 {priv: in};
// // State Variables
// var u_sas: J_per_m3 {pub: in};
// var u_ra: J_per_m3 {pub: in};
// var u_corepi: J_per_m3 {init: 0.0, pub: out};
// var u_corint: J_per_m3 {init: 0.0, pub: out};
// var u_corcap: J_per_m3 {init: 0.0, pub: out};
// var u_corvns: J_per_m3 {init: 0.0, pub: out};
// var v_corepi: m3_per_s {init: 0.0, pub: out};
// var v_corint: m3_per_s {init: 0.0, pub: out};
// var v_corcap: m3_per_s {init: 0.0, pub: out};
// var v_corvns_in: m3_per_s {init: 0.0, pub: out};
// var v_corvns_out: m3_per_s {init: 0.0, pub: out};
// // Constitutive Relations
// ode(v_corepi, t) = (u_sas-u_corepi-v_corepi*R_corepi)/I_corepi;
// ode(v_corint, t) = (u_corepi-u_corint-v_corint*R_corint)/I_corint;
// ode(v_corcap, t) = (u_corint-u_corcap-v_corcap*R_corcap)/I_corcap;
// ode(v_corvns_in, t) = (u_corcap-u_corvns-v_corvns_in*R_corvns/2{dimensionless})/(I_corvns/2{dimensionless});
// ode(v_corvns_out, t) = (u_corvns-u_ra-v_corvns_out*R_corvns/2{dimensionless})/(I_corvns/2{dimensionless});
// ode(u_corepi, t) = (v_corepi-v_corint)/C_corepi;
// ode(u_corint, t) = (v_corint-v_corcap)/C_corint;
// ode(u_corcap, t) = (v_corcap-v_corvns_in)/C_corcap;
// ode(u_corvns, t) = (v_corvns_in-v_corvns_out)/C_corvns;
// enddef;
//------------ Heart ------------//
// def comp Heart as
// var t: second {pub: in};
// // Constitutive Parameters
// var T: second {priv: in};
// var T_sys_0: second {priv: in};
// var k_sys: s2 {priv: in};
// var T_sys: second;
// var phi: dimensionless;
// var uu: dimensionless;
// var xi: dimensionless {init: 0.0};
// var C_la: m6_per_J {priv: in};
// var C_ra: m6_per_J {priv: in};
// var R_la: J_s_per_m6 {priv: in};
// var R_ra: J_s_per_m6 {priv: in};
// var R_lv: J_s_per_m6;
// var R_rv: J_s_per_m6;
// var E_max_lv: J_per_m6 {priv: in};
// var E_max_rv: J_per_m6 {priv: in};
// var k_E_lv: per_m3 {priv: in};
// var k_E_rv: per_m3 {priv: in};
// var k_R_lv: s_per_m3 {priv: in};
// var k_R_rv: s_per_m3 {priv: in};
// var u_0_lv: J_per_m3 {priv: in};
// var u_0_rv: J_per_m3 {priv: in};
// var q_0_lv: m3 {priv: in};
// var q_0_rv: m3 {priv: in};
// // State Variables
// var u_sas: J_per_m3 {pub: in};
// var u_la: J_per_m3 {init: 0.0, pub: out};
// var u_ra: J_per_m3 {init: 0.0, pub: out};
// var u_par: J_per_m3 {pub: in};
// var u_lv: J_per_m3 {pub: out};
// var u_rv: J_per_m3 {pub: out};
// var u_max_lv: J_per_m3;
// var u_max_rv: J_per_m3;
// var q_rv: m3 {init: 172.0e-6};
// var q_lv: m3 {init: 136.0e-6};
// var q_ra: m3 {init: 67.0e-6};
// var q_la: m3 {init: 71.0e-6};
// var v_miv: m3_per_s {pub: out};
// var v_aov: m3_per_s {pub: out};
// var v_trv: m3_per_s {pub: out};
// var v_puv: m3_per_s {pub: out};
// var v_sup_venacava: m3_per_s {pub: in};
// var v_inf_venacava: m3_per_s {pub: in};
// // var v_coronary_sinus: m3_per_s {pub: in};
// var v_pvn: m3_per_s {pub: in};
// // Blood Flow
// v_miv = sel
// case u_la > u_lv:
// (u_la-u_lv)/R_la;
// case u_la <= u_lv:
// 0{m3_per_s};
// endsel;
// v_aov = sel
// case u_max_lv > u_sas:
// (u_max_lv-u_sas)/R_lv;
// case u_max_lv <= u_sas:
// 0{m3_per_s};
// endsel;
// v_trv = sel
// case u_ra > u_rv:
// (u_ra-u_rv)/R_ra;
// case u_ra <= u_rv:
// 0{m3_per_s};
// endsel;
// v_puv = sel
// case u_max_rv >= u_par:
// (u_max_rv-u_par)/R_rv;
// case u_max_rv < u_par:
// 0{m3_per_s};
// endsel;
// // Pressure of the Atria
// ode(u_la, t) = (v_pvn-v_miv)/C_la;
// ode(u_ra, t) = (v_sup_venacava+v_inf_venacava-v_trv)/C_ra;
// // Blood Volume of the Ventricles
// ode(q_lv, t) = v_miv-v_aov;
// ode(q_rv, t) = v_trv-v_puv;
// ode(q_la, t) = v_pvn-v_miv;
// ode(q_ra, t) = v_sup_venacava+v_inf_venacava-v_trv;
// // Isometric pressure of the Ventricles
// u_max_lv = phi*E_max_lv*(q_lv-q_0_lv)+(1{dimensionless}-phi)*u_0_lv*(exp(k_E_lv*q_lv)-1{dimensionless});
// u_max_rv = phi*E_max_rv*(q_rv-q_0_rv)+(1{dimensionless}-phi)*u_0_rv*(exp(k_E_rv*q_rv)-1{dimensionless});
// // Viscous Resistance of the Ventricles
// R_lv = k_R_lv*u_max_lv;
// R_rv = k_R_rv*u_max_rv;
// // Instantaneous Pressure of the Ventricles
// u_lv = u_max_lv-R_lv*v_aov;
// u_rv = u_max_rv-R_rv*v_puv;
// // Ventricle Activation Function
// phi = sel
// case (uu >= 0{dimensionless}) and (uu <= T_sys/T):
// sqr(sin(pi*T*uu/T_sys));
// case (uu > T_sys/T) and (uu <= 1{dimensionless}):
// 0{dimensionless};
// endsel;
// // Fraction of Cardiac Cycle
// uu = xi-floor(xi);
// ode(xi, t) = 1{dimensionless}/T;
// // Systole Duration
// T_sys = T_sys_0-k_sys/T;
// enddef;
def comp Heart as
var t: second {pub: in};
// Constitutive Parameters
var mt: second;
var e_a: dimensionless;
var e_v: dimensionless;
var T: second {priv: in};
var t_ac: dimensionless {priv: in};
var t_ar: dimensionless {priv: in};
var T_ac: dimensionless {priv: in};
var T_ar: dimensionless {priv: in};
var T_vc: dimensionless {priv: in};
var T_vr: dimensionless {priv: in};
var R_trv: J_s_per_m6 {priv: in};
var R_puv: J_s_per_m6 {priv: in};
var R_miv: J_s_per_m6 {priv: in};
var R_aov: J_s_per_m6 {priv: in};
var B_trv: J_s2_per_m9 {priv: in};
var B_puv: J_s2_per_m9 {priv: in};
var B_miv: J_s2_per_m9 {priv: in};
var B_aov: J_s2_per_m9 {priv: in};
var L_trv: J_s2_per_m6 {priv: in};
var L_puv: J_s2_per_m6 {priv: in};
var L_miv: J_s2_per_m6 {priv: in};
var L_aov: J_s2_per_m6 {priv: in};
var E_lv_A: J_per_m6 {priv: in};
var E_lv_B: J_per_m6 {priv: in};
var E_la_A: J_per_m6 {priv: in};
var E_la_B: J_per_m6 {priv: in};
var E_rv_A: J_per_m6 {priv: in};
var E_rv_B: J_per_m6 {priv: in};
var E_ra_A: J_per_m6 {priv: in};
var E_ra_B: J_per_m6 {priv: in};
var q_ra_0: m3 {priv: in};
var q_rv_0: m3 {priv: in};
var q_la_0: m3 {priv: in};
var q_lv_0: m3 {priv: in};
// State Variables
var u_ra: J_per_m3 {pub: out};
var u_rv: J_per_m3;
var u_la: J_per_m3 {pub: out};
var u_lv: J_per_m3;
var u_sas: J_per_m3 {pub: in};
var u_par: J_per_m3 {pub: in};
var v_trv: m3_per_s {init: 0};
var v_puv: m3_per_s {init: 0, pub: out};
var v_miv: m3_per_s {init: 0};
var v_aov: m3_per_s {init: 0, pub: out};
var v_sup_venacava: m3_per_s {pub: in};
var v_inf_venacava: m3_per_s {pub: in};
var v_pvn: m3_per_s {pub: in};
var q_ra: m3 {init: 20.0e-6};
var q_rv: m3 {init: 500.0e-6};
var q_la: m3 {init: 20.0e-6};
var q_lv: m3 {init: 500.0e-6};
// Elastance Functions
mt = t-T*floor(t/T);
e_a = sel
case (mt >= 0{second}) and (mt <= (t_ar+T_ar)*T-T):
0.5{dimensionless}*(1{dimensionless}+cos(pi*(mt+T-t_ar*T)/(T_ar*T)));
case (mt > (t_ar+T_ar)*T-T) and (mt <= t_ac*T):
0{dimensionless};
case (mt > t_ac*T) and (mt <= (t_ac+T_ac)*T):
0.5{dimensionless}*(1{dimensionless}-cos(pi*(mt-t_ac*T)/(T_ac*T)));
case (mt > (t_ac+T_ac)*T) and (mt <= T):
0.5{dimensionless}*(1{dimensionless}+cos(pi*(mt-t_ar*T)/(T_ar*T)));
endsel;
e_v = sel
case (mt >= 0{second}) and (mt <= T_vc*T):
0.5{dimensionless}*(1{dimensionless}-cos(pi*mt/(T_vc*T)));
case (mt > T_vc*T) and (mt <= (T_vc+T_vr)*T):
0.5{dimensionless}*(1{dimensionless}+cos(pi*(mt-T_vc*T)/(T_vr*T)));
case (mt > (T_vc+T_vr)*T) and (mt < T):
0{dimensionless};
endsel;
// Constitutive Relations
ode(v_trv, t) = sel
case u_ra >= u_rv:
(u_ra-u_rv-(R_trv+B_trv*abs(v_trv))*v_trv)/L_trv;
case u_ra < u_rv:
-(R_trv+B_trv*abs(v_trv))*v_trv/L_trv;
endsel;
ode(v_puv, t) = sel
case u_rv >= u_par:
(u_rv-u_par-(R_puv+B_puv*abs(v_puv))*v_puv)/L_puv;
case u_rv < u_par:
-(R_puv+B_puv*abs(v_puv))*v_puv/L_puv;
endsel;
ode(v_miv, t) = sel
case u_la >= u_lv:
(u_la-u_lv-(R_miv+B_miv*abs(v_miv))*v_miv)/L_miv;
case u_la < u_lv:
-(R_miv+B_miv*abs(v_miv))*v_miv/L_miv;
endsel;
ode(v_aov, t) = sel
case u_lv >= u_sas:
(u_lv-u_sas-(R_aov+B_aov*abs(v_aov))*v_aov)/L_aov;
case u_lv < u_sas:
-(R_aov+B_aov*abs(v_aov))*v_aov/L_aov;
endsel;
u_ra = (e_a*E_ra_A+E_ra_B)*(q_ra-q_ra_0);
u_rv = (e_v*E_rv_A+E_rv_B)*(q_rv-q_rv_0);
u_la = (e_a*E_la_A+E_la_B)*(q_la-q_la_0);
u_lv = (e_v*E_lv_A+E_lv_B)*(q_lv-q_lv_0);
// Conservation Laws
ode(q_ra, t) = v_sup_venacava+v_inf_venacava-v_trv;
ode(q_rv, t) = v_trv-v_puv;
ode(q_la, t) = v_pvn-v_miv;
ode(q_lv, t) = v_miv-v_aov;
enddef;
//------------ Mapping ------------//
def map between Heart and Systemic for
vars v_sup_venacava and v_sup_venacava;
vars v_inf_venacava and v_inf_venacava;
vars u_ra and u_ra;
vars u_sas and u_sas;
vars v_aov and v_aov;
enddef;
def map between Heart and Pulmonary for
vars u_par and u_pas;
vars v_pvn and v_pvn;
vars u_la and u_la;
vars v_puv and v_puv;
enddef;
// def map between Heart and Coronary for
// vars v_coronary_sinus and v_corvns_out;
// vars u_ra and u_ra;
// enddef;
// def map between Systemic and Coronary for
// vars v_coronary and v_corepi;
// vars u_sas and u_sas;
// enddef;
// def map between Parameters_Coronary and Coronary for
// vars C_corepi and C_corepi;
// vars C_corint and C_corint;
// vars C_corcap and C_corcap;
// vars C_corvns and C_corvns;
// vars R_corepi and R_corepi;
// vars R_corint and R_corint;
// vars R_corcap and R_corcap;
// vars R_corvns and R_corvns;
// vars I_corepi and I_corepi;
// vars I_corint and I_corint;
// vars I_corcap and I_corcap;
// vars I_corvns and I_corvns;
// enddef;
// def map between environment and Coronary for
// vars time and t;
// enddef;
def map between Parameters_Systemic and Systemic for
vars C_svl and C_svl;
vars C_svn and C_svn;
vars C_svc and C_svc;
vars C_ivl and C_ivl;
vars C_ivn and C_ivn;
vars C_ivc and C_ivc;
vars R_svl and R_svl;
vars R_svn and R_svn;
vars R_svc and R_svc;
vars R_ivl and R_ivl;
vars R_ivn and R_ivn;
vars R_ivc and R_ivc;
vars I_svl and I_svl;
vars I_svn and I_svn;
vars I_svc and I_svc;
vars I_ivl and I_ivl;
vars I_ivn and I_ivn;
vars I_ivc and I_ivc;
enddef;
def map between Parameters_Pulmonary and Pulmonary for
vars C_pas and C_pas;
vars C_pat and C_pat;
vars C_par and C_par;
vars C_pcp and C_pcp;
vars C_pvn and C_pvn;
vars C_pvc and C_pvc;
vars R_pas and R_pas;
vars R_pat and R_pat;
vars R_par and R_par;
vars R_pcp and R_pcp;
vars R_psh and R_psh;
vars R_pvn and R_pvn;
vars I_pas and I_pas;
vars I_pat and I_pat;
vars I_par and I_par;
vars I_pcp and I_pcp;
vars I_pvn and I_pvn;
vars I_psh and I_psh;
enddef;
// def map between Parameters_Heart and Heart for
// vars T and T;
// vars T_sys_0 and T_sys_0;
// vars k_sys and k_sys;
// vars C_la and C_la;
// vars C_ra and C_ra;
// vars R_la and R_la;
// vars R_ra and R_ra;
// vars E_max_lv and E_max_lv;
// vars E_max_rv and E_max_rv;
// vars k_E_lv and k_E_lv;
// vars k_E_rv and k_E_rv;
// vars k_R_lv and k_R_lv;
// vars k_R_rv and k_R_rv;
// vars u_0_lv and u_0_lv;
// vars u_0_rv and u_0_rv;
// vars q_0_rv and q_0_rv;
// vars q_0_lv and q_0_lv;
// enddef;
def map between Parameters_Heart and Heart for
vars T and T;
vars t_ac and t_ac;
vars t_ar and t_ar;
vars T_ac and T_ac;
vars T_ar and T_ar;
vars T_vc and T_vc;
vars T_vr and T_vr;
vars R_trv and R_trv;
vars R_puv and R_puv;
vars R_miv and R_miv;
vars R_aov and R_aov;
vars B_trv and B_trv;
vars B_puv and B_puv;
vars B_miv and B_miv;
vars B_aov and B_aov;
vars L_trv and L_trv;
vars L_puv and L_puv;
vars L_miv and L_miv;
vars L_aov and L_aov;
vars E_lv_A and E_lv_A;
vars E_lv_B and E_lv_B;
vars E_la_A and E_la_A;
vars E_la_B and E_la_B;
vars E_rv_A and E_rv_A;
vars E_rv_B and E_rv_B;
vars E_ra_A and E_ra_A;
vars E_ra_B and E_ra_B;
vars q_ra_0 and q_ra_0;
vars q_rv_0 and q_rv_0;
vars q_la_0 and q_la_0;
vars q_lv_0 and q_lv_0;
enddef;
def map between environment and Heart for
vars time and t;
enddef;
def map between environment and Systemic for
vars time and t;
enddef;
def map between environment and Pulmonary for
vars time and t;
enddef;
// Terminal
def map between Systemic and posterior_intercostal_T1_R98_module for
vars u_ivl and u_out;
vars v_posterior_intercostal_T1_R98 and v_T;
vars t and t;
enddef;
def map between Systemic and posterior_intercostal_T1_L102_module for
vars u_ivl and u_out;
vars v_posterior_intercostal_T1_L102 and v_T;
vars t and t;
enddef;
def map between Systemic and posterior_intercostal_T2_R106_module for
vars u_ivl and u_out;
vars v_posterior_intercostal_T2_R106 and v_T;
vars t and t;
enddef;
def map between Systemic and posterior_intercostal_T2_L110_module for
vars u_ivl and u_out;
vars v_posterior_intercostal_T2_L110 and v_T;
vars t and t;
enddef;
def map between Systemic and left_gastric_T3_C120_module for
vars u_ivl and u_out;
vars v_left_gastric_T3_C120 and v_T;
vars t and t;
enddef;
def map between Systemic and dorsal_pancreatic_T1_C124_module for
vars u_ivl and u_out;
vars v_dorsal_pancreatic_T1_C124 and v_T;
vars t and t;
enddef;
def map between Systemic and splenic_T2_C126_module for
vars u_ivl and u_out;
vars v_splenic_T2_C126 and v_T;
vars t and t;
enddef;
def map between Systemic and hepatic_artery_proper_left_branch_C132_module for
vars u_ivl and u_out;
vars v_hepatic_artery_proper_left_branch_C132 and v_T;
vars t and t;
enddef;
def map between Systemic and hepatic_artery_proper_right_branch_C134_module for
vars u_ivl and u_out;
vars v_hepatic_artery_proper_right_branch_C134 and v_T;
vars t and t;
enddef;
def map between Systemic and middle_colic_T8_C140_module for
vars u_ivl and u_out;
vars v_middle_colic_T8_C140 and v_T;
vars t and t;
enddef;
def map between Systemic and jejunal_3_T10_C144_module for
vars u_ivl and u_out;
vars v_jejunal_3_T10_C144 and v_T;
vars t and t;
enddef;
def map between Systemic and jejunal_6_T11_C148_module for
vars u_ivl and u_out;
vars v_jejunal_6_T11_C148 and v_T;
vars t and t;
enddef;
def map between Systemic and ileocolic_T9_C152_module for
vars u_ivl and u_out;
vars v_ileocolic_T9_C152 and v_T;
vars t and t;
enddef;
def map between Systemic and ileal_4_T12_C156_module for
vars u_ivl and u_out;
vars v_ileal_4_T12_C156 and v_T;
vars t and t;
enddef;
def map between Systemic and ileal_6_T13_C160_module for
vars u_ivl and u_out;
vars v_ileal_6_T13_C160 and v_T;
vars t and t;
enddef;
def map between Systemic and superior_mesenteric_T4_C162_module for
vars u_ivl and u_out;
vars v_superior_mesenteric_T4_C162 and v_T;
vars t and t;
enddef;
def map between Systemic and inferior_segmental_T5_L170_module for
vars u_ivl and u_out;
vars v_inferior_segmental_T5_L170 and v_T;
vars t and t;
enddef;
def map between Systemic and superior_segmental_T4_L172_module for
vars u_ivl and u_out;
vars v_superior_segmental_T4_L172 and v_T;
vars t and t;
enddef;
def map between Systemic and renal_posterior_branch_T3_L174_module for
vars u_ivl and u_out;
vars v_renal_posterior_branch_T3_L174 and v_T;
vars t and t;
enddef;
def map between Systemic and superior_segmental_T4_R182_module for
vars u_ivl and u_out;
vars v_superior_segmental_T4_R182 and v_T;
vars t and t;
enddef;
def map between Systemic and inferior_segmental_T5_R184_module for
vars u_ivl and u_out;
vars v_inferior_segmental_T5_R184 and v_T;
vars t and t;
enddef;
def map between Systemic and renal_posterior_branch_T3_R186_module for
vars u_ivl and u_out;
vars v_renal_posterior_branch_T3_R186 and v_T;
vars t and t;
enddef;
def map between Systemic and inferior_mesenteric_T5_C190_module for
vars u_ivl and u_out;
vars v_inferior_mesenteric_T5_C190 and v_T;
vars t and t;
enddef;
def map between Systemic and internal_iliac_T1_R218_module for
vars u_ivl and u_out;
vars v_internal_iliac_T1_R218 and v_T;
vars t and t;
enddef;
def map between Systemic and profundus_T2_R224_module for
vars u_ivl and u_out;
vars v_profundus_T2_R224 and v_T;
vars t and t;
enddef;
def map between Systemic and anterior_tibial_T3_R230_module for
vars u_ivl and u_out;
vars v_anterior_tibial_T3_R230 and v_T;
vars t and t;
enddef;
def map between Systemic and posterior_tibial_T4_R236_module for
vars u_ivl and u_out;
vars v_posterior_tibial_T4_R236 and v_T;
vars t and t;
enddef;
def map between Systemic and internal_iliac_T1_L196_module for
vars u_ivl and u_out;
vars v_internal_iliac_T1_L196 and v_T;
vars t and t;
enddef;
def map between Systemic and profundus_T2_L202_module for
vars u_ivl and u_out;
vars v_profundus_T2_L202 and v_T;
vars t and t;
enddef;
def map between Systemic and anterior_tibial_T3_L208_module for
vars u_ivl and u_out;
vars v_anterior_tibial_T3_L208 and v_T;
vars t and t;
enddef;
def map between Systemic and posterior_tibial_T4_L214_module for
vars u_ivl and u_out;
vars v_posterior_tibial_T4_L214 and v_T;
vars t and t;
enddef;
def map between Systemic and posterior_interosseous_T3_R40_module for
vars u_svl and u_out;
vars v_posterior_interosseous_T3_R40 and v_T;
vars t and t;
enddef;
def map between Systemic and ulnar_T2_R42_module for
vars u_svl and u_out;
vars v_ulnar_T2_R42 and v_T;
vars t and t;
enddef;
def map between Systemic and radial_T1_R44_module for
vars u_svl and u_out;
vars v_radial_T1_R44 and v_T;
vars t and t;
enddef;
def map between Systemic and posterior_interosseous_T3_L88_module for
vars u_svl and u_out;
vars v_posterior_interosseous_T3_L88 and v_T;
vars t and t;
enddef;
def map between Systemic and ulnar_T2_L90_module for
vars u_svl and u_out;
vars v_ulnar_T2_L90 and v_T;
vars t and t;
enddef;
def map between Systemic and radial_T1_L92_module for
vars u_svl and u_out;
vars v_radial_T1_L92 and v_T;
vars t and t;
enddef;
def map between Systemic and external_carotid_T2_R26_module for
vars u_svl and u_out;
vars v_external_carotid_T2_R26 and v_T;
vars t and t;
enddef;
def map between Systemic and external_carotid_T2_L62_module for
vars u_svl and u_out;
vars v_external_carotid_T2_L62 and v_T;
vars t and t;
enddef;
def map between Systemic and internal_carotid_R8_C_module for
vars u_svl and u_out;
vars v_internal_carotid_R8_C and v_T;
vars t and t;
enddef;
def map between Systemic and internal_carotid_L50_C_module for
vars u_svl and u_out;
vars v_internal_carotid_L50_C and v_T;
vars t and t;
enddef;
def map between Systemic and vertebral_R272_module for
vars u_svl and u_out;
vars v_vertebral_R272 and v_T;
vars t and t;
enddef;
def map between Systemic and vertebral_L2_module for
vars u_svl and u_out;
vars v_vertebral_L2 and v_T;
vars t and t;
enddef;
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
////////////////////////////////////////////// Importing Vessels //////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
// Aorta
def import using "BG_Modules.cellml" for
comp ascending_aorta_A_module using comp vv_type;
enddef;
def import using "BG_Modules.cellml" for
comp ascending_aorta_B_module using comp pv_type;
enddef;
def import using "BG_Modules.cellml" for
comp ascending_aorta_C_module using comp pv_type;
enddef;
def import using "BG_Modules.cellml" for
comp ascending_aorta_D_module using comp pv_type;
enddef;
def import using "BG_Modules.cellml" for
comp aortic_arch_C2_module using comp pv_jII_type;
enddef;
def import using "BG_Modules.cellml" for
comp brachiocephalic_trunk_C4_module using comp pv_jII_type;
enddef;
def import using "BG_Modules.cellml" for
comp aortic_arch_C46_module using comp pv_jII_type;
enddef;
def import using "BG_Modules.cellml" for
comp aortic_arch_C64_module using comp pv_jII_type;
enddef;
def import using "BG_Modules.cellml" for
comp aortic_arch_C94_module using comp pv_type;
enddef;
def import using "BG_Modules.cellml" for
comp thoracic_aorta_C96_module using comp pv_jII_type;
enddef;
def import using "BG_Modules.cellml" for
comp thoracic_aorta_C100_module using comp pv_jII_type;
enddef;