Pith. sign in

REVIEW 2 major objections 5 minor 52 references

A hierarchy of blood vessel models, Part II: 3D-3D to 3D-1D and 1D

T0 review · 2 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Pith's one-line read The pressure fields of the full 3D-3D perfusion model converge to the reduced 3D-1D model at rate C ε^{1/6}|log ε| |p0| as ε → 0, all the way to the tapering free tip.

desk verdict The 3D-3D-to-1D convergence program is novel and the asymptotic analysis is serious, but the proof of the 3D-3D existence theorem rests on a non-H^1 vector field and the main results currently do not close. read the letter →

arxiv 2507.13330 v1 pith:TNAM74HD submitted 2025-07-17 math.AP physics.bio-phphysics.flu-dynq-bio.TO

classification math.APphysics.bio-phphysics.flu-dynq-bio.TO MSC 35B2576D0776S0576Z05
keywords bloodperfusionmodellingDarcy-StokesflowDarcy-Poiseuillereduction1DGreen'sfunctionmodelslenderbodyapproximationdegeneratefreetipconvergencerateporousmedium
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper is the second half of an argument that three standard descriptions of blood perfusion through tissue around a thin vessel — a full three-dimensional Darcy–Stokes system, a reduced model with a one-dimensional vessel coupled to a three-dimensional porous tissue, and a one-dimensional Green's-function model — describe the same physics in a quantitatively controlled way. Its central result is a convergence theorem: for a vessel whose radius decays to zero at the free tip, the pressures computed by the most detailed 3D-3D model differ from those of the reduced 3D-1D model by at most $C \epsilon^{1/6}|\log\epsilon|\,|p_0|$, where $\epsilon$ is the maximum vessel radius and $p_0$ is the incoming pressure. The proof is indirect: it constructs an explicit approximate velocity field from the 1D model's pressure, proves the 3D-3D system converges directly to that 1D solution at the $\epsilon^{1/6}|\log\epsilon|$ rate, and then borrows the 3D-1D-to-1D convergence of Part I, together with the 1D pressure bounds proved there, to close the remaining gap. If correct, this gives the much cheaper 1D and 3D-1D models used in microvascular simulation a rigorous error guarantee that extends up to the degenerate endpoint where the vessel becomes indistinguishable from a capillary.

What carries the argument

The load-bearing object is the 1D 'slender-body' pressure $p_{\mathrm{SB}}(s)$, the solution of the integrodifferential equation (1.21), in which the exterior pressure is expressed through the Neumann Green's function of the half-space as $q_{\mathrm{SB}} = \tfrac{\pi}{8\zeta\mu} S_N\big[\tfrac{d}{ds}(a^4 \tfrac{dp_{\mathrm{SB}}}{ds})\big]$. From the companion Part I the paper imports three ingredients: weighted $L^2$ and $L^\infty$ bounds for $p_{\mathrm{SB}}$ whose endpoint weights gain one or two powers of the radius $a(s)$ near the tip (Lemma 5.1), a near-$\theta$-independence estimate for $q_{\mathrm{SB}}$ on the vessel surface (Lemma 5.2), and residual bounds showing that $(p_{\mathrm{SB}}, q_{\mathrm{SB}})$ almost satisfies the 3D-1D weak form (Lemma 5.3). The comparison object is the velocity ansatz $U$ of (5.17), built from the variable-radius Poiseuille law and multiplied by a cutoff $\phi_\epsilon(s)$ that vanishes on the final tip segment of length $\epsilon^{4/3}$; the cutoff length is the optimized balance between the ansatz failing near the curved tip and the geometric smallness of the tapering end, and that balance produces the $\epsilon^{1/6}|\log\epsilon|$ rate. Supporting this core, Appendix A supplies $\epsilon$-scaled Korn, pressure-operator, and trace inequalities for the slender vessel with a degenerate free end.

What would settle it

