放置内部流动的入口和出口边界
CFD 仿真通常需要大量计算,我们自然会尽量减少仿真中的自由度。如果将自由度减少到极致,我们可能会得到一个入口边界和出口边界相交的几何结构。让我们设想一个横截面为半圆形的管道中的 90°管道弯头。
弯头为 90°、横截面为半圆形的管道。如果使用上图显示的几何结构建立仿真,则入口和出口边界共享一条边。在许多情况下,单单这个现象就可能导致严重的收敛问题。不过,在这个特定示例中,经过几次迭代后,解会收敛。我们还考虑建立一个恰当的仿真,其中入口和出口管道延伸到十倍半径的长度(如下图所示)。
带延伸的入口和出口管道的 QA 90°弯头。
基于水力直径,在雷诺数为 120 的情况下执行仿真,其中 是横截面面积, 是横截面周长。在入口处应用均匀的速度分布,在出口处施加 0 法向应力。下图显示两个仿真中弯头处的压力,左侧是带延伸的入口和出口管道的情况,右侧是不带延伸管道的情况。对于带延伸的入口和出口管道的情况,从压力中减去下游边界上的平均压力,使两个结果在此边界上的平均压力都为 0。
横截面为半圆形的管道中 90°弯头上的压力变化,其中绘制了上游边界上的压力表面图和管道壁上的等压线。上侧的绘图显示带延伸的入口和出口管道情况下的结果,下侧的绘图显示不带延伸管道情况下的结果。
仿真结果显示,不带延伸的入口和出口管道的情况下,压力变化要大得多。由于所应用的均匀速度分布和壁面无滑移 边界条件之间的不相容性,靠近入口的壁上存在急剧的压力梯度。左侧的绘图显示弯头上游侧的压力更加均匀,这表明当流动到达弯头时已充分发展。不过,压力并不是完全均匀的,在靠近锐角的位置,压力看起来略低,这表明弯头的上游效应。我们还看到上游边界对面的管壁上有一个滞流点。弯头的损耗系数定义为
不带延伸的入口和出口管道的情况下,值为 2.3,带延伸的入口和出口管道的情况下,值为 0.60。通过观察速度场可以了解更多信息。
横截面为半圆形的管道中 90°弯头的速度分布和流线。
上图显示了弯头上游四个位置和弯头下游四个位置的速度分布以及中心平面的流线。在上游区域,我们可以看到均匀的速度分布如何演变成充分发展的速度分布。在弯头处,我们看到管壁上的滞流点朝向入口管道和相关的再循环区域。在急转弯头的下游还有一个再循环区域,我们可以看到,充分发展的速度分布首先在出口管道末端获得。所有这些在简单的几何结构(只包含 90°弯头)中都不存在,我们得到错误的压降也就不足为奇了。
“CFD 模块”是 COMSOL Multiphysics® 软件的附加模块,在入口和出口边界特征中提供了充分发展的流动选项,以避免入口和出口管道过长。前面两个图中的结果充分表明,我们应该在距离弯头一定距离处应用这些条件以获得好的结果。那么我们需要在上游和下游多远的位置应用充分发展的流动 选项?如果我们将入口管道和出口管道分别从弯头处延伸 1 个半径的长度,弯头的损耗系数将变为 0.54,而如果每个方向延伸 2 个半径的长度,则损耗系数变为 0.58。从此时开始,收敛到 0.60 的速度会变慢。因此,在本例中,每个方向延伸 2 个半径的长度是不错的权衡方法。
随着雷诺数的增大,弯头下游的再循环区域长度将增大,最终变得不稳定。对于雷诺数为 1200 的情况,如果管道末端应用充分发展的流动 选项,当出口管道延伸超过 20 个半径长度时,损耗系数不会有明显变化。根据管道入口长度的相关性,
适用于层流