4. 管路系统几何建模

管路系统几何建模的目标,不是把每个管件都做成可制造级别的三维实体,而是为流动与传热仿真提供一套几何上足够真实、参数上可追溯、网格上可控制的计算域。直管、弯头、三通、变径管和阀门在仿真中的处理方式差异很大:直管主要影响沿程阻力与换热面积,弯头决定二次流与局部损失,三通控制分流与汇流,变径管引入加速/减速与压力恢复,阀门则往往被简化为局部阻力元件。理解这些差异,是参数化建模的起点。

一、从真实管路到计算域:几何简化的基本约定

几何简化必须遵循一条原则:保留对流动与传热有主导影响的几何特征,忽略对结果影响小于网格与模型误差的特征。典型做法包括:忽略法兰螺栓、焊缝余高、外壁保温层细节;将螺纹连接简化为等径圆柱段;将阀门内腔简化为具有等效流通面积的收缩—扩张结构或局部阻力段。简化后应记录每一项简化的依据,例如“忽略法兰,因其轴向长度小于管径的5%,对压降贡献可忽略”。

参数化建模的核心是建立几何参数与仿真参数之间的映射。常用几何参数包括:内径 D、壁厚 t、弯曲半径 R、分支夹角 θ、锥角 α、阀门开度 h/D。仿真参数则包括:水力直径 D_h、局部阻力系数 ζ、换热面积 A、网格第一层高度 y₁。参数化建模要求这些量能由几何参数直接计算或查表得到,而不是每次重新画图。

管路管件几何简化与参数化建模要素 几何参数 → 仿真参数映射关系(示例:D=50 mm,R/D=1.5,θ=45°) 直管段 参数:D、L、t、ε(粗糙度) 输出:D_h=D,A=πDL,ζ=f·L/D 弯头 参数:D、R、θ(90°/45°) 输出:R/D、ζ弯头、二次流强度 三通 参数:D主、D支、θ、流量比 输出:ζ直通、ζ支路、分流系数 变径管 参数:D₁、D₂、α(锥角) 输出:ζ渐扩/渐缩、压力恢复 阀门 参数:D、开度h/D、类型 输出:ζ(h/D)、等效流通面积 参数化映射:几何参数 → 仿真输入 D_h = 4A/P | ζ = f(Re, 几何) | y₁ = D_h·√(f/8)·(Re·√(f/8))⁻¹ 估算 网格:弯头/三通局部加密,直管段可粗化,阀门前后设过渡段 参数化建模流程与关键控制点 1. 提取几何参数 D、R、θ、α、h/D 2. 几何简化 去细节、保主导特征 3. 参数化脚本 尺寸驱动、批量生成 4. 网格划分 边界层、局部加密 5. 参数校验与迭代 ζ、压降、网格无关性 关键控制点: • 弯头 R/D ≥ 1.5 时二次流显著,需加密周向网格;R/D < 1.0 时按局部阻力处理更经济 • 三通支管与主管夹角 θ ≤ 45° 时分流损失较小,θ > 60° 时 ζ支路急剧增大 • 阀门开度 h/D < 0.2 时流通面积骤减,建议用等效收缩段替代真实阀芯几何

二、直管与弯头的参数化处理

直管是最简单的管件,但参数化时仍需明确三件事:内径 D、长度 L、壁面粗糙度 ε。水力直径直接取 D_h = D,换热面积 A = πDL。沿程阻力系数由 f 决定,层流时 f = 64/Re,湍流时可用 Colebrook 公式迭代或查 Moody 图。在参数化脚本中,直管通常用“圆柱段”模板生成,只需输入 D 和 L,网格沿轴向可设置渐变,入口段加密以捕捉发展中的边界层。

