{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"institution":[{"name":"bioRxiv"}],"indexed":{"date-parts":[[2026,1,15]],"date-time":"2026-01-15T19:47:47Z","timestamp":1768506467818,"version":"3.49.0"},"posted":{"date-parts":[[2023,6,1]]},"group-title":"Bioengineering","reference-count":62,"publisher":"openRxiv","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"accepted":{"date-parts":[[2023,6,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                <jats:p>Clot formation is a crucial process that prevents bleeding, but can lead to severe disorders when imbalanced. This process is regulated by the coagulation cascade, a biochemical network that controls the enzyme thrombin, which converts soluble fibrinogen into the fibrin fibers that constitute clots. Coagulation cascade models are typically complex and involve dozens of partial differential equations (PDEs) representing various chemical species\u2019 transport, reaction kinetics, and diffusion. Solving these PDE systems computationally is challenging, due to their large size and multi-scale nature.<\/jats:p>\n                <jats:p>\n                  We propose a multi-fidelity strategy to increase the efficiency of coagulation cascade simulations. Leveraging the slower dynamics of molecular diffusion, we transform the governing PDEs into ordinary differential equations (ODEs) representing the evolution of species concentrations versus blood residence time. We then Taylor-expand the ODE solution around the zero-diffusivity limit to obtain spatiotemporal maps of species concentrations in terms of the statistical moments of residence time,\n                  <jats:inline-formula>\n                    <jats:alternatives>\n                      <jats:inline-graphic xmlns:xlink=\"http:\/\/www.w3.org\/1999\/xlink\" xlink:href=\"542763v1_inline1.gif\"\/>\n                    <\/jats:alternatives>\n                  <\/jats:inline-formula>\n                  , and provide the governing PDEs for\n                  <jats:inline-formula>\n                    <jats:alternatives>\n                      <jats:inline-graphic xmlns:xlink=\"http:\/\/www.w3.org\/1999\/xlink\" xlink:href=\"542763v1_inline2.gif\"\/>\n                    <\/jats:alternatives>\n                  <\/jats:inline-formula>\n                  . This strategy replaces a high-fidelity system of\n                  <jats:italic>N<\/jats:italic>\n                  PDEs representing the coagulation cascade of\n                  <jats:italic>N<\/jats:italic>\n                  chemical species by\n                  <jats:italic>N<\/jats:italic>\n                  ODEs and\n                  <jats:italic>p<\/jats:italic>\n                  PDEs governing the residence time statistical moments. The multi-fidelity order(\n                  <jats:italic>p<\/jats:italic>\n                  ) allows balancing accuracy and computational cost, providing a speedup of over\n                  <jats:italic>N\/p<\/jats:italic>\n                  compared to high-fidelity models.\n                <\/jats:p>\n                <jats:p>\n                  Using a simplified coagulation network and an idealized aneurysm geometry with a pulsatile flow as a benchmark, we demonstrate favorable accuracy for low-order models of\n                  <jats:italic>p<\/jats:italic>\n                  = 1 and\n                  <jats:italic>p<\/jats:italic>\n                  = 2. These models depart from the high-fidelity solution by under 16% (\n                  <jats:italic>p<\/jats:italic>\n                  = 1) and 5% (\n                  <jats:italic>p<\/jats:italic>\n                  = 2) after 20 cardiac cycles.\n                <\/jats:p>\n                <jats:p>The favorable accuracy and low computational cost of multi-fidelity models could enable unprecedented coagulation analyses in complex flow scenarios and extensive reaction networks. Furthermore, it can be generalized to advance our understanding of other systems biology networks affected by blood flow.<\/jats:p>","DOI":"10.1101\/2023.05.29.542763","type":"posted-content","created":{"date-parts":[[2023,6,1]],"date-time":"2023-06-01T13:25:18Z","timestamp":1685625918000},"source":"Crossref","is-referenced-by-count":1,"title":["Efficient multi-fidelity computation of blood coagulation under flow"],"prefix":"10.64898","author":[{"given":"Manuel","family":"Guerrero-Hurtado","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-6953-2270","authenticated-orcid":false,"given":"Manuel","family":"Garcia-Villalba","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7492-5741","authenticated-orcid":false,"given":"Alejandro","family":"Gonzalo","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-0511-7120","authenticated-orcid":false,"given":"Pablo","family":"Martinez-Legazpi","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Andy M.","family":"Kahn","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-1684-2432","authenticated-orcid":false,"given":"Elliot","family":"McVeigh","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5241-2758","authenticated-orcid":false,"given":"J.","family":"Bermejo","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-5683-0239","authenticated-orcid":false,"given":"Juan C.","family":"del Alamo","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2365-0738","authenticated-orcid":false,"given":"Oscar","family":"Flores","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"54368","reference":[{"key":"2023060508000474000_2023.05.29.542763v1.1","doi-asserted-by":"publisher","DOI":"10.1136\/hrt.82.5.547"},{"key":"2023060508000474000_2023.05.29.542763v1.2","doi-asserted-by":"publisher","DOI":"10.1016\/j.jtbi.2008.04.015"},{"key":"2023060508000474000_2023.05.29.542763v1.3","doi-asserted-by":"publisher","DOI":"10.1152\/ajpheart.00461.2014"},{"issue":"09","key":"2023060508000474000_2023.05.29.542763v1.4","doi-asserted-by":"crossref","first-page":"1969","DOI":"10.1142\/S0218127402005637","article-title":"Spatio-temporal dynamics of blood coagulation and pattern formation: an experimental study","volume":"12","year":"2002","journal-title":"Int. J. Bifurcat. Chaos"},{"issue":"6","key":"2023060508000474000_2023.05.29.542763v1.5","doi-asserted-by":"crossref","first-page":"619","DOI":"10.1070\/PU2002v045n06ABEH001090","article-title":"A new class of stopping self-sustained waves: a factor determining the spatial dynamics of blood coagulation","volume":"45","year":"2002","journal-title":"Phys-Usp+"},{"key":"2023060508000474000_2023.05.29.542763v1.6","first-page":"266","article-title":"An integrated fluid-chemical model toward modeling the formation of intra-luminal thrombus in abdominal aortic aneurysms","volume":"3","year":"2012","journal-title":"Front. Physiol"},{"issue":"10","key":"2023060508000474000_2023.05.29.542763v1.7","doi-asserted-by":"crossref","first-page":"1010","DOI":"10.1016\/j.jacc.2022.01.011","article-title":"Left ventricular thrombus following acute myocardial infarction: Jacc state-of-the-art review","volume":"79","year":"2022","journal-title":"J. Am. Coll. Cardiol"},{"issue":"25","key":"2023060508000474000_2023.05.29.542763v1.8","first-page":"2619","article-title":"Antithrombotic therapy in abdominal aortic aneurysm: beneficial or detrimental?","volume":"132","year":"2018","journal-title":"Am. J. Hematol"},{"key":"2023060508000474000_2023.05.29.542763v1.9","doi-asserted-by":"publisher","DOI":"10.1016\/j.jns.2006.12.008"},{"key":"2023060508000474000_2023.05.29.542763v1.10","doi-asserted-by":"publisher","DOI":"10.1126\/science.145.3638.1310"},{"issue":"716","key":"2023060508000474000_2023.05.29.542763v1.11","first-page":"519","volume":"10","year":"2020","journal-title":"Mapping and quantifying blood stasis and thrombus risk in the heart"},{"key":"2023060508000474000_2023.05.29.542763v1.12","doi-asserted-by":"publisher","DOI":"10.1098\/rspa.2014.0163"},{"key":"2023060508000474000_2023.05.29.542763v1.13","doi-asserted-by":"publisher","DOI":"10.1136\/heartjnl-2020-316977"},{"key":"2023060508000474000_2023.05.29.542763v1.14","doi-asserted-by":"publisher","DOI":"10.1159\/000089933"},{"key":"2023060508000474000_2023.05.29.542763v1.15","doi-asserted-by":"crossref","first-page":"2","DOI":"10.1016\/j.euromechflu.2012.01.019","article-title":"High-speed particle image velocimetry to assess cardiac fluid dynamics in vitro: From performance to validation","volume":"35","year":"2012","journal-title":"Eur. J. Mech. B\/Fluids"},{"key":"2023060508000474000_2023.05.29.542763v1.16","doi-asserted-by":"crossref","unstructured":"O. Flores , L. Rossini , A Gonzalo , D. Vigneault , J. Bermejo , AM Kahn , E. McVeigh , M. Garcia-Villalba , and J C. del Alamo . Evaluation of blood stasis in the left atrium using patient-specific direct numerical simulations. In ERCOFTAC Workshop Direct and Large Eddy Simulation, pages 485\u2013490. Springer, 2019.","DOI":"10.1007\/978-3-030-42822-8_64"},{"key":"2023060508000474000_2023.05.29.542763v1.17","doi-asserted-by":"publisher","DOI":"10.1109\/TMI.2005.859204"},{"key":"2023060508000474000_2023.05.29.542763v1.18","doi-asserted-by":"crossref","first-page":"596596","DOI":"10.3389\/fphys.2021.596596","article-title":"Demonstration of patient-specific simulations to assess left atrial appendage thrombogenesis risk","volume":"12","year":"2021","journal-title":"Frontiers Physiol"},{"key":"2023060508000474000_2023.05.29.542763v1.19","doi-asserted-by":"publisher","DOI":"10.1016\/S0894-7317(99)70105-7"},{"issue":"6","key":"2023060508000474000_2023.05.29.542763v1.20","doi-asserted-by":"crossref","first-page":"e3597","DOI":"10.1002\/cnm.3597","article-title":"Non-newtonian blood rheology impacts left atrial stasis in patient-specific simulations","volume":"38","year":"2022","journal-title":"Int. J. Numer. Method Biomed. Eng"},{"issue":"9","key":"2023060508000474000_2023.05.29.542763v1.21","doi-asserted-by":"crossref","first-page":"e1009331","DOI":"10.1371\/journal.pcbi.1009331","article-title":"Computational modeling of blood component transport related to coronary artery thrombosis in kawasaki disease","volume":"17","year":"2021","journal-title":"PLoS Comput. Biol"},{"key":"2023060508000474000_2023.05.29.542763v1.22","doi-asserted-by":"publisher","DOI":"10.1111\/j.1542-4758.2006.00119.x"},{"key":"2023060508000474000_2023.05.29.542763v1.23","doi-asserted-by":"publisher","DOI":"10.1007\/s10439-013-0853-z"},{"key":"2023060508000474000_2023.05.29.542763v1.24","unstructured":"C. Hirsch . Numerical computation of internal and external flows. (2nd edition). Elsevier, 2007."},{"key":"2023060508000474000_2023.05.29.542763v1.25","doi-asserted-by":"publisher","DOI":"10.1074\/jbc.M201173200"},{"key":"2023060508000474000_2023.05.29.542763v1.26","doi-asserted-by":"crossref","unstructured":"K. It\u00f4 , P. Henry Jr , et al. Diffusion processes and their sample paths: Reprint of the 1974 edition. Springer Science & Business Media, 1996.","DOI":"10.1007\/978-3-642-62025-6"},{"issue":"1","key":"2023060508000474000_2023.05.29.542763v1.27","doi-asserted-by":"crossref","first-page":"202","DOI":"10.1006\/jcph.1996.0130","article-title":"Efficient implementation of weighted ENO schemes","volume":"126","year":"1996","journal-title":"J. Comput. Phy"},{"issue":"37","key":"2023060508000474000_2023.05.29.542763v1.28","doi-asserted-by":"crossref","first-page":"23367","DOI":"10.1016\/S0021-9258(17)31662-9","article-title":"A model for the tissue factor pathway to thrombin. II. A mathematical simulation","volume":"269","year":"1994","journal-title":"J. Biol. Chem"},{"key":"2023060508000474000_2023.05.29.542763v1.29","doi-asserted-by":"publisher","DOI":"10.1114\/1.1349704"},{"key":"2023060508000474000_2023.05.29.542763v1.30","doi-asserted-by":"crossref","unstructured":"JP Keener and James Sneyd . Mathematical physiology 1: Cellular physiology. Springer New York, NY, USA, 2009.","DOI":"10.1007\/978-0-387-75847-3_1"},{"issue":"5","key":"2023060508000474000_2023.05.29.542763v1.31","doi-asserted-by":"crossref","first-page":"e2821","DOI":"10.1002\/cnm.2821","article-title":"A comparison of hemodynamic metrics and intraluminal thrombus burden in a common iliac artery aneurysm","volume":"33","year":"2017","journal-title":"Int. J. Numer. Method Biomed. Eng"},{"key":"2023060508000474000_2023.05.29.542763v1.32","doi-asserted-by":"publisher","DOI":"10.1146\/annurev.fluid.29.1.399"},{"key":"2023060508000474000_2023.05.29.542763v1.33","doi-asserted-by":"publisher","DOI":"10.1016\/S0006-3495(01)76085-7"},{"issue":"37","key":"2023060508000474000_2023.05.29.542763v1.34","doi-asserted-by":"crossref","first-page":"23357","DOI":"10.1016\/S0021-9258(17)31661-7","article-title":"A model for the tissue factor pathway to thrombin. i. an empirical study","volume":"269","year":"1994","journal-title":"J Biol Chem"},{"key":"2023060508000474000_2023.05.29.542763v1.35","doi-asserted-by":"publisher","DOI":"10.1093\/imammb\/dqq005"},{"key":"2023060508000474000_2023.05.29.542763v1.36","doi-asserted-by":"crossref","first-page":"S12","DOI":"10.1016\/j.thromres.2014.03.005","article-title":"An overview of mathematical modeling of thrombus formation under flow","volume":"133","year":"2014","journal-title":"Thromb. Res"},{"key":"2023060508000474000_2023.05.29.542763v1.37","doi-asserted-by":"crossref","first-page":"111211","DOI":"10.1016\/j.jbiomech.2022.111211","article-title":"Flow residence time in intracranial aneurysms evaluated by in vitro 4d flow mri","volume":"141","year":"2022","journal-title":"J. Biomech"},{"issue":"2-3","key":"2023060508000474000_2023.05.29.542763v1.38","doi-asserted-by":"crossref","first-page":"121","DOI":"10.1159\/000089932","article-title":"The effect of convective flows on blood coagulation processes","volume":"34","year":"2005","journal-title":"Pathophysiol. Haemos. Thromb"},{"key":"2023060508000474000_2023.05.29.542763v1.39","doi-asserted-by":"publisher","DOI":"10.1038\/202498a0"},{"key":"2023060508000474000_2023.05.29.542763v1.40","doi-asserted-by":"publisher","DOI":"10.1097\/MAT.0b013e3180683b7c"},{"key":"2023060508000474000_2023.05.29.542763v1.41","doi-asserted-by":"crossref","first-page":"395","DOI":"10.1017\/jfm.2017.508","article-title":"On the aerodynamic forces on heaving and pitching airfoils at low reynolds number","volume":"828","year":"2017","journal-title":"J. Fluid Mech"},{"key":"2023060508000474000_2023.05.29.542763v1.42","doi-asserted-by":"crossref","first-page":"306","DOI":"10.3389\/fphys.2018.00306","article-title":"Thrombosis in cerebral aneurysms and the computational modeling thereof: a review","volume":"9","year":"2018","journal-title":"Front. Physiol"},{"key":"2023060508000474000_2023.05.29.542763v1.43","doi-asserted-by":"publisher","DOI":"10.1177\/0954411912444080"},{"key":"2023060508000474000_2023.05.29.542763v1.44","doi-asserted-by":"publisher","DOI":"10.1529\/biophysj.105.069062"},{"key":"2023060508000474000_2023.05.29.542763v1.45","doi-asserted-by":"publisher","DOI":"10.23919\/CinC53138.2021.9662687"},{"key":"2023060508000474000_2023.05.29.542763v1.46","doi-asserted-by":"publisher","DOI":"10.1161\/ATVBAHA.114.304488"},{"key":"2023060508000474000_2023.05.29.542763v1.47","doi-asserted-by":"publisher","DOI":"10.1007\/s10441-019-09372-w"},{"key":"2023060508000474000_2023.05.29.542763v1.48","doi-asserted-by":"publisher","DOI":"10.1007\/s10439-010-0065-8"},{"key":"2023060508000474000_2023.05.29.542763v1.49","doi-asserted-by":"publisher","DOI":"10.1115\/1.2910291"},{"key":"2023060508000474000_2023.05.29.542763v1.50","unstructured":"Stanley L Robbins and Ramzi S Cotran . Pathologic basis of disease. Saunders, 1979."},{"issue":"11","key":"2023060508000474000_2023.05.29.542763v1.51","doi-asserted-by":"crossref","first-page":"2152","DOI":"10.1016\/j.jbiomech.2015.11.049","article-title":"A clinical method for mapping and quantifying blood stasis in the left ventricle","volume":"49","year":"2016","journal-title":"J. Biomech"},{"issue":"12","key":"2023060508000474000_2023.05.29.542763v1.52","doi-asserted-by":"crossref","first-page":"1557","DOI":"10.1002\/ange.200353428","article-title":"Minimal functional model of hemostasis in a biomimetic microfluidic system","volume":"116","year":"2004","journal-title":"Angew. Chem"},{"issue":"11","key":"2023060508000474000_2023.05.29.542763v1.53","doi-asserted-by":"crossref","first-page":"H1567","DOI":"10.1152\/ajpheart.00855.2015","article-title":"A coupled chemo-fluidic computational model for thrombogenesis in infarcted left ventricles","volume":"310","year":"2016","journal-title":"Amer. J. Physiol.-Heart Circul. Physiol"},{"key":"2023060508000474000_2023.05.29.542763v1.54","doi-asserted-by":"publisher","DOI":"10.1023\/B:PHAM.0000008049.40370.5a"},{"issue":"2","key":"2023060508000474000_2023.05.29.542763v1.55","doi-asserted-by":"crossref","first-page":"740","DOI":"10.1016\/j.bbe.2020.02.010","article-title":"Experimental and numerical flow analysis through arteries with stent using particle image velocimetry and computational fluid dynamics method","volume":"40","year":"2020","journal-title":"Biocybern. Biomed. Eng"},{"issue":"2","key":"2023060508000474000_2023.05.29.542763v1.56","doi-asserted-by":"crossref","first-page":"448","DOI":"10.1016\/j.jcp.2005.03.017","article-title":"An immersed boundary method with direct forcing for the simulation of particulate flows","volume":"209","year":"2005","journal-title":"J. Comput. Phys"},{"issue":"2","key":"2023060508000474000_2023.05.29.542763v1.57","first-page":"263","article-title":"Broad-based intracranial aneurysms: thrombosis induced by stent placement","volume":"24","year":"2003","journal-title":"Am. J. Neuroradiol"},{"key":"2023060508000474000_2023.05.29.542763v1.58","doi-asserted-by":"crossref","unstructured":"C. Vlachopoulos , M. O\u2019Rourke , and W. W. Nichols . McDonald\u2019s blood flow in arteries: theoretical, experimental and clinical principles. CRC press, 2011.","DOI":"10.1201\/b13568"},{"key":"2023060508000474000_2023.05.29.542763v1.59","doi-asserted-by":"publisher","DOI":"10.1161\/CIRCRESAHA.115.306841"},{"key":"2023060508000474000_2023.05.29.542763v1.60","doi-asserted-by":"publisher","DOI":"10.1146\/annurev.bioeng.1.1.299"},{"issue":"1","key":"2023060508000474000_2023.05.29.542763v1.61","doi-asserted-by":"crossref","first-page":"e1005291","DOI":"10.1371\/journal.pcbi.1005291","article-title":"A general shear-dependent model for thrombus formation","volume":"13","year":"2017","journal-title":"PLoS Comput. Biol"},{"key":"2023060508000474000_2023.05.29.542763v1.62","doi-asserted-by":"publisher","DOI":"10.1063\/1.1345728"}],"container-title":[],"original-title":[],"link":[{"URL":"https:\/\/syndication.highwire.org\/content\/doi\/10.1101\/2023.05.29.542763","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,1,15]],"date-time":"2026-01-15T04:43:14Z","timestamp":1768452194000},"score":1,"resource":{"primary":{"URL":"http:\/\/biorxiv.org\/lookup\/doi\/10.1101\/2023.05.29.542763"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,6,1]]},"references-count":62,"URL":"https:\/\/doi.org\/10.1101\/2023.05.29.542763","relation":{"is-preprint-of":[{"id-type":"doi","id":"10.1371\/journal.pcbi.1011583","asserted-by":"subject"}]},"subject":[],"published":{"date-parts":[[2023,6,1]]},"subtype":"preprint"}}