也许这是个需要指出的问题_湍流_船舶_汽车_水工

也许这是个需要指出的问题

7月前浏览462

正文共： 2058 字 5 图预计阅读时间： 6 分钟

今天我们讨论一个可能对少数读者而言会忽略的问题，考虑如下的流道，假设液态水从母管右侧向左流动，希望通过 CFD 方法计算两个向下的支管流量分别为多少。

可能有少数读者直接就创建几何模型，然后将介质设置成液态水，估计几分钟就完成计算得出结果。但似乎这里有个重要的因素缺乏讨论，那就是重力影响要不要考虑。从生活常识就能看出，当右侧入口的水流量较小时，母管是不会被充满的，也即非满管流动。此时必须采用两相流模型才可以正确地模拟该流动问题。

而当右侧入口的水流量足够大时，母管才会被充满，也即满管流动。此时采用前文所说的几分钟就可以完成的单相流模拟方法即可获得正确的结果。为什么这个简单的问题今天拿出来讨论呢，因为确实在笔者周围出现了“粗心的” CFDer 开展了错误的模拟工作……

实际上，该问题的核心就是流体力学一个非常重要的无量纲数 —— 弗劳德数（ Fr ），即惯性力和重力之比。

当 Fr >> 1 时：惯性力占主导，重力影响可忽略。

例如：高速喷气、子弹飞行、空气中的高速气流。

此时流体的行为主要受惯性力和压力支配，重力效应相对微小。

当 Fr ≈ 1 或 Fr << 1 时：重力影响显著，必须考虑。

例如：明渠流动、船舶兴波、自由表面流动（如水闸泄流）、密度分层流等。

当流动存在自由表面或由重力驱动时（如浮力对流），重力是关键因素。

显然，重力加速度和特征长度确定后，速度就是是否需要考虑重力影响的唯一因素了，也就是上文所述的两种情况。这里需要强调的是特征长度的选取问题，特征长度 L 的选择没有绝对统一的标准，但它必须能代表 “ 重力影响流动的关键尺度 ” 。选择的核心原则是： L 应取在重力作用下，流动发生显著变化的那个方向上的特征尺度。选择不当会导致弗劳德数失去物理意义，从而做出错误判断。对于上述问题，重力影响水流在重力方向上的分布，因此特征长度可考虑选择母管的高度。

以下是让 Deepseek 帮忙整理的关于特征长度 L 的选取方法。

PART 01

按流动类型选择特征长度

1. 自由表面流动（最具代表性）

明渠流动（河流、水渠）：

通常选择水深 h ，因为重力沿水深方向驱动流动，自由表面的波动高度也与水深直接相关。

例如：分析河水流速是否会造成急流（ Fr>1 ）或缓流（ Fr<1 ）时，用平均水深。

船舶兴波、水工结构物（如桥墩、闸门）：

通常选择物体的特征长度，如船长、桥墩宽度。因为物体扰动了自由表面，产生的重力波尺度与物体尺寸相关。

波浪传播：

对于深水波，选择波长；对于浅水波，选择水深。

2. 物体在流体中运动（部分浸没或有自由表面）

船舶航行：常用水线长。

潜艇近水面航行：可用潜艇长度或潜深。

鱼雷入水：可用鱼雷直径或入水深度。

3. 内部流动（无自由表面，但重力可能影响）

管道中的气液两相流：若存在分层流，可选管道直径 D

自然对流（由温差引起的流动）：通常选择发生换热的垂直高度（如加热板的高度）。此时更多用格拉晓夫数，但特征长度选择原则相同。

4. 其他复杂流动

流体越过障碍物（如坝顶溢流）：可选择溢流头部的水深或坝高。

射流进入静止流体：可选择射流出口的孔径或宽度。

PART 02

核心选择原则总结

指向重力方向的原则：优先选择沿重力方向的尺度。对于水平明渠流，重力虽竖直，但通过压力梯度驱动水平流动，水深正是这个驱动力的关键尺度。

