[论文解读] Pore-Scale Modeling of Navier-Stokes Flow in Distensible Networks and Porous Media

Taha Sochi|arXiv (Cornell University)|Sep 29, 2013

Lattice Boltzmann Simulation Studies参考文献 41被引用 10

一句话总结

该论文提出了一种基于流量连续性残差的牛顿-拉夫森迭代方法，用于模拟可变形管道和多孔介质中的一维纳维-斯托克斯流。通过利用弹性管道的解析压力-流量关系，该方法在复杂非线性流态下实现了比传统有限元方法更优的收敛性，同时具备高精度、高鲁棒性和计算高效性。

ABSTRACT

In this paper, a pore-scale network modeling method, based on the flow continuity residual in conjunction with a Newton-Raphson non-linear iterative solving technique, is proposed and used to obtain the pressure and flow fields in a network of interconnected distensible ducts representing, for instance, blood vasculature or deformable porous media. A previously derived analytical expression correlating boundary pressures to volumetric flow rate in compliant tubes for a pressure-area constitutive elastic relation has been used to represent the underlying flow model. Comparison to a preceding equivalent method, the one-dimensional Navier-Stokes finite element, was made and the results were analyzed. The advantages of the new method have been highlighted and practical computational issues, related mainly to the rate and speed of convergence, have been discussed.

研究动机与目标

• 开发一种计算高效且鲁棒的孔尺度网络建模方法，用于模拟可变形管道和顺应性多孔介质中的流体流动。
• 克服传统有限元方法在可变形结构中求解一维纳维-斯托克斯流时计算成本高、复杂度大的问题。
• 利用弹性管道的解析压力-面积本构关系，实现在非线性、时间无关系统中的高精度流量预测。
• 从收敛性、精度和计算效率三个方面，将所提出的孔尺度方法与传统的基于一维纳维-斯托克斯方程的有限元方法进行比较。
• 通过收敛性度量验证该方法，并借助系统化的数值测试与参数调优，证明其可靠性。

提出的方法

• 该方法基于网络节点处体积流量连续性条件，采用残差驱动的公式化方法，以在相互连接的可变形管道间强制实现质量守恒。
• 采用牛顿-拉夫森迭代格式求解所得的非线性方程组，每次迭代中更新各节点的压力和流量变量。
• 采用一维纳维-斯托克斯流在顺应性管道中的解析解，通过压力-面积本构关系将边界压力与体积流量关联。
• 该方法在由相互连接的管道构成的网络上实现，将每个管道视为具有由杨氏模量和泊松比定义的弹性壁属性的顺应性单元。
• 通过使用泊肃叶流作为初始猜测、压力渐变法以及自适应降低误差容差等策略，提升收敛性。
• 通过残差范数和收敛性度量对方法进行验证，并探索动量通量校正因子α对稳定性的影响。

实验结果

研究问题

• RQ1 与基于一维纳维-斯托克斯方程的有限元方法相比，孔尺度残差驱动牛顿-拉夫森方法在求解可变形网络流动时表现如何？
• RQ2 在高度非线性流动状态下，所提出的孔尺度方法的收敛特性与计算效率如何？
• RQ3 能否在网络建模框架中有效整合弹性管道的解析压力-流量关系，以确保精度与鲁棒性？
• RQ4 在高压梯度条件下，哪些实用策略可提升孔尺度方法的收敛速度与稳定性？
• RQ5 当通过解析基准和收敛性准则进行验证时，孔尺度方法在多大程度上保持了精度与可靠性？

主要发现

• 在高压梯度的非线性区域，该孔尺度方法的收敛速度和稳定性显著优于有限元方法。
• 该方法的计算结果与顺应性管道流动的解析模型高度一致，表现出极高的精度与可靠性。
• 通过使用泊肃叶流解作为初始猜测，并采用压力渐变技术逐步增加压力梯度，可显著提升收敛性。
• 调整动量通量校正因子α可增强收敛性与解的质量，表明该参数在非线性系统中具有显著影响。
• 与有限元方法相比，该方法在计算成本上更低、鲁棒性更强，尤其在复杂、大规模可变形网络中表现更优。
• 通过残差范数和收敛性度量的验证，证实该方法能够可靠地检测并排除不准确解，确保仿真的高保真度。