📚 课程目录第1章二维管道层流流动的Fluent求解与速度剖面验证

1. 二维管道层流流动的Fluent求解与速度剖面验证

案例目标与物理背景

本案例模拟二维直管内不可压缩层流流动。管道半径 R = 0.01 m（直径 D = 0.02 m），长度 L = 0.5 m（L/D = 25），流体为空气，密度 ρ = 1.225 kg/m³，动力粘度 μ = 1.7894×10⁻⁵ Pa·s。入口采用充分发展速度剖面（抛物线分布），对应雷诺数 Re = ρ·U_avg·D/μ = 200（层流）。

Poiseuille 解析解给出：

轴向速度分布：u(r) = 2·U_avg·[1 − (r/R)²]
轴线最大速度：u_max = 2·U_avg
沿程压降：Δp = 8·μ·L·U_avg / R²

取 U_avg = 0.146 m/s，则 u_max = 0.292 m/s，Δp ≈ 0.0836 Pa。求解后提取 x = 0.4 m 截面上的速度剖面，与解析解对比。

几何与网格策略

在 SpaceClaim 或 DesignModeler 中创建矩形面：x 从 0 到 0.5 m，y 从 −0.01 到 0.01 m。采用结构化四边形网格，轴向 200 等分，径向 40 等分（壁面附近加密）。网格总数约 8000 单元，满足层流计算精度要求。

边界命名：inlet（x=0）、outlet（x=0.5）、wall_top（y=0.01）、wall_bottom（y=−0.01）。

Fluent 求解设置

启动 Fluent 2D 双精度求解器，导入网格并检查最小正交质量（应 > 0.7）。
求解器设为 Pressure-Based，Absolute 速度格式，Steady 稳态。
Viscous 模型选择 Laminar。
材料：Air（ρ=1.225，μ=1.7894e-5）。
入口边界：Velocity Inlet，速度大小采用 UDF 或表达式 0.292*(1-(y/0.01)^2)，方向沿 +x。
出口边界：Pressure Outlet，表压 0 Pa。
壁面：无滑移静止壁面。
压力-速度耦合：SIMPLE；压力离散：Second Order；动量离散：Second Order Upwind。
初始化后设置收敛残差 1e-6，迭代 500 步。

入口速度剖面也可用 UDF 定义：

DEFINE_PROFILE(inlet_vel, thread, position){ real y[ND_ND]; face_t f; begin_f_loop(f,thread){ F_CENTROID(y,f,thread); F_PROFILE(f,thread,position)=0.292*(1-(y[1]/0.01)*(y[1]/0.01)); } end_f_loop(f,thread); }

速度剖面提取与验证

求解收敛后，在 Fluent 中创建 x = 0.4 m 的等值线（Iso-Surface），导出该面上 y 方向各点的 x 速度。将数据与解析解 u(y) = 0.292·[1 − (y/0.01)²] 对比。

预期结果：中心线速度约 0.292 m/s，壁面速度为 0，剖面呈标准抛物线。若中心速度偏差超过 2%，需检查网格分辨率或入口剖面是否准确施加。

结果判读与常见问题

现象	可能原因	处理方式
中心速度低于解析解	入口剖面未正确施加或网格过粗	检查 UDF/表达式，加密径向网格
剖面不对称	网格不对称或边界条件不一致	检查上下壁面设置，确保几何居中
残差不下降	入口速度过大导致非物理	确认 Re < 2300，降低松弛因子
出口出现回流	管道长度不足	延长管道至 L/D > 20

层流充分发展段长度约为 0.05·Re·D = 0.2 m。本案例管道长 0.5 m，x = 0.4 m 处已充分发展，可直接与 Poiseuille 解对比。若在 x < 0.2 m 处提取，剖面尚未完全发展，偏差会偏大。