Parameter-robust multiphysics algorithms for Biot model with application in brain edema simulation
Parameter-robust multiphysics algorithms for Biot model with application in brain edema simulation
复制标题
Biot模型参数鲁棒多物理场算法在脑水肿模拟中的应用
DOI：
10.1016/j.matcom.2020.04.027
复制
发表时间：
2020
影响因子：
4.6
通讯作者：
Tian, Jing
中科院分区：
文献类型：
--
作者：
Ju, Guoliang;Cai, Mingchao;Li, Jingzhi;Tian, Jing
本文提出了Biot模型的参数稳健数值算法，并将其应用于脑水肿模拟。通过引入一个中间变量，我们得到了Biot模型的多物理重新表述。在重新表述的基础上，将Biot模型看作是一个结合了反应扩散子问题的广义Stokes子问题。同时或分开解决这两个子问题会产生耦合或解耦算法。我们进行了大量的数值实验，证明了这两种算法对于关键的物理参数是稳健的。将该算法应用于脑损伤区域脑脊液异常积聚引起的脑肿胀的研究。仔细研究了关键物理参数对脑肿胀的影响。结果表明，通透性对颅内压和组织变形影响最大，杨氏模数和泊松比对最大值影响不大，但对组织变形和脑肿胀发展速度有较大影响。
In this paper, we develop parameter-robust numerical algorithms for Biot model and apply the algorithms in brain edema simulations. By introducing an intermediate variable, we derive a multiphysics reformulation of the Biot model. Based on the reformulation, the Biot model is viewed as a generalized Stokes subproblem combining with a reaction–diffusion subproblem. Solving the two subproblems together or separately leads to a coupled or a decoupled algorithm. We conduct extensive numerical experiments to show that the two algorithms are robust with respect to the key physical parameters. The algorithms are applied to study the brain swelling caused by abnormal accumulation of cerebrospinal fluid in injured areas. The effects of the key physical parameters on brain swelling are carefully investigated. It is observed that the permeability has the biggest influence on intracranial pressure (ICP) and tissue deformation; the Young’s modulus and the Poisson ratio do not affect the maximum value of ICP too much but have big influence on the tissue deformation and the developing speed of brain swelling.