[论文解读] Comparing Poiseuille with 1D Navier-Stokes Flow in Rigid and Distensible Tubes and Networks

Taha Sochi|arXiv (Cornell University)|May 11, 2013

Enhanced Oil Recovery Techniques参考文献 24被引用 11

一句话总结

本研究通过解析关系和有限元模型，对比了刚性管道及网络中的泊肃叶流与可变形管道及网络中的1D纳维-斯托克斯流。结果表明，两种模型在流量预测上存在显著差异，尤其是在可变形系统中，由于压力依赖性、非线性和物理保真度的根本差异，二者不应互换使用。

ABSTRACT

A comparison is made between the Hagen-Poiseuille flow in rigid tubes and networks on one side and the time-independent one-dimensional Navier-Stokes flow in elastic tubes and networks on the other. Analytical relations, a Poiseuille network flow model and two finite element Navier-Stokes one-dimensional flow models have been developed and used in this investigation. The comparison highlights the differences between Poiseuille and one-dimensional Navier-Stokes flow models which may have been unjustifiably treated as equivalent in some studies.

研究动机与目标

• 严格比较刚性管道/网络中的泊肃叶流与可变形管道/网络中的1D纳维-斯托克斯流。
• 在不同物理和几何条件下，识别并量化两种模型之间的差异。
• 评估在血流动力学和多孔介质建模中，将泊肃叶流与1D纳维-斯托克斯模型视为等价模型的合理性。
• 评估每种模型在计算和数值方面的影响，包括内存成本和求解器收敛性。
• 确定在不同流动情景下，哪种模型更能真实反映物理现实，尤其在压力依赖性和瞬态行为方面。

提出的方法

• 基于恒定导率和节点处流量守恒，开发了基于解析的泊肃叶网络流模型。
• 为弹性管道中的1D纳维-斯托克斯方程构建了两种有限元实现，采用两种不同的压力-面积（p-A）本构关系。
• 将模型应用于单根管道和具有不同代数（5至15）、分支比和几何参数的分形状管道网络。
• 对泊肃叶模型使用直接线性系统求解器（内存需求为O(N²)），对非线性的1D纳维-斯托克斯模型使用迭代求解器（内存需求为O(4N²)）。
• 在一系列参数范围内比较流量（QP 与 QN）：α、β、μ、压力幅值、压力极限、ρ 和网络规模。
• 通过独立指标评估模型性能，如与物理原理的一致性以及实验测得的顺应性。

实验结果

研究问题

• RQ1 在单根刚性和可变形管道中，泊肃叶流与1D纳维-斯托克斯流在体积流量预测上存在哪些差异？
• RQ2 管道可变形性及压力-面积本构关系对泊肃叶流与1D纳维-斯托克斯流预测差异的影响如何？
• RQ3 在管网系统中，α、β、μ等参数及网络拓扑如何影响两种模型之间的差异？
• RQ4 计算成本和求解器收敛性问题如何影响每种模型的可靠性与实用性？
• RQ5 泊肃叶流在多大程度上能准确表征瞬态或空间变化的流场，相较于1D纳维-斯托克斯流？

主要发现

• 由于压力依赖性不同，泊肃叶流与1D纳维-斯托克斯流在可变形管道中产生显著不同的流量预测：泊肃叶流依赖于压力差，而1D纳维-斯托克斯流依赖于入口和出口的绝对压力。
• 1D纳维-斯托克斯模型具有非线性特性，需使用迭代求解器，导致更高的计算成本（内存需求O(4N²)）并可能引发收敛问题，而泊肃叶模型为线性，可在单次迭代内求解。
• 在管网系统中，压力-面积关系中的参数α和β对两种模型差异的影响最大，而μ和压力幅值等其他参数在所研究案例中仅产生中等或可忽略的影响。
• 泊肃叶模型提供线性压力场，可直接计算内部点的压力；而1D纳维-斯托克斯模型的非线性压力场无法实现这种直接计算。
• 1D纳维-斯托克斯模型可模拟瞬态流动，而泊肃叶模型本质上是时间无关的，仅能通过顺序稳态帧近似时间依赖行为。
• 本研究警告不应将两种模型视为可互换，尤其是在血流动力学背景下，强调模型选择应基于物理保真度，而非计算便利性。