Solve the 1D integrodifferential equation (1.21) numerically for a straight vessel with a spheroidal tip and check the weighted $L^\infty$ bounds of Lemma 5.1 directly: if $\|p_{\mathrm{SB}}\|_{L^\infty}$ or $\|a\,p_{\mathrm{SB},s}\|_{L^\infty}$ grows faster than $\epsilon^{-1/2}|p_0|$ as $\epsilon \to 0$, the tip-cutoff estimates of Section 5 fail and Theorem 1.4 cannot hold. A complementary check is to solve the 3D-3D Darcy–Stokes and 3D-1D Darcy–Poiseuille systems in the same geometry and confirm that $(1/|V_\epsilon|^{1/2})\|p^\epsilon-p\|_{L^2(V_\epsilon)} + \|q^\epsilon-q\|_{D^{1,2}(\Omega_\epsilon)}$ decays like $\epsilon^{1/6}|\log\epsilon|$ rather than saturating at a slower rate.

Watch

Extended reading notes

Core claim

The central claim, Theorem 1.4, is that for sufficiently small $\epsilon$ the pressure difference between the 3D-3D Darcy–Stokes solution $(u^\epsilon, p^\epsilon, q^\epsilon)$ and the 3D-1D Darcy–Poiseuille solution $(p,q)$ satisfies $(1/|V_\epsilon|^{1/2})\|p^\epsilon-p\|_{L^2(V_\epsilon)} + \|q^\epsilon-q\|_{D^{1,2}(\Omega_\epsilon)} \le C \epsilon^{1/6}|\log\epsilon|\,|p_0|$, with $C$ independent of $\epsilon$, and the estimate holds up to the free tip $s=1$ where the vessel radius decays spheroidally to zero and the reduced models become degenerate. Because the angle-averaged Robin condition in the 3D-1D model does not easily control the exterior pressure pointwise on the vessel surface, the proof detours through the 1D slender-body model: a velocity ansatz built from the variable-radius Poiseuille law is cut off in a tip region of length $\epsilon^{4/3}$, and the 3D-3D system is shown to converge to this explicit 1D solution (Theorem 1.6) at the same rate; the 3D-3D-to-3D-1D result then follows from the Part I estimate of Theorem 1.5. The $\epsilon^{1/6}$ rate is entirely a consequence of the free-end behavior; away from the degenerate endpoint the limiting factor would instead be the faster Part I rate $\epsilon^{1/2}|\log\epsilon|$.

Load-bearing premise

The convergence proof stands on endpoint-weighted estimates for the one-dimensional model's pressure near the vessel tip that are proved in the companion Part I paper and assumed rather than reproved here; if those estimates fail, the $\epsilon^{1/6}|\log\epsilon|$ rate does not close.

Editorial extensions