主导物理机制原则：问自己 “ 在这个问题中，重力主要通过哪种方式影响流动？ ”

如果是驱动流动（如河道），选水深。

如果是产生波浪（如船行波），选物体长度。

如果是引起分层（如大气），选层厚或高度。

一致性原则：在模型实验（如水工模型）中，必须保证模型与原型的弗劳德数相等。此时，特征长度的选择必须在模型和原型中对应同一个物理尺度（如都是水深或都是船长）。

PART 03

重要注意事项与技巧

当存在多个候选尺度时：如果水深 h 和渠宽 b

量级不同，需要根据主导作用选择。通常对于宽浅河道，选 h ；对于窄深河道，可能需要综合分析。

无量纲分析的本质：特征长度的选择本质上是为了使弗劳德数能够正确地表征惯性力与重力的比值关系。同一个流动，如果关心的问题不同（比如是整体流动还是局部冲刷），特征长度可能不同。

CFD/ 模拟中的设置：在软件中设置参考值时，特征长度必须与你选择的尺度一致，否则计算出的 Fr 数无意义。

与雷诺数选择的区别：雷诺数中的特征长度通常取流动方向的尺度（如管径、机翼弦长），而弗劳德数更常取与重力作用方向相关的尺度。

PART 04

实例辨析

案例 1 ：圆形管道满管湍流

是否需要考虑重力？通常不需要。因为无自由表面，且重力仅贡献于静压梯度（已包含在压力项中），不影响湍流结构。此时雷诺数主导，特征长度取管径 D 。

案例 2 ：圆柱绕流（空气中）

高速流动（如汽车风洞实验）： Fr 很大，重力影响可忽略。

但若圆柱是部分浸没在水中的桥墩，则必须考虑重力引起的自由表面波动和兴波阻力，特征长度应取圆柱直径或吃水深度。

PS ： AI 这样发展真是要替代好多人工了，但是怎么说呢， AI 本质上应该是供人类使用的工具，至少要懂得怎么使用吧，比如今天讨论的这个问题，至少要懂得想问的是什么吧？从这个角度看，基本的专业知识门槛还是有的吧。在 AI 面前，愿上帝保佑每个人类 ~~

来源：仿真与工程

FLUENT用pid算法模拟温度调节

正文共：2696字8图预计阅读时间：7分钟1前言今天我们做一个比较有意思的案例，之所以说有意思，是因为通常不会在设备层级做这个模拟。有如下的二维问题，左侧为热空气入口，右侧为冷空气入口，两者汇合后从下侧流出。假设热空气流量一定，但是温度呈现323.15~363.15K做正弦规律波动（任何规律都可以），拟采用pid算法调节冷空气入口流量（冷空气温度一定），使得出口的平均温度控制在314.5K。显然，需要用UDF来赋予冷空气入口边界条件，本案例的核心在于pid算法的实现。2边界条件与求解设置基本的设置这里不再赘述，仅仅说明一下关键设置，注意，本案例忽略了温度对介质热物性的影响。将热空气入口设置为速度入口，速度恒定在5m/s，温度采用UDF设定。将冷空气入口设置为速度入口，温度恒定在293.15K，速度值采用UDF设定。出口设置为压力出口，表压0Pa。本案例需要用Adjust挂载UDF实现pid调节。按如下值进行初始化，瞬态求解，时间步长0.01s。3计算结果先看一下热空气入口的温度曲线，按照正弦规律波动。再看一下冷空气入口的速度调节结果，符合预期，冷空气流量将随着热空气温度做同向波动，即热空气温度越高，需要越多的冷空气来降温。最后看一下出口的温度曲线，基本上稳定在314.5K附近，误差不大。本案例的UDF代码如下，需通过编译来调用，注意，pid三个参数（比例、积分、微分）需要自行整定，本案例提供了一组不错的设置。来源：仿真与工程