面向多孔组织流体与固体应力建模的计算框架：应用于小梁化胚胎心脏
A computational framework for modeling fluid and solid stressesin porous tissues with application to the trabeculated embryonic heart
一个针对多孔组织流体-固体相互作用的高效计算框架被提出，并用于模拟小梁化阶段小鼠胚胎心脏心室壁的血流。研究基于 E12.5 小鼠心室的共聚焦显微成像与三维重建，先构建显式孔隙尺度有限元模型，再用理想化与图像衍生的单元胞估算水力渗透率，并与 Kozeny-Carman 固体-孔隙与纤维理论的刚性情形趋势对照，同时用超弹性单元胞材料量化组织弹性的影响。随后以完整流体与固体方程求解心室壁的多孔表示，并用基于 Brinkman 方程的连续介质模型替代显式孔隙结构。
来源：PLOS Computational Biology · journals.plos.org
by Christine Miller Buffinton, James W. Baish
This study provides a framework for efficient solutions of computational fluid-solid interaction models of porous tissue, with a focus on blood flow through the developing heart wall in the critical chamber and valve formation stages, when trabeculation fills much of the ventricles. This porous, trabeculated tissue is the site of complex mechanotransduction processes, but accurately estimating solid and fluid stresses within such intricate geometries is difficult. We addressed this with three levels of modeling and verification. Images of E12.5 mouse ventricles were first obtained through confocal microscopy and 3D reconstruction. These guided creation of explicit pore-scale finite element simulations with fluid-solid interaction. Idealized and image-derived unit cells were analyzed to determine hydraulic permeability. Next, a multi-pore representation of the ventricular wall was solved with full fluid and solid equations. A continuum model based on Brinkman equations then replaced the explicit pore structure, and sensitivity to permeability was evaluated in terms of fluid flow in the open channel and porous region. Permeabilities obtained from unit cell simulations followed trends from Kozeny-Carman solid-pore and fiber theory for the rigid case, while simulations for hyperelastic unit-cell material quantified the effects of tissue elasticity. Comparison of continuum and explicit model results showed the Brinkman equation, which incorporates the diffusion of shear stress from flow in the free fluid into the porous media, but not explicitly inertial effects, could significantly improve the accuracy of continuum models compared to Darcy's law. A 3D model of blood flow through the ventricles of the developing mouse heart demonstrated the application of the technique to more realistic geometry, showed locations of peak solid and shear stress, and produced physiologically reasonable flow. Furthermore, we developed methodologies to estimate local wall shear and solid stress from continuum-based calculations. This approach can enable efficient, whole-organ parametric analyses of fluid flow, circumventing the need for computationally expensive explicit pore modeling.