## 5 Data Availability  
The model source code is made openly available via GitHub using which all data in this work can be reproduced. **Supplementary Materials:** The following supporting information can be downloaded at [[https://www.mdpi.com/article/10.3390/app12010307/s1](https://www.mdpi.com/article/10.3390/app12010307/s1)]([https://www.mdpi.com/article/10.3390/app12010307/s1](https://www.mdpi.com/article/10.3390/app12010307/s1)), Table S1: Model resistance, compliance, and cardiac elastance parameters used to generate populations; Figure S1: Schematic diagram of the detailed right-left kidneys' complex. Each kidney consists of six microvascular components, implemented as three element Windkessel models. Each component consists of two microvascular beds that is assigned two resistance values and one compliance value. All resistances and com-pliances under baseline conditions were assumed to have the same value (see Table S1 for pa-rameter values); Figure S2: Schematic of the dialysis unit adapted from the literature. Further details including baseline parameter values used in this work; Figure S3: Schematic of the baroreflex sub-model; Figure S4: Representative parameter Gaussian distributions in the baseline model population (see Figure 2, main manuscript) where mean values are taken from the unperturbed model, and a coefficient of variation of 0.25 is used to generate the population. Left panel shows distribution of aortic compliance, and right panel shows distribution of left kidney inlet resistance; Figure S5: Distribution of flows in model organs under baseline (grey), hemodialysis (red), therapeutic hypothermia (blue), and simultaneous hemodialysis with therapeutic hypothermia (green). Top row: upper body; second row: splanchic; third row: lower body; bottom row: two kidney com-plex. This illustration accompanies Figures 2 and 3 in the main manuscript; Figure S6: Baroreflex control of heart rate by the parasympathetic and sympathetic tone mechanisms. The control value is shown by the black circle in the bottom left hand corner of the panel. Both parasympathetic (horizontal axis) and sympathetic (vertical axis) rate values are shown after subtracting the base-line value. See Table S1 for parameter values; Figure S7: Blood flow under baseline (left column), AF (middle column), and exercise (rightcolumn) conditions. All data are normalized to the baseline values of the left panel; Figure S8: Population distributions of intracellular volume (left panel), interstitial volume (middle panel), and plasma volume (right panel). See Figures 2 and 3 and related text for details. Study design: J.J.J., S.R.K., D.G. and C.W.M. ; model development, J.J.J., T.J.H., S.R.K. and C.S. simulations: J.J.J. and S.R.K. ; data analysis and figures, J.J.J., T.J.H. and C.S. ; draft writing, J.J.J., C.S. and S.R.K. ; draft completion, S.R.K. ; final manuscript, S.R.K., J.J.J., D.G. and C.W.M. ; funding secured by, S.R.K., C.W.M. and D.G. All authors have read and agreed to the published version of the manuscript. This research was funded by Canada Canarie Inc. grant number [RS3-111], Canada Heart and Stroke Foundation grant number [G-20-0028717] and the APC was funded by the publishing MDPI journal. This study did not use sensitive data, patient information, or any other identity data. This computational study did not require ethics approvals of any nature. This computational study did not require consent from patients or otherwise. The MATLAB source code to reproduce all data in this work is openly available in our lab's github repository with URL: [[https://github.com/mccssk2/HumanCirculationModel202_MDPI](https://github.com/mccssk2/HumanCirculationModel202_MDPI)]([https://github.com/mccssk2/HumanCirculationModel202_MDPI](https://github.com/mccssk2/HumanCirculationModel202_MDPI)) (accessed on 21 December 2021). Computationally efficient source code can be obtained from the corresponding authors upon reasonable request.