阀件与集成流道内流 CFD：ζ/Kv-Cv 特性曲线族的生成与交付

课程大纲

• 能对阀体与集成壳体生成随阀位成族的内流负体积域，按阀位真值表定出端口对矩阵，并说明端口延伸段长度与取压面位置的选取理由（判据＝出口面无回流单元、延伸段加倍后 Δp 变化落进声明带内）
• 能对密封面做几何愈合，并用层流膜间隙泄漏量级式给出每一处微间隙「留还是封」的决策，说清该决策把交付物变成含内漏还是不含内漏
• 能用 h/Δx_n ≥ n_layer 与间隙雷诺数 Re_h 两条判据定间隙网格与流态处理，并在最恶阀位上做面向压降（目标量＝端口对 Δp，不是温度）的网格无关性
• 能先按 N_cases = N_pos × N_pair × N_Q 算出算例预算与机时量级再开算，把几何—边界—后处理三段全脚本化，输出 Δp-Q 与 ζ-Re 曲线族并显式标注外插边界
• 能用总压沿程分账把 ζ 偏高定位到死区 / 弯段与分配腔内的二次流 / 分配腔进口冲击 / 急转损失四类病灶之一，提出一条改型并回到整条曲线族复核而不是只看一个工况点
• 能按参考截面、Δp 类型、介质与 SG 折算、取压位置、阀位定义五项写出一维侧与台架能同时消费的交付口径，并用「网格未收敛 / 间隙流 regime 用错 / 台架直管段不足」三条归因清单排 CFD 与台架 Cv 的差值

1. 无热源流体件的问题定义与内流域 / 负体积抽取：随阀位成族的域、端口对矩阵与微间隙留封决策
   • 先定响应量再定域：本课的响应量是压降与端口流量分配（最终折成 ζ / Kv-Cv），不是壁温、不是换热量——这一句决定后面全部取舍（加密加在哪、y+ 目标定多少、网格无关性拿什么当目标量）。共轭传热与带热源的部件级问题归 J2-04 / J3-01 / J3-02，本课是单相、无热源、无相变的内流。
   • 负体积（negative volume）有两条生成路线：从水密外壳做布尔减、或对内表面直接包面。验收判据是三条——闭合、单一连通（除非本次就是要算多路并行）、端口面为平面且法向明确。通用脏数模清理、破面缝合与水密性检查归 J2-06，本课只补它 boundary 第⑥条让出来的阀件专属部分。
   • 阀件负体积的独有难点是「域的拓扑随阀位变」：同一份几何在不同阀位得到的根本不是同一个域，有的端口对连通、有的封断。所以域抽取必须以阀位为参数、一次生成一族域；抽出一个域只够算一个阀位，这是本课与 J2-06 通用流程的第一个分岔点。
   • 端口对矩阵的定义：N 端口阀的几何可能组合是 N(N−1)/2，但只有阀位真值表允许的那些组合有工程意义。真值表本体归 G8-06，本课不重新定义模式逻辑，只把它翻译成算例矩阵，并标出哪些组合是「按真值表排除」而不是「忘了算」。
   • 端口延伸段不是可选项：阀口出流强旋、出口面上会出现回流单元，取压面若贴在阀口上，测到的压差里混着动压畸变与回流噪声。做法是端口面外接直管延伸段（工程量级为数倍水力直径），取压面另设在延伸段内、位置与台架取压位置对齐（台架侧口径承 L1-02）。判据：出口面无回流单元，且延伸段长度加倍后 Δp 变化落在声明带内。
   • 密封面处的几何愈合：CAD 里密封件常缺省或已被简化，负体积会在密封面处出现零厚度接触、面重叠或虚假贯通。三种处置——按压缩后形态补出密封体、直接封断该路径、留一条名义间隙——选哪一种必须跟着「本次交付要不要含内漏」走，而不是跟着哪种最好画走。
   • 微间隙「留还是封」是本课最贵的一个决定：只交主通路流阻时封（可省一到两个数量级的网格量），要算内漏与模式串扰时必须留，且留出的 h 要取装配态 / 热态而非名义图纸值。开算前先用层流膜泄漏量级式估内漏占主流量的比例，低于本次要分辨的量级就封。内漏率目标本身由 G4-03 给，本课不设目标、只做建模决策。
   • 域交付前先声明「这份域代表哪一个装配态」——冷态名义、热态变形后、还是磨损末期。同一把阀这三个态给出的不是一条曲线；变形量与间隙链属结构侧（承 G4-03 的密封与间隙设计口径），本课不自算，但必须把所用假设写进第 5 讲的适用性声明卡。