弯头的关键参数是弯曲半径 R 与内径之比 R/D。工程中常见 R/D = 1.0(短半径)和 R/D = 1.5(长半径)。R/D 越小,内侧分离越强,局部阻力系数 ζ 越大。以 90° 弯头为例,R/D = 1.0 时 ζ ≈ 0.9,R/D = 1.5 时 ζ ≈ 0.3,R/D = 3.0 时 ζ ≈ 0.15(示例值,实际随 Re 和表面粗糙度变化)。参数化建模时,弯头应沿弯曲方向生成扫掠网格,并在内侧壁面附近加密,以捕捉二次流和分离区。

三、三通、变径管与阀门的等效建模

三通的几何参数包括主管内径 D主、支管内径 D支、分支夹角 θ 和流量比。仿真中三通有两种处理方式:显式建模和局部阻力等效。显式建模适合需要观察分流流场、温度混合或支管回流的情况;等效建模则用直通和支路两个阻力系数 ζ直通、ζ支路 替代,适合系统级管网计算。参数化时,三通可用“主管+支管”布尔运算生成,支管与主管相交处需局部加密,避免网格畸变。

变径管的锥角 α 决定压力恢复能力。渐扩管 α ≤ 7° 时流动不易分离,压力恢复好;α > 15° 时出口附近出现分离,ζ 显著增大。渐缩管则相对稳定,ζ 通常小于渐扩管。参数化建模时,变径管用“圆台”模板,输入 D₁、D₂、L,锥角由 α = 2·arctan((D₂-D₁)/(2L)) 计算。网格沿轴向应保证锥面附近有足够分辨率。

阀门是管路中几何最复杂的管件,但仿真中通常不建真实阀芯。常见简化方式有三种:等效流通面积法,将阀门简化为一个收缩—扩张段,收缩面积等于阀门开度对应的流通面积;局部阻力系数法,直接给定 ζ(h/D) 曲线,用多孔介质或阻力源项替代;多孔介质法,将阀门区域设为各向异性多孔介质,通过渗透率控制压降。参数化时,阀门应记录开度 h/D、类型(闸阀、球阀、蝶阀、截止阀)和对应的 ζ 查表数据。

四、参数化建模的落地方式与校验

实际工程中,参数化建模通常借助脚本或专用工具完成。以下是一个简化的参数化建模伪代码示例,用于批量生成不同 R/D 的弯头几何:

# 弯头参数化生成示例(伪代码)
D = 50e-3          # 内径 50 mm
R_ratio = 1.5      # R/D
R = R_ratio * D    # 弯曲半径
theta = 90         # 弯曲角度
# 生成扫掠路径:圆弧,半径 R,角度 theta
path = arc(radius=R, angle=theta)
# 沿路径扫掠圆形截面,直径 D
bend = sweep(circle(diameter=D), path)
# 输出几何并标记边界:入口、出口、壁面
export(bend, boundaries=["inlet", "outlet", "wall"])

生成几何后,必须进行参数校验。校验内容包括:几何尺寸是否与设计值一致(如 D、R、θ);水力直径 D_h 是否按定义计算;局部阻力系数 ζ 是否与经验值或查表值吻合;网格质量指标(如正交性、长宽比、边界层 y⁺)是否满足求解器要求。对于弯头和三通,还应检查弯曲段和分支处的网格是否出现负体积或过度畸变。

实用建议:建立一份“管件参数表”,每行记录管件类型、几何参数、简化方式、仿真参数(D_h、ζ、A)和校验结果。这样在后续章节进行系统级仿真时,可以直接引用,避免重复建模。
注意:几何简化不是随意删减。任何简化都应评估其对压降、流量分配和换热的影响。例如,忽略阀门内腔细节可能使局部压降偏差超过 20%,此时应改用等效阻力元件而非直接删除阀门几何。

总结来说,管路系统几何建模的关键在于:用参数描述几何,用几何推导仿真输入,用校验保证可信度。直管、弯头、三通、变径管和阀门各有其主导参数和简化策略,参数化建模的价值在于让这些策略可复用、可批量执行、可追溯。掌握这一节的方法后,后续的网格划分、边界条件设置和求解器配置才能建立在可靠的几何基础之上。