If this is right

  • The full 3D-3D Darcy–Stokes description of perfusion can be replaced, with a certified error $C\epsilon^{1/6}|\log\epsilon|\,|p_0|$, by the reduced 3D-1D Darcy–Poiseuille system for a thin vessel whose radius decays at the free tip.
  • Through Part I's Theorem 1.5, the same guarantee descends to the fully 1D Green's-function model, closing the three-level hierarchy at a single rate.
  • The estimate is valid all the way to the endpoint $s=1$, where the vessel radius vanishes and neither reduced model admits a classical boundary value, so the degenerate capillary-scale tip does not spoil the approximation.
  • The vessel-interior velocity field is also approximated by the explicit Poiseuille-type ansatz $U$ at the same $\epsilon^{1/6}|\log\epsilon|$ rate, in the sense of Theorem 1.6.
  • The $\epsilon^{1/6}$ rate is imposed entirely by the free end; refined asymptotics for the tip boundary layer would give more accurate but more complicated reduced models.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the paper's route goes 3D-3D to 1D and then re-uses the 3D-1D-to-1D result, the convergence of the full model hinges far more on the quality of the 1D pressure bounds (Lemma 5.1) than on the details of the 3D-1D coupling; this suggests the 3D-1D model's role in the hierarchy is largely motivational and structural.
  • The $\epsilon^{1/6}$ exponent is probably not sharp: a boundary-layer ansatz that correctly describes the flow near the spheroidal tip should push the rate toward the faster $\epsilon^{1/2}|\log\epsilon|$ of Part I, since away from the tip the present proof is limited only by that faster rate.
  • The half-space geometry is used purely for the explicit Neumann Green's function; for a bounded domain with a vessel reconnecting to the boundary at both ends, one would expect an analogous and faster convergence with no degenerate tip, a setting the paper leaves as future work.
  • A numerical check of Lemma 5.1 on the 1D equation (1.21) would independently test the assumed Part I input, since the weighted $L^\infty$ bounds can be computed without ever solving the 3D systems.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The paper proposes a hierarchy of three blood-perfusion models and proves convergence between the 3D-3D Darcy--Stokes model and the 3D-1D Darcy--Poiseuille model, as well as to a 1D Green's-function model. The main results are well-posedness for the 3D-3D system (Theorem 1.2), well-posedness for the 3D-1D system (Theorem 1.3), a 3D-3D to 1D convergence estimate (Theorem 1.6), and, by chaining Theorem 1.6 with the companion result Theorem 1.5 from Part I [42], a 3D-3D to 3D-1D error estimate at rate O(epsilon^{1/6}|log epsilon|) (Theorem 1.4). The proof uses a slender-body velocity ansatz built from the 1D pressure, a cutoff isolating the degenerate tip, and a series of epsilon-dependent Korn, pressure-operator, and trace estimates proved in the appendices.

Significance. If the results are correct and the dependency on Part I is resolved, the paper would provide a rare quantitative link among three widely used modeling levels of tissue perfusion, with an explicit convergence rate that is valid up to the degenerate vessel tip. The construction of the velocity ansatz, the optimized epsilon^{4/3} tip cutoff, and the careful epsilon tracking in Sections 5 and Appendix A are genuine strengths, as is the explicit treatment of the free-tip geometry. However, the main theorems are conditional on unproved estimates from the companion paper [42], and one load-bearing existence proof in Section 2 contains a gap that must be repaired.

major comments (2)
  1. [§2.2, Eqs. (2.17)–(2.19)] The surjectivity proof for the divergence operator on V_sigma^perp is not valid as written. The constructed field v_1 has e_theta-component proportional to c r theta (Eq. (2.17)); on the full tubular domain V_epsilon, the angle theta is a periodic coordinate and is not single-valued, so any representative of v_1 has a jump discontinuity across a branch cut. Hence v_1 is not in H^1(V_epsilon), the identity div v_1 = c fails in the weak sense (a surface distribution on the cut appears), and the decomposition in (2.18) does not produce the claimed zero-mean divergence. Consequently Corollary 2.3, the inf-sup condition, and the recovery of p_epsilon in Theorem 1.2 are not established by the given argument. This gap is local and likely repairable by a different construction, but the proof must be rewritten.
  2. [§1.4, §5.1, §5.3] The central results are conditional on the companion paper [42]. Theorem 1.4 is obtained by chaining Theorem 1.6 with Theorem 1.5 [42], and even the proof of Theorem 1.6 invokes, without proof in this manuscript, Lemma 5.1 (weighted L^2 and L^infinity bounds (5.2)–(5.3)), Lemma 5.2 (near-theta-independence (5.4)), and Lemma 5.3 (residual estimates (5.7)–(5.8)). These imported bounds enter at load-bearing points, for example in (5.21), (5.38), and the E_9 estimate (5.58). If Part I is not independently available in refereed form or included as an appendix, the manuscript does not by itself prove Theorem 1.4. The authors should either supply the needed proofs or explicitly state that Theorems 1.4 and 1.6 are conditional on [42].