2. 阀腔 / 急转 / 突扩 / 间隙的网格与流态策略：面向压降而非温度
   • 本讲承载本课最硬的独有块：缝隙尺度的网格判据。J2-01 的网格节只到边界层棱柱层与偏斜度 / 纵横比 / 正交质量 / 体积突变，没有「跨一条缝要几层」这一条；本判据不限于阀件，可泛化到齿顶 / 端面缝隙类问题（G1-09 的齿顶与齿侧缝隙、G5-07 的端面缝隙与困油腔与本问题同构，那两门现把缝隙网格难点前指 J2-01，实际落空）。
   • 跨间隙法向层数判据：h / Δx_n ≥ n_layer，工程量级 3~10 层——只要总泄漏量取下限，要解析膜内速度剖面或估局部冲蚀速度取上限。层数不够时数值扩散会给出一个与 h 弱相关的假值，而且误差方向随离散格式而变、事后无法修正，所以只能靠层数保证，不能靠「看着够密」或后处理修正。
   • 间隙内先判 regime 再定模型：以水力直径 2h 为特征尺度的 Re_h 常落在层流甚至爬流区，而同一算例的主通路同时是湍流——两个 regime 并存是阀件内流的常态。判据是分区看 Re_h，不是全域一个湍流模型定生死；全域套一个高 Re 模型会把间隙里的黏性主导流动算成湍流耗散。
   • regime 与泄漏定律的对应承 G4-03 的结论（层流膜区用 h³·Δp 式、冲蚀成孔口后才用 √Δp）；本课把这条判据反过来用——先用两式估泄漏量级，据此判该不该把缝建进域、要不要为它加密。G4-03 给判据本体与用途边界，本课给它在网格与建模决策里的用法，两课不重讲。
   • 阀腔急转、突扩与射流处的加密要跟着剪切层走，不是均匀加密全域：分离区、射流剪切层与再附着点决定 ζ 的大头。判据有两条——剪切层厚度方向要有若干层单元横跨，且再附着点位置在粗细两套网格间不跳变；位置一跳变，ζ 的绝对值就还没进入网格无关区。
   • 湍流模型选型、壁面律、y+ 与壁面处理的本体全部前指 J2-01，本课只声明选型口径：响应量是压降而非壁温时，对壁面处理的敏感度低于换热问题，y+ 落进壁面函数合法区间即可，不必为 ζ 把全域压到 y+ 约 1；例外是间隙内与射流剪切层——那里的黏性剪切直接进 Δp，必须解析。
   • 面向压降的网格无关性：目标量取端口对 Δp（及由它折出的 ζ / Kv），不取温度、不取壁面热流；除了 Δp 的相对变化要落进声明带内，还要看第二个量——间隙泄漏占主流量的比例是否随网格漂移，漂移就说明缝隙网格还没到位。GCI 与 Richardson 外推的本体见 J2-01，本课只换目标量与验收对象。
   • 网格无关性必须在最恶阀位（最小节流面积 / 最小间隙）上做，不能只在全开位做：全开位收敛的网格搬到中间开度与小开度常常不成立，因为那里的梯度最陡、间隙最小。这条顺序错了，整族曲线会在小开度段系统性偏。
