pith. sign in

arxiv: 2512.03446 · v3 · submitted 2025-12-03 · ⚛️ physics.bio-ph · physics.flu-dyn

A novel multiscale modelling for the hemodynamics in retinal microcirculation with an analytic solution for the capillary-tissue coupled system

Chang Lin , Zilong Song , Robert Eisenberg , Shixin Xu , Huaxiong Huang This is my paper

Pith reviewed 2026-05-17 02:43 UTC · model grok-4.3

classification ⚛️ physics.bio-ph physics.flu-dyn
keywords retinal microcirculationhemodynamicsmultiscale modelingDarcy equationsanalytic solutioncapillary-tissue couplingblood flow simulation
0
0 comments X

The pith

A multiscale retinal blood flow model uses an analytic solution to couple capillaries and tissue with one-dimensional vessel flows.

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

The paper develops a mathematical model to simulate blood movement through the retina's vessels and surrounding tissue, addressing the challenge of connecting different scales in the vasculature for studying eye diseases. It pairs one-dimensional models for larger arterioles and venules with Darcy equations for the capillary network and tissue, but solves the capillary-tissue part exactly in closed form. This analytic solution speeds up calculations, gives clearer physical insight into the flows, and creates a dynamic link between the capillary bed and the incoming and outgoing vessel flows. A reader would care because such models can help explain how changes in blood pressure or vessel properties affect retinal health without heavy computation. The approach is checked for mathematical consistency and compared to data from experiments and prior simulations.

Core claim

The central claim is that an analytic solution exists for the coupled Darcy system describing flow in the capillary bed and interstitial tissue. This solution supplies explicit expressions that replace numerical solution of the tissue-capillary equations, yielding faster computation and a direct interpretation of pressure and velocity fields. It also generates a dynamic coupling condition that transmits information between the capillary level and the one-dimensional models of upstream arterioles and downstream venules. The resulting multiscale framework is shown to be mathematically robust by truncation-error and convergence analysis, and its outputs are compared with experimental retinal-he

What carries the argument

The analytic solution of the capillary-tissue Darcy system, which supplies closed-form expressions for pressure and velocity that enforce mass conservation and enable direct coupling to 1D vessel equations.

If this is right

  • The model supplies a practical tool for exploring how changes in arteriolar or venular resistance alter tissue perfusion.
  • Parameter studies become feasible at low computational cost because the capillary-tissue part is solved exactly.
  • The dynamic coupling condition provides a mathematically consistent way to embed capillary beds inside larger vascular network models.
  • Validation against existing data sets supports use for interpreting clinical observations of retinal blood flow.

Where Pith is reading between the lines

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

  • The same analytic coupling could be tested in other microvascular beds such as the brain or kidney where similar scale separation exists.
  • If vessel-wall compliance were added to the 1D segments, the analytic tissue solution might still close the system and allow study of pulse-wave effects on capillary exchange.
  • Patient-specific retinal geometries obtained from imaging could be inserted into the 1D network while retaining the same capillary-tissue solution.

Load-bearing premise

The capillary bed behaves as a porous medium governed by Darcy flow whose parameters and boundary conditions remain representative across the retinal conditions examined.

What would settle it

Direct comparison of the model's predicted capillary pressures or flow rates against high-resolution intravital measurements in living retinas would show large discrepancies under parameter values taken from the paper.

Figures

Figures reproduced from arXiv: 2512.03446 by Chang Lin, Huaxiong Huang, Robert Eisenberg, Shixin Xu, Zilong Song.