minor comments (5)
  1. [§1.4, final paragraph] The phrase 'the 3D-3D to 1D convergence result of Theorem 1.5 in section 5' should refer to Theorem 1.6, not Theorem 1.5.
  2. [§5.3, estimate (5.58)] The inequality 'a^{-1}(s) >= C epsilon^{-2/3}' should read 'a^{-1}(s) <= C epsilon^{-2/3}' (equivalently, a(s) >= C epsilon^{2/3}); the displayed direction is reversed and, as written, does not yield the stated epsilon^{2/3} bound.
  3. [§2.2, Eq. (2.22)] There is a typographical error in the norm on the last line: the expression 'V x D^{1,2}(Omega_epsilon ||(w,w)||...' is missing a closing parenthesis and should read '||(u,q)||_{V x D^{1,2}(Omega_epsilon)} ||(w,w)||_{V x D^{1,2}(Omega_epsilon)}'.
  4. [§2.2, after Eq. (2.23)] In the displayed coercivity computation, '+(zeta_epsilon nabla q, nabla q)' contains a doubled plus sign; it should be a single plus.
  5. [§1.3, after Theorem 1.2] The sentence 'The proof of Theorem 1.2 appears in section 2, including justification for the choice of epsilon-scaling of the coefficients zeta epsilon^4 and kappa epsilon^3 appearing in (1.11c)' should refer to (1.10c), not (1.11c), since (1.11c) belongs to the 3D-1D system.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the 3D-3D-to-1D estimate is proven by residual comparison against the independent 1D model, and the imported Part I results are parameter-free inputs, not fitted/predicted outputs.

full rationale

The paper's central estimate (Theorem 1.4) is obtained by transitivity: Theorem 1.6 compares the 3D-3D solution to the slender-body/1D pair (U, pSB, qSB) defined by (1.20)-(1.21), and Theorem 1.5 from Part I [42] compares that same 1D pair to the 3D-1D solution. The 1D model is not fitted to the 3D-3D solution; pSB and qSB are defined by an independent integrodifferential equation, and the proof of Theorem 1.6 inserts (U, pSB, qSB) into the 3D-3D weak form and bounds the residuals. The residual estimates use weighted L∞ and near-theta-independence bounds (Lemmas 5.1-5.2) and residual bounds (Lemma 5.3) imported from Part I. These are self-citations by the same authors, but they are load-bearing only in the sense of being unproved in this submission; they are parameter-free statements whose assumptions (admissible radius, p0 data) do not include the target 3D-3D-to-3D-1D convergence, so under the rubric they count as independent support rather than circularity. The lack of a reproduced proof of Part I is a completeness/reproducibility concern, not a circular reduction. The flagged Lemma 2.2 objection (the vector field v1 contains a multi-valued angle theta and may fail to lie in H1(Vepsilon)) is a potential correctness gap in the existence proof, not an instance of a conclusion reducing to its inputs, so it does not affect the circularity score.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The central proof rests primarily on modeling assumptions about coefficient scalings and the degenerate endpoint geometry, plus imported estimates from the companion paper Part I. No new physical entities are introduced; the only hand-chosen proof parameter is the tip cutoff exponent.

free parameters (1)
  • Tip cutoff exponent 4/3 = 4/3 (length scale of the tip region epsilon^{4/3})
    Chosen by hand in (1.24)/(5.18) to balance error terms near the free end; this choice determines the final O(epsilon^{1/6}|log epsilon|) rate. It is a proof construction parameter, not a physical constant.
assumptions (5)
  • domain assumption The coupling coefficients in the 3D-3D system scale as zeta^epsilon = zeta epsilon^4 and kappa^epsilon = kappa epsilon^3 with zeta, kappa, mu of order one.
    Introduced in (1.10) and justified only by energy considerations in Section 2.4; the convergence result depends on this precise scaling.
  • domain assumption The vessel has one end embedded in the plane z=0, one free end, and radius function a(s) with spheroidal decay (1.4) up to error C epsilon^2 sqrt(1-s^2).
    The upper half-space and degenerate free end are needed for the explicit Green's function and for the endpoint trace estimates; Section 1.1, Definition 1.1.
  • domain assumption The exterior tissue is a porous medium governed by Darcy's law with hydraulic conductivity zeta epsilon^4, and the vessel wall is semi-permeable with no-slip tangential velocity.
    These are the modeling choices in (1.10a)-(1.10f), stated in Section 1.2 and Section 3.
  • ad hoc to paper The 1D slender-body solution pSB satisfies the weighted bounds of Lemma 5.1, near-theta-independence of Lemma 5.2, and the weak-form residuals of Lemma 5.3, as proved in Part I [42].
    These imported estimates are used throughout Section 5 and are not proved in this submission.
  • standard math Standard PDE tools: 3D Sobolev inequality in D1,2(R^3_+), Korn inequality in slender domains, Bogovskii-type pressure operators, and epsilon-dependent trace inequalities.
    Used in Sections 2 and 5; the epsilon-dependent variants are proved in Appendix A, relying on standard results [15,18,30].