3. 阀位与端口对扫描的参数化、批量后处理与 ζ / Kv-Cv 曲线族生成
   • 扫描规模先算再开算：N_cases = N_pos × N_pair × N_Q（必要时再乘温度点数 N_T），一次曲线族动辄上百个算例。这是本课区别于单点 CFD 的第一个工程约束——网格量级必须在算例预算的框架下定，而不是先把网格做到极致再发现跑不完。
   • 参数化的落点是「几何 + 边界 + 后处理」三段全脚本化：几何按阀位重生成、边界按流量点改、后处理按同一模板出数与出图。任何一段留手工，几十上百个算例既无法复现，也无法在几何改版后重跑——而集成模块的几何在样件阶段必然改版。
   • 阀位变化本课默认走「多位置分算」：每个阀位一份静态几何、一个独立算例。动网格 / 重叠网格与多位置分算的取舍已划归 J4-04，本课引用不重讲，只声明取分算的理由——本课交付的是稳态曲线族，不需要解切换过程本身。
   • 每个流量点的收敛判据落在监测量而不是残差：监测端口对 Δp 与各端口质量流量，残差平了而 Δp 还在漂就是假收敛（识别方法承 J2-01）。批量跑时这条必须写成自动判定并逐算例留痕，否则上百个算例里混进几个未收敛的没人看得出来。
   • ζ-Re 曲线的形态与外插：充分湍流段 ζ 趋于常数、低 Re 段随黏性作用抬升；拟合式只是工程拟合、不是物理定律，外插到扫描区间之外必须在交付物上显式标注。这条最容易出事的地方是低温——低温黏度陡升会把工作点推到扫描区间的低 Re 端（低温回路后果见 A2-01），常温三点扫出来的曲线族到冬季工况全靠外插。
   • Kv / Cv 的提取口径：定义式与选型用法承 G4-01 / G4-04 / G4-05，本课只管「从 CFD 场里怎么取 Q 与 Δp」——取压面位置、取静压差还是全压差、介质与 SG 折算，三项必须与台架口径（L1-02）逐条对齐，否则第 5 讲的互校差值全是口径差。
   • 交付物是「曲线族 + 适用性声明卡」，不是一叠曲线：卡上至少写清几何版本与装配态、介质与温度（即黏度）、扫描区间与外插边界、取压面与参考截面定义、含不含内漏、网格与收敛证据。A2-01 那条「有整体 Δp-Q 曲线就当黑箱」的红线，其成立前提正是这张卡。
   • 批量后处理要顺带产出下一讲要用的固定视图：同一套截面、云图与沿程总压曲线随算例自动生成。在上百算例规模下靠人工挑图不可执行，病灶定位能不能做，取决于这一步有没有把模板固化下来。
