pith. sign in

arxiv: 2510.26138 · v2 · submitted 2025-10-30 · ⚛️ physics.flu-dyn

Diffuse interface approach to oxygen transport and metabolism under blood flow dynamics in microcirculations

Naoki Takeishi , Junya Kobayashi , Shigeo Wada , Satoshi Ii This is my paper

Pith reviewed 2026-05-18 03:47 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords diffuse interface methodoxygen transportred blood cellsmicrocirculationcapillary networksimmersed boundary methodtissue oxygenationblood flow dynamics
0
0 comments X

The pith

Red blood cells autonomously adjust oxygen supply to produce uniform tissue oxygenation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a computational approach to simulate how oxygen moves and gets consumed while red blood cells move and deform inside capillary networks. It rewrites the oxygen transport equations as mixture forms that fold all interface rules directly into the equations, so moving boundaries do not need separate tracking. This setup is paired with an immersed-boundary treatment of the fluid flow to run everything on one fixed mesh that covers cells, plasma, and tissue. The resulting calculations for straight vessels and branching networks indicate that each red blood cell senses nearby tissue oxygen and alters its release accordingly. A reader would care because the outcome points to a built-in mechanism that keeps oxygen levels steady across tissue even when blood flow changes.

Core claim

The paper proposes a diffuse interface approach for oxygen transport using a mixture formulation. Oxygen transport is formulated as an advection-diffusion-reaction equation, and all governing equations are rewritten in mixture forms using phase indicator functions so that every interface condition is carried inside the equations themselves. This removes the need to treat discontinuities on largely moving and deforming surfaces even when red blood cells are packed densely. Cellular flow is handled as a fluid-membrane interaction problem with the immersed boundary method. The combined scheme runs the coupled flow and oxygen problems across cytoplasm, plasma, and tissue on a single fixed mesh.

What carries the argument

Mixture formulation with phase indicator functions that folds every interface condition for oxygen transport and metabolism into a single advection-diffusion-reaction equation on a fixed mesh.

If this is right

  • Simulations of oxygen delivery become practical for dense, deforming red blood cells inside branching capillary networks without explicit interface tracking.
  • Each red blood cell alters its oxygen release in response to the local tissue oxygenation level it encounters.
  • This local adjustment produces spatially homogeneous oxygen levels throughout the surrounding tissue.
  • The same fixed-mesh framework couples fluid dynamics and reaction-diffusion transport across internal cell fluid, external plasma, and tissue.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same mixture approach could be tested on other diffusion-reaction systems that involve many moving deformable objects.
  • If the self-regulation holds, tissue oxygenation may stay stable even under irregular or pulsatile flow patterns not yet simulated.
  • Adding spatial variation in tissue metabolism rate would provide a direct test of how robust the homogeneity result remains.

Load-bearing premise

The mixture formulation using phase indicator functions correctly embeds all interface conditions for oxygen transport and metabolism without loss of accuracy even under highly dense RBC conditions in capillary networks.

What would settle it

High-resolution measurements of oxygen concentration at many points inside living tissue surrounding a capillary network, taken while flow rate and local oxygen demand are varied, would show whether the predicted uniform distribution actually occurs.

Figures

Figures reproduced from arXiv: 2510.26138 by Junya Kobayashi, Naoki Takeishi, Satoshi Ii, Shigeo Wada.