how reviews work

0 comments
Cite this review

Pith. "Pith review of A hierarchy of blood vessel models, Part II: 3D-3D to 3D-1D and 1D." pith.science (2026). https://pith.science/paper/TNAM74HD

@misc{pith2026250713330,
  author       = {Pith},
  title        = {Pith review of: A hierarchy of blood vessel models, Part II: 3D-3D to 3D-1D and 1D},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TNAM74HD}},
  note         = {Machine review of arXiv:2507.13330}
}
abstract

We propose and analyze a hierarchy of three models of blood perfusion through a tissue surrounding a thin arteriole or venule. Our goal is to rigorously link 3D-3D Darcy--Stokes, 3D-1D Darcy--Poiseuille, and 1D Green's function methods commonly used to model this process. Here in Part II, we consider the most detailed level, a 3D-3D Darcy-Stokes system coupled across the permeable vessel surface by mass conservation and pressure/stress balance conditions. We derive a convergence result between the 3D-3D model and both the 3D-1D Darcy--Poiseuille model and 1D Green's function model proposed in Part I [Ohm \& Strikwerda, arXiv preprint July 2025] at a rate proportional to $\epsilon^{1/6}|\log\epsilon|$, where $\epsilon$ is the maximum vessel radius. The rate is limited by the inclusion of a degenerate endpoint where the vessel radius vanishes, i.e. becomes indistinguishable from a capillary. Key to our proof are \emph{a priori} estimates for the 1D integrodifferential model obtained in Part I.

Figures

Figures reproduced from arXiv: 2507.13330 by the authors.

