A Lattice Boltzmann Method for Non-Newtonian Blood Flow in Coiled Intracranial Aneurysms
Pith reviewed 2026-05-24 03:14 UTC · model grok-4.3
The pith
A lattice Boltzmann method with volume-averaged Navier-Stokes equations models non-Newtonian blood flow through coiled aneurysms by treating the coil as an inhomogeneous porous medium.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the volume-averaged Navier-Stokes equations for non-Newtonian blood, closed by an inhomogeneous porous-medium representation of the coil, can be solved with a lattice Boltzmann method on patient-derived geometries to produce flow predictions whose essential features match those obtained from fully resolved simulations, thereby enabling patient-specific assessment of flow changes induced by coil insertion.
What carries the argument
Volume-averaged Navier-Stokes equations closed by an inhomogeneous porous-medium model for the coil and discretized by a problem-adapted lattice Boltzmann method.
If this is right
- Flow predictions for coiled aneurysms become feasible on patient-specific meshes without explicit resolution of individual coil filaments.
- The same workflow can be used to compare alternative coil placements or packing densities before intervention.
- Non-Newtonian rheology is retained throughout the porous-medium region, so shear-thinning effects remain present after treatment.
- The lattice Boltzmann discretization handles the transition between free-flow and porous regions without additional interface conditions.
Where Pith is reading between the lines
- The method could be combined with time-resolved imaging to update permeability parameters from post-treatment angiograms and thereby refine predictions for follow-up interventions.
- Because the porous-medium description is local, the framework extends in principle to other endovascular devices such as flow-diverter stents whose effective permeability can be calibrated similarly.
- If the permeability maps can be derived from coil packing density alone, the workflow becomes a candidate for rapid in-clinic estimation of hemodynamic risk after embolization.
Load-bearing premise
Representing the coil as an inhomogeneous porous medium with suitable permeability parameters, together with the volume-averaged Navier-Stokes equations, sufficiently captures the essential flow changes caused by coil insertion.
What would settle it
A side-by-side comparison in the same patient geometry in which the volume-averaged model produces statistically significant differences in wall shear stress, intra-aneurysmal velocity, or recirculation patterns relative to both fully resolved simulations and available experimental data would falsify the validity statement.
read the original abstract
Intracranial aneurysms are the leading cause of hemorrhagic stroke. One of the established treatment approaches is the embolization induced by coil insertion. However, the prediction of treatment and subsequent changed flow characteristics in the aneurysm is still an open problem. In this work, we present an approach based on a patient-specific geometry and parameters including a coil representation as inhomogeneous porous medium. The model consists of the volume-averaged Navier-Stokes equations for a non-Newtonian blood rheology. We solve these equations using a problem-adapted lattice Boltzmann method and present a comparison between fully-resolved and volume-averaged simulations. The results indicate the validity of the model. Overall, this workflow allows for patient specific assessment of the flow due to potential treatment.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a lattice Boltzmann discretization of the volume-averaged Navier-Stokes equations with a non-Newtonian blood rheology for patient-specific intracranial aneurysm geometries. Coils are represented as an inhomogeneous porous medium whose permeability parameters are chosen to match a fully-resolved coil geometry; the authors perform a comparison between the fully-resolved and volume-averaged simulations and conclude that the results indicate model validity, thereby enabling patient-specific post-treatment flow assessment.
Significance. If the porous-medium approximation is shown to reproduce the essential flow diversion and intra-aneurysmal stasis with quantifiable accuracy, the method would offer a computationally tractable route to treatment planning that avoids explicit resolution of individual coil wires. The use of a problem-adapted LBM and non-Newtonian rheology are standard strengths in this domain.
major comments (2)
- [Abstract] Abstract: the statement that 'the results indicate the validity of the model' is unsupported because the manuscript reports no quantitative comparison metrics (relative L2 errors on velocity or wall-shear-stress fields, integrated flow-rate differences, or sensitivity of results to the chosen permeability values) between the fully-resolved and volume-averaged simulations.
- [Model description / Results] Model description and results sections: the central claim that the inhomogeneous porous-medium representation 'sufficiently captures the essential flow changes for treatment prediction' rests on an unquantified agreement; without reported error thresholds or parameter-selection procedure, it is impossible to judge whether the approximation meets the accuracy needed for the stated clinical workflow.
minor comments (1)
- The permeability tensor is described as 'inhomogeneous' but the precise spatial variation and the fitting procedure used to obtain its components from the fully-resolved geometry are not stated explicitly.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback on our manuscript. The comments correctly identify that the comparison between fully-resolved and volume-averaged simulations is presented primarily through qualitative figures without accompanying quantitative error metrics. We address each point below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Abstract] Abstract: the statement that 'the results indicate the validity of the model' is unsupported because the manuscript reports no quantitative comparison metrics (relative L2 errors on velocity or wall-shear-stress fields, integrated flow-rate differences, or sensitivity of results to the chosen permeability values) between the fully-resolved and volume-averaged simulations.
Authors: We agree that the abstract claim would be strengthened by quantitative support. The current manuscript presents the comparison via side-by-side visualizations of velocity and flow patterns that demonstrate agreement in flow diversion and stasis, but does not include the listed numerical metrics. In revision we will add relative L2 errors on velocity, wall-shear-stress differences, integrated flow-rate discrepancies, and a brief sensitivity study with respect to permeability values, then update the abstract to reflect these results. revision: yes
-
Referee: [Model description / Results] Model description and results sections: the central claim that the inhomogeneous porous-medium representation 'sufficiently captures the essential flow changes for treatment prediction' rests on an unquantified agreement; without reported error thresholds or parameter-selection procedure, it is impossible to judge whether the approximation meets the accuracy needed for the stated clinical workflow.
Authors: The model description section explains the permeability calibration by matching the effective hydraulic resistance of the fully-resolved coil geometry via separate simulations on the same domain. The results section then shows that the volume-averaged model reproduces the main flow features. We acknowledge, however, that no explicit error thresholds or quantitative assessment of clinical relevance are supplied. We will therefore include error thresholds, report the parameter-selection procedure more explicitly, and discuss the observed accuracy level relative to the requirements of post-treatment flow assessment. revision: yes
Circularity Check
No circularity; validity rests on independent fully-resolved comparison
full rationale
The paper sets up volume-averaged Navier-Stokes with non-Newtonian rheology and an inhomogeneous porous-medium coil representation, then solves via adapted LBM. The central validity claim is explicitly tied to a comparison against separate fully-resolved simulations rather than any fitted parameter, self-citation chain, or definitional equivalence. No step reduces a prediction to its own inputs by construction, and the derivation chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- porous medium permeability parameters
axioms (1)
- domain assumption Volume-averaged Navier-Stokes equations with non-Newtonian rheology accurately represent flow inside and around the coiled region when coils are treated as porous medium
Reference graph
Works this paper leans on
-
[1]
Ashrafizaadeh, M., Bakhshaei, H.: A comparison of non-Newtonian models for lattice Boltz- mann blood flow simulations. Comput. Math. Appl. 58(5), 1045–1054 (2009). DOI 10.1016/j.camwa.2009.02.021
-
[2]
Ataei, M., Salehipour, H.: XLB: A differentiable massively parallel lattice Boltzmann library in Python. Comput. Phys. Commun. 300, 109187 (2024). DOI 10.1016/j.cpc.2024.109187
-
[3]
Boyd, J., Buick, J.M., Green, S.: Analysis of the Casson and Carreau-Yasuda non-Newtonian blood models in steady and oscillatory flows using the lattice Boltzmann method. Phys. Fluids 19(9), 093103 (2007). DOI 10.1063/1.2772250
-
[4]
Eddleman, C.S., Welch, B.G., Vance, A.Z., Rickert, K.L., White, J.A., Pride, G.L., Purdy, P.D.: Endovascular coils: properties, technical complications and salvage techniques. J. Neurointerv. Surg. 5(2), 104–109 (2013). DOI 10.1136/neurintsurg-2012-010263
-
[5]
Ergun, S.: Fluid flow through packed columns. Chem. Eng. Prog. 48, 89–94 (1952) 10 M. Horvat, S. B. Lunowa, D. Sytnyk, and B. Wohlmuth
work page 1952
-
[6]
Fritz, M., K ¨oppl, T., Oden, J.T., Wagner, A., Wohlmuth, B., Wu, C.: A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue. Int. J. Numer. Meth. Bio. 38(7), e3612 (2022). DOI 10.1002/cnm.3612
-
[7]
Geier, M., Sch ¨onherr, M., Pasquali, A., Krafczyk, M.: The cumulant lattice Boltzmann equation in three dimensions: Theory and validation. Comput. Math. Appl.70(4), 507–547 (2015). DOI 10.1016/j.camwa.2015.05.001
-
[8]
Guo, Z., Zheng, C., Shi, B.: Discrete lattice effects on the forcing term in the lattice Boltzmann method. Phys. Rev. E 65(4), 046308 (2002). DOI 10.1103/physreve.65.046308
-
[9]
He, Y.L., Liu, Q., Li, Q., Tao, W.Q.: Lattice Boltzmann methods for single-phase and solid- liquid phase-change heat transfer in porous media: A review. Int. J. Heat Mass Tran. 129, 160–197 (2019). DOI j.ijheatmasstransfer.2018.08.135
work page 2019
-
[10]
H ¨ocker, S.B., Trunk, R., D¨orfler, W., Krause, M.J.: Towards the simulations of inertial dense particulate flows with a volume-averaged lattice Boltzmann method. Comput. Fluids 166, 152–162 (2018). DOI 10.1016/j.compfluid.2018.02.011
-
[11]
Holzberger, F., Muhr, M., Wohlmuth, B.: A comprehensive numerical approach to coil place- ment in cerebral aneurysms: Mathematical modeling and in silico occlusion classification (2024). ArXiv:2402.02798. To appear in Biomech. Model. Mechan
-
[12]
Hosseini, S.A., Huang, F., Th ´evenin, D.: Lattice Boltzmann model for simulation of flow in intracranial aneurysms considering non-Newtonian effects. Phys. Fluids34(7), 073105 (2022). DOI 10.1063/5.0098383
-
[13]
Kr¨ uger, T., Kusumaatmaja, H., Kuzmin, A., Shardt, O., Silva, G., Viggen, E.M.: The Lattice Boltzmann Method. Springer (2017). DOI 10.1007/978-3-319-44649-3
-
[14]
Maier, M.L.: Coupled lattice Boltzmann and discrete element method for reactive particle fluid flows with applications in process engineering. Ph.D. thesis, Karlsruher Institut f¨ ur Technologie (KIT) (2021). DOI 10.5445/IR/1000132643
-
[15]
Matyka, M., Koza, Z., Lukasz Miros law: Wall orientation and shear stress in the lattice Boltzmann model. Comput. Fluids73, 115–123 (2013). DOI 10.1016/j.compfluid.2012.12.018
-
[16]
https://doi.org/10.1007/978-1- 4684-3384-5_11
Nield, D.A., Bejan, A.: Convection in Porous Media. Springer (2013). DOI 10.1007/978-1- 4614-5541-7
-
[17]
Otani, T., Ii, S., Shigematsu, T., Fujinaka, T., Hirata, M., Ozaki, T., Wada, S.: Computational study for the effects of coil configuration on blood flow characteristics in coil-embolized cerebral aneurysm. Med. Biol. Eng. Comput. 55, 697–710 (2017). DOI 10.1007/s11517-016- 1541-6
-
[18]
Ouared, R., Chopard, B.: Lattice Boltzmann simulations of blood flow: Non-Newtonian rhe- ology and clotting processes. J. Stat. Phys. 121, 209–221 (2005). DOI 10.1007/s10955-005- 8415-x
-
[19]
Pierot, L., Wakhloo, A.K.: Endovascular treatment of intracranial aneurysms. Stroke 44, 2046–2054 (2013). DOI 10.1161/STROKEAHA.113.000733
-
[20]
Romero Bhathal, J., Chassagne, F., Marsh, L., Levitt, M.R., Geindreau, C., Aliseda, A.: Mod- eling flow in cerebral aneurysm after coils embolization treatment: a realistic patient-specific porous model approach. Cardiovasc. Eng. Tech.14(1), 115–128 (2023). DOI 10.1007/s13239- 022-00639-x
-
[21]
Sequeira, A.: Hemorheology: Non-Newtonian Constitutive Models for Blood Flow Simula- tions, pp. 1–44. Springer (2018). DOI 10.1007/978-3-319-74796-5 1
-
[22]
Succi, S.: The Lattice Boltzmann Equation: For Complex States of Flowing Matter. Oxford University Press (2018). DOI 10.1093/oso/9780199592357.001.0001
-
[23]
Yadollahi-Farsani, H., Herrmann, M., Frakes, D., Chong, B.: A new method for simulating embolic coils as heterogeneous porous media. Cardiovasc. Eng. Techn. 10, 32–45 (2019). DOI 10.1007/s13239-018-00383-1
-
[24]
Zhao, J., Lin, H., Summers, R., Yang, M., Cousins, B.G., Tsui, J.: Current treatment strate- gies for intracranial aneurysms: an overview. Angiology 69(1), 17–30 (2018). DOI 10.1177/0003319717700503
-
[25]
Zou, Q., He, X.: On pressure and velocity boundary conditions for the lattice Boltzmann BGK model. Phys. Fluids 9(6), 1591–1598 (1997). DOI 10.1063/1.869307
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.