Figure 1
Figure 1. Figure 1: (a) Schematic of the oxygen transport system in the pre [PITH_FULL_IMAGE:figures/full_fig_p027_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) Schematic of an finite interface region [PITH_FULL_IMAGE:figures/full_fig_p028_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Definition of the domains for a spherical diffusion proble [PITH_FULL_IMAGE:figures/full_fig_p028_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Numerical validation of oxygen transport with a moving (bu [PITH_FULL_IMAGE:figures/full_fig_p029_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) Computational domains for oxygen transport with RB [PITH_FULL_IMAGE:figures/full_fig_p030_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Microvascular network model, where the middle area consis [PITH_FULL_IMAGE:figures/full_fig_p031_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: (a–c) Snapshots of the distribution of RBCs ( [PITH_FULL_IMAGE:figures/full_fig_p032_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (a–c) Snapshots of oxygen saturation s and oxygen concentration ci at t = 0.6 s for Hcta = 0.112 (a), Hcta = 0.224 (b), and Hcta = 0.336 (c). (d) Ratio of oxygen concentration Rhc3i and RBC phase Rhψ1i at t = 1 s. (e and f) Time-averaged oxygen extraction rates ∆s and the number of RBCs (nRBC ) passing through each Area for Hcta = 0.244 and 0.336, where time averaging is performed over 0.2 s between t = 0.… view at source ↗
Figure 9
Figure 9. Figure 9: (a) Snapshot of the analysis area in the modified microvasc [PITH_FULL_IMAGE:figures/full_fig_p034_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: (a and b) Relationships between the apparent viscosity [PITH_FULL_IMAGE:figures/full_fig_p035_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: (a and b) Snapshots of the distribution of oxygen satur [PITH_FULL_IMAGE:figures/full_fig_p036_11.png] view at source ↗
read the original abstract

The relationship between the spatiotemporal distribution of oxygen transport and blood flow dynamics, accounting for the motion and deformation of individual red blood cells (RBCs), is of fundamental importance for understanding microcirculation systems. Three-dimensional (3D) modeling is indispensable for addressing complex oxygen transport and cellular behaviors in capillary networks; however, the computational approach is formidable for enforcing interface (or jump) conditions on largely moving and deforming interfaces. In this paper, we propose a diffuse interface approach for the oxygen transport using a mixture formulation. We formulate oxygen transport using an advection-diffusion-reaction equation and rewrite all governing equations in mixture forms using phase indicator functions, where all the interface conditions are included in the governing equations. This innovation avoids the complexity associated with discontinuities for largely moving interfaces in highly dense RBC conditions. We model cellular flow as a fluid-membrane interaction problem using the immersed boundary method (IBM). The method allows the seamless calculation of coupling problems for cellular flows and oxygen transports in the cytoplasm (internal fluid) of the RBC, plasma (external fluid), and tissue regions using a fixed Cartesian coordinate mesh. The proposed method accurately captures the analytical solution for spherically symmetric diffusion, and successfully demonstrates oxygen transport in both straight capillaries and their networks. These rigorous analyses suggest that RBCs can autonomously regulate the oxygen supply to tissues in response to the local tissue oxygenation level, resulting in the establishment of homogeneous tissue oxygenation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper develops a diffuse-interface method for oxygen transport and metabolism coupled to RBC motion and deformation in microcirculatory networks. The advection-diffusion-reaction equations are rewritten in mixture form via phase-indicator functions so that all interface conditions are embedded in the governing equations on a fixed Cartesian mesh; RBC flow is handled by the immersed-boundary method. The approach is shown to recover the analytical solution for spherically symmetric diffusion and is then applied to straight capillaries and capillary networks, from which the authors conclude that RBCs autonomously regulate oxygen supply in response to local tissue oxygenation, producing homogeneous tissue oxygen levels.

Significance. If the mixture formulation is shown to preserve interface conditions to sufficient accuracy under dense RBC packing, the method would supply a practical tool for three-dimensional simulations of coupled flow and oxygen transport without explicit interface tracking. The reported match to the spherical-diffusion solution and the demonstration in network geometries constitute concrete strengths; the autonomous-regulation observation, if robust, would be a noteworthy mechanistic insight for microcirculation physiology.

major comments (2)
  1. [Abstract] Abstract (paragraph on governing-equation rewrite): the central claim that RBCs autonomously regulate oxygen supply to yield homogeneous tissue oxygenation rests on the assertion that the mixture formulation with phase-indicator functions correctly embeds all jump conditions for oxygen flux and metabolism. No derivation of the mixture equations, no explicit verification that normal-flux continuity and reaction terms remain accurate when RBCs approach contact or walls, and no mesh-convergence or error-bar data are supplied; these omissions directly affect the load-bearing interpretation of the network results.
  2. [Abstract] Abstract (validation paragraph): the statement that the method 'accurately captures the analytical solution for spherically symmetric diffusion' is presented without quantitative error norms, grid-resolution details, or comparison against a sharp-interface reference solution at comparable hematocrit; without such metrics it is unclear whether the observed regulation effect survives the O(ε) smoothing errors inherent to the diffuse interface.
minor comments (1)
  1. [Abstract] The abstract does not specify the precise form of the phase-indicator functions or the value of the smoothing length ε used in the reported simulations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful and constructive review of our manuscript. The comments on the need for explicit derivation of the mixture equations and quantitative validation metrics are well taken. We address each major comment below and will incorporate revisions to strengthen the presentation and support for our conclusions.

read point-by-point responses

  1. Referee: [Abstract] Abstract (paragraph on governing-equation rewrite): the central claim that RBCs autonomously regulate oxygen supply to yield homogeneous tissue oxygenation rests on the assertion that the mixture formulation with phase-indicator functions correctly embeds all jump conditions for oxygen flux and metabolism. No derivation of the mixture equations, no explicit verification that normal-flux continuity and reaction terms remain accurate when RBCs approach contact or walls, and no mesh-convergence or error-bar data are supplied; these omissions directly affect the load-bearing interpretation of the network results.

    Authors: We appreciate the referee highlighting the importance of a transparent derivation. The manuscript introduces the mixture formulation via phase-indicator functions to embed interface conditions on a fixed mesh, but we acknowledge that an expanded, step-by-step derivation would improve clarity. In the revised manuscript we will add a dedicated subsection deriving the mixture advection-diffusion-reaction equations from the sharp-interface jump conditions, explicitly showing how normal-flux continuity and volumetric reaction terms are preserved. We will also include targeted numerical tests for RBCs in near-contact and near-wall configurations, together with mesh-convergence studies and error bars on the network results, to confirm that the autonomous-regulation observation remains robust. revision: yes

  2. Referee: [Abstract] Abstract (validation paragraph): the statement that the method 'accurately captures the analytical solution for spherically symmetric diffusion' is presented without quantitative error norms, grid-resolution details, or comparison against a sharp-interface reference solution at comparable hematocrit; without such metrics it is unclear whether the observed regulation effect survives the O(ε) smoothing errors inherent to the diffuse interface.

    Authors: We agree that quantitative error metrics are essential to substantiate the validation claim. The revised manuscript will report L2 and L∞ error norms for the spherically symmetric diffusion test, specify the grid resolutions employed, and add a side-by-side comparison against a sharp-interface reference solution at comparable hematocrit. These additions will directly address whether the homogeneous oxygenation pattern persists within the O(ε) smoothing tolerance of the diffuse interface. revision: yes

Circularity Check

0 steps flagged

No circularity: forward simulation of mixture formulation yields regulation result as output

full rationale

The paper introduces a diffuse-interface mixture formulation by rewriting the advection-diffusion-reaction equation with phase indicator functions so that interface conditions are embedded in the governing equations. This modeling step is presented as a direct rewriting that avoids explicit discontinuity handling; it is not derived from or fitted to the target homogeneous-oxygenation outcome. The method is validated against an independent analytical solution for spherically symmetric diffusion, then applied via IBM-coupled simulations to capillary networks. The observed autonomous regulation and homogeneous tissue oxygenation emerge as simulation outputs under the stated assumptions, without reduction to a fitted parameter, self-citation chain, or definitional equivalence. The central claim therefore remains independent of its own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard transport equations plus the unproven assumption that the phase-indicator mixture form preserves interface conditions exactly for moving dense RBCs; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • standard math Oxygen transport obeys an advection-diffusion-reaction equation in each region
    Invoked when the authors formulate oxygen transport and then rewrite it in mixture form.
  • domain assumption Phase indicator functions allow all interface conditions to be included inside the governing equations without explicit enforcement
    This is the key modeling step stated in the abstract as the innovation that avoids discontinuity handling.

pith-pipeline@v0.9.0 · 5793 in / 1401 out tokens · 50003 ms · 2026-05-18T03:47:37.481985+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

79 extracted references · 79 canonical work pages

  1. [1]

    Krogh, The progress of physiology, Science 70 (1929) 200–2 04

    A. Krogh, The progress of physiology, Science 70 (1929) 200–2 04. doi:10.1152/ajplegacy.1929.90.2.243

    work page doi:10.1152/ajplegacy.1929.90.2.243 1929

  2. [2]

    A. H. Ramsook, P. B. Dominelli, S. A. Angus, J. W. Senefeld, C. C. W igginsa, M. J. Joyner, The oxygen transport cascade and exercise: Less ons from compara- tive physiology, Comp. Biochem. Physol. A Mol. Integr. Physiol. 282 ( 2023) 111442. doi:10.1016/j.cbpa.2023.111442

    work page doi:10.1016/j.cbpa.2023.111442 2023

  3. [3]

    Manohar, Right heart pressures and blood-gas tensions in po nies during exercise and laryngeal hemiplegia, Am

    O. Manohar, Right heart pressures and blood-gas tensions in po nies during exercise and laryngeal hemiplegia, Am. J. Physiol. Heart Circ. Physiol. 251 (19 86) H121–H126. doi:10.1152/ajpheart.1986.251.1.H121

    work page doi:10.1152/ajpheart.1986.251.1.h121 1986

  4. [4]

    E. K. Weir, J. L´ opez-Barneo, K. J. Buckler, S. L. Archer, Mec hanisms of dis- ease acute oxygen-sensing mechanisms, N. Engl. J. Med. 353 (200 5) 2042–2055. doi:10.1056/NEJMc053396

    work page doi:10.1056/nejmc053396 2042

  5. [5]

    Krogh, The number and distribution of capillaries in muscles with c alculations of the oxygen pressure head necessary for supplying the tissue, J

    A. Krogh, The number and distribution of capillaries in muscles with c alculations of the oxygen pressure head necessary for supplying the tissue, J. Phy siol. 52 (1919) 409–415. doi:10.1113/jphysiol.1919.sp001839

    work page doi:10.1113/jphysiol.1919.sp001839 1919

  6. [6]

    J. D. Hellums, The resistance to oxygen transport in the capillarie s rela- tive to that in the surrounding tissue, Microvasc. Res. 13 (1977) 1 31–136. doi:10.1016/0026-2862(77)90122-4

    work page doi:10.1016/0026-2862(77)90122-4 1977

  7. [7]

    W. J. Federspiel, A. S.Popel, A theoretical analysis of the effect o f the particulate nature of blood on oxygen release in capillaries, Microvasc. Res. 32 ( 1986) 164–189. doi:10.1016/0026-2862(86)90052-x

    work page doi:10.1016/0026-2862(86)90052-x 1986

  8. [8]

    J. D. Hellums, P. K. Nair, N. S. Huang, N. Ohshima, Simulation of intr aluminal gas transport processes in the microcirculation, Ann. Biomed. Eng 24 (1995) 1–24. doi:10.1007/BF02770991

    work page doi:10.1007/bf02770991 1995

  9. [9]

    C. H. Wang, A. S. Popel, Effect of red blood cell shape on oxygen t ransport in capillaries, Math Biosci. 116 (1993) 89–110. doi:10.1016/0025-5564(93)90062-f

    work page doi:10.1016/0025-5564(93)90062-f 1993

  10. [10]

    C. D. Eggleton, A. Vadapalli, T. K. Roy, A. S. Popel, Calculations o f intracap- illary oxygen tension distributions in muscle, Math. Biosci. 167 (2000) 123–143. doi:10.1016/s0025-5564(00)00038-9

    work page doi:10.1016/s0025-5564(00)00038-9 2000

  11. [11]

    L¨ ucker, B

    A. L¨ ucker, B. Weber, P. Jenny, A dynamic model of oxygen tr ansport from capillaries to tissue with moving red blood cells, Am. J. Physiol. Heart. Circ. Phys iol. 308 (2015) H206–H216. doi:10.1152/ajpheart.00447.2014

    work page doi:10.1152/ajpheart.00447.2014 2015

  12. [12]

    L¨ ucker, T

    A. L¨ ucker, T. W. Secomb, B. Weber, P. Jenny, The relative infl uence of hematocrit and red blood cell velocity on oxygen transport from capillaries to tissue , Microcirculation 24 (2017) e12337. doi:10.1111/micc.12337. 37

    work page doi:10.1111/micc.12337 2017

  13. [13]

    C. B. Schaffer, B. Friedman, N. Nishimura, L. F. Schroeder, P. S. Tsai, F. F. Ebner, P. D. Lyden, D. Kleinfeld, Two-photon imaging of cortical surface m icrovessels reveals a robust redistributionin blood flow after vascular occlusion, PLoS B iol. 4 (2006) e22. doi:10.1371/journal.pbio.0040022

    work page doi:10.1371/journal.pbio.0040022 2006

  14. [14]

    Lecoq, A

    J. Lecoq, A. Parpaleix, E. Roussakis, M. Ducros, Y. G. Housse n, S. A. Vinogradov, S. Charpak, Simultaneous two-photon imaging of oxygen and blood fl ow in deep cerebral vessels, Nat. Med. 17 (2011) 893–898. doi:10.1038/nm.2394

    work page doi:10.1038/nm.2394 2011

  15. [15]

    D. G. Lyons, A. Parpaleix, M. Roche, S. Charpak, Mapping oxyg en concentration in the awake mouse brain, eLIFE 5 (2016) e12024. doi:10.7554/eLife.12024

    work page doi:10.7554/elife.12024 2016

  16. [16]

    Moeini, X

    M. Moeini, X. Lu, P. K. Avti, R. Damseh, S. B. F. Picard, D. Boas, A. Kakkar, F. Lesage, Compromised microvascular oxygen delivery increases b rain tissue vulnera- bility with age, Sci. Rep. 8 (2018) 8219. doi:10.1038/s41598-018-26543-w

    work page doi:10.1038/s41598-018-26543-w 2018

  17. [17]

    Gagnon, A

    L. Gagnon, A. F. Smith, D. A. Boas, A. Devor, T. W. Secomb, S. Sakadˇ zi´ c, Modeling of cerebral oxygen transport based on in vivo microscopic imaging o f microvascular network structure, blood flow, and oxygenation, Front. Comput . Neurosci. 10 (2016)

    work page 2016

  18. [18]

    doi:10.3389/fncom.2016.00082

    work page doi:10.3389/fncom.2016.00082 2016

  19. [19]

    I. G. Gould, P. Tsai, D. Kleinfeld, A. Linninger, The capillary bed off ers the largest hemodynamic resistance to the cortical blood supply, J. Cereb. Blo od Flow Metab. 37 (2017) 52––68. doi:10.1177/0271678X16671146

    work page doi:10.1177/0271678x16671146 2017

  20. [20]

    Hartung, C

    G. Hartung, C. Vesel, R. Morley, A. Alaraj, J. Sled, D. Kleinfeld, A. Linninger, Simulations of blood as a suspension predicts a depth dependent hem atocrit in the circulation throughout the cerebral cortex, PLoS Comput. Biol. 1 4 (2018) e1006549. doi:10.1371/journal.pcbi.1006549

    work page doi:10.1371/journal.pcbi.1006549 2018

  21. [21]

    P. W. Sweeney, S. Walker-Samuel, R. J. Shipley, Insights into ce rebral haemodynamics and oxygenation utilising in vivo mural cell imaging and mathematical mo delling, Sci. Rep. 8 (2018) 2373. doi:10.1038/s41598-017-19086-z

    work page doi:10.1038/s41598-017-19086-z 2018

  22. [22]

    Goldman, A

    D. Goldman, A. S. Popel, A computational study of the effect of capillary network anastomoses and tortuosity on oxygen transport, J. Theor. Bio l. 206 (2000) 181–194. doi:10.1006/jtbi.2000.2113

    work page doi:10.1006/jtbi.2000.2113 2000

  23. [23]

    K¨ oppl, E

    T. K¨ oppl, E. Vidotto, B. Wohlmuth, A 3D-1D coupled blood flow an d oxygen transport model to generate microvascular networks, Int. J. Numer. Meth . Biomed. Engng. 36 (2020) e3386. doi:10.1002/cnm.3386

    work page doi:10.1002/cnm.3386 2020

  24. [24]

    Takeishi, Y

    N. Takeishi, Y. Imai, K. Nakaaki, T. Yamaguchi, T. Ishikawa, Leu kocyte margination at arteriole shear rate, Physiol. Rep. 2 (2014) e12037. doi:10.14814/phy2.12037

    work page doi:10.14814/phy2.12037 2014

  25. [25]

    Takeishi, Y

    N. Takeishi, Y. Imai, T. Yamaguchi, T. Ishikawa, Flow of a circulat ing tu- mor cell and red blood cells in microvessels, Phys. Rev. E. 92 (2015) 0 63011. doi:10.1103/PhysRevE.92.063011. 38

    work page doi:10.1103/physreve.92.063011 2015

  26. [26]

    S. Ii, K. Shimizu, K. Sugiyama, S. Takagi, Continuum and stochas tic approach for cell adhesion process based on eulerian fluid-capsule coupling with lagran gian markers, J. Comput. Phys. 374 (2018) 769–786. doi:10.1016/j.jcp.2018.08.002

    work page doi:10.1016/j.jcp.2018.08.002 2018

  27. [27]

    Balogh, P

    P. Balogh, P. Bagchi, Acomputationalapproachtomodeling cellula r-scale blood flow in complex geometry, J. Comput. Phys. 334 (2018) 280–3 07. doi:10.1016/j.jcp.2017.01.007

    work page doi:10.1016/j.jcp.2017.01.007 2018

  28. [28]

    Balogh, P

    P. Balogh, P. Bagchi, Direct numerical simulation of cellular-scale blood flow in 3D microvascular networks, Biophys. J. 113 (2017) 2815–28 26. doi:10.1016/j.bpj.2017.10.020

    work page doi:10.1016/j.bpj.2017.10.020 2017

  29. [29]

    Gibou, L

    F. Gibou, L. Chen, D. N. adn S. Banerjee, A level-set based sha rp interface method for the multiphase incompressible navier–stokes equations with phase c hange, J. Comput. Phys. 222 (2007) 536–555. doi:10.1016/j.jcp.2006.07.035

    work page doi:10.1016/j.jcp.2006.07.035 2007

  30. [30]

    Gibou, R

    F. Gibou, R. P. Fedkiw, L.-T. Cheng, M. Kang, A second-order- accurate symmetric discretization of the poisson equation on irregular domains, Journa l of Computational Physics 176 (1) (2002) 205–227. doi:10.1006/jcph.2001.6977

    work page doi:10.1006/jcph.2001.6977 2002

  31. [31]

    R. J. LeVeque, Z. Li, The immersed interface method for elliptic e quations with discon- tinuous coefficients and singular sources, SIAM J. Numer. Anal. 31 ( 4) (1994) 1019–1044. doi:10.1137/0731054

    work page doi:10.1137/0731054 1994

  32. [32]

    A. T. Layton, Modeling water transport across elastic bounda ries using an explicit jump method, SIAM J. Numer. Anal. 28 (6) (2006) 2189–2207. doi:10.1137/050642198

    work page doi:10.1137/050642198 2006

  33. [33]

    P. G. Jayathilake, B. C. Khoo, Z. Than, N. E. Wijeysundera, De formation and os- motic swelling of an elastic membrane capsule in stokes flows by the imme rsed interface method, Chem. Eng. Sci. 65 (2010) 1237–1252. doi:10.1016/j.ces.2009.09.078

    work page doi:10.1016/j.ces.2009.09.078 2010

  34. [34]

    P. G. Jayathilake, B. C. Khoo, Z. Than, N. E. Wijeysundera, Ef ect of membrane perme- ability on capsule substrate adhesion: Computation using immersed in terface method, Chem. Eng. Sci. 65 (2010) 3565–3578. doi:10.1016/j.ces.2010.02.050

    work page doi:10.1016/j.ces.2010.02.050 2010

  35. [35]

    P. G. Jayathilake, G. Liu, Z. Tan, B. C. Khoo, Numerical study o f a permeable capsule under stokes flows by the immersed interface method, Chem. Eng. Sci. 66 (2011) 2080–

    work page 2011

  36. [36]

    doi:10.1016/j.ces.2011.02.005

    work page doi:10.1016/j.ces.2011.02.005 2011

  37. [37]

    Miyauchi, S

    S. Miyauchi, S. Takeuchi, T. Kajishima, A numerical method for m ass transfer by a thin moving membrane with selective permeabilities, J. Comput. Phys. 284 (2015) 490–504. doi:10.1016/j.jcp.2014.12.048

    work page doi:10.1016/j.jcp.2014.12.048 2015

  38. [38]

    Z. Wang, F. Colin, G. Le, J. Zhang, Counter-extrapolation met hod for conjugate heat and mass transfer with interfacial discontinuity, Int. J. Numer. M ethods Heat Fluid Flow 27 (10) (2017) 2231–2258. doi:10.1108/HFF-10-2016-0422 . 39

    work page doi:10.1108/hff-10-2016-0422 2017

  39. [39]

    F. A. Amiri, J. Zhang, Oxygen transport across tank-treadin g red blood cell: Individual and joint roles of flow convection and oxygen-hemoglobin reaction, Microvasc. Res. 145 (2023) 104447. doi:10.1016/j.mvr.2022.104447

    work page doi:10.1016/j.mvr.2022.104447 2023

  40. [40]

    F. A. Amiri, J. Zhang, Tissue oxygenation around capillaries: Effe cts of hematocrit and arteriole oxygen condition, Bull. Math. Biol. 85 (2023 ) 50. doi:10.1007/s11538-023-01155-2

    work page doi:10.1007/s11538-023-01155-2 2023

  41. [41]

    Z. Zhao, J. Yan, Enriched immersed boundary method (eibm) fo r interface-coupled multi-physics and applications to convective conjugate heat trans fer, Comput. Meth. Appl. Mech. Eng. 401 (2022) 115667. doi:10.1016/j.cma.2022.115667

    work page doi:10.1016/j.cma.2022.115667 2022

  42. [42]

    Henneaux, P

    D. Henneaux, P. Schrooyen, P. Chatelain, T. Magin, High-orde r enforcement of jumps conditions between compressible viscous phases: An extended inte rior penalty discon- tinuous galerkin method for sharp interface simulation, Comput. Me th. Appl. Mech. Eng. 415 (2023) 116215. doi:10.1016/j.cma.2023.116215

    work page doi:10.1016/j.cma.2023.116215 2023

  43. [43]

    Huang, K

    H. Huang, K. Sugiyama, S. Takagi, An immersed boundary metho d for restricted diffusion with permeable interfaces, J. Comput. Phys. 228 (15) (20 09) 5317–5322. doi:10.1016/j.jcp.2009.04.040

    work page doi:10.1016/j.jcp.2009.04.040 2009

  44. [44]

    X. Gong, Z. Gong, H. Huang, An immersed boundary method for mass transfer across permeable moving interfaces, J. Comput. Phys. 278 (2014) 148–1 68

    work page 2014

  45. [45]

    X. Wang, X. Gong, K. Sugiyama, S. Takagi, H. Huang, An immerse d boundary method for mass transfer through porous biomembranes under large def ormations, J. Comput. Phys. 413 (2020) 109444. doi:10.1016/j.jcp.2020.109444

    work page doi:10.1016/j.jcp.2020.109444 2020

  46. [46]

    F. A. Amiri, J. Zhang, An immersed membrane method for mass tr ansfer across flexible semipermeable membranes in flows, Int. Commun. Heat Mass Transf . 128 (2021) 105601. doi:10.1016/j.icheatmasstransfer.2021.105601

    work page doi:10.1016/j.icheatmasstransfer.2021.105601 2021

  47. [47]

    Reder, A

    M. Reder, A. Prahs, D. Schneider, B. Nestler, Viscous stress approximations in diffuse interface methods for two-phase flow based on mechanical jump c onditions, Comput. Meth. Appl. Mech. Eng. 432 (2024) 117341. doi:10.1016/j.cma.2024.117341

    work page doi:10.1016/j.cma.2024.117341 2024

  48. [48]

    S. K. Stoter, P. M¨ uller, L. Cicalese, M. Tuveri, D. Schillinger, T. J. Hughes, A diffuse interface method for the navier–stokes/darcy equations: Perf usion profile for a patient- specific human liver based on mri scans, Comput. Meth. Appl. Mech. Eng. 321 (2017) 70–102. doi:10.1016/j.cma.2017.04.002

    work page doi:10.1016/j.cma.2017.04.002 2017

  49. [49]

    C. S. Peskin, The immersed boundary method, Acta Numer. 11 ( 2002) 479–517. doi:10.1017/S0962492902000077

    work page doi:10.1017/s0962492902000077 2002

  50. [50]

    Clark Jr., W

    A. Clark Jr., W. J. Federspiel, P. A. A. Clark, G. R. Cokelet, Oxyg en delivery from red cells, Biophys. J. 47 (1985) 171–181. doi:10.1016/s0006-3495(85)83890-x

    work page doi:10.1016/s0006-3495(85)83890-x 1985

  51. [51]

    Goldman, Theoretical models of microvascular oxygen trans port to tissue, Microcir- culation 15 (2008) 795–811

    D. Goldman, Theoretical models of microvascular oxygen trans port to tissue, Microcir- culation 15 (2008) 795–811. doi:10.1080/10739680801938289. 40

    work page doi:10.1080/10739680801938289 2008

  52. [52]

    Nicoli, M

    M. Nicoli, M. Plapp, H. Henry, Tensorial mobilities for accurate so lution of trans- port problems in models with diffuse interfaces, Phys. Rev. E 84 (201 1) 046707. doi:10.1103/PhysRevE.84.046707

    work page doi:10.1103/physreve.84.046707

  53. [53]

    Skalak, A

    R. Skalak, A. Tozeren, R. P. Zarda, S. Chien, Strain energy fu nction of red blood cell membranes, Biophys. J. 13 (1973) 245–264. doi:10.1016/S0006-3495(73)85983-1

    work page doi:10.1016/s0006-3495(73)85983-1 1973

  54. [54]

    Zhong-can, W

    O.-Y. Zhong-can, W. Helfrich, Bending energy of vesicle membra nes: Gen- eral expressions for the first, second, and third variation of the shape energy and applications to spheres and cylinders, Phys. Rev. A 39 (1989) 5 280–5288. doi:10.1103/PhysRevA.39.5280

    work page doi:10.1103/physreva.39.5280 1989

  55. [55]

    Xu, Convergence of discrete laplace-beltrami operators o ver surfaces, Comput

    G. Xu, Convergence of discrete laplace-beltrami operators o ver surfaces, Comput. Math. Appl. 48 (3) (2004) 347–360. doi:10.1016/j.camwa.2004.05.001

    work page doi:10.1016/j.camwa.2004.05.001 2004

  56. [56]

    Reuter, S

    M. Reuter, S. Biasotti, D. Giorgi, G. Patan` e, M. Spagnuolo, Dis crete laplace–beltrami operators for shape analysis and segmentation, Comput. Graph. 33 (3) (2009) 381–390. doi:10.1016/j.cag.2009.03.005

    work page doi:10.1016/j.cag.2009.03.005 2009

  57. [57]

    Yazdani, P

    A. Yazdani, P. Bagchi, Three-dimensional numerical simulation o f vesicle dy- namics using a front-tracking method, Phys. Rev. E 85 (2012) 056 308. doi:10.1103/PhysRevE.85.056308

    work page doi:10.1103/physreve.85.056308 2012

  58. [58]

    S. O. Unverdi, G. Tryggvason, A front-tracking method for v iscous, in- compressible, multi-fluid flows, J. Comput. Phys. 100 (1992) 25–37 . doi:10.1016/0021-9991(92)90307-K

    work page doi:10.1016/0021-9991(92)90307-k 1992

  59. [59]

    S. Ii, K. Sugiyama, S. Takeuchi, S. Takagi, Y. Matsumoto, F. Xia o, An interface capturing method with a continuous function: the THINC method with multi-dimensional reconstruction, J. Comput. Phys. 231 (201 2) 2328–2358. doi:10.1016/j.jcp.2011.11.038

    work page doi:10.1016/j.jcp.2011.11.038 2011

  60. [60]

    Weymouth, D

    G. Weymouth, D. K. Yue, Boundary data immersion method for c artesian-grid simu- lations of fluid-body interaction problems, Journal of Computation al Physics 230 (16) (2011) 6233–6247. doi:10.1016/j.jcp.2011.04.022

    work page doi:10.1016/j.jcp.2011.04.022 2011

  61. [61]

    L. Fu, X. Y. Hu, N. A. Adams, A family of high-order targeted en o schemes for compressible-fluid simulations, J. Comput. Phys. 305 ( 2016) 333–359. doi:10.1016/j.jcp.2015.10.037

    work page doi:10.1016/j.jcp.2015.10.037 2016

  62. [62]

    H. F. Harlow, J. E. Welch, Numerical calculation of time–depende nt viscous incompress- ible flow of fluid with free surface, Phys. Fluids 8 (1965) 2182. doi:10.1063/1.1761178

    work page doi:10.1063/1.1761178 1965

  63. [63]

    S. Ii, M. A. H. M. Adib, Y. Watanabe, S. Wada, Physically consiste nt data assimilation method based on feedback control for patient-specific blood flow analysis, Int. J. Numer. Methods Biomed Eng. 34 (1) (2018) e2910. doi:10.1002/cnm.2910

    work page doi:10.1002/cnm.2910 2018

  64. [64]

    Jiang and C.-W

    G.-S. Jiang, C.-W. Shu, Efficient implementation of weighted ENO sc hemes, J. Comput. Phys. 126 (1996) 202–228. doi:10.1006/jcph.1996.0130. 41

    work page doi:10.1006/jcph.1996.0130 1996

  65. [65]

    F. Xiao, Y. Honma, T. kono, A simple algebraic interface capturin g scheme using hy- perbolic tangent function, Int. J. Numer. Meth. Fluid. 48 (2005) 1 023–1040

    work page 2005

  66. [66]

    Mohandas, P

    N. Mohandas, P. G. Gallagher, Red cell membrane: past, prese nt, and future, Blood 112 (2008) 3939–3948

    work page 2008

  67. [67]

    Harkness, R

    J. Harkness, R. B. Whittington, Blood-plasma viscosity: an app roximate temperature- invariant arising from generalised concepts, Biorheology 6 (1970) 1 69–187

    work page 1970

  68. [68]

    Evans, Y

    E. Evans, Y. C. Fung, Improved measurements of the erythr ocyte geometry, Microvasc. Res. 4 (1972) 335–347

    work page 1972

  69. [69]

    Takeishi, Y

    N. Takeishi, Y. Imai, Capture of microparticles by bolus flow of re d blood cells in capillaries, Sci Rep 7 (2017) 5381. doi:10.1038/s41598-017-05924-7

    work page doi:10.1038/s41598-017-05924-7 2017

  70. [70]

    A. R. Pries, D. Neuhaus, P. Gaehtgens, Blood viscosity in tube fl ow: dependence on diameter and hematocrit, Am. J. Physiol. Heart Circ. Physiol. 263 (6 ) (1992) H1770– H1778. doi:10.1152/ajpheart.1992.263.6.H1770

    work page doi:10.1152/ajpheart.1992.263.6.h1770 1992

  71. [71]

    Hartridge, F

    H. Hartridge, F. J. W. Roughton, The rate of distribution of dis solved gases between the red blood corpuscle and its fluid environment, J. Physiol. 62 (192 7) 232–242

    work page

  72. [72]

    Z. Sun, S. Inaba, F. Xiao, Boundary variation diminishing (bvd) r econstruction: A new approach to improve godunov schemes, J. Comput. Phys. 322 (2016) 309–325. doi:10.1016/j.jcp.2016.06.051

    work page doi:10.1016/j.jcp.2016.06.051 2016

  73. [73]

    Wakimura, T

    H. Wakimura, T. Li, K.-M. Shyue, T. Aoki, F. Xiao, High-resolution boundary variation diminishing scheme for two-phase compressible flow with cavitation an d evaporation, J. Comput. Phys. 513 (2024) 113164. doi:10.1016/j.jcp.2024.113164

    work page doi:10.1016/j.jcp.2024.113164 2024

  74. [74]

    H. G. Augustin, G. Y. Koh, Organotypic vascularture: from de scriptive heterogeneity to functional pathophysiology, Science 357 (2017) 771. doi:10.1126/science.aal2379

    work page doi:10.1126/science.aal2379 2017

  75. [75]

    Popel, Theory of oxygen transport to tissue, Crit

    A. Popel, Theory of oxygen transport to tissue, Crit. Rev. Bio med. Eng. 17 (1989) 257–321

    work page 1989

  76. [76]

    B. Li, T. V. Esipova, I. Sencan, K. Kilic, B. Fu, M. Desjardins, M. Moeini, S. Kura, M. A. Yaseen, F. Lesage, L. Østergaard, A. Devor, D. A. Boas, S. A. V inogradov, S. Sakadˇ zi´ c, More homogeneous capillary flow and oxygenation in deeper cortical layers correlate with increased oxygen extraction, eLife 8 (2019) e42299. doi:10.7554/eLife.42299

    work page doi:10.7554/elife.42299 2019

  77. [77]

    T. J. McMahon, Red blood cell deformability, vasoactive mediato rs, and adhesion, Front. Physiol. 10 (2019) 1417. doi:10.3389/fphys.2019.01417

    work page doi:10.3389/fphys.2019.01417 2019

  78. [78]

    D. G. Buerk, Can we model nitric oxide biotransport? A survey o f mathematical models for a simple diatomic molecule with surprisingly complex biological activitie s, Annu. Rev. Biomed. Eng. 3 (2001) 109–143. doi:10.1146/annurev.bioeng.3.1.109

    work page doi:10.1146/annurev.bioeng.3.1.109 2001

  79. [79]

    Carmeliet, R

    P. Carmeliet, R. K. Jain, Principles and mechanisms of vessel nor malization for cancer and other angiogenic diseases, Nat. Rev. Drug Discov. 10 ( 2011) 417–427. doi:10.1038/nrd3455. 42

    work page doi:10.1038/nrd3455 2011