4. 流场病灶定位与改型：总压沿程分账定位四类病灶，含低于 p_v 区域的补充证据
   • 本讲不新起体系，只回答一个问题：曲线族已经算出来，ζ 偏高的那部分损失长在哪里。工具就是上一讲批处理产出的固定视图与沿程总压降落曲线，不额外引入新的求解方法。
   • 四类常见病灶及其判读特征：死区（低速回流区、总压几乎不降却占着体积，后果多在滞留与排气）、弯段与分配腔内的二次流旋涡对（把动压转成损失）、分配腔进口冲击（射流直冲对面壁，进口即掉一大段总压）、急转损失（弯转半径不足、内侧分离）。术语防撞——「二次流」在本体系另有两处同形异义（L1-01 风洞、O2-03 氢气引射器），交付文件须写全「弯段 / 分配腔内的二次流旋涡」。
   • 定位判据用沿流向的总压分账，不是看云图好不好看：把主通路切成若干截面，逐段算总压降占比，占比最大的那一两段才是改型对象。静压 / 动压 / 全压的口径与取压方式承 A2-06，本课只用其结论。
   • 改型手段与它们各自动哪个占比：加大弯转半径、加导流筋（治急转与二次流），改分配腔进口方向或加缓冲腔（治进口冲击），削死区（常同时改善排气与清洗）。流道几何准则本体（截面渐变、弯转半径、分流三通与汇流点、高低点与放净、流速许用区间）归 G8-02，本课只负责指出改哪一处、预期动多少占比。
   • 改型必须回到整条曲线族复核，不能只看一个点：局部改动常在大流量点见效、在小开度反向；判据是改前改后的整条 ζ-Re 与各端口对 Δp-Q 一起比，任一端口对变坏都要摆到台面上。
   • NPSHa 的补充证据：从 CFD 压力场里找出局部静压低于该温度饱和蒸汽压 p_v 的区域，作为泵入口与阀口汽蚀风险的位置性证据——它回答「哪里先掉到 p_v」，不替代 NPSHa 手算。汽蚀机理与 NPSH 判据归 G5-04，泵工作点与吸入条件归 A2-02，去气与定压归 G10-02，模块内进口去气与最低液位归 G8-02，本课一律引用不重讲。
   • 这条证据有两条硬边界，必须随结论一起说出来：稳态 RANS 给的是时均压力，瞬时低压谷比时均更低，所以判据偏乐观；且单相求解在真正发生汽蚀之后已不成立（含相变的两相求解见 J3-04）。结论只能用来给风险位置排序，不能拿来宣布「不会汽蚀」。
5. 三方交付与互校：喂一维、对台架、把瞬态边界交出去
   • 喂一维的口径统一：交给 A2-01 / J1-03 的是 ζ（或 K）与 Kv / Cv 两种表达，必须随数交清参考截面 A_ref 与 Δp 的类型（全压差还是静压差）——ζ 虽是无量纲数却依赖参考面，参考面选错整条曲线整体缩放。一维侧的叠加 ζ_total = Σζ_i、当量长度换算与网络求解归 A2-01 与 G8-02，本课只生产那一项数据。
   • 阀内过流槽这类局部阻力在一维侧要作为独立一项计入、不可按直管当量长度折算（G8-02 已立的规矩），且欠流校核要按该支路最恶阀位取值——本课交付的曲线族正是这两条规矩的数据源，所以最恶阀位必须在扫描矩阵里，不能只给全开位。
   • 与 L1-02 台架实测 Cv / Kv 的互校是同一个量的两条独立获取路径：比对前必须先对三项口径——介质与 SG 折算、取压位置、阀位定义（这三项 L1-02 已给，本课只接不重讲）。口径没对齐时，差值全部是口径差，不是仿真误差。
   • 偏差归因清单三条（CFD 侧两条 + 试验侧一条）：① 网格未收敛——回第 2 讲重查间隙层数与剪切层加密，并确认网格无关性是在最恶阀位做的；② 间隙流 regime 用错——层流膜区套了孔口式或反过来，表现是小开度段整族偏，判据回 Re_h；③ 台架直管段不足——取压面处流动未恢复、测到的压差含动压畸变，此条属试验侧，本课只列条目并前指 L1-02，不在此展开。
   • 归因顺序是先口径、再试验侧、最后才动模型：口径差与试验侧问题都会伪装成模型误差，先动模型就会把网格或湍流设置调歪来凑一个错的目标。通用相关性框架（误差口径、分层排查顺序、不达标时的可信度声明）归 J7-01，本课只留上面这几条阀件专属的。
   • 交给瞬态的是边界条件、不是瞬态解：阀位切换过程中的水锤与流量超调属一维压力波问题，由 J1-03 的水力模型承接，设计侧规避与阀位组合逻辑见 G8-06 / G4-05。本课不自建瞬态章节，交出去的是各阀位的稳态 ζ / Cv 与切换前后的端口对连通表。
   • 交付包的最小集是四件：曲线族数据文件、适用性声明卡、网格与收敛证据、互校与归因记录。通用的仿真任务输入输出与模型资产管理归 J7-07、不确定度量化的方法体系归 J7-06，本课只定这四件与它们的口径。
   • 要知道自己在给谁交货：A2-01（K 值手册覆盖不到的非标件）、J1-03（一维网络的 ζ / Cv 输入）、G4-03（内漏与流阻评估的方法层）、G8-02（集成后按实际流道重算模块内阻曲线）、B3-01 与 G8-06（拓扑与模式切换下的支路阻力）、L2-03（回路级流量分配反标）全是这份曲线族的消费者。交付口径一旦改动，这些下游的假设都要跟着复核。

