流线 — CAE术语解说
流线
在CFD结果中会显示「流线」，但流线究竟精确地表示什么呢？
流线的理论基础
流线的定义与物理意义
「流线」就是CFD结果中经常看到的那条线对吧。不过流线、迹线、流脉线有什么不同呢？我一直以为它们是一样的。
很好的观察角度。这些概念在时间处理方式上有根本性的差异。流线是指某个瞬间，根据该时刻的速度向量所作的切线。而迹线是指某个流体质点随时间推进所经过的轨迹，流脉线是指在某一瞬间通过特定点的全部质点所形成的线。在定常流中，这三者是重合的，但在非定常流中则完全不同。比如，飞机翼周围的定常分析看流线，而直升机旋翼转动的非定常分析则应该看迹线。
那么，流线的数学定义是什么？能用公式表示吗？
设某个固定时刻 $t_0$ 的速度场为 $\mathbf{u}(x, y, z, t_0)$。流线是满足以下微分方程的曲线 $\mathbf{r}(s) = (x(s), y(s), z(s))$：
那Navier-Stokes方程与流线是什么关系？是方程求解速度场，然后再从中描绘流线吗？
完全正确。控制方程非压缩性Navier-Stokes方程是：
流线的数值计算方法
流线的计算算法
CFD软件从离散化的单元中心或节点速度数据出发，是怎样计算出光滑的流线的呢？
这是一个很核心的问题。主要有两种算法。第一种是Euler法的应用，最为简单。从种子点 $\mathbf{r}_0$ 开始，用 $\mathbf{r}_{n+1} = \mathbf{r}_n + \mathbf{u}(\mathbf{r}_n) \Delta s$ 求下一个点。但精度低，误差会不断积累，实务中几乎不用。
那实际应用中使用的方法是什么？
第二种也是标准使用的方法，是4阶/5阶Runge-Kutta法（RK45）。它在一个步长内多次评估速度，达到高阶精度。Ansys Fluent和OpenFOAM的 `streamLine` 函数默认采用这种RK法。步长 $\Delta s$ 通常设为单元大小的10%到50%，有些求解器还有自动控制功能。
离散数据「之间」的速度怎样确定？只有单元中心的值啊。
这是很重要的一点，需要插值。基于FVM（有限体积法）的软件，将单元中心的值插值到节点，再通过形状函数对任意点进行内插。Ansys CFX在单元内进行二次插值。插值的精度直接关系到流线的光滑性和物理正确性。在粗网格上尝试追踪细小涡时，插值误差会导致流线呈锯齿状，或出现实际不存在的结构。
流线计算的停止条件是什么？总不能一直计算下去吧。
主要有4个终止条件。1) 计算域的边界；2) 进入极低速区域（$|\mathbf{u}| < 10^{-6} \, \text{m/s}$ 等）；3) 流线长度超过用户定义的最大值（比如100m）；4) 流线形成与自身相交的闭合曲线（循环）。Siemens Star-CCM+ 可以对这些条件进行详细设置。
流线的实务应用
有效的可视化与解释
看分析结果时，流线的种子点应该放在哪里最有效？随意放置的话会很乱。
目的不同，策略就不同。比如要查看流入条件的均匀性，就在流入边界面上均匀间隔放置数十个点。要可视化剥离涡，就在翼或车身产生剥离的边缘后面线性放置。这就是Ansys Fluent的「Line/Rake」种子功能。要寻找死区（循环区），可以在整个区域随机放置数百个点，用「Point Cloud」功能。核心原则是「在要观察的物理现象的上游放置种子」。
给流线着色时，通常用什么物理量来区分颜色？速度大小吗？
速度大小虽然基础，但仅靠这个不够。实务中常用压力（特别是全压损失）、涡度（涡的强度）、温度（热的传输路径）、湍流能（混合的程度）来着色。比如，对于泵内部流动，按压力给流线着色，可以一眼看出流线是否通过易发生空化的低压区。颜色范围应设置为有物理意义的值。比如温度就设在流入温度和壁面温度之间。
有时流线看起来像是穿过了网格。这是计算出错的标志吗？
不一定是计算误差。可视化软件不显示网格，只显示流线时就经常出现这种情况。但流线穿过壁面就是明显的问题。主要原因有两个。1) 壁面边界条件（通常是无滑移）没有正确设置；2) 网格太粗，无法解析壁面附近的速度梯度，插值时壁面处的速度不为零。后者需要检查壁面的剪应力或 y+ 分布。
流线能用来定量评估，比如「流的均匀性」吗？
直接从流线本身提取数值很难。流线是定性评估的专家。定量评估需要用其基础的速度场。比如评估热交换器翅片间流动的均匀性，就计算下游特定截面的速度分布标准差或最大速度/平均速度比。流线的作用是帮助你从视觉上理解不均匀在哪里、为什么发生。可视化（流线）和定量数据（图表、等高线）的结合才是金科玉律。
流线的软件比较
各软件中的实现与特点
Ansys Fluent和Siemens Star-CCM+ 在流线计算功能和显示方式上有区别吗？
差异很大。首先是种子设置。Fluent传统上提供「点」「线」「面」「体」等多种方式，但设置界面不太直观。Star-CCM+ 用「零件」（几何要素：点、线、面、体）直接设置种子，界面很直观，特别是用「Field Function」可以用函数定义复杂的种子分布，这是很大的优势。
计算算法的性能有差异吗？
两者都实现了高阶Runge-Kutta法，但并行计算支持有差异。Fluent流线计算本身已经并行化，但生成极多数量（数万条）的流线时，后处理阶段会变慢。Star-CCM+ 模拟和后处理集成在一起，流线生成也能高效地进行并行分布处理。在1000万单元模型上生成数万条流线时，Star-CCM+ 通常更快。
像ParaView这样的免费可视化软件呢？
ParaView非常强大。可以读入Fluent或OpenFOAM数据，用「Stream Tracer」过滤器。算法可以选RK2、RK4、RK45，种子设置也很多样。在功能上不亚于商业软件，但学习曲线陡，超大数据（数十GB以上）的内存管理很头疼。VisIt也类似，但不如ParaView流行。在有经费限制的研究机构或项目中，ParaView是首选。
像COMSOL Multiphysics这样的多物理场软件呢？
COMSOL的特点是把流线视为「粒子追踪」模块的一部分。这样可以很容易地沿流线计算任意变量（温度、化学种浓度、电位等）的积分值。比如在微流控芯片内沿特定路径的反应物平均停留时间，可以用流线很容易算出。这是纯CFD软件没有的优势，特别是针对多物理场分析。但对大规模湍流分析的性能比不上专业CFD软件。
流线的故障排除
常见问题与对策
流线在计算域边缘处突然切断，不再延伸。原因是什么？
最常见的原因有两个。第一是出口边界条件。在流出（Outflow）或压力出口（Pressure Outlet），流线在该面处终止，这是正确的行为。问题是流线在途中消失的情况，可能是该位置的单元数据损坏（NaN或异常值）。第二是网格质量。极度细长或扭曲的单元会导致插值计算失败，无法获得速度向量。解决办法是先检查问题区域的网格，改进单元质量（Skewness < 0.95）。
流线呈锯齿状，不光滑。细化网格能解决吗？
细化网格是根本解决办法，但计算代价大。首先要试的是调整流线计算求解器设置。把步长 $\Delta s$ 从默认值（比如单元大小的20%）减小到5%或10%。其次，把算法从一阶Euler法改为四阶Runge-Kutta法。在ParaView中，把「Integration Step Unit」从「Cell Length」改为「Length」，指定绝对步长（比如0.001 m）能更稳定。都不行的话，对问题区域进行局部网格细化最保险。
用同一组数据，Fluent里看流线和ParaView读入数据后看的流线形状略有差异。哪个可信？
这很常见。原因主要是数据插值方法和种子点精确位置的区别。Fluent用自有内部插值法，ParaView通过VTK库读数据用VTK的插值法。种子点坐标也可能存在舍入误差。可信的是大的流动模式（是否剥离、大涡位置）。细节差异不必过分关注。但如果流线上的压力或速度极值相差很大，要怀疑数据导出/导入过程中有损失或转换误差。
非定常分析里，某个时刻的流线做成动画后，流线瞬间变断断续续。这是正确的吗？
这是完全正确的。非定常流中，流线始终基于「瞬间的速度场」。时间推进时，速度场本身就变化了，流线形态当然也不连续。如果看到流线平顺地变形的动画，那可能是显示的迹线而非流线，或者可视化本身有问题。要观察非定常流的时间演化，应该用粒子追踪（Particle Tracking）功能生成迹线动画，或者把多个时刻的流线并排比较。Ansys Ensight和Techplot在这类时序可视化上很强。
相关主题
价值
更详细
错误