Figure 1
Figure 1. Figure 1: Spatial discretization of vessel. where U n+1 i = U(si , tn+1), and similar notations are used for flux and source vectors. Thus FVM equation becomes U n+1 i = U n i − ∆t ∆s  F n+1/2 i+1/2 − F n+1/2 i−1/2  + ∆t 2  S n+1/2 i+1/2 + S n+1/2 i−1/2  . (17) The unknown interface variables F n+1/2 i+1/2 , F n+1/2 i−1/2 , S n+1/2 i+1/2 and S n+1/2 i−1/2 depend on U n+1/2 m with m = i ± 1/2, which can be determ… view at source ↗
Figure 2
Figure 2. Figure 2: Synthetic and real retinal vasculature. (a) The arterial (red) and venous (blue) trees [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The dependence of relative ratio Em,m+1 on the order m for rn = 0.25Rt. Solid lines denote the values of relative ratio while dashed lines are the upper bound. Different colours correspond to various values of radial coordinate. (a) Value of drainage rate is of α = 2 × 10−9 cm · s/g. (b) Value of drainage rate is of α = 2 × 10−8 cm · s/g. (c) Value of drainage rate is of α = 2 × 10−7 cm · s/g. Other parame… view at source ↗
Figure 4
Figure 4. Figure 4: The dependence of relative ratio Em,m+1 on the order m for rn = 0.25Rt. Solid lines denote the values of relative ratio while dashed lines are the upper bound. Different colours correspond to various values of radial coordinate. The drainage rate is α = 2 × 10−6 cm · s/g. (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Two different bifurcation boundaries. obtain Rm,n(r) ≤ R0,n(r)(min{r, rn}/ max{r, rn}) m and following estimate for the truncation error: EM ≤ 1 π X∞ m=M+1 |Rm,n(r)| ≤ 1 π |R0,n(r)| X∞ m=M+1  min{r, rn} max{r, rn} m . (86) For any r ̸= rn the series P∞ m=M+1(min{r, rn}/ max{r, rn}) m converges, then the truncation error satisfies EM ≤ 1 π |R0,n(r)| (min{r, rn}/ max{r, rn})M 1 − (min{r, rn}/ max{r, rn}) .… view at source ↗
Figure 6
Figure 6. Figure 6: The dependence of radial component Rm,n on the order m for rn = 0.25Rt. Circles are the exact values of radial component computed by equation (68) while lines denote the approximation (126). Different colours correspond to various values of radial coordinate. (a) Value of drainage rate is of α = 2 × 10−9 cm · s/g. (b) Value of drainage rate is of α = 2 × 10−8 cm · s/g. (c) Value of drainage rate is of α = … view at source ↗
Figure 7
Figure 7. Figure 7: The dependence of radial component Rm,n on the order m for rn = 0.25Rt. Circles are the exact values of radial component computed by equation (68) while lines denote the approximation (126). Different colours correspond to various values of radial coordinate. The drainage rate is α = 2 × 10−6 cm · s/g. and Ak and qk denote the area and flow rate of the terminal vessels connected to arterioles or veins, res… view at source ↗
Figure 8
Figure 8. Figure 8: The behaviour of f2 on r/Rt with rn = 0.25Rt and θn = 3π/2. Circles are results of the truncation (78) with orderM = 26, while lines are the values of approximation (128). Different colours correspond to various values of θ. The tissue permeability is set as kt = 2 × 10−12 cm2 . (a) Value of drainage rate is of α = 2 × 10−9 cm · s/g. (b) Value of drainage rate is of α = 2 × 10−8 cm · s/g. (c) Value of drai… view at source ↗
Figure 9
Figure 9. Figure 9: The dependence of total flow on truncation order [PITH_FULL_IMAGE:figures/full_fig_p025_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Flow rate and pressure difference distributions in arteries and veins. Lines denote [PITH_FULL_IMAGE:figures/full_fig_p025_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Capillary pressure, tissue pressure and exchange pressure. (a) Capillary pressure. [PITH_FULL_IMAGE:figures/full_fig_p026_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The dependence of total flow on capillary and tissue permeabilities. Other param [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Domain-averaged pressure as the function of permeabilities. (a) Averaged pressure [PITH_FULL_IMAGE:figures/full_fig_p027_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Domain-averaged exchange pressure as the function of permeabilities. (a) Aver [PITH_FULL_IMAGE:figures/full_fig_p028_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Total flow as a function of drainage rate for various permeability ratios. Line [PITH_FULL_IMAGE:figures/full_fig_p029_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: The pulsatile pressure profile at CRA inlet. [PITH_FULL_IMAGE:figures/full_fig_p030_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: The temporal evolution of flow rate at bifurcation boundaries of arterial and venous [PITH_FULL_IMAGE:figures/full_fig_p031_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: The evolution of pressure at bifurcation boundaries of arterial and venous trees. [PITH_FULL_IMAGE:figures/full_fig_p032_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Pulsatile flow rates and pressures as functions of diameter. (a) Averaged pulsatile [PITH_FULL_IMAGE:figures/full_fig_p032_19.png] view at source ↗
read the original abstract

Mathematical modelling of the microcirculatory hemodynamics in the retina is an essential tool for understanding various diseases of the retina, yet remains challenging due to the multiscale nature of the retinal vasculature and its coupling to surrounding tissue. To address this, we develop a multiscale model that couples retinal vasculature across scales with interstitial tissue. Our model combines the one-dimensional (1D) model for arterioles and venules with the coupled Darcy equations for capillaries and tissue. The model uses an analytic solution for capillary-tissue coupled system that provides a simple interpretation of the results along with much faster computation. The analytic solution implies a dynamic coupling condition that links the capillary bed with upstream arteriolar and downstream venular flows. The model is mathematically robust, demonstrated through analysis of the solution's truncation error and convergence. Its predictive accuracy is validated against experimental data and other models, making it useful in interpreting experimental results. Finally, the role of various parameters in controlling retinal hemodynamics is explored.

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 / 2 minor

Summary. The manuscript develops a multiscale hemodynamic model for retinal microcirculation that couples one-dimensional representations of arterioles and venules to a capillary-tissue subsystem governed by Darcy equations. An analytic solution is derived for the coupled capillary-tissue system; this solution is asserted to yield a dynamic coupling condition linking the capillary bed to upstream and downstream flows, to reduce computational cost relative to fully discrete approaches, and to admit straightforward interpretation. The model is reported to be mathematically robust on the basis of truncation-error analysis and convergence checks, and its predictions are validated against experimental data and other models. Parameter studies are used to explore control of retinal hemodynamics.

Significance. If the analytic solution and the implied dynamic coupling remain valid under physiological conditions, the approach would supply a computationally efficient framework for multiscale retinal flow modeling. The combination of an explicit analytic form, truncation-error bounds, and experimental validation would constitute a concrete advance over purely numerical multiscale schemes, enabling rapid parameter exploration and mechanistic interpretation of retinal hemodynamics in health and disease.

major comments (2)
  1. [Model formulation and analytic solution (abstract and §3)] The central claim that the analytic solution of the coupled Darcy system correctly implies a dynamic coupling condition to the 1D arteriolar/venular models rests on the adequacy of the Darcy continuum approximation for the discrete retinal capillary network. The manuscript does not supply a quantitative error bound or direct comparison against a discrete vessel network under realistic distributions of capillary diameter, tortuosity, and pressure drop; this comparison is load-bearing for the asserted computational speed-up and interpretive simplicity.
  2. [Analytic solution and error analysis] The truncation-error analysis and convergence checks cited in the abstract are presented without the explicit form of the analytic solution or the precise boundary conditions used to derive the dynamic coupling. Without these, it is not possible to verify that the coupling condition remains representative when permeability or source terms vary across the physiological range examined in the parameter study.
minor comments (2)
  1. [Validation against experimental data] Clarify the precise experimental datasets used for validation and any exclusion criteria applied; this would strengthen the claim of predictive accuracy.
  2. [Parameter selection] Ensure that all permeability and resistance parameters introduced in the Darcy subsystem are listed with their sources or fitting procedures so that readers can assess dependence on data.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their detailed and constructive report. We address the major comments point by point below. We agree that greater transparency regarding the analytic solution and the justification for the Darcy approximation is warranted, and we will revise the manuscript accordingly to incorporate these improvements.

read point-by-point responses

  1. Referee: [Model formulation and analytic solution (abstract and §3)] The central claim that the analytic solution of the coupled Darcy system correctly implies a dynamic coupling condition to the 1D arteriolar/venular models rests on the adequacy of the Darcy continuum approximation for the discrete retinal capillary network. The manuscript does not supply a quantitative error bound or direct comparison against a discrete vessel network under realistic distributions of capillary diameter, tortuosity, and pressure drop; this comparison is load-bearing for the asserted computational speed-up and interpretive simplicity.

    Authors: We acknowledge this valid point. The Darcy continuum approximation is a standard approach in retinal microcirculation modeling, supported by our truncation error analysis which bounds the discretization and continuum errors. To address the request for quantitative validation, we will add in the revised version a discussion of error estimates derived from literature on capillary network homogenization and provide a sample calculation of the approximation error for typical retinal capillary parameters (diameter ~5-10 μm, tortuosity effects averaged). A full discrete-to-continuum comparison simulation is not feasible within the current computational framework but is noted as a direction for future validation. This revision will clarify the robustness of the dynamic coupling condition. revision: partial

  2. Referee: [Analytic solution and error analysis] The truncation-error analysis and convergence checks cited in the abstract are presented without the explicit form of the analytic solution or the precise boundary conditions used to derive the dynamic coupling. Without these, it is not possible to verify that the coupling condition remains representative when permeability or source terms vary across the physiological range examined in the parameter study.

    Authors: We thank the referee for highlighting this presentation issue. The explicit analytic solution for the capillary-tissue system (solving the coupled Darcy equations with appropriate source terms for oxygen consumption) and the boundary conditions (Dirichlet-type matching to 1D vessel pressures at the capillary interfaces and no-flux conditions on the outer tissue boundaries) are derived in Section 3. In the revision, we will include the closed-form expressions in the main text or a new appendix for clarity. We will also augment the error analysis with additional results showing that the derived dynamic coupling condition (relating effective capillary pressure to flow rates) remains accurate for permeability variations of ±50% and source term changes within physiological ranges, supported by extended convergence studies. revision: yes

standing simulated objections not resolved
  • Performing a direct, quantitative comparison to a fully resolved discrete capillary network model incorporating realistic tortuosity and diameter distributions, which would require developing and running a separate high-fidelity discrete model not part of the present study.

Circularity Check

0 steps flagged

Analytic solution derived directly from governing equations; no circularity detected

full rationale

The paper derives its central analytic solution for the capillary-tissue coupled Darcy system from the standard governing PDEs (Darcy flow in capillaries and tissue with appropriate source/sink terms). This closed-form solution is obtained by mathematical solution of the boundary-value problem rather than by fitting to target outputs or by self-referential definition. The implied dynamic coupling condition to the 1D arteriolar/venular models follows directly from continuity of pressure and flux at the interfaces once the analytic form is substituted; it is not presupposed. Validation against external experimental data and other models occurs after the derivation and does not enter the solution procedure itself. No load-bearing self-citation, ansatz smuggling, or renaming of known results is present in the derivation chain. The model therefore remains self-contained against its own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on standard continuum assumptions for blood flow and tissue perfusion; no new particles or forces are introduced. Specific free parameters such as permeability coefficients or resistance values are likely fitted or chosen to match data, but exact values are not visible in the abstract.

axioms (2)
  • domain assumption Darcy's law governs flow in the capillary and tissue domains
    Invoked when the capillary-tissue system is modeled with coupled Darcy equations.
  • domain assumption One-dimensional Poiseuille-type flow applies to arterioles and venules
    Used for the 1D model of larger vessels.

pith-pipeline@v0.9.0 · 5488 in / 1306 out tokens · 62361 ms · 2026-05-17T02:43:11.505036+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

64 extracted references · 64 canonical work pages · 2 internal anchors

  1. [1]

    T. E. Kornfield, E. A. Newman, Measurement of Retinal Blood Flow Using Fluorescently Labeled Red Blood Cells, eNeuro 2 (2) (2015) ENEURO.0005–15.2015.doi:10.1523/ENEURO.0005-15.2015. URLhttps://www.eneuro.org/lookup/doi/10.1523/ENEURO.0005-15. 2015 39

    work page doi:10.1523/eneuro.0005-15.2015 2015

  2. [2]

    Garh¨ ofer, C

    G. Garh¨ ofer, C. Zawinka, H. Resch, K. H. Huemer, L. Schmetterer, G. T. Dorner, Response of Retinal Vessel Diameters to Flicker Stimulation in Patients with Early Open Angle Glaucoma:, J. Glaucoma 13 (4) (2004) 340–344.doi:10.1097/00061198-200408000-00013. URLhttp://journals.lww.com/00061198-200408000-00013

    work page doi:10.1097/00061198-200408000-00013 2004

  3. [3]

    Caprioli, A

    J. Caprioli, A. L. Coleman, Blood Pressure, Perfusion Pres- sure, and Glaucoma, Am. J. Ophthalmol. 149 (5) (2010) 704–712. doi:10.1016/j.ajo.2010.01.018. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0002939410000346

    work page doi:10.1016/j.ajo.2010.01.018 2010

  4. [4]

    T. W. Ch’ng, C. Y. Chua, M. A. Ummi Kalsom, Y. Azhany, V. H. M. Gong, A. H. G. Rasool, A. T. Liza-Sharmini, Ocular Perfusion Pressure and Severity of Glaucoma: Is There a Link?, J. Curr. Glaucoma Pract. 15 (2) (2021) 78–85.doi:10.5005/jp-journals-10078-1305. URLhttps://www.jocgp.com/doi/10.5005/jp-journals-10078-1305

    work page doi:10.5005/jp-journals-10078-1305 2021

  5. [5]

    N. N. Osborne, R. J. Casson, J. P. M. Wood, G. Chidlow, M. Graham, J. Melena, Retinal ischemia: mechanisms of damage and potential therapeutic strategies, Prog. Retin. Eye Res. 23 (1) (2004) 91–147. doi:10.1016/j.preteyeres.2003.12.001. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S1350946203000788

    work page doi:10.1016/j.preteyeres.2003.12.001 2004

  6. [6]

    Ruffini, M

    A. Ruffini, M. Dvoriashyna, A. Govetto, M. R. Romano, R. Repetto, A Mathematical Model of Interstitial Fluid Flow and Retinal Tissue Deformation in Macular Edema, Invest. Ophthalmol. Vis. Sci. 65 (11) (2024) 19.doi:10.1167/iovs.65.11.19. URLhttps://iovs.arvojournals.org/article.aspx?articleid= 2800791

    work page doi:10.1167/iovs.65.11.19 2024

  7. [7]

    J. H. Siggers, C. R. Ethier, Fluid Mechanics of the Eye, Annu. Rev. Fluid Mech. 44 (1) (2012) 347–372.doi:10.1146/ annurev-fluid-120710-101058. URLhttps://www.annualreviews.org/doi/10.1146/ annurev-fluid-120710-101058

    work page 2012

  8. [8]

    Flores, J

    J. Flores, J. Alastruey, E. Corvera Poir´ e, A Novel Analytical Approach to Pulsatile Blood Flow in the Arterial Network, Ann. Biomed. Eng. 44 (10) (2016) 3047–3068.doi:10.1007/s10439-016-1625-3. URLhttp://link.springer.com/10.1007/s10439-016-1625-3

    work page doi:10.1007/s10439-016-1625-3 2016

  9. [9]

    Guidoboni, A

    G. Guidoboni, A. Harris, S. Cassani, J. Arciero, B. Siesky, A. Amireskan- dari, L. Tobe, P. Egan, I. Januleviciene, J. Park, Intraocular Pressure, Blood Pressure, and Retinal Blood Flow Autoregulation: A Mathematical Model to Clarify Their Relationship and Clinical Relevance, Invest. Ophthalmol. Vis. Sci. 55 (7) (2014) 4105.doi:10.1167/iovs.13-13611. 40 U...

    work page doi:10.1167/iovs.13-13611 2014

  10. [10]

    Ganesan, S

    P. Ganesan, S. He, H. Xu, Analysis of retinal circulation using an image- based network model of retinal vasculature, Microvasc. Res. 80 (1) (2010) 99–109.doi:10.1016/j.mvr.2010.02.005. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0026286210000373

    work page doi:10.1016/j.mvr.2010.02.005 2010

  11. [11]

    B. C. Fry, E. B. Coburn, S. Whiteman, A. Harris, B. Siesky, J. Arciero, Predicting retinal tissue oxygenation using an image-based theoretical model, Math. Biosci. 305 (2018) 1–9.doi:10.1016/j.mbs.2018.08.005. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0025556417306831

    work page doi:10.1016/j.mbs.2018.08.005 2018

  12. [12]

    D. Liu, N. B. Wood, N. Witt, A. D. Hughes, S. A. Thom, X. Y. Xu, Com- putational Analysis of Oxygen Transport in the Retinal Arterial Network, Curr. Eye Res. 34 (11) (2009) 945–956.doi:10.3109/02713680903230079. URLhttp://www.tandfonline.com/doi/full/10.3109/ 02713680903230079

    work page doi:10.3109/02713680903230079 2009

  13. [13]

    Causin, G

    P. Causin, G. Guidoboni, F. Malgaroli, R. Sacco, A. Harris, Blood flow mechanics and oxygen transport and delivery in the retinal microcir- culation: multiscale mathematical modeling and numerical simulation, Biomech. Model. Mechanobiol. 15 (3) (2016) 525–542.doi:10.1007/ s10237-015-0708-7. URLhttp://link.springer.com/10.1007/s10237-015-0708-7

    work page doi:10.1007/s10237-015-0708-7 2016

  14. [14]

    Julien, S

    L. Julien, S. Bonnin, M. Paques, J.-M. Fullana, One-dimensional modeling of microvascular hemodynamics in the retina using multimodal imaging, Phys. Fluids 35 (6) (2023) 061901.doi:10.1063/5.0152499. URLhttps://pubs.aip.org/pof/article/35/6/061901/2894457/ One-dimensional-modeling-of-microvascular

    work page doi:10.1063/5.0152499 2023

  15. [15]

    Rebhan, L

    J. Rebhan, L. P. Parker, L. J. Kelsey, F. K. Chen, B. J. Doyle, A computational framework to investigate retinal haemodynamics and tis- sue stress, Biomech. Model. Mechanobiol. 18 (6) (2019) 1745–1757.doi: 10.1007/s10237-019-01172-y. URLhttp://link.springer.com/10.1007/s10237-019-01172-y

    work page doi:10.1007/s10237-019-01172-y 2019

  16. [16]

    Malek, A

    J. Malek, A. T. Azar, B. Nasralli, M. Tekari, H. Kamoun, R. Tourki, Computational analysis of blood flow in the retinal arteries and veins using fundus image, Comput. Math. Appl. 69 (2) (2015) 101–116. doi:10.1016/j.camwa.2014.11.017. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0898122114005859

    work page doi:10.1016/j.camwa.2014.11.017 2015

  17. [17]

    K. C. Tripathy, A. Siddharth, A. Bhandari, Image-based insil- ico investigation of hemodynamics and biomechanics in healthy 41 and diabetic human retinas, Microvasc. Res. 150 (2023) 104594. doi:10.1016/j.mvr.2023.104594. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0026286223001206

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

  18. [18]

    Eisenberg, Structural analysis of fluid flow in complex biological systems, MAIO 4 (1) (Feb

    R. Eisenberg, Structural analysis of fluid flow in complex biological systems, MAIO 4 (1) (Feb. 2023).doi:10.35119/maio.v4i1.126. URLhttps://www.maio-journal.com/index.php/MAIO/article/view/ 126

    work page doi:10.35119/maio.v4i1.126 2023

  19. [19]

    R. J. Shipley, A. F. Smith, P. W. Sweeney, A. R. Pries, T. W. Secomb, A hybrid discrete–continuum approach for modelling microcirculatory blood flow, Math. Med. Biol. (Mar. 2019).doi:10.1093/imammb/dqz006. URLhttps://academic.oup.com/imammb/advance-article/doi/10. 1093/imammb/dqz006/5393647

    work page doi:10.1093/imammb/dqz006 2019

  20. [20]

    U. N. A. Qohar, A. Z. Munthe-Kaas, J. M. Nordbotten, E. A. Hanson, A nonlinear multi-scale model for blood circulation in a realistic vascular sys- tem, R. Soc. open sci. 8 (12) (2021) 201949.doi:10.1098/rsos.201949. URLhttps://royalsocietypublishing.org/doi/10.1098/rsos. 201949

    work page doi:10.1098/rsos.201949 2021

  21. [21]

    Fritz, T

    M. Fritz, T. K¨ oppl, J. Oden, A. Wagner, B. Wohlmuth, C. Wu, A 1D–0D–3D coupled model for simulating blood flow and transport pro- cesses in breast tissue, Int. J. Numer. Meth. Biomed. Engng. 38 (7) (2022) e3612.doi:10.1002/cnm.3612. URLhttps://onlinelibrary.wiley.com/doi/10.1002/cnm.3612

    work page doi:10.1002/cnm.3612 2022

  22. [22]

    P. W. Sweeney, C. Walsh, S. Walker-Samuel, R. J. Shipley, A three- dimensional, discrete-continuum model of blood pressure in microvascu- lar networks, Int. J. Numer. Methods Biomed. Eng. 40 (8) (2024) e3832. doi:10.1002/cnm.3832. URLhttps://onlinelibrary.wiley.com/doi/10.1002/cnm.3832

    work page doi:10.1002/cnm.3832 2024

  23. [23]

    d’Esposito, P

    A. d’Esposito, P. W. Sweeney, M. Ali, M. Saleh, R. Ramasawmy, T. A. Roberts, G. Agliardi, A. Desjardins, M. F. Lythgoe, R. B. Pedley, R. Ship- ley, S. Walker-Samuel, Computational fluid dynamics with imaging of cleared tissue and of in vivo perfusion predicts drug uptake and treat- ment responses in tumours, Nat. Biomed. Eng. 2 (10) (2018) 773–787. doi:10...

    work page doi:10.1038/s41551-018-0306-y 2018

  24. [24]

    Z. Song, S. Xu, R. Eisenberg, H. Huang, A multi-domain model for microcirculation in optic nerve: Blood flow and oxygen transport, Physica D 467 (2024) 134272.doi:10.1016/j.physd.2024.134272. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0167278924002239 42

    work page doi:10.1016/j.physd.2024.134272 2024

  25. [25]

    X.-F. Wang, S. Nishi, M. Matsukawa, A. Ghigo, P.-Y. Lagr´ ee, J.-M. Fullana, Fluid friction and wall viscosity of the 1D blood flow model, J. Biomech. 49 (4) (2016) 565–571.doi:10.1016/j.jbiomech.2016.01.010. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0021929016000191

    work page doi:10.1016/j.jbiomech.2016.01.010 2016

  26. [26]

    W. Song, Q. Wei, W. Liu, T. Liu, J. Yi, N. Sheibani, A. A. Fawzi, R. A. Linsenmeier, S. Jiao, H. F. Zhang, A combined method to quantify the retinal metabolic rate of oxygen using photoacoustic ophthalmoscopy and optical coherence tomography, Sci. Rep. 4 (1) (2014) 6525.doi:10.1038/ srep06525. URLhttps://www.nature.com/articles/srep06525

    work page 2014

  27. [27]

    J. Liu, H. Suito, A multicore parallel algorithm for multiscale modelling of an entire human blood circulation network (2024).doi:10.48550/ARXIV. 2404.01680. URLhttps://arxiv.org/abs/2404.01680

    work page internal anchor Pith review doi:10.48550/arxiv 2024

  28. [28]

    L. O. M¨ uller, E. F. Toro, A global multiscale mathematical model for the human circulation with emphasis on the venous system: A GLOBAL MULTISCALE MATHEMATICAL MODEL FOR THE HUMAN CIRCU- LATION, Int. J. Numer. Meth. Biomed. Engng. 30 (7) (2014) 681–725. doi:10.1002/cnm.2622. URLhttps://onlinelibrary.wiley.com/doi/10.1002/cnm.2622

    work page doi:10.1002/cnm.2622 2014

  29. [29]

    A. R. Pries, T. W. Secomb, T. Gessner, M. B. Sperandio, J. F. Gross, P. Gaehtgens, Resistance to blood flow in microvessels in vivo., Circ. Res. 75 (5) (1994) 904–915.doi:10.1161/01.RES.75.5.904. URLhttps://www.ahajournals.org/doi/10.1161/01.RES.75.5.904

    work page doi:10.1161/01.res.75.5.904 1994

  30. [30]

    T. W. Secomb, A. R. Pries, Blood viscosity in microvessels: Experiment and theory, C. R. Phys. 14 (6) (2013) 470–478. doi:10.1016/j.crhy.2013.04.002. URLhttps://comptes-rendus.academie-sciences.fr/physique/ articles/10.1016/j.crhy.2013.04.002/

    work page doi:10.1016/j.crhy.2013.04.002 2013

  31. [31]

    Bappoo, L

    N. Bappoo, L. J. Kelsey, L. Parker, T. Crough, C. M. Moran, A. Thomson, M. C. Holmes, C. S. Wyrwoll, B. J. Doyle, Viscosity and haemodynamics in a late gestation rat feto-placental arterial network, Biomech. Model. Mechanobiol. 16 (4) (2017) 1361–1372.doi:10.1007/s10237-017-0892-8. URLhttp://link.springer.com/10.1007/s10237-017-0892-8

    work page doi:10.1007/s10237-017-0892-8 2017

  32. [32]

    V. B. Kolachalama, N. W. Bressloff, P. B. Nair, C. P. Shearman, Predictive Haemodynamics in a One-Dimensional Human Carotid Artery Bifurcation. Part I: Application to Stent Design, IEEE Trans. Biomed. Eng. 54 (5) (2007) 802–812.doi:10.1109/TBME.2006.889188. URLhttps://ieeexplore.ieee.org/document/4155000/ 43

    work page doi:10.1109/tbme.2006.889188 2007

  33. [33]

    M. S. Olufsen, C. S. Peskin, W. Y. Kim, E. M. Pedersen, A. Nadim, J. Larsen, Numerical Simulation and Experimental Validation of Blood Flow in Arteries with Structured-Tree Outflow Conditions, Ann. Biomed. Eng. 28 (11) (2000) 1281–1299.doi:10.1114/1.1326031. URLhttp://link.springer.com/10.1114/1.1326031

    work page doi:10.1114/1.1326031 2000

  34. [34]

    E. E. Brown, A. A. Guy, N. A. Holroyd, P. W. Sweeney, L. Gourmet, H. Coleman, C. Walsh, A. E. Markaki, R. Shipley, R. Rajendram, S. Walker-Samuel, Physics-informed deep generative learning for quan- titative assessment of the retina, Nat. Commun. 15 (1) (2024) 6859. doi:10.1038/s41467-024-50911-y. URLhttps://www.nature.com/articles/s41467-024-50911-y

    work page doi:10.1038/s41467-024-50911-y 2024

  35. [35]

    Zamir, Arterial Branching within the Confines of Fractal L- System Formalism, J

    M. Zamir, Arterial Branching within the Confines of Fractal L- System Formalism, J. Gen. Physiol. 118 (3) (2001) 267–276. doi:10.1085/jgp.118.3.267. URLhttps://rupress.org/jgp/article/118/3/267/42446/ Arterial-Branching-within-the-Confines-of-Fractal

    work page doi:10.1085/jgp.118.3.267 2001

  36. [36]

    IEEE Transactions on Medical Imaging44(10), 4049–4062 (2025).https://doi.org/10.1109/TMI

    J. Staal, M. D. Abramoff, M. Niemeijer, M. A. Viergever, B. van Gin- neken, Ridge-Based Vessel Segmentation in Color Images of the Retina, IEEE Trans. Med. Imaging 23 (4) (2004) 501–509.doi:10.1109/TMI. 2004.825627. URLhttp://ieeexplore.ieee.org/document/1282003/

    work page doi:10.1109/tmi 2004

  37. [37]

    T. I. J´ ozsa, R. M. Padmos, N. Samuels, W. K. El-Bouri, A. G. Hoekstra, S. J. Payne, A porous circulation model of the human brain forin silico clinical trials in ischaemic stroke, Interface Focus 11 (1) (2021) 20190127. doi:10.1098/rsfs.2019.0127. URLhttps://royalsocietypublishing.org/doi/10.1098/rsfs.2019. 0127

    work page doi:10.1098/rsfs.2019.0127 2021

  38. [38]

    Vidotto, T

    E. Vidotto, T. Koch, T. K¨ oppl, R. Helmig, B. Wohlmuth, Hybrid Models for Simulating Blood Flow in Microvascular Networks, Multiscale Model. Simul. 17 (3) (2019) 1076–1102.doi:10.1137/18M1228712. URLhttps://epubs.siam.org/doi/10.1137/18M1228712

    work page doi:10.1137/18m1228712 2019

  39. [39]

    S. Xu, B. Eisenberg, Z. Song, H. Huang, Osmosis through a Semi-permeable Membrane: a Consistent Approach to Interactions, version Number: 3 (2018).doi:10.48550/ARXIV.1806.00646. URLhttps://arxiv.org/abs/1806.00646

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1806.00646 2018

  40. [40]

    J. R. Levick, C. C. Michel, Microvascular fluid exchange and the revised Starling principle, Cardiovasc. Res. 87 (2) (2010) 198–210. doi:10.1093/cvr/cvq062. URLhttps://academic.oup.com/cardiovascres/article-lookup/ doi/10.1093/cvr/cvq062 44

    work page doi:10.1093/cvr/cvq062 2010

  41. [41]

    K. J. Kim, J. D. Jackson, Proof that the Neumann Green’s function in electrostatics can be symmetrized, Am. J. Phys. 61 (12) (1993) 1144–1146. doi:10.1119/1.17311. URLhttps://pubs.aip.org/ajp/article/61/12/1144/1054244/ Proof-that-the-Neumann-Green-s-function-in

    work page doi:10.1119/1.17311 1993

  42. [42]

    K. F. Riley, M. P. Hobson, S. J. Bence, Mathematical Methods for Physics and Engineering, 3rd Edition, Cambridge University Press, Cambridge, 2006.doi:10.1017/CBO9780511810763. URLhttps://www.cambridge.org/highereducation/books/ mathematical-methods-for-physics-and-engineering/ FC466374D5B94E86D969100070CA6483

    work page doi:10.1017/cbo9780511810763 2006

  43. [43]

    A new type of sharp bounds for ratios of modified Bessel functions

    D. Ruiz-Antol´ ın, J. Segura, A new type of sharp bounds for ratios of modified Bessel functions, J. Math. Anal. Appl. 443 (2) (2016) 1232–1246. doi:10.1016/j.jmaa.2016.06.011. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0022247X16302402

    work page doi:10.1016/j.jmaa.2016.06.011 2016

  44. [44]

    Y. He, H. Liu, R. Himeno, A one-dimensional thermo-fluid model of blood circulation in the human upper limb, Int. J. Heat Mass Transf. 47 (12-13) (2004) 2735–2745.doi:10.1016/j.ijheatmasstransfer.2003.10.041. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S001793100400002X

    work page doi:10.1016/j.ijheatmasstransfer.2003.10.041 2004

  45. [45]

    K¨ oppl, R

    T. K¨ oppl, R. Helmig, Dimension Reduced Modeling of Blood Flow in Large Arteries: An Introduction for Master Students and First Year Doctoral Students, Mathematical Engineering, Springer International Publishing, Cham, 2023.doi:10.1007/978-3-031-33087-2. URLhttps://link.springer.com/10.1007/978-3-031-33087-2

    work page doi:10.1007/978-3-031-33087-2 2023

  46. [46]

    Guidoboni, A

    G. Guidoboni, A. Harris, L. Carichino, Y. Arieli, B. A. Siesky, Effect of intraocular pressure on the hemodynamics of the central retinal artery: A mathematical model, Mathematical Biosciences and Engineering 11 (3) (2014) 523–546.doi:10.3934/mbe.2014.11.523. URLhttp://www.aimspress.com/article/10.3934/mbe.2014.11.523

    work page doi:10.3934/mbe.2014.11.523 2014

  47. [47]

    Ertop, Y

    M. Ertop, Y. T. Fulya, Y. K. Arıba¸ s, H. B. ¨Ozdemir, G. G¨ urelik, M. E. Atılgan, S ¸. ¨Ozdek, Evaluation of Central and Peripheral Reti- nal Vascular Changes in the Fellow Eyes of Patients with Unilateral Retinal Vein Occlusions, Turk. J. Ophthalmol. 53 (6) (2023) 349–355. doi:10.4274/tjo.galenos.2023.94389. URLhttps://oftalmoloji.org/articles/doi/tjo....

    work page doi:10.4274/tjo.galenos.2023.94389 2023

  48. [48]

    G. T. Dorner, E. Polska, G. Garh¨ ofer, C. Zawinka, B. Frank, L. Schmetterer, Calculation of the diameter of the central retinal artery from noninvasive measurements in humans, Curr. Eye Res. 25 (6) 45 (2002) 341–345.doi:10.1076/ceyr.25.6.341.14231. URLhttp://www.tandfonline.com/doi/full/10.1076/ceyr.25.6. 341.14231

    work page doi:10.1076/ceyr.25.6.341.14231 2002

  49. [49]

    R. J. Hernandez, S. Madhusudhan, Y. Zheng, W. K. El-Bouri, Link- ing Vascular Structure and Function: Image-Based Virtual Popula- tions of the Retina, Invest. Ophthalmol. Vis. Sci. 65 (4) (2024) 40. doi:10.1167/iovs.65.4.40. URLhttps://iovs.arvojournals.org/article.aspx?articleid= 2793616

    work page doi:10.1167/iovs.65.4.40 2024

  50. [50]

    C. E. Riva, J. E. Grunwald, S. H. Sinclair, B. L. Petrig, Blood velocity and volumetric flow rate in human retinal vessels, Invest. Ophthalmol. Vis. Sci. 26 (8) (1985) 1124–1132

    work page 1985

  51. [51]

    Abramowitz, I

    M. Abramowitz, I. A. Stegun (Eds.), Handbook of mathematical functions: with formulas, graphs, and mathematical tables, 9th Edition, Dover books on mathematics, Dover Publ, New York, NY, 2013

    work page 2013

  52. [52]

    R. M. Werkmeister, N. Dragostinoff, S. Palkovits, R. Told, A. Boltz, R. A. Leitgeb, M. Gr¨ oschl, G. Garh¨ ofer, L. Schmetterer, Measurement of Absolute Blood Flow Velocity and Blood Flow in the Human Retina by Dual-Beam Bidirectional Doppler Fourier-Domain Optical Coher- ence Tomography, Invest. Ophthalmol. Vis. Sci. 53 (10) (2012) 6062. doi:10.1167/iovs...

    work page doi:10.1167/iovs.12-9514 2012

  53. [53]

    Kristen, L

    A. Kristen, L. Kelsey, E. Wintermantel, B. Doyle, Fundus Image Based Blood Flow Simulation of the Retinal Arteries, in: G. R. Joldes, B. Doyle, A. Wittek, P. M. Nielsen, K. Miller (Eds.), Computational Biomechanics for Medicine, Springer International Publishing, Cham, 2016, pp. 143–154. doi:10.1007/978-3-319-28329-6_13. URLhttp://link.springer.com/10.100...

    work page doi:10.1007/978-3-319-28329-6_13 2016

  54. [54]

    I. Fatt, K. Shantinath, Flow conductivity of retina and its role in retinal adhesion, Exp. Eye Res. 12 (2) (1971) 218–226. doi:10.1016/0014-4835(71)90094-7. URLhttps://linkinghub.elsevier.com/retrieve/pii/ 0014483571900947

    work page doi:10.1016/0014-4835(71)90094-7 1971

  55. [55]

    R. J. Antcliff, Hydraulic Conductivity of Fixed Retinal Tissue After Se- quential Excimer Laser Ablation: Barriers Limiting Fluid Distribution and Implications for Cystoid Macular Edema, Arch. Ophthalmol. 119 (4) (2001) 539.doi:10.1001/archopht.119.4.539. URLhttp://archopht.jamanetwork.com/article.aspx?doi=10.1001/ archopht.119.4.539 46

    work page doi:10.1001/archopht.119.4.539 2001

  56. [56]

    Y. Zhu, S. Xu, R. S. Eisenberg, H. Huang, A Bidomain Model for Lens Microcirculation, Biophysical Journal 116 (6) (2019) 1171–1184. doi:10.1016/j.bpj.2019.02.007. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0006349519301341

    work page doi:10.1016/j.bpj.2019.02.007 2019

  57. [57]

    Y. Zhu, S. Xu, R. S. Eisenberg, H. Huang, Membranes in Optic Nerve Mod- els, arXiv:2105.14411 [math] (Jun. 2021).doi:10.48550/arXiv.2105. 14411. URLhttp://arxiv.org/abs/2105.14411

    work page doi:10.48550/arxiv.2105 2021

  58. [58]

    S. Xiao, H. Huang, R. Eisenberg, Z. Song, S. Xu, Potassium Clearance in Optic Nerve: A Multidomain Model, Front. Biosci. (Landmark Ed) 30 (7) (2025) 39722.doi:10.31083/FBL39722. URLhttps://www.imrpress.com/journal/FBL/30/7/10.31083/ FBL39722

    work page doi:10.31083/fbl39722 2025

  59. [59]

    S. Xiao, H. Huang, R. Eisenberg, Z. Song, S. Xu, Glymphatic Clearance in the Optic Nerve: A Multidomain Electro-osmostic Model, version Number: 1 (2025).doi:10.48550/ARXIV.2510.22271. URLhttps://arxiv.org/abs/2510.22271

    work page doi:10.48550/arxiv.2510.22271 2025

  60. [60]

    Y. Zhu, S. Xu, R. S. Eisenberg, H. Huang, A tridomain model for potassium clearance in optic nerve of Necturus, Biophys. J. 120 (15) (2021) 3008–3027.doi:10.1016/j.bpj.2021.06.020. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0006349521005075

    work page doi:10.1016/j.bpj.2021.06.020 2021

  61. [61]

    Y. Zhu, S. Xu, R. S. Eisenberg, H. Huang, Optic nerve microcirculation: Fluid flow and electrodiffusion, Phys. Fluids 33 (4) (2021) 041906. doi:10.1063/5.0046323. URLhttps://pubs.aip.org/pof/article/33/4/041906/994351/ Optic-nerve-microcirculation-Fluid-flow-and

    work page doi:10.1063/5.0046323 2021

  62. [62]

    Baricz, D

    ´A. Baricz, D. J. Maˇ sirevi´ c, S. Ponnusamy, S. Singh, Bounds for the product of modified Bessel functions, Aequat. Math. 90 (4) (2016) 859–870.doi: 10.1007/s00010-016-0414-2. URLhttp://link.springer.com/10.1007/s00010-016-0414-2

    work page doi:10.1007/s00010-016-0414-2 2016

  63. [63]

    Bounds for modified Bessel functions of the first and second kinds

    ´A. Baricz, Bounds for modified Bessel functions of the first and second kinds, Proc. Edinb. Math. Soc. 53 (3) (2010) 575–599. doi:10.1017/S0013091508001016. URLhttps://www.cambridge.org/core/product/identifier/ S0013091508001016/type/journal_article

    work page doi:10.1017/s0013091508001016 2010

  64. [64]

    C. M. Joshi, S. K. Bissu, Inequalities for some special functions, J. Comput. Appl. Math. 69 (2) (1996) 251–259.doi:10.1016/0377-0427(95)00042-9. URLhttps://linkinghub.elsevier.com/retrieve/pii/ 0377042795000429 47

    work page doi:10.1016/0377-0427(95)00042-9 1996