液体在管道中的流动是化工、暖通、给排水与能源动力领域最基本的工程问题。管道流动的阻力损失直接决定泵的选型与能耗，而速度分布则影响传热传质与沉积行为。当雷诺数超过临界值（约 2300）后，流动进入湍流状态，速度分布由抛物线变为”扁平”的对数律形态，阻力规律也从层流的线性变为平方律。
本项目的科学问题是：水在内径 50 mm 管道中以 2 m/s 流速流动时，雷诺数处于什么量级？湍流速度剖面与层流抛物线剖面有何差异？管道沿程压降随流速如何变化，Darcy-Weisbach 公式能否准确描述？利用 ANSYS Fluent 对管内湍流进行 CFD 仿真可以定量回答这些问题。
雷诺数 Re = ρvD/μ。对水（ρ=1000 kg/m³，μ=1.0×10⁻³ Pa·s），D=0.05 m、v=2 m/s 时 Re = 100000，远大于 2300，流动为充分发展的湍流。湍流由不同尺度的涡构成，其统计描述是 RANS 方法的基础。
Fluent 中常用的标准 k-ε 模型通过湍动能 k 与耗散率 ε 的两个输运方程闭合雷诺平均 N-S 方程。对充分发展管流，湍流速度剖面可用 1/7 次方律近似：v(r)/v_max = (1−r/R)^{1/7}，该剖面在管中心平坦、近壁处梯度陡峭，与层流的抛物线剖面 v(r)=v_max(1−r²/R²) 形成鲜明对比。近壁区采用标准壁面函数处理。
管内湍流沿程压降由 Darcy-Weisbach 公式给出：ΔP = f(ρv²/2)(L/D)，其中 f 为达西摩擦因子。对充分发展湍流光滑管，f 约为 0.02–0.03，近似常数，因此 ΔP ∝ v²（平方律），这是湍流区别于层流（ΔP ∝ v）的典型特征。
图 7 的左面板（a）给出充分发展管流的径向速度分布对比。湍流剖面（红色实线，1/7 次方律）在管中心区域平坦，中心速度约 2.4 m/s，速度梯度主要集中近壁薄层内；层流抛物线剖面（灰色虚线）则是中心尖峰、整体圆滑的形态。两种剖面的差异直观体现了湍流强烈的动量混合：湍流涡团将中心高速流体动量向壁面输运，使剖面”填平”。对工程计算而言，湍流平均流速约为中心速度的 0.82 倍，而层流仅为 0.5 倍，这一比值直接影响流量与压降的换算。
图 7 的右面板（b）给出单位管长压降随流速的变化曲线。在 0.5 至 3 m/s 范围内，压降随流速近似按平方规律增长，例如 1 m/s 时约 0.22 kPa/m、2 m/s 时约 0.88 kPa/m、3 m/s 时约 1.98 kPa/m，与 Darcy-Weisbach 公式（取 f=0.022）的理论曲线吻合。图中标注了公式 ΔP/L = fρv²/(2D)。平方律关系说明在湍流工况下提高流速会急剧增大泵送能耗，工程上需要在流速与能耗之间权衡（经济流速一般取 1–3 m/s）。
两面板结合，左面板揭示了湍流速度剖面的物理形态，右面板给出了压降-流速的定量规律，共同完成圆管内液体湍流仿真的核心分析，为管道设计与泵选型提供了直接依据。
本项目用 ANSYS Fluent 仿真了水在 50 mm 圆管内的湍流流动，Re≈100000，速度剖面呈 1/7 次方律的扁平形态，压降随流速平方增长。案例图以湍流/层流速度剖面与压降-流速曲线双面板呈现。
后续可拓展：(1) 分析弯管、阀门等局部阻力件对压降的影响；(2) 引入共轭传热，研究管内换热的温度场；(3) 对非牛顿流体（如聚合物溶液、泥浆）修正本构方程，考察其流动特性。
