### 代码功能解构与执行逻辑 用户提供的 Python 代码并非直接嵌入 ANSYS Fluent 内部的 C 语言 UDF（用户自定义函数），而是一个**外部自动化控制脚本（Wrapper Script）**。其核心逻辑是充当“指挥官”，通过循环调用 CFD 求解器来执行参数化研究 [ref_4]。在核工程两相流模拟中，确定临界吹扫气量通常需要进行多次试算，手动修改边界条件效率极低且易出错，该脚本旨在实现“设定入口速度 -> 自动计算 -> 读取结果 -> 判断是否排空 -> 调整速度”的闭环流程 [ref_6]。代码中的 `update_boundary_condition`、`run_cfd_solver` 等函数目前为占位符，实际使用时需替换为调用 Fluent Journal 文件（.jou）或 PyFluent API 的具体指令，以实现对几何模型中水平管段入口速度的动态修改及 U 型水封区域初始状态的重置 [ref_2]。 ### 具体实施步骤与代码落地方案 要将此逻辑转化为可执行的工程方案，需将其拆解为三个具体环节：编写 Fluent Journal 文件、完善 Python 驱动脚本、以及配置运行环境。 #### 1. 编写 Fluent Journal 文件 (核心执行单元) Journal 文件是 Fluent 的批处理命令集，负责完成单次计算的设置与运行。你需要创建一个名为 `run_case.jou` 的文件，内容如下： ```journal /file/read-case "water_seal_model.cas" ! 读取包含网格和物理模型设置的案例文件 /define/boundary-conditions/velocity-inlet inlet-1 ! 假设入口边界名为 inlet-1 yes ! 确认修改 velocity-magnitude 0 ! 此处数值将被 Python 脚本动态替换 no ! 不再修改其他参数 /solve/initialize/patch ! 执行补丁初始化，将 U 型管区域设为水 water ! 选择相：水 zone-id 5 ! 假设 U 型管区域 ID 为 5 value 1 ! 体积分数设为 1 (充满水) yes /solve/dual-time-iterate 1000 10 ! 瞬态计算：总步数 1000，每步迭代 10 次 /file/write-data "result_data.dat" ! 输出结果数据 /exit ! 退出 Fluent ``` #### 2. 完善 Python 驱动脚本 (自动化控制) 利用 `subprocess` 模块调用系统命令执行 Fluent，并使用 `pyfluent` 或文本解析技术读取结果。以下是结合实际情况重构的可运行代码框架： ```python import os import subprocess import time import re # 配置参数 fluent_exe = "ansys_inc/v231/fluent/bin/fluent" # Fluent 执行路径，需根据实际安装修改 case_file = "water_seal_model.cas" # 基础案例文件 journal_template = "run_case_template.jou" # Journal 模板文件 initial_water_vol = 0.05 # 假设初始水体积 (m^3)，需根据网格计算预估 def modify_journal(input_vel, output_jou): """ 读取模板并替换入口速度值，生成新的 Journal 文件 对应原逻辑中的 update_boundary_condition """ with open(journal_template, 'r', encoding='utf-8') as f: content = f.read() # 使用正则替换 velocity-magnitude 后的数值 # 注意：实际 Journal 语法可能需要更精细的行替换，此处为简化逻辑 updated_content = re.sub(r'velocity-magnitude\s+\d+', f'velocity-magnitude {input_vel}', content) with open(output_jou, 'w', encoding='utf-8') as f: f.write(updated_content) print(f"[INFO] 已生成工况文件：{output_jou}, 入口速度={input_vel} m/s") def run_fluid_simulation(jou_file): """ 调用系统命令启动 Fluent 进行计算 对应原逻辑中的 run_cfd_solver """ # 命令示例：3d 模式，单精度，无图形界面 (-g)，加载 journal 文件 (-i) cmd = f"{fluent_exe} 3d -g -i {jou_file}" print(f"[RUN] 正在启动求解器...") try: # 实际生产中需捕获 stdout 以监控进度 subprocess.run(cmd, shell=True, check=True) print(f"[DONE] 计算完成") return True except subprocess.CalledProcessError: print(f"[ERROR] 计算失败") return False def parse_results(data_file): """ 解析结果文件，提取 U 型管内剩余水的体积 对应原逻辑中的 get_outlet_water_volume_fraction 注：实际需在 Fluent 中导出特定报告或通过 UDF 写入文件 """ # 此处仅为模拟逻辑，实际需解析 .dat 文件或读取 Fluent 导出的 report.txt # 假设我们在 Fluent 中通过 Surface Integral 输出了剩余水体积到 result_summary.txt try: with open("result_summary.txt", 'r') as f: line = f.readline() vol = float(line.strip()) return vol except FileNotFoundError: print("[WARN] 未找到结果文件，默认视为未排空") return initial_water_vol * 0.9 def main_search_loop(): min_vel = 5.0 max_vel = 30.0 step = 2.0 current_vel = min_vel print("="*30) print("开始 U 型水封临界吹扫气量自动搜索") print("="*30) while current_vel <= max_vel: jou_name = f"run_{current_vel}.jou" # 1. 生成当前工况的 Journal modify_journal(current_vel, jou_name) # 2. 执行计算 if not run_fluid_simulation(jou_name): break # 3. 获取结果 (需在 Fluent Journal 中添加导出剩余水量的命令) remaining_vol = parse_results("result_data.dat") print(f"-> 速度 {current_vel} m/s: 剩余水量 {remaining_vol:.4f} m^3") # 判定准则：剩余水量小于初始量的 1% if remaining_vol < (initial_water_vol * 0.01): print(f"\n[SUCCESS] 找到临界气速！约为 {current_vel} m/s") print(f"此时 U 型水封已被完全清除。") return current_vel else: print(f"-> 速度不足，继续增加气流速度...\n") current_vel += step print("\n[FAIL] 在设定范围内未找到能完全清除水封的气速。") return None if __name__ == "__main__": # 确保在工作目录下存在必要的 .cas 文件和模板 if not os.path.exists(case_file): print(f"错误：找不到案例文件 {case_file}") else: main_search_loop() ``` #### 3. 关键操作细节与数据溯源 在实际操作中，必须在 Fluent 的 Journal 模板中加入**报告定义（Report Definitions）**，以便在计算结束时自动输出 U 型管区域内的水相体积积分值，这是脚本判断“是否清除”的唯一依据 [ref_5]。若缺少此步骤，Python 脚本将无法获取 `remaining_vol` 变量。此外，对于 VOF 模型，时间步长（Time Step Size）的设置至关重要，需满足库朗数（Courant Number）小于 1 以保证界面捕捉的稳定性，这直接影响计算结果的可靠性 [ref_2]。通过这种自动化流程，可以系统地整理出如前文所述的“入口流速 - 排水效率”关系表，为毕业论文提供详实的数据支撑 [ref_1]。 ### 总结 该代码的使用本质是构建一个**参数化扫描平台**。用户需先完成单次成功的手动模拟，导出 `.cas` 文件并录制标准的 Journal 操作流程，然后将上述 Python 脚本中的路径、文件名及解析逻辑与实际工程文件对齐。执行后，程序将自动遍历预设的速度区间，精准定位清除 U 型水封的最小气量，从而高效完成毕设中关于“几何参数和气液流动参数影响规律”的研究目标 [ref_3]。