• 以为把阀的负体积抽出来就能开算，其实同一份几何在不同阀位得到的根本不是同一个域——端口对的连通关系随阀位改变，域必须以阀位为参数成族生成；抽一个域只够算一个阀位。
• 以为端口面直接贴在阀口上省事，其实阀口出流强旋、出口面上会有回流单元，取压面落在未恢复的流场里，测到的压差混着动压畸变；必须外接延伸段并把取压面放到与台架同样的位置（台架侧口径见 L1-02）。
• 以为 ζ 是纯粹的无量纲数、报一个数就行，其实它依赖参考截面——换参考面 ζ 按面积比的平方缩放，同一算例能报出相差数倍的两个值；交付必须随数给出 A_ref 与「取全压差还是静压差」。
• 以为微间隙一律封死最省事，其实封不封取决于本次交付含不含内漏；封了却不写进声明卡，下游会拿一条不含内漏的曲线去评模式串扰。反过来全都留着会把网格量抬一个数量级，而绝大部分算例根本用不上。
• 以为间隙只要加密就能算准，其实跨间隙方向层数不够时数值扩散给出的是与 h 弱相关的假值，误差方向随离散格式而变、事后无法修正；判据是 h/Δx_n ≥ n_layer，不是「看着够密」。
• 以为全域一个湍流模型就够，其实阀件内流常常主通路湍流、间隙层流甚至爬流并存；regime 要按 Re_h 分区判，层流膜区套孔口式或反过来，都会让小开度段整族曲线系统性偏（regime 判据本体见 G4-03）。
• 以为网格无关性照搬 J2-01 那套就行，其实目标量换了——本课要的是端口对 Δp 收敛而非温度收敛，还要看间隙泄漏占比是否随网格漂移；而且必须在最恶阀位上做，全开位收敛的网格搬到小开度常常不成立。
• 以为曲线族算完就算交付，其实没有适用性声明卡的曲线族在一维侧必然被外插误用——低温高黏工况的工作点常落在扫描区间之外的低 Re 端，卡上不写外插边界，下游看不出自己已经出界。
• 以为 CFD 的 Cv 与台架的 Cv 对不上就是仿真不准，其实先要对三项口径（介质与 SG 折算、取压位置、阀位定义），不一致时差值全是口径差；口径对齐后剩下的差值才允许进归因清单，归因顺序与通用相关性框架见 J7-01。
• 以为从 CFD 压力场看到局部低于 p_v 就等于判定汽蚀，其实稳态 RANS 给的是时均压力、瞬时低压谷更低，且单相求解在真正发生汽蚀之后已不成立；这条证据只能给风险位置排序，NPSH 判据与机理仍归 G5-04 / A2-02 / G10-02。
• 以为阀位切换的水锤与流量超调也该由这门课算，其实那是一维压力波问题、归 J1-03 的水力模型，设计侧规避与阀位组合逻辑见 G8-06 / G4-05；本课只交各阀位的稳态 ζ / Cv 与切换前后的端口对连通表。
• 以为一次曲线族和一次单点 CFD 的工作量差不多，其实 N_pos × N_pair × N_Q（× N_T）动辄上百个算例，任何一段没脚本化都会让几何改版后无法重跑；算力不够时该砍的是算例覆盖面并写进声明卡，不是砍网格。

需求区 · 想听众筹

想听/点名均为需求登记，不涉及任何付款（意向金仅登记不收款）。满 5 人点名 → 自动向平台专家发邀约。