Figure 1
Figure 1. An example geometry for the blood vessel Vϵ. 1.2. The models. We begin with the 3D Darcy–3D Stokes description of perfusion outside of the vessel and flow within. Within the vessel, we consider velocity and pressure fields (u ϵ , pϵ ) : Vϵ × Vϵ → R 3 × R. Letting E(v) = 1 2 (∇v + ∇v T) denote the symmetric gradient, we also define the stress tensor Σϵ = 2µ E(u ϵ ) − p ϵ I, where µ is the blood viscosity and is indep… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages

  1. [42]

    Ohm and S

    L. Ohm and S. Strikwerda. A hierarchy of blood vessel models, part I: 3D-1D to 1D. arXiv preprint, 2025

  2. [1]

    Berrone, C

    S. Berrone, C. Giverso, D. Grappein, L. Preziosi, and S. Scial` o. An optimization based 3d-1d coupling strategy for tissue perfusion and chemical transport during tumor-induced angiogenesis. Computers & Mathematics with Applications, 151:252–270, 2023

  3. [2]

    R. L. Bishop. There is more than one way to frame a curve. The American Mathematical Monthly, 82(3):246–251, 1975

  4. [3]

    Blanco, R

    P. Blanco, R. Feij´ oo, and S. Urquiza. A unified variational approach for coupling 3d–1d models and its blood flow applications. Computer Methods in Applied Mechanics and Engineering , 196(41-44):4391–4410, 2007

  5. [4]

    P. J. Blanco, M. Pivello, S. Urquiza, and R. A. Feij´ oo. On the potentialities of 3d–1d coupled models in hemo- dynamics simulations. Journal of biomechanics , 42(7):919–930, 2009

  6. [5]

    Bogovskiı

    M. Bogovskiı. Solutions of some problems of vector analysis, associated with the operators div and grad. Theory of cubature formulas and the application of functional analysis to problems of mathematical physics , 1980:5–40, 1980

  7. [6]

    Boulakia, C

    M. Boulakia, C. Grandmont, F. Lespagnol, and P. Zunino. Mathematical and numerical analysis of reduced order interface conditions and augmented finite elements for mixed dimensional problems. Computers & Mathematics with Applications, 175:536–569, 2024

  8. [7]

    Canic and A

    S. Canic and A. Mikelic. Effective equations modeling the flow of a viscous incompressible fluid through a long elastic tube arising in the study of blood flow through small arteries. SIAM Journal on Applied Dynamical Systems, 2(3):431–463, 2003

Show all 52 references
  1. [8]

    Cassot, F

    F. Cassot, F. Lauwers, S. Lorthois, P. Puwanarajah, V. Cances-Lauwers, and H. Duvernoy. Branching patterns for arterioles and venules of the human cerebral cortex. Brain research, 1313:62–78, 2010

  2. [9]

    Castineira, E

    G. Castineira, E. Maruˇ si´ c-Paloka, I. Paˇ zanin, and J. M. Rodriguez. Rigorous justification of the asymptotic model describing a curved-pipe flow in a time-dependent domain, 2019

  3. [10]

    D’Angelo

    C. D’Angelo. Multiscale modelling of metabolism and transport phenomena in living tissues . PhD thesis, EPFL, 2007

  4. [11]

    D’Angelo

    C. D’Angelo. Finite element approximation of elliptic problems with dirac measure terms in weighted spaces: applications to one-and three-dimensional coupled problems. SIAM Journal on Numerical Analysis , 50(1):194– 215, 2012

  5. [12]

    D’angelo and A

    C. D’angelo and A. Quarteroni. On the coupling of 1d and 3d diffusion-reaction equations: application to tissue perfusion problems. Mathematical Models and Methods in Applied Sciences , 18(08):1481–1504, 2008. 40 LAUREL OHM AND SARAH STRIKWERDA

  6. [13]

    D’Angelo and P

    C. D’Angelo and P. Zunino. Robust numerical approximation of coupled stokes’ and darcy’s flows applied to vascular hemodynamics and biochemical transport. ESAIM: Mathematical Modelling and Numerical Analysis , 45(3):447–476, 2011

  7. [14]

    Discacciati, A

    M. Discacciati, A. Quarteroni, et al. Navier-stokes/darcy coupling: modeling, analysis, and numerical approxi- mation. Rev. Mat. Complut , 22(2):315–426, 2009

  8. [15]

    Fabricius, E

    J. Fabricius, E. Miroshnikova, A. Tsandzana, and P. Wall. Pressure-driven flow in thin domains. Asymptotic Analysis, 116(1):1–26, 2019

  9. [16]

    Formaggia, A

    L. Formaggia, A. Moura, and F. Nobile. Coupling 3d and 1d fluid-structure interaction models for blood flow simulations. In PAMM: Proceedings in Applied Mathematics and Mechanics, volume 6, pages 27–30. Wiley Online Library, 2006

  10. [17]

    Formaggia, A

    L. Formaggia, A. Moura, and F. Nobile. On the stability of the coupling of 3d and 1d fluid-structure interaction models for blood flow simulations. ESAIM: Mathematical Modelling and Numerical Analysis , 41(4):743–769, 2007

  11. [18]

    G. P. Galdi. An introduction to the mathematical theory of the Navier-Stokes equations: Steady-state problems . Springer Science & Business Media, 2011

  12. [19]

    Ghosh, V

    A. Ghosh, V. Kozlov, and S. Nazarov. Modified reynolds equation for steady flow through a curved pipe. Journal of Mathematical Fluid Mechanics , 23:1–22, 2021

  13. [20]

    Heltai and P

    L. Heltai and P. Zunino. Reduced lagrange multiplier approach for non-matching coupling of mixed-dimensional domains. Mathematical Models and Methods in Applied Sciences , 33(12):2425–2462, 2023

  14. [21]

    Hsu and T

    R. Hsu and T. W. Secomb. A green’s function method for analysis of oxygen delivery to tissue by microvascular networks. Mathematical biosciences, 96(1):61–78, 1989

  15. [22]

    E. R. Hyde, C. Michler, J. Lee, A. N. Cookson, R. Chabiniok, D. A. Nordsletten, and N. P. Smith. Parameterisa- tion of multi-scale continuum perfusion models from discrete vascular networks. Medical & biological engineering & computing, 51:557–570, 2013

  16. [23]

    S. C. Johnston, D. K. Wallace, S. F. Freedman, T. L. Yanovitch, and Z. Zhao. Tortuosity of arterioles and venules in quantifying plus disease.Journal of American Association for Pediatric Ophthalmology and Strabismus, 13(2):181–185, 2009

  17. [24]

    T. Koch, M. Schneider, R. Helmig, and P. Jenny. Modeling tissue perfusion in terms of 1d-3d embedded mixed- dimension coupled problems with distributed sources. Journal of Computational Physics , 410:109370, 2020

  18. [25]

    K¨ oppl, E

    T. K¨ oppl, E. Vidotto, and B. Wohlmuth. A 3d-1d coupled blood flow and oxygen transport model to generate microvascular networks. International journal for numerical methods in biomedical engineering , 36(10):e3386, 2020

  19. [26]

    Kuchta, F

    M. Kuchta, F. Laurino, K.-A. Mardal, and P. Zunino. Analysis and approximation of mixed-dimensional pdes on 3d-1d domains coupled with lagrange multipliers. SIAM Journal on Numerical Analysis , 59(1):558–582, 2021

  20. [27]

    Laurino and P

    F. Laurino and P. Zunino. Derivation and analysis of coupled pdes on manifolds with high dimensionality gap arising from topological model reduction. ESAIM: Mathematical Modelling and Numerical Analysis , 53(6):2047– 2080, 2019

  21. [28]

    Lorthois, F

    S. Lorthois, F. Lauwers, and F. Cassot. Tortuosity and other vessel attributes for arterioles and venules of the human cerebral cortex. Microvascular research, 91:99–109, 2014

  22. [29]

    S. R. Manjate. Mathematical modelling of pressure-driven flow in thin domains.Doctoral thesis, Luleøa University of Technology, 2024

  23. [30]

    Maruˇ si´ c and E

    S. Maruˇ si´ c and E. Maruˇ si´ c-Paloka. Two-scale convergence for thin domains and its applications to some lower- dimensional models in fluid mechanics. Asymptotic Analysis, 23(1):23–57, 2000

  24. [31]

    Maruˇ si´ c-Paloka

    E. Maruˇ si´ c-Paloka. The effects of flexion and torsion on a fluid flow through a curved pipe.Applied mathematics and optimization , 44:245–272, 2001

  25. [32]

    Mironescu

    P. Mironescu. The role of the hardy type inequalities in the theory of function spaces. Revue roumaine de math´ ematiques pures et appliqu´ ees, 63(4):447–525, 2018

  26. [33]

    Miroshnikova

    E. Miroshnikova. Pressure-driven flow in a thin pipe with rough boundary.Zeitschrift f¨ ur angewandte Mathematik und Physik , 71(4):138, 2020

  27. [34]

    Mori and L

    Y. Mori and L. Ohm. An error bound for the slender body approximation of a thin, rigid fiber sedimenting in stokes flow. Research in the Mathematical Sciences , 7(2):1–27, 2020

  28. [35]

    Mori and L

    Y. Mori and L. Ohm. Accuracy of slender body theory in approximating force exerted by thin fiber on viscous fluid. Studies in Applied Mathematics , 2021

  29. [36]

    Y. Mori, L. Ohm, and D. Spirn. Theoretical justification and error analysis for slender body theory. Communi- cations on Pure and Applied Mathematics , 73(6):1245–1314, 2020

  30. [37]

    Y. Mori, L. Ohm, and D. Spirn. Theoretical justification and error analysis for slender body theory with free ends. Archive for Rational Mechanics and Analysis , 235(3):1905–1978, 2020

  31. [38]

    Nishimura, C

    N. Nishimura, C. B. Schaffer, B. Friedman, P. D. Lyden, and D. Kleinfeld. Penetrating arterioles are a bottleneck in the perfusion of neocortex. Proceedings of the National Academy of Sciences , 104(1):365–370, 2007. A HIERARCHY OF BLOOD VESSEL MODELS, PART II: 3D-3D TO 3D-1D ...

  32. [39]

    F. Nobile. Coupling strategies for the numerical simulation of blood flow in deformable arteries by 3d and 1d models. Mathematical and Computer Modelling , 49(11-12):2152–2160, 2009

  33. [40]

    Notaro, L

    D. Notaro, L. Cattaneo, L. Formaggia, A. Scotti, and P. Zunino. A mixed finite element method for model- ing the fluid exchange between microcirculation and tissue interstitium. Advances in discretization methods: discontinuities, virtual elements, fictitious domain methods , ...

  34. [41]

    L. Ohm. On an angle-averaged Neumann-to-Dirichlet map for thin filaments. arXiv preprint arXiv:2308.06592 , 2023

  35. [43]

    R. N. Pittman. Influence of microvascular architecture on oxygen exchange in skeletal muscle. Microcirculation, 2(1):1–18, 1995

  36. [44]

    Possenti, A

    L. Possenti, A. Cicchetti, R. Rosati, D. Cerroni, M. L. Costantino, T. Rancati, and P. Zunino. A mesoscale computational model for microvascular oxygen transfer. Annals of Biomedical Engineering , 49:3356–3373, 2021

  37. [45]

    A. R. Pries and T. W. Secomb. Blood flow in microvascular networks. In Microcirculation, pages 3–36. Elsevier, 2008

  38. [46]

    Qi, S.-S

    Y. Qi, S.-S. Chang, Y. Wang, C. Chen, K. I. Baek, T. Hsiai, and M. Roper. Hemodynamic regulation allows stable growth of microvascular networks. Proceedings of the National Academy of Sciences , 121(9):e2310993121, 2024

  39. [47]

    Qi and M

    Y. Qi and M. Roper. Control of low flow regions in the cortical vasculature determines optimal arterio-venous ratios. Proceedings of the National Academy of Sciences , 118(34):e2021840118, 2021

  40. [48]

    T. W. Secomb, R. Hsu, E. Y. Park, and M. W. Dewhirst. Green’s function methods for analysis of oxygen delivery to tissue by microvascular networks. Annals of biomedical engineering , 32:1519–1529, 2004

  41. [49]

    A. Y. Shih, C. R¨ uhlmann, P. Blinder, A. Devor, P. J. Drew, B. Friedman, P. M. Knutsen, P. D. Lyden, C. Mateo, L. Mellander, et al. Robust and fragile aspects of cortical blood flow in relation to the underlying angioarchitecture. Microcirculation, 22(3):204–218, 2015

  42. [50]

    R. J. Shipley, A. F. Smith, P. W. Sweeney, A. R. Pries, and T. W. Secomb. A hybrid discrete–continuum approach for modelling microcirculatory blood flow. Mathematical medicine and biology: a journal of the IMA , 37(1):40–57, 2020

  43. [51]

    P. W. Sweeney, C. Walsh, S. Walker-Samuel, and R. J. Shipley. A three-dimensional, discrete-continuum model of blood pressure in microvascular networks. International Journal for Numerical Methods in Biomedical Engi- neering, 40(8):e3832, 2024

  44. [52]

    P. Zunino. Mathematical and numerical modeling of mass transfer in the vascular system . PhD thesis, Ecole Polytechnique F´ ed´ erale de Lausanne, Switzerland, 2002. Department of Mathematics, University of Wisconsin - Madison, Madison, WI 53706 Email address : lohm2@wisc.edu ...

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.