ANSYS Fluent 流体仿真全流程详解与自检表

一、仿真流程概述

ANSYS Fluent 是一款功能强大的流体力学仿真软件,广泛应用于航空航天、汽车、能源、建筑等多个领域。本文将详细介绍使用 SpaceClaim 建模、Fluent Meshing 划分网格,以及 Fluent 进行求解的完整流程,并提供一份全面的自检表,帮助快速掌握流体仿真的关键步骤和注意事项。

1.1 仿真流程基本框架

Fluent 仿真的基本流程主要包括三个关键环节:

  1. 几何建模:使用 SpaceClaim 创建或处理几何模型

  2. 网格划分:使用 Fluent Meshing 生成计算网格

  3. 求解计算:在 Fluent 中设置物理模型、边界条件并进行求解

这三个环节相互关联,任何一个环节出现问题都可能导致整个仿真失败或结果不准确。因此,每个步骤都需要严格按照规范操作,并进行必要的检查。

1.2 软件版本说明

本文基于 ANSYS 2025R1 版本进行讲解,该版本中已经将 SCDM (ANSYS SpaceClaim Direct Modeler) 的功能集成到了 Discovery 中,但基本命令操作没有大的变化,只是界面有所调整,需要一定时间适应。

二、SpaceClaim 几何建模详解

2.1 启动与基本设置

第一步:启动 ANSYS Workbench

在 Windows 开始菜单中找到 ANSYS 2025R1 文件夹,点击ANSYS Workbench 图标启动软件。Workbench 是 ANSYS 的项目管理平台,可以方便地管理多个仿真项目。

第二步:创建新项目并添加 Geometry 模块

  1. 在 Workbench 界面左侧的 Toolbox 中找到 “Component Systems” 下的 “Geometry” 模块

  2. 点击并拖动 “Geometry” 模块到右侧的 Project Schematic 区域

  3. 右键点击 Geometry 模块,选择 “Properties” 可以设置项目名称和保存路径

  4. 双击 Geometry 模块的 "Setup"项,启动 SpaceClaim 进行几何建模

第三步:设置单位系统

在 SpaceClaim 中,正确的单位设置至关重要:

  1. 选择 “Units” 菜单(通常位于界面右下角)
  2. 根据模型实际尺寸选择合适的单位,如毫米 (mm)、米 (m) 或英寸 (in)
  3. 确保所有后续操作使用统一的单位系统,避免单位转换错误

设置单位:文件→SpaceClaim选项→单位→选择合适的长度单位(通常为mm或m) SpaceClaimSpaceClaim

2.2 几何模型创建与编辑

基本几何体创建

SpaceClaim 提供了多种创建基本几何体的方法:

  1. 立方体创建:选择 “Create” 选项卡下的 “Box” 命令,在工作区拖动鼠标即可创建立方体。也可以通过输入精确的尺寸参数来创建。

  2. 圆柱体创建:选择 “Create” 选项卡下的"Cylinder" 命令,设置圆柱的轴线、半径和高度。同样支持精确尺寸输入。

  3. 球体创建:选择 “Create” 选项卡下的 “Sphere” 命令,定义中心点和半径即可创建球体。

几何编辑与变换

创建基本几何体后,需要进行编辑和变换以满足实际需求:

  1. 移动操作:选择几何体,激活 “Move” 工具,然后拖动几何体到新的位置。也可以通过输入精确的坐标值来移动。

  2. 旋转操作:选择几何体,激活 “Rotate” 工具,拖动旋转轴到期望的角度。同样支持输入精确的旋转角度。

  3. 缩放操作:选择几何体,激活 “Scale” 工具,拖动缩放把手来改变大小。可以按比例缩放或自定义各方向的缩放比例。

核心建模工具与技巧

四个基础工具(处理80%的建模任务):

拉伸工具(Pull Tool - P键)

  • 拉伸2D草图创建3D实体或曲面
  • 偏移面以增厚或减薄几何体
  • 旋转创建旋转特征
  • 沿路径扫掠
  • 添加拔模角度
  • 创建平滑过渡 SpaceClaim

操作流程:

  1. 选择面、边或草图
  2. 激活拉伸工具(P键或功能区) spaceclaim
  3. 拖动手柄或输入精确值
  4. 使用工具指南进行特定操作(倒角、圆角等) SpaceClaim

移动工具(Move Tool - M键)

  • 沿轴平移几何体
  • 围绕轴或点旋转
  • 复制移动(Ctrl+拖动)
  • 精确坐标定位

填充工具(Fill Tool - F键)

  • 移除不需要的特征(孔、圆角、倒角)
  • 修复曲面间隙
  • 简化复杂几何体以便网格划分
  • 清理导入的CAD数据

组合工具(Combine Tool - I键)

  • 添加:合并多个实体
  • 减去:从一个实体切除另一个
  • 相交:仅保留交集
  • 分割:切分实体并保留两部分
高级建模技术

对于复杂模型,需要使用高级建模技术:

  1. 布尔运算:包括并集 (Union)、交集 (Intersect) 和差集 (Subtract)。通过布尔运算可以将多个简单几何体组合成复杂形状。

  2. 曲面建模:使用 "Surface"工具创建和编辑曲面,适用于创建有机形态或复杂的自由形态设计。曲面编辑包括曲面融合、修剪和延伸等操作。

  3. 参数化设计:使用 “Variables” 表定义关键参数,如长度、半径等,然后将这些参数应用于几何体的尺寸定义。参数化设计便于后续修改和优化。

主要界面元素:
  • 设计窗口:显示3D模型的主工作区,包含2D草图网格
  • 功能区(Ribbon):包含所有建模工具,按设计(Design)、准备(Prepare)等选项卡组织
  • 结构树(Structure Panel):层次化显示所有设计对象,支持可见性控制
  • 图层面板:管理对象分组和视觉属性
  • 选择面板:提供高级几何选择工具
  • 组面板(Groups Panel):创建CFD边界条件的命名选择 SpaceClaimspaceclaim

2.3 几何模型简化与修复

几何简化

对于从 CAD 软件导入的复杂模型,通常需要进行简化:

  1. 选择 “准备” 菜单中的 “简化” 功能

  2. 根据模型复杂程度选择简化级别:

  • 级别 0 (边界框):最简单的简化方式

  • 级别 1 (长方体 / 圆柱体拟合):适用于简单几何

  • 级别 2 (多边形拟合):适用于中等复杂程度几何

  • 级别 3 (CAD 对象):最复杂但最准确的简化方式,适用于复杂模型

几何简化建议:

  • 移除不必要的小特征(文字、标志、小圆角)
  • 填充不影响流动的小孔
  • 使用填充工具去特征化以提高网格质量
  • 关闭间隙和重叠

SpaceClaim中的共享拓扑(Share Topology)是一个重要的几何处理功能,主要作用包括:

主要功能

确保几何连续性 共享拓扑能够识别并合并相邻实体之间的公共边界和顶点,确保模型各部分之间保持几何连续性。这对于后续的网格划分和仿真分析至关重要。

简化网格生成 当启用共享拓扑后,相邻实体在公共界面上会自动共享节点和单元,避免了网格不匹配的问题。这样可以确保载荷和边界条件在界面上正确传递。

提高仿真精度 通过共享拓扑,可以消除模型中的小间隙、重叠和不连续区域,这些缺陷可能会导致仿真结果不准确或求解失败。共享的界面确保了应力、热量等物理量能够正确地在部件之间传递。

典型应用场景

装配体分析 在分析包含多个零件的装配体时,共享拓扑可以自动处理零件之间的接触面,无需手动定义接触关系,大大简化了前处理工作。

多体零件处理 对于由多个实体组成的零件,共享拓扑能够确保这些实体被视为一个整体进行分析,而不是独立的部分。

CAD模型修复 导入的CAD模型经常存在微小的间隙或重叠,共享拓扑功能可以自动检测并修复这些问题,使模型适合进行有限元分析。

使用共享拓扑时,SpaceClaim会根据指定的容差值自动识别应该共享的几何特征,这个容差值可以根据模型的尺度和精度要求进行调整。

几何修复

修复几何模型中的常见问题:

  1. 面片修复:使用 " 修复"菜单中的" 面片修复 " 功能,修复模型中的破面、裂缝等问题。

  2. 间隙检测:选择 “检查” 菜单中的 “间隙检测”,找出模型中的小间隙,并设置合适的间隙公差进行修复。

  3. 小型面处理:使用 “修复” 菜单中的 “小型面” 功能,删除或合并模型中的小面,提高网格划分成功率。

2.4 开口与流域创建

创建流体流域

方法1:封闭工具(Enclosure Tool)

  1. 进入准备(Prepare)选项卡
  2. 选择封闭(Enclosure)工具
  3. 点击并拖动覆盖整个几何体创建边界框 Verus Engineeringverus-engineering
  4. 封闭自动包围所有可见几何体 Verus Engineering
  5. 隐藏原始几何体,仅处理流体域 Verus Engineering
  6. 使用拉伸工具根据需要扩展域边界 Ansys +3

方法2:体积提取工具(Volume Extract Tool) 用于内部流动分析:

  1. 打开准备选项卡→体积提取 SpaceClaimspaceclaim
  2. 点击选择边缘工具指南 SpaceClaimspaceclaim
  3. 选择包围流体区域的边缘环 SpaceClaimspaceclaim
  4. 双击进行连续边缘选择 SpaceClaimspaceclaim
  5. 如需要,使用选择种子面识别内部体积 SpaceClaimspaceclaim
  6. 点击完成创建流体体 SpaceClaimspaceclaim
  7. 在结构树中重命名体积 CFD.NINJA +2
设置进出水口

在流体域上创建开口作为边界条件:

  1. 选择 “开口” 工具(位于 “准备” 菜单中)
  2. 选择需要设置为开口的面
  3. 输入开口名称,如 “inlet” 或 “outlet”
  4. 对于复杂模型,可以使用 “按平面分割” 工具创建内部流道

创建命名选择(Named Selections):

  1. 组面板→创建NS(命名选择) Accessengineeringlibrary
  2. 选择边界条件的面(入口、出口、壁面)
  3. 适当命名(如"inlet"、“outlet”、“wall”、“symmetry”) Accessengineeringlibrary
  4. 这些将直接转换为Fluent边界条件 Ansys Innovation SpaceAnsys Learning Forum

2.5 几何模型检查与导出

几何检查清单

在完成几何建模后,需要进行以下检查:

  1. 模型尺寸检查:确保所有尺寸与设计要求一致,可以通过 “测量” 工具进行验证。

  2. 模型连续性检查:检查模型是否有间隙、重叠或不连续的部分,这些问题会导致网格划分失败。

  3. 流体域完整性检查:确保流体域是封闭的,没有遗漏的开口或孔洞。

模型导出

完成检查后,将几何模型导出为 Fluent Meshing 可以识别的格式:

  1. 选择 “文件” 菜单中的 “另存为”

  2. 选择合适的文件格式,推荐使用 ANSYS 默认的.agdb 格式

  3. 保存模型到项目文件夹中,便于后续操作访问

三、Fluent Meshing 网格划分详解

3.1 网格划分前的准备工作

启动 Fluent Meshing

在 ANSYS Workbench 中启动 Fluent Meshing 有两种方法:

  1. 在 Workbench 的 Project Schematic 中,右键点击 Geometry 模块,选择 “Transfer Data To New”->“Fluent Meshing”
  2. 双击新生成的 Fluent Meshing 模块的 “Setup” 项,启动 Fluent Meshing

导入几何模型

在 Fluent Meshing 中导入已创建的几何模型:

  1. 选择 “Import Geometry” 任务
  2. 在 “File Format” 中选择几何文件格式(如.agdb)
  3. 点击 “Browse” 找到之前保存的几何模型文件
  4. 设置单位与 SpaceClaim 中保持一致
  5. 点击 “Update” 加载几何模型

3.2 网格划分基本流程

Fluent Meshing 提供了多种网格划分流程,推荐使用 “Watertight Geometry” 工作流程,该流程适用于大多数情况,且不易遗漏关键步骤。

网格划分工作流选择

水密几何工作流(Watertight Geometry Workflow):

  • 适用于清洁、水密的CAD模型 Nimbix
  • 自动CAD分面,可自定义公差
  • 支持体标签识别,便于高级几何组织
  • 几何缩放和单位转换选项 Ansys +4

容错工作流(Fault-Tolerant Workflow):

  • 用于复杂或有问题的几何体 Nimbix
  • 包含包裹技术处理脏CAD装配体 Ansys
  • 高级几何修复功能
  • 处理非水密几何体 AnsysNimbix
Watertight Geometry 工作流程
  1. 导入几何体:已在上一步完成
  2. 添加局部尺寸控制:设置全局和局部网格尺寸
  3. 生成表面网格:先生成表面网格,为体网格做准备
  4. 描述几何体:指定几何体类型(流体、固体或混合)
  5. 封闭流体区域:自动封闭开口,创建流体区域
  6. 更新区域:确认和更新流体区域和固体区域
  7. 添加边界层控制:对壁面附近添加边界层网格
  8. 生成体网格:最终生成三维体网格

3.3 网格尺寸控制

全局网格尺寸设置

全局网格尺寸决定了整体网格密度:

在 “Add Local Sizing” 任务中选择 “Yes” 添加局部尺寸控制

在 “Size Control Type” 中选择 “Global”

设置 “Target Mesh Size” 为合适的全局网格尺寸

通常,全局网格尺寸应根据模型特征尺寸和计算精度要求来确定,初始值可设为特征尺寸的 1/10 到 1/20

局部网格加密

对于关键区域,需要进行局部网格加密:

  1. 边尺寸控制:选择 “Edge Size” 控制类型,对指定的棱线进行加密。适用于需要精确捕捉流动变化的边,如模型的锐边或曲率大的区域。
  2. 面尺寸控制:选择 “Face Size” 控制类型,对指定的面进行加密。适用于边界层、狭窄通道或流动分离区域。
  3. 体尺寸控制:选择 “Body Size” 控制类型,对指定的体进行加密。适用于整个流体域或特定的子区域。
  4. 基于影响的尺寸控制
    • Body of Influence (BOI):通过其他实体表征加密区域,BOI 是封闭的三维体,规定了其边界内的网格区域的最大单元尺寸。
    • Face of Influence (FOI):通过面来控制附近的网格尺寸,FOI 规定了其表面附近的最大单元尺寸。
网格尺寸控制参数

在设置网格尺寸时,需要考虑以下参数:

  1. 单元尺寸:最大允许全局单元尺寸
  2. Growth Rate:网格尺寸的增长率,默认为 1.2,值越大,网格增长越快,但可能影响计算精度。
  3. Min Cell LengthMax Cell* Length(最大/最小尺寸):设置网格单元的最小和最大尺寸限制,避免出现过小或过大的单元。
  4. 高级尺寸函数(ASF)
    1. 基于曲率的细化:使用"曲率法向角"(典型值:6-12°)
    2. 基于邻近度的细化:每个间隙的单元参数(典型值:2-4) Computationalfluiddynamics

3.4 网格类型选择与设置

网格类型概述

Fluent Meshing 支持多种网格类型,选择合适的网格类型对仿真结果至关重要:

四面体网格(Tetrahedral)

  • 标准非结构化方法
  • 适合复杂几何体
  • 单元数量较高但自动生成容易

适用于复杂几何,易于生成,但计算精度相对较低。


六面体/六核网格(Hexcore)

  • 体积区域采用结构化核心
  • 最佳计算效率
  • 复杂几何体具有挑战性

适用于规则几何,计算精度高,但生成难度较大。


混合网格 (Hybrid)

结合四面体和六面体网格的优点,在复杂区域使用四面体网格,在规则区域使用六面体网格。


多面体网格(Polyhedral)

  • 相比四面体减少3-5倍单元数量 Ansys
  • 更好的收敛特性
  • 可从四面体网格自动转换

适用于复杂几何,单元数量少且计算精度高,是近年来推荐使用的网格类型。


Poly-Hexcore网格(Mosaic技术)

  • 创新混合方法:
    • 体积区域使用八叉树六面体
    • 多面体单元提供几何灵活性
    • 高质量棱柱层用于边界层 ansys +2
  • 优势:

薄体网格划分

对于薄体结构(如泡棉、隔板等),Fluent Meshing 2025R1 引入了新的薄体网格划分功能:

  1. 在完成面网格划分后,右键选择 “Add Thin Volume Meshing Controls”
  2. 设置薄体的层数 (Number of Intervals) 和增长率 (Growth Rate)
  3. 选择需要划分薄体网格的区域
  4. 选择源面和目标面,确保网格正确对齐
  5. 点击 “Update” 生成薄体网格
边界层网格设置

在壁面附近添加边界层网格以捕捉粘性效应:

  1. 选择 “Boundary Layers” 任务
  2. 选择需要添加边界层的壁面
  3. 设置边界层的层数、总厚度和增长率
  4. 对于湍流流动,通常需要至少 3 层边界层网格,以确保 y + 值在合理范围内

第一层高度:用于直接y+控制

  • 壁面函数:y+ = 30-100
  • 增强壁面处理:y+ < 5(理想y+ ≈ 1)
  • 低雷诺数模型:y+ < 1

层数:最少10-20层用于边界层解析 Ansys Learning Forum

增长率:典型值1.1-1.3(高雷诺数流动最大1.2)

过渡比:控制到核心网格的体积变化

3.5 网格质量检查

网格质量指标

生成网格后,需要检查以下质量指标:

  1. 正交性 (Orthogonal Quality**):衡量网格单元各边之间的正交程度,理想值为 1,一般要求大于 0.2。
  2. 偏斜度 (Skewness):衡量网格单元偏离理想形状的程度,理想值为 0,一般要求小于 0.9。
  3. 纵横比 (Aspect Ratio):网格单元最长边与最短边的比值,一般要求不超过 5:1,在边界层区域可适当放宽至 10:1。
  4. 雅可比比率 (Jacobian Ratio)**:衡量网格单元变形程度的指标,一般要求大于 0.6。
质量指标优秀良好可接受
偏斜度(Skewness)0-0.250.25-0.50.5-0.8>0.95
正交质量(Orthogonal Quality)0.95-1.00.7-0.950.1-0.7<0.1
纵横比(Aspect Ratio)<1010-5050-100>100
网格质量检查方法

在 Fluent Meshing 中检查网格质量:

  1. 选择 “Check Mesh” 工具,查看网格统计信息,包括单元总数、节点数、单元类型分布等

  2. 使用 “Mesh Statistics” 查看详细的网格质量指标分布

  3. 使用 “Mesh Quality” 工具可视化网格质量分布,识别低质量区域

  4. 对于质量较差的区域,需要重新调整网格尺寸参数或进行局部网格优化

网格优化

改善网格质量的方法:

  1. 局部网格细化:对低质量区域进行局部网格细化

  2. 网格平滑:使用 “Smooth” 工具对网格进行平滑处理

  3. 网格交换:使用 “Swap” 工具交换边或面,改善网格质量

  4. 重新划分:对于严重低质量区域,考虑重新划分网格

3.6 网格导出与保存

保存网格

在完成网格划分并通过质量检查后,需要保存网格:

  1. 选择 “File” 菜单中的 “Save As”

  2. 输入网格文件名,推荐使用.msh 格式

  3. 选择保存路径,建议与项目文件放在同一目录下

  4. 点击 “Save” 保存网格文件

导出到 Fluent

在 Fluent Meshing 中直接导出网格到 Fluent 求解模块:

  1. 选择 “File” 菜单中的 “Export to Fluent”

  2. 选择需要导出的网格区域

  3. 设置导出选项,如单位、精度等

  4. 点击 “OK” 导出网格到 Fluent

四、Fluent 求解设置与计算

4.1 启动 Fluent 并导入网格

启动 Fluent

在 ANSYS Workbench 中启动 Fluent 求解模块:

  1. 在 Project Schematic 中,右键点击 Fluent Meshing模块,选择 “Transfer Data To New”->“Fluent”

  2. 双击新生成的 Fluent 模块的 “Setup” 项,启动 Fluent 求解器

导入网格

在 Fluent 中导入已划分好的网格:

  1. 选择 “File” 菜单中的 “Read”->“Mesh”

  2. 找到之前保存的网格文件(.msh 格式)

  3. 点击 “OK” 加载网格

  4. 网格加载后,Fluent 会自动显示网格信息,包括节点数、单元数和区域数等

检查网格

在开始求解前,需要检查网格:

  1. 选择 “General” 面板中的 “Check” 按钮,检查网格完整性和一致性

  2. 查看网格统计信息,确保最小体积和面积为正值

  3. 检查单位是否正确,必要时使用 “Scale” 功能调整单位

  4. 使用 “Display” 功能可视化网格,确认模型正确无误

4.2 基本求解设置

求解器选择

Fluent 提供两种主要的求解器类型:

  1. 压力基求解器 (Pressure-Based):适用于不可压缩或低马赫数可压缩流动,求解效率高,是默认选择。

  2. 密度基求解器 (Density-Based):适用于高马赫数可压缩流动,如超音速或高超音速流动。

选择求解器的基本原则:

  • 对于不可压缩流动或 Ma<0.3 的可压缩流动,选择压力基求解器

  • 对于 Ma≥0.3 的可压缩流动,选择密度基求解器

  • 对于冲击射流等特殊情况,可选择密度基显式求解器

时间类型设置

根据问题特性选择时间类型:

  1. 稳态 (Steady):适用于流场不随时间变化的情况,计算资源需求较少。

  2. 瞬态 (Transient):适用于流场随时间变化的情况,如非定常流动、周期性流动等,需要设置时间步长和总时间。

速度公式设置

在压力基求解器中,需要设置速度公式:

  1. 绝对速度 (Absolute):速度相对于全局坐标系

  2. 相对速度 (Relative):速度相对于邻近单元区域,适用于旋转坐标系或移动参考系问题

4.3 物理模型设置

湍流模型选择

根据流动特性选择合适的湍流模型:

  1. 层流模型 (Laminar):适用于低雷诺数流动或精确已知为层流的情况。

  2. k-ε 模型(有限使用):标准 k-ε 模型适用于高雷诺数的完全湍流流动,计算效率高,但对复杂流动预测能力有限。

    • 标准k-ε:仅用于向后兼容,需要激活生产限制器
    • Realizable k-ε:可实现性限制器效果有限
    • RNG k-ε:仅用于向后兼容
  3. k-ω 模型(大多数应用推荐):标准 k-ω 模型适用于边界层和近壁区域,但对自由剪切流预测效果较差。

    • SST模型:推荐用于大多数工业应用,对逆压梯度和分离流具有高精度

    • BSL模型:SST的替代方案,当SST预测过度分离时使用

    • GEKO模型:最通用,设计用于替代所有其他两方程模型

  4. SST k-ω 模型:结合了 k-ε 和 k-ω 模型的优点,适用于边界层和分离流动,是工程应用中最常用的湍流模型之一。

  5. 雷诺应力模型 (RSM):考虑了雷诺应力的各向异性,适用于复杂流动和强旋流,但计算资源需求大。

[!TIP]

选择指南:

  • 分离流、逆压梯度使用SST或GEKO
  • 外部空气动力学使用Spalart-Allmaras
  • GEKO可模拟其他模型:CSEP=1.0模拟k-ε行为,默认设置模拟SST行为
多相流模型设置

如果涉及多相流,需要选择合适的多相流模型:

  1. VOF (Volume of Fluid) 模型:适用于不混溶流体的界面追踪,如气液两相流。

  2. Mixture 模型:适用于相之间有较强耦合但不需要精确追踪界面的情况。

  3. Eulerian 模型:适用于多相流中各相均为连续相的情况,计算复杂度较高。

能量方程设置

对于涉及热量传递的问题,需要启用能量方程:

  1. 选择 “Models”->“Energy”

  2. 勾选 “Energy Equation” 启用能量方程

  3. 设置相关参数,如导热系数、比热容等

4.4 材料属性设置

定义材料

在 Fluent 中定义流体和固体材料:

  1. 选择 “Materials” 面板

  2. 点击 “Create/Edit” 按钮打开材料编辑器

  3. 可以选择 Fluent 内置材料库中的材料,也可以自定义新材料

  4. 对于自定义材料,需要设置密度、粘度、导热系数、比热容等属性

设置流体材料

以空气为例,设置流体材料属性:

  1. 在材料编辑器中选择 “air”

  2. 设置密度为理想气体 (ideal-gas),适用于可压缩流动

  3. 设置粘度为常数或使用 Sutherland 定律(适用于温度变化较大的情况)

  4. 设置导热系数和比热容为常数或温度函数

  5. 点击 “Change/Create” 保存设置

设置固体材料

如果模型包含固体区域,需要设置固体材料属性:

  1. 在材料编辑器中选择或创建固体材料

  2. 设置密度、导热系数和比热容等属性

  3. 对于多孔介质,还需要设置渗透率和粘性阻力系数等参数

4.5 边界条件设置

入口边界条件

常见的入口边界条件包括:

速度入口 (Velocity Inlet):指定入口速度和湍流参数,适用于已知入口速度的情况。

设置参数:

  • 速度大小和方向
  • 湍流强度(典型1-10%)
  • 湍流粘度比(典型1-10)
  • 水力直径 (Hydraulic Diameter)
  • 对于热流问题,还需设置入口温度或热流密度
  • 如适用,指定组分质量分数

压力入口 (****Pressure Inlet):指定入口总压和静压,适用于入口压力已知的情况。

设置参数:

  • 总压和静压
  • 湍流参数
  • 流动方向规格
  • 对于可压缩流动,需设置总温
  • 当质量流量随迭代变化时使用

质量流量入口 (Mass-Flow Inlet):指定入口质量流量,适用于质量流量已知的情况。

  • 直接质量流量规格
  • 入口均匀分布
  • 湍流参数如上
出口边界条件

常见的出口边界条件包括:

压力出口 (Pressure Outlet**)**:指定出口静压,适用于出口压力已知的情况。

设置参数:

  • 静压
  • 回流条件(当流动反向时的参数)

自由出口 (Outflow):适用于出口流动充分发展,且出口压力未知的情况。

  • 零扩散通量边界
  • 应放置在流动充分发展的位置
  • 出口压力未知时使用

质量流量出口 (Mass-F****low Outlet):指定出口质量流量,适用于出口质量流量已知的情况。

壁面边界条件

壁面边界条件设置:

  1. 无滑移条件 (No Slip):默认设置,壁面处流体速度为零
  2. 指定剪切应力 (Specified Shear):适用于需要指定壁面剪切应力的特殊情况
  3. 绝热壁面 (Adiabatic Wall):默认设置,壁面热通量为零
  4. 指定热通量 (Specified Heat Flux):适用于已知壁面热通量的情况
  5. 指定壁面温度 (Specified Wall Temperature):适用于已知壁面温度的情况
特殊边界条件

其他特殊边界条件:

  1. 对称边界 (Symmetry):适用于几何和流动对称的情况,可以减少计算量
  2. 周期性边界 (Periodic):适用于周期性流动,如旋转机械中的叶片通道
  3. 内部界面 (Internal):用于连接两个相邻的流体区域或固体区域

4.6 求解控制设置

求解方法设置

在 “Solution Methods” 面板中设置求解方法:

压力 - 速度耦合算法:在压力基求解器中,常用的耦合算法有:

  • SIMPLE 算法:稳定性好,收敛速度适中,是默认选择
  • SIMPLEC 算法:收敛速度比 SIMPLE 快,适用于复杂几何
  • PISO 算法:适用于瞬态问题或需要快速收敛的情况
方案特点应用场景亚松弛因子
SIMPLE稳健但收敛较慢稳态流动默认压力0.3,动量0.7
SIMPLEC比SIMPLE收敛更好简单问题可使用更高因子
PISO推荐用于瞬态计算高度倾斜网格因子接近1.0
Coupled稳态最佳收敛需要快速收敛Courant数200
空间离散格式

压力格式:

  • 标准:大多数应用的默认值
  • PRESTO!:旋流、高瑞利数、VOF
  • 二阶:需要更高精度时

动量和标量:

  • 一阶迎风:最稳定,精度较低
  • 二阶迎风:默认,精度/稳定性平衡良好
  • QUICK:更高精度,需要有界值

[!TIP]

  • 梯度计算:通常选择最小二乘单元法 (Least Squares Cell Based)
  • 流动项:推荐使用二阶迎风格式 (Second Order Upwind),精度高于一阶
  • 湍流项:同样推荐使用二阶迎风格式
求解控制参数

在 “Solution Controls” 面板中设置求解控制参数:

松弛因子 (Under-Relaxation Factors):控制迭代过程的收敛稳定性

  • 压力 (Pressure):默认值 0.3,可根据收敛情况调整

  • 密度 (Density):对于可压缩流动需要调整

  • 动量 (Momentum):默认值 0.7,可适当增大以加快收敛

  • 其他方程(如 k、ω 等):默认值 0.8-1.0


库朗数 (Courant Number):仅在密度基求解器或瞬态问题中出现

  • 库朗数越高,收敛速度越快,但稳定性可能降低

  • 初始计算时可设置较低库朗数(如 5-10),待收敛稳定后逐步提高至 50-100

4.7 初始化与求解

初始化方法

标准初始化 (Standard Initialization)

使用简单的初始条件填充整个计算域,适用于大多数情况。


混合初始化 (Hybrid Initialization)

基于网格信息生成更合理的初始场,通常能更快收敛,推荐使用。

  • 求解压力和速度的拉普拉斯方程

  • 比标准初始化更真实的初始场

  • 基于域平均值自动修补变量

  • 推荐用于大多数应用


完全多重网格(FMG)初始化

  • 使用多重网格技术获得优越的初始解

  • 对大型问题计算成本低 ANSYS +2

  • 过程:

    1. 切换到基于密度的显式求解器
    2. 在多个网格级别上求解无粘欧拉方程 ANSYS
    3. 从粗网格开始,插值到细网格
    4. 返回原始求解器设置 ANSYS

分区技术(Patching):

  • 用不同值初始化特定区域
  • 对复杂流动有用(热点、物种浓度)
  • 可按区域、单元寄存器或几何区域分区

初始化步骤

  1. 选择 “Solution Initialization” 面板
  2. 选择初始化方法(推荐混合初始化)
  3. 设置初始化参数,如参考压力、温度等
  4. 点击 “Initialize” 进行初始化
  5. 初始化完成后,可以通过 “Patch” 功能对特定区域进行初始值调整
残差监控设置
路径方式:
  • 主菜单:Solution → Monitors → Residual
  • 功能区:Solution标签 → Monitors组 → Residual → Edit...
  • 树形菜单:Solution → Monitors → Residual → 右键Edit...
  • TUI命令:/solve/monitors/residual
残差监控对话框详细设置

方程列表区域(Equations):

方程类型默认监控状态典型收敛标准物理含义
Continuity(连续性)10⁻³质量守恒
x-velocity10⁻³X方向动量守恒
y-velocity10⁻³Y方向动量守恒
z-velocity10⁻³Z方向动量守恒(3D)
Energy10⁻⁶能量守恒
k (turbulent kinetic energy)10⁻³湍动能
epsilon/omega10⁻³湍流耗散率
Species-n10⁻⁶组分守恒
DO Intensity10⁻⁶辐射强度

每个方程的控制选项:

  1. Monitor(监控)复选框

    • 功能:控制是否监控该方程的残差
    • 默认:所有求解的方程均开启
    • 操作:取消勾选将停止该方程的残差监控(但仍然求解)
    • 应用场景:
      • 某些辅助方程不需要严格收敛时可关闭
      • 减少残差图的混乱程度
  2. Check Convergence(检查收敛)复选框

    • 功能:控制是否将该方程纳入自动收敛判断
    • 默认:所有监控的方程均开启
    • 操作:取消勾选后该方程不参与收敛判断
    • 应用场景:
      • 某些方程波动较大但不影响整体解时
      • 想要手动控制收敛判断时
  3. Absolute Criteria(绝对收敛准则)

    • 功能:设置每个方程的收敛阈值

    • 输入范围:通常10⁻³到10⁻⁸

    • 判断逻辑:当残差降至该值以下时认为收敛

    • 推荐设置:

      [!TIP]

      连续性:1e-3 到 1e-4
      动量: 1e-3 到 1e-5
      能量: 1e-6 到 1e-8
      湍流: 1e-3 到 1e-5
      组分: 1e-6 到 1e-8
      辐射: 1e-6

高级选项设置(Show Advanced Options)
  1. Scale(缩放)选项

    • 功能:控制是否使用缩放残差
    • 默认:开启
    • 缩放方式:
      • Global Scaling(全局缩放):默认方式
        • 公式:R_scaled = R_raw / (a_p × φ)
        • 适用:大多数问题
      • Local Scaling(局部缩放):需额外开启
        • 每个单元独立缩放
        • 更严格的收敛判断
    • 何时关闭:
      • 需要查看原始残差值时
      • 对比不同算例时
  2. Normalize(归一化)选项

    • 功能:将残差归一化到初始值
    • 默认:关闭
    • 作用:使所有方程初始残差为1
    • 归一化公式:R_normalized = R_current / R_reference
    • 应用场景:
      • 判断残差下降了多少个数量级
      • 不同方程间的对比
  3. Iterations(迭代次数)

    • 功能:设置用于计算归一化因子的迭代次数
    • 默认值:5
    • 范围:1-100
    • 说明:取前N次迭代的最大残差作为归一化分母
    • 调整建议:
      • 稳定问题:5-10次
      • 振荡问题:20-50次
  4. Normalization Factor(归一化因子)

    • 功能:显示/手动设置每个方程的归一化分母
    • 默认:自动计算(前N次迭代的最大值)
    • 手动设置场景:
      • 重启计算时保持一致性
      • 特定的收敛研究
  5. Compute Local Scale(计算局部缩放)

    • 功能:启用基于单元的局部缩放
    • 影响:
      • 收敛准则自动调整为1e-5(稳态)或1e-4(瞬态)
      • 计算和存储开销增加
    • 报告选项(Reporting Option):
      • Global:显示全局缩放残差
      • Local:显示局部缩放残差
收敛准则类型(Convergence Criterion)

下拉菜单选项:

  1. absolute(绝对准则)

    • 默认选项(稳态)
    • 判断:残差 < 绝对值
    • 适用:大多数稳态问题
  2. relative(相对准则)

    • 判断:(R_current - R_min) / R_initial < 相对值
    • 适用:瞬态计算
    • 设置值:通常0.01-0.001
  3. iterate(迭代准则)

    • 基于迭代次数
    • 固定迭代次数后停止
    • 适用:参数研究
  4. none(无准则)

    • 不自动停止
    • 手动控制求解
    • 适用:复杂非线性问题
操作按钮功能
  1. Renormalize(重新归一化)
    • 功能:重新计算归一化因子
    • 使用时机:
      • 残差最大值出现在5次迭代之后
      • 切换缩放/非缩放模式后
      • 重启计算时
  2. Reset(重置)
    • 功能:清除残差历史
    • 影响:重新开始残差记录
    • 使用场景:
      • 更改物理模型后
      • 显著修改边界条件后
  3. Options…(选项)
    • 打开绘图选项对话框
    • 设置残差图显示属性
同问题类型的推荐设置

简单稳态层流问题:

连续性:1e-3
动量:  1e-3
能量:  1e-6(如果求解)
缩放:  开启(全局)
归一化:关闭

复杂湍流问题(RANS):

连续性:1e-4
动量:  1e-4
能量:  1e-6
k:     1e-4
ε/ω:   1e-4
缩放:  开启(全局)
归一化:开启(便于判断下降量级)

传热主导问题:

连续性:1e-3
动量:  1e-3
能量:  1e-8(严格)
缩放:  开启
归一化:开启
监控额外的传热率

多相流(VOF):

连续性:1e-3
动量:  1e-4
体积分数:1e-6
缩放:开启(局部)
相对准则:适用于瞬态

化学反应流:

连续性:1e-4
动量:  1e-4
能量:  1e-8
组分:  1e-8(每个组分)
缩放:  开启
监控关键组分生成率
求解执行
访问求解控制

路径:

  • 功能区:Solution → Solution → Methods/Controls
  • 树形菜单:Solution → Solution MethodsSolution → Controls
  • TUI:/solve/set/
压力基求解器控制(Pressure-Based Solver)

亚松弛因子设置(Under-Relaxation Factors):

变量默认值稳定范围激进范围调整策略
Pressure0.30.1-0.30.3-0.5发散时降低0.1
Momentum0.70.5-0.70.7-0.9与压力成反比
Turbulent Kinetic Energy0.80.5-0.80.8-1.0稳定后可提高
Turbulent Dissipation Rate0.80.5-0.80.8-1.0与k保持一致
Energy1.00.8-1.01.0传热问题可降至0.9
Species1.00.8-1.01.0反应流降至0.8
Density1.00.5-1.01.0可压缩流降低
Body Forces1.00.5-1.01.0自然对流降低

压力-速度耦合方案详细说明:

  1. SIMPLE(Semi-Implicit Method)

    • 特点:最稳健,收敛较慢
    • 迭代/时间步:50-200次
    • 亚松弛因子:压力0.3,动量0.7
    • 适用:
      • 首次运行的新问题
      • 网格质量一般
      • 稳态计算默认
  2. SIMPLEC(SIMPLE-Consistent)

    • 特点:比SIMPLE快30-50%
    • 迭代/时间步:30-150次
    • 亚松弛因子:压力0.3-0.5,动量0.7-0.9
    • 适用:
      • 网格质量好
      • 简单几何
      • 层流或简单湍流
  3. PISO(Pressure-Implicit Split-Operator)

    • 特点:瞬态计算优选
    • 迭代/时间步:20-50次
    • 邻域修正(Neighbor Correction):1-2
    • 偏斜修正(Skewness Correction):1(自动)
    • 适用:
      • 瞬态流动
      • 高度倾斜网格
      • VOF多相流
  4. Coupled(耦合求解器)

    • 特点:收敛最快,内存需求大(1.5-2倍)
    • 迭代次数:减少50-75%
    • 伪瞬态设置:
      • Courant Number:初始5-10,逐步增至200
      • 时间尺度因子:1.0
    • 适用:
      • 稳态问题需要快速收敛
      • 强耦合问题(高速流动)
      • 充足的计算资源

Pseudo Transient(伪瞬态)选项:

  • 功能:为稳态问题添加伪时间项
  • 优势:提高收敛稳定性
  • 设置:
    • 自动时间步长:推荐
    • 时间尺度因子:0.1-1.0
    • Length Scale Method:
      • Conservative(保守):0.1
      • Aggressive(激进):1.0
      • User-Specified(用户指定)
密度基求解器控制(Density-Based Solver)

Courant Number(库朗数)设置:

求解器类型初始值稳定值最大值调整策略
Explicit0.1-0.50.5-1.01.0保持<1.0
Implicit1.05.050-200指数增长

时间步进策略:

初始阶段(0-100迭代):CFL = 0.1-1.0
发展阶段(100-500):CFL = 1.0-10.0
收敛阶段(500+):CFL = 10.0-200.0
运行计算设置(Run Calculation)
稳态计算设置

基本参数:

  1. Number of Iterations(迭代次数)

    • 范围:100-10000

    • 初始运行:500-1000

    • 精细收敛:2000-5000

    • 设置策略:

      简单问题:500-1000
      中等复杂:1000-3000
      复杂问题:3000-5000
      使用监控点判断,不依赖固定次数
      
  2. Reporting Interval(报告间隔)

    • 默认:1(每次迭代)
    • 建议设置:
      • 快速计算:5-10
      • 长时间计算:10-50
    • 影响:控制台输出频率、监控更新频率
  3. Profile Update Interval(型线更新间隔)

    • 默认:1
    • 用途:更新随时间变化的边界条件
    • 设置:通常与报告间隔一致
瞬态计算设置

时间步设置决策树:

物理时间尺度计算
├── 特征长度 L
├── 特征速度 V
└── 物理时间尺度 τ = L/V

数值稳定性要求
├── 最小网格尺寸 Δx_min
├── 最大速度 V_max
└── CFL时间步 Δt_CFL = CFL × Δx_min/V_max

最终时间步
└── Δt = min(τ/10, Δt_CFL)

详细参数设置:

  1. Time Step Size(时间步长)

    • 计算公式:

      • 对流限制:Δt < Δx/V
      • 扩散限制:Δt < Δx²/ν
      • 源项限制:Δt < 1/S
    • 典型值:

      低速流(<50 m/s):1e-3 到 1e-2 秒
      中速流(50-200 m/s):1e-4 到 1e-3 秒
      高速流(>200 m/s):1e-5 到 1e-4 秒
      VOF自由面:1e-5 到 1e-4 秒
      化学反应:1e-6 到 1e-5 秒
      
  2. Number of Time Steps(时间步数)

    • 计算:总时间 / 时间步长

    • 示例:

      模拟1秒,Δt=1e-3:需要1000步
      模拟0.1秒,Δt=1e-4:需要1000步
      
  3. Max Iterations/Time Step(每时间步最大迭代)

    • SIMPLE/SIMPLEC:20-40

    • PISO:5-20

    • Coupled:10-30

    • 判断准则:

      理想:15-25次收敛
      过快(<10):可能时间步太小
      过慢(>30):时间步可能太大
      
  4. Time Stepping Method(时间推进方法)

    • Fixed(固定):恒定时间步

    • Adaptive(自适应)

      • 设置CFL目标值

      • 自动调整时间步

      • 参数:

        Truncation Error Tolerance: 0.01
        Ending Time: 物理总时间
        Minimum Time Step Size: 1e-8
        Maximum Time Step Size: 1e-2
        Minimum Step Change Factor: 0.5
        Maximum Step Change Factor: 1.5
        
求解过程监控

在求解过程中需要监控以下内容:

  1. 残差曲线:观察各方程的残差是否下降并趋于稳定
  2. 物理量监测:设置关键物理量(如进出口流量、压力降等)的监测,确保其趋于稳定
  3. 计算资源使用情况:监控 CPU 和内存使用情况,确保计算正常进行

求解过程监控(Solution Monitoring)

Report Definitions(报告定义)
创建报告定义

访问方法:

  1. 功能区:Solution → Reports → Definitions → New
  2. 树形菜单:右键Report Definitions
  3. 图形窗口:选择面后右键菜单
Surface Report(表面报告)类型详解

1. Area-Weighted Average(面积加权平均)

定义:Σ(φᵢ × Aᵢ) / Σ(Aᵢ)
用途:温度、压力、速度的平均值
设置示例:
  Name: outlet-temperature-avg
  Field Variable: Temperature
  Surfaces: outlet
  Average Over: 1(瞬时)或 10(移动平均)

2. Mass-Weighted Average(质量加权平均)

定义:Σ(φᵢ × ṁᵢ) / Σ(ṁᵢ)
用途:出口平均温度、组分浓度
设置示例:
  Name: outlet-temp-mass-avg
  Field Variable: Temperature
  Surfaces: outlet
  适用:需要考虑流量分布的平均

3. Facet Average(面元平均)

定义:Σ(φᵢ) / n
用途:几何平均(不考虑面积)
注意:网格依赖性强

4. Flow Rate(流率)

类型选择:
  - Mass Flow Rate:kg/s
  - Volume Flow Rate:m³/s
设置:
  Name: inlet-mass-flow
  Surfaces: inlet
  用途:质量守恒检查

5. Integral(积分)

定义:Σ(φᵢ × Aᵢ)
用途:总传热率、总力
示例:
  Name: wall-heat-transfer
  Field Variable: Total Surface Heat Flux
  Surfaces: heated-wall

6. Standard Deviation(标准偏差)

用途:均匀性评估
设置:
  Name: outlet-velocity-std
  Field Variable: Velocity Magnitude
  Surfaces: outlet
  判断:值越小越均匀
Volume Report(体积报告)类型

1. Volume Average(体积平均)

设置示例:
  Name: fluid-temp-avg
  Field Variable: Temperature
  Cell Zones: fluid
  用途:整体温度监控

2. Volume Integral(体积积分)

示例:热源总功率
  Name: total-heat-generation
  Field Variable: Energy Source
  Cell Zones: heat-source

3. Min/Max(最小/最大值)

监控极值:
  Name: max-temperature
  Report Type: Max
  Field Variable: Temperature
  Cell Zones: all-zones
  用途:安全限制检查
Force & Moment Reports(力和力矩报告)

Drag Coefficient(阻力系数)设置:

Name: drag-coefficient
Force Vector: [1, 0, 0](X方向)
Wall Zones: [选择所有壁面]
Reference Values设置:
  Area: 参考面积
  Density: 参考密度
  Velocity: 参考速度

Lift Coefficient(升力系数)设置:

Name: lift-coefficient  
Force Vector: [0, 1, 0](Y方向)
Wall Zones: [翼型表面]

Moment Coefficient(力矩系数)设置:

Name: moment-coefficient
Moment Center: [0.25, 0, 0](1/4弦长)
Moment Axis: [0, 0, 1](Z轴)
Wall Zones: [选择表面]
Report Files & Plots(报告文件和图表)
创建Report File(报告文件)

设置步骤:

  1. Solution → Reports → Files → New

  2. 选择Report Definitions

  3. 设置输出参数:

    Name: convergence-history
    Report Definitions: [选择多个]
    Frequency: 
      - Get Data Every: 1 Iteration/Time-step
      - Report File Frequency: 10(每10次写入)
    File Name: monitor.out
    

文件格式:

("变量名1" "变量名2" "变量名3")
迭代/时间  值1  值2  值3
1          0.5  300  0.001
2          0.48 301  0.0008
创建Report Plot(报告图表)

设置流程:

  1. Solution → Reports → Plots → New

  2. 配置选项:

    Name: residual-plot
    Report Definitions: [选择监控量]
    X-Axis:
      - Iteration(稳态)
      - Flow-Time(瞬态)
      - Time-Step(瞬态步数)
    Y-Axis:
      - Log Scale(对数坐标):残差用
      - Auto Scale(自动缩放):推荐
    Plot Window: 0-10(窗口ID)
    Frequency: 1(更新频率)
    

多图显示配置:

窗口0:残差图
窗口1:升阻力系数
窗口2:出口参数
窗口3:能量平衡
Convergence Conditions(收敛条件)
基于报告定义的收敛判断

设置方法:

  1. 创建关键参数的Report Definition

  2. Solution → Reports → Convergence Conditions

  3. 添加条件:

    Report Definition: outlet-mass-flow
    Condition Type: 
      - stop-if-all:所有条件满足时停止
      - stop-if-any:任一条件满足时停止
    Convergence Criteria:
      - Absolute: 值变化 < 1e-6
      - Relative: 相对变化 < 0.001
      - Iterate: 固定迭代次数
    Previous Iterations: 50(比较最近50次)
    

多条件组合示例:

条件1:残差全部 < 1e-4
条件2:出口流量变化 < 0.1%
条件3:阻力系数变化 < 0.01%
逻辑:所有条件同时满足
高级监控技术
Execute Commands(执行命令)

自动保存设置:

Command-1: 自动保存
  When: Every 100 iterations
  Command: /file/write-case-data backup-%i.cas.h5
  
Command-2: 导出图片
  When: Every 500 iterations
  Command: /display/save-picture contour-%i.png
Animation(动画)设置

Solution Animation设置:

  1. Solution → Calculation Activities → Solution Animations

  2. 创建动画序列:

    Animation-1:
      Type: Contour
      Variable: Velocity Magnitude
      Surfaces: symmetry-plane
      Storage:
        - Memory(内存):快速回放
        - HSF File(文件):节省内存
      Frequency: Every 10 time-steps
    
监控数据分析与判断
残差曲线分析

正常收敛模式:

特征:单调下降
斜率:-45°到-60°(对数坐标)
判断:良好收敛

振荡收敛:

特征:围绕趋势线振荡
原因:
  - 网格质量
  - 非稳态物理
  - 亚松弛因子过高
处理:
  - 降低亚松弛因子
  - 检查是否需要瞬态

平台期:

特征:残差不再下降
可能原因:
  - 达到机器精度(单精度~1e-6)
  - 达到离散误差极限
  - 物理振荡
判断方法:
  - 检查监控点
  - 检查质量守恒

发散:

特征:残差上升
紧急处理:
  1. 立即停止(Ctrl+C)
  2. 降低亚松弛因子50%
  3. 减小时间步/CFL数
  4. 检查边界条件
监控量判断标准

质量守恒检查:

Report → Fluxes
Net Mass Flow Rate: < 1% of inlet
理想值: < 0.1% of inlet

能量平衡(传热问题):

Total Heat Input - Total Heat Output
误差容限: < 1% of total

力系数稳定性:

Cd变化: < 0.001 (最近100次迭代)
Cl振幅: < 0.01 (除非物理振荡)
故障诊断流程图
残差不收敛?
├── 检查类型
│   ├── 持续振荡
│   │   ├── 降低亚松弛因子
│   │   ├── 检查网格质量
│   │   └── 考虑瞬态计算
│   ├── 缓慢下降
│   │   ├── 提高亚松弛因子
│   │   ├── 改用Coupled求解器
│   │   └── 检查初始化
│   └── 停滞不降
│       ├── 检查是否达到精度极限
│       ├── 监控其他物理量
│       └── 网格独立性研究
│
监控量异常?
├── 质量不守恒
│   ├── 检查边界条件设置
│   ├── 提高连续性收敛标准
│   └── 检查是否有泄漏
├── 能量不平衡
│   ├── 检查所有传热边界
│   ├── 降低能量方程亚松弛
│   └── 使用双精度求解器
└── 力系数振荡
    ├── 检查y+值
    ├── 时间步是否太大
    └── 是否存在涡脱落
关键成功要素
  1. 分阶段监控策略

    • 初始阶段:关注残差快速下降
    • 中期阶段:监控关键物理量稳定性
    • 收敛阶段:验证守恒性和精度
  2. 多重判断准则

    • 不仅依赖残差
    • 必须监控物理量
    • 检查质量能量守恒
  3. 自适应调整

    • 根据收敛行为调整参数
    • 问题特定的收敛标准
    • 经验积累很重要
  4. 文档记录

    • 记录成功的设置
    • 建立问题类型模板
    • 形成标准化流程

通过系统掌握这些监控和控制技术,可以显著提高CFD仿真的效率和可靠性,确保获得准确的计算结果。

4.8 结果后处理与分析

一、定量分析技术

定量分析技术总览

技术类型计算内容输出类型典型应用精度考虑注意事项验证方法
表面积分面积、流率、平均值数值出口温度、质量流量、压降网格依赖选择正确表面网格独立性
体积积分体积、总和、平均数值热源功率、平均温度、总质量单元类型影响2D需×2π守恒检查
通量报告质量/能量平衡平衡表守恒性检查、泄漏检测高精度包含所有边界Net≈0
力报告压力/粘性力力/系数升阻力、扭矩、压力中心壁面y+影响参考值设置实验对比
XY图沿线/面分布曲线图压力分布、速度剖面、温度曲线插值影响采样密度理论对比
直方图统计分布柱状图粒径分布、温度概率、速度统计分组数影响Bins=20-50正态性检验
监控点点值历史时间序列收敛判断、振荡检测、瞬态响应位置敏感避免壁面多点对比
自定义函数用户定义派生量效率、无量纲数、特殊参数公式正确性单位一致解析解验证
2.1 表面积分(Surface Integrals)
2.1.1 积分类型详解

访问路径:
Results → Reports → Surface Integrals

可用积分类型:

积分类型计算公式用途单位示例
Area∑Ai表面积
Area-Weighted Average∑(φi×Ai)/∑Ai平均温度、压力K, Pa
Mass-Weighted Average∑(φi×ṁi)/∑ṁi出口平均-
Mass Flow Rate∑ṁi质量流量kg/s
Volume Flow Rate∑(V·n×A)体积流量m³/s
Integral∑(φi×Ai)总传热率W
Standard Deviationσ = √(∑(φi-φ̄)²/n)均匀性-
Flow Rate∑(φi×V·n×A)物质流率varies
Sum∑φi简单求和-
Facet Average∑φi/n几何平均-
Facet Min/Maxmin/max(φi)极值-
Vertex Average节点平均光滑场-
Uniformity Index1-σ/φ̄均匀度0-1
2.1.2 典型应用案例

传热分析:

Wall Heat Transfer Rate:
- Report Type: Integral
- Field Variable: Total Surface Heat Flux
- Surfaces: heated-walls
- 结果单位: W

压降计算:

Pressure Drop:
- Inlet: Area-Weighted Average of Total Pressure
- Outlet: Area-Weighted Average of Total Pressure
- ΔP = P_inlet - P_outlet

混合效率:

Mixing Uniformity:
- Report Type: Uniformity Index (Area)
- Field Variable: Species Mass Fraction
- Surface: outlet
- 值域: 0(完全不均)到 1(完全均匀)
2.2 体积积分(Volume Integrals)
2.2.1 积分选项

访问路径:
Results → Reports → Volume Integrals

积分类型:

Volume:总体积
Sum:求和(2D轴对称×2π)
Volume Integral:体积积分
Volume Average:体积平均
Mass Integral:质量积分
Mass Average:质量平均
Min/Max:极值
2.2.2 应用示例

热源功率:

Total Heat Generation:
- Report Type: Volume Integral
- Field Variable: Energy Source
- Cell Zones: heat-source-zone
- 结果: 总发热功率 [W]

平均温度:

Bulk Temperature:
- Report Type: Mass-Weighted Average
- Field Variable: Temperature
- Cell Zones: fluid
2.3 通量报告(Flux Reports)
2.3.1 质量流量平衡

访问路径:
Results → Reports → Fluxes → Mass Flow Rate

输出信息:

各边界流量:
├── inlet: +X kg/s(流入为正)
├── outlet: -Y kg/s(流出为负)
├── User Sources: Z kg/s
└── Net: 应接近0(<0.1%)
2.3.2 能量平衡

Total Heat Transfer Rate包含:

能量组成:
├── Convection(对流)
├── Conduction(导热)
├── Radiation(辐射)
├── Species(组分焓)
├── Viscous Work(粘性功)
└── Pressure Work(压力功)

典型检查:

能量守恒检查:
Input - Output = Generation
误差容限: < 1% of total
2.4 力和力矩报告(Forces & Moments)
2.4.1 力系数计算

设置参数:

Force Report:
├── Wall Zones:选择所有壁面
├── Force Vector:[1,0,0](阻力方向)
├── Reference Values:
│   ├── Area:参考面积
│   ├── Density:参考密度
│   ├── Velocity:参考速度
│   └── Length:参考长度
└── 输出:
    ├── Pressure Force
    ├── Viscous Force
    └── Total Force

无量纲系数:

Cd = Fd/(0.5×ρ×V²×A)
Cl = Fl/(0.5×ρ×V²×A)
Cm = M/(0.5×ρ×V²×A×L)
2.4.2 力矩中心设置
Moment Center:
├── 翼型:1/4弦长
├── 汽车:前轴位置
├── 建筑:基础中心
└── 旋转机械:旋转轴
2.5 XY图表(XY Plots)
2.5.1 创建XY Plot

访问路径:
Results → Plots → XY Plot → New

设置选项:

Plot Direction:
├── X方向:[1,0,0]
├── Y方向:[0,1,0]
├── 自定义:[a,b,c]

Y Axis Function:
├── 选择绘图变量
└── 可多选(多曲线)

X Axis Function:
├── Direction:沿方向距离
├── Arc Length:弧长
└── Index:点索引
2.5.2 典型应用

边界层速度分布:

设置:
- 创建垂直壁面的线
- Y轴:Velocity Magnitude
- X轴:Wall Distance
- 用途:验证边界层理论

温度分布曲线:

管道中心线温度:
- Surface:centerline
- Y轴:Temperature
- X轴:Axial Coordinate
2.6 直方图(Histograms)
2.6.1 创建方法

Histogram设置:

参数:
├── Field Variable:统计变量
├── Surfaces/Zones:计算域
├── Number of Bins:分组数(10-100)
└── Range:自动或手动
2.6.2 应用场景

粒径分布:

  • DPM粒子直径统计
  • 喷雾特性分析

温度分布统计:

  • 燃烧室温度分布
  • 热点概率分析
2.7 自定义场函数(Custom Field Functions)
2.7.1 创建自定义函数

访问路径:
User Defined → Field Functions → Custom

函数定义示例:

压力系数:
Cp = (pressure - p_inf)/(0.5*density*velocity-mag^2)

总压损失系数:
K = (pt_inlet - pt_outlet)/(0.5*density*velocity^2)

传热系数:
h = wall-heat-flux/(T_wall - T_ref)
2.7.2 使用Report Definitions

创建复杂计算:

效率计算:
1. 创建入口总压报告:pt_in
2. 创建出口总压报告:pt_out
3. 定义效率:η = pt_out/pt_in
4. 监控效率变化
二、定性分析技术

定性分析技术总览

技术类型主要功能适用变量典型应用场景优势局限性推荐设置
云图(Contours)标量场分布可视化压力、温度、速度大小、浓度温度场、压力分布、浓度场直观、连续分布仅显示标量Filled+Node Values
矢量图(Vectors)矢量场方向和大小速度、热通量、应力流场结构、回流区、边界层显示方向和大小密集时混乱Scale=10, Skip=2
流线(Streamlines)瞬时流动路径速度场稳态流场结构、涡流识别流动模式清晰仅适用稳态从关键位置释放
路径线(Pathlines)质点运动轨迹速度场(瞬态)混合过程、停留时间、瞬态流动真实粒子路径计算量大Steps=1000-5000
等值面(Iso-surface)三维等值面任意标量涡核(Q准则)、温度分层、浓度边界3D结构清晰需选择合适值透明度50-70%
剖面(Profiles)变量沿法向分布任意变量边界层剖面、射流剖面定量分布需定义方向配合云图使用
场景(Scenes)多图形叠加组合显示综合分析、报告展示信息丰富可能混乱2-3层叠加为宜
动画(Animation)时间/空间变化任意显示瞬态过程、参数研究动态展示文件较大20-30 fps
2.1 云图显示(Contours)

云图是显示标量场分布的最常用可视化方法,通过颜色梯度表示变量在空间中的变化。

操作细节:

  • 路径:Results → Graphics → Contours
  • 选择显示变量(压力、速度、温度等)
  • 设置色标范围(自动或手动)
  • 选择显示表面或截面
  • 调整Levels数量(通常20-40)
  • 开启/关闭Node Values获得更平滑显示

应用场景:

  • 识别高低压区域
  • 分析温度分布
  • 查找分离区和回流区
  • 评估混合均匀性

温度场分析:

  • 变量:Static Temperature
  • 色彩:bgr(冷到热)
  • 等级:50-100(光滑)
  • 用途:传热分析、热点识别

压力分布:

  • 变量:Static Pressure
  • 色彩:rgb
  • 等级:20-30
  • 用途:压力梯度、分离区识别

速度分布:

  • 变量:Velocity Magnitude
  • 色彩:spectrum
  • 等级:30-50
  • 用途:流场结构、滞止区

基本设置

Contours of:选择显示变量
├── Pressure(压力)
│ ├── Static Pressure(静压)
│ ├── Total Pressure(总压)
│ └── Dynamic Pressure(动压)
├── Velocity(速度)
│ ├── Velocity Magnitude(速度大小)
│ ├── X/Y/Z Velocity(分量速度)
│ └── Tangential Velocity(切向速度)
├── Temperature(温度)
│ ├── Static Temperature(静温)
│ └── Total Temperature(总温)
└── Turbulence(湍流)
├── Turbulent Kinetic Energy
└── Turbulent Viscosity

显示选项详解

选项功能描述推荐设置应用场景
Filled填充色彩开启光滑渐变显示
Node Values节点插值开启平滑显示
Global Range全局范围视情况多面对比
Auto Range自动范围默认开启单面分析
Clip to Range裁剪范围关闭特定值范围
Draw Mesh显示网格关闭检查网格时开启
Contour Lines等值线视需求精确值定位

高级设置:

Levels(等级数):5-100
├── 默认:20
├── 光滑显示:50-100
└── 离散显示:5-10

Color Map(色彩映射):
├── bgr(蓝-绿-红):温度场
├── rgb(红-绿-蓝):速度场
├── gray(灰度):打印友好
└── custom(自定义):特殊需求
2.2 矢量图(Vectors)

矢量图用箭头表示流场的方向和大小,是分析流动模式的重要工具。

操作细节:

  • 路径:Results → Graphics → Vectors
  • 设置矢量密度(Skip值)
  • 调整缩放因子(Scale)
  • 选择颜色映射变量
  • 设置矢量样式(箭头、锥形等)
  • 可叠加在云图上显示

应用场景:

  • 分析流动方向和模式
  • 识别涡旋和分离点
  • 评估二次流
  • 研究边界层发展

矢量显示控制:

参数功能推荐值说明
Scale箭头缩放1-20控制箭头长度
Skip跳过间隔1-10减少矢量密度
Fixed Length固定长度关闭所有箭头等长
Normalize归一化关闭仅显示方向
Style箭头样式arrow可选cone、harpoon
Color by着色变量Velocity Magnitude可选其他变量

矢量过滤选项:

Display限制:
├── In Plane:平面内分量
├── X Component:仅X分量
├── Y Component:仅Y分量
└── Z Component:仅Z分量(3D)

范围控制:
├── Min:最小显示值
├── Max:最大显示值
└── Skip:每n个显示一个

特殊应用技巧

  • 边界层分析:
设置组合:
- Scale: 5-10(较小)
- Skip: 0(全显示)
- 局部放大关键区域
- 配合y+云图
  • 涡流识别:
设置组合:
- In Plane显示
- 适当Skip值(2-5)
- 配合涡量云图
- 关键截面分析
2.3 流线图(Pathlines/Streamlines)

概念区分

类型定义适用特点
Streamline瞬时切线稳态空间线
Pathline质点轨迹瞬态时间线
Streakline染料线实验对比连续释放
Time Line时间线瞬态同时释放

操作细节:

  • 路径:Results → Graphics → Pathlines
  • 选择释放表面(inlet、特定点或线)
  • 设置粒子数量(通常100-500)
  • 选择着色变量(速度、压力、停留时间)
  • 调整积分步长和最大步数
  • 可制作动画显示瞬态过程

应用场景:

  • 追踪流体质点轨迹
  • 分析混合过程
  • 识别死区和短路流
  • 优化流道设计

Pathlines详细设置

  • 释放源设置:
Release from Surfaces:
├── inlet(入口释放)
├── interior points(内部点)
├── line surface(线释放)
└── point surface(点释放)


释放参数:
├── Number of Paths:10-1000
├── Step Length Factor:0.01-1.0
├── Max Steps:500-10000
└── Path Skip:1-10

显示控制:

Style选项:
├── line(线型):默认
├── ribbon(带状):显示旋转
├── cylinder(圆柱):3D效果
└── sphere(球形):离散显示

着色选项:
├── Velocity Magnitude
├── Pressure
├── Temperature
├── Residence Time
└── Path ID

高级应用

停留时间分析:

  • 变量:Residence Time
  • 用途:混合效率、死区识别
  • 设置:大量路径线(>500)

粒子追踪:

  • DPM耦合显示
  • 质量/热量传递可视化
  • 多相流分析
2.4 等值面(Iso-surfaces)

操作细节:

  • 路径:Results → Surface → Iso-Surface
  • 选择等值变量
  • 设置等值大小
  • 可在等值面上显示其他变量云图
  • 用于三维流场可视化

应用场景:

  • 显示涡核(Q准则、λ2准则)
  • 追踪火焰面
  • 分析激波位置
  • 显示浓度分布
2.5 动画制作(Animations)

操作细节:

  • 路径:Results → Animations → Solution Animation
  • 选择动画类型(时间序列或参数扫描)
  • 设置帧数和时间间隔
  • 选择显示对象(云图、矢量等)
  • 导出格式(AVI、MPEG、图片序列)

应用场景:

  • 展示瞬态过程
  • 分析周期性流动
  • 演示参数影响
  • 结果展示汇报

后处理工具快速选择指南

分析目的首选工具次选工具辅助工具输出格式典型设置
流场结构流线+矢量云图(速度)涡量等值面PNG/矢量图多截面组合
温度分布温度云图等温线XY曲线高分辨率图50-100 levels
压力分析压力云图Cp曲线压力积分数据表格全局范围
混合评估浓度云图粒子轨迹统计分布动画时间序列
性能评估数值报告XY曲线监控历史Excel表自动输出
分离检测壁面剪切力流线回流区标记矢量+云图近壁放大
涡识别Q/λ2准则涡量云图流线3D渲染等值面透明
瞬态演化动画相位平均FFT分析MP4/GIF适当帧率
三、结果验证技术

不同物理问题的后处理策略

问题类型关键定性分析关键定量分析重点关注参数常见错误检查报告必备内容
内部流动速度云图、流线、压力分布压降、流量分配、摩擦系数Re数、压降系数、速度分布均匀性入口发展段、出口回流压降-流量曲线、速度剖面
外部绕流压力系数、流线、涡量升阻力、St数、压力积分Cd、Cl、分离点、尾流区远场边界、阻塞比力系数vs攻角、压力系数分布
传热问题温度云图、热通量矢量Nu数、传热率、热阻温差、h分布、热流密度能量不平衡、Bi数能量平衡表、Nu关联式
燃烧反应温度、组分、反应率NOx生成、燃烧效率火焰温度、当量比、停留时间组分守恒、反应速率温度分布、污染物排放
多相流相分布、界面、速度差含气率、滑移比、压降体积分数、界面面积、气泡尺寸质量守恒、界面捕捉流型图、压降关联式
旋转机械相对速度、压力、流线扭矩、功率、效率压比、效率、流量系数周期性、MRF设置性能曲线、速度三角形
瞬态问题时间动画、相位平均频谱分析、RMS值频率、振幅、相位差时间步、采样率FFT谱、时间历程
传质问题浓度云图、扩散通量Sh数、传质率、混合效率浓度梯度、Pe数、混合时间组分守恒、Schmidt数混合曲线、传质系数
3.1 残差监测(Residuals)

操作细节:

  • 路径:Results → Plots → Residuals
  • 检查所有方程残差趋势
  • 目标:连续性<1e-3,其他<1e-6
  • 观察振荡情况判断收敛性
  • 可保存残差历史用于后处理

验证内容:

  • 数值收敛性
  • 求解稳定性
  • 迭代充分性
3.2 监测点设置(Monitor Points)

操作细节:

  • 路径:Solve → Monitors → Surface Monitors
  • 创建点、线或面监测器
  • 选择监测变量
  • 设置输出频率和文件
  • 实时绘图显示

验证内容:

  • 关键位置参数稳定性
  • 周期性流动特征
  • 瞬态响应
3.3 网格独立性验证

操作细节:

  • 准备3套不同密度网格(粗、中、细)
  • 使用相同边界条件求解
  • 对比关键参数(压降、传热量等)
  • 计算Richardson外推误差
  • 绘制网格收敛性曲线

验证内容:

  • 数值误差评估
  • 网格质量影响
  • 计算精度保证
3.4 Y+值检查

操作细节:

  • 路径:Results → Plots → XY Plot
  • 选择Y+ (Wall Yplus)
  • 选择所有壁面
  • 导出数据进行统计分析
  • 标准:Y+<1(低雷诺数)或30<Y+<300(壁面函数)

验证内容:

  • 近壁面网格质量
  • 湍流模型适用性
  • 边界层解析度
3.5 实验数据对比

操作细节:

  • 导入实验数据(File → Read → Profile)
  • 创建XY Plot对比曲线
  • 计算相对误差和均方根误差
  • 使用自定义场函数(Custom Field Functions)
  • 生成验证报告

验证内容:

  • 模型准确性
  • 物理真实性
  • 参数可靠性
四、数据导出与第三方软件接口
导出类型文件格式适用软件用途导出命令/路径注意事项
网格数据.msh/.casTecplot/ParaView高级可视化File→Export→Mesh包含结果数据
场数据.dat/.csvExcel/Matlab数据分析File→Export→Solution DataASCII格式
图片.png/.jpg/.epsPPT/论文报告展示Display→Save Picture高分辨率
动画.mpeg/.avi媒体播放器动态展示Solve→Animation→Playback压缩设置
XY数据.xy/.csvOrigin/Excel曲线拟合Plot→File→Write包含标题行
报告.txt/.html文本编辑器结果记录Report→Summary自动格式化
EnSight.case/.encasEnSight专业后处理File→Export→EnSight时间序列
CGNS.cgns多种CFD软件数据交换File→Export→CGNS标准格式
五、后处理技术应用场景总结表
技术类别具体技术主要用途适用场景输出形式
定量分析表面积分计算通量、流量进出口平衡、传热分析数值报告
体积积分平均值计算混合效率、能量分布数值报告
通量报告守恒性检查质量/能量平衡验证表格数据
力和力矩受力分析气动设计、结构载荷数值+分量
定性分析云图场变量分布温度场、压力场分析彩色图像
矢量图流动方向速度场、涡旋识别箭头图
流线图流动轨迹流型分析、混合评估线条图
等值面3D特征提取涡结构、激波定位3D曲面
动画动态过程瞬态演化、参数研究视频文件
结果验证残差监测收敛性判断所有稳态计算对数曲线
监测点局部追踪关键位置参数变化时间历程
网格独立性数值误差精度验证对比曲线
Y+检查网格质量湍流模拟验证分布曲线
实验对比模型验证准确性评估对比图表

后处理常见问题与解决方案

问题现象可能原因诊断方法解决方案预防措施
云图不连续Cell Values模式检查Node Values选项开启Node Values默认使用节点值
矢量太密Skip值太小查看Vector设置增加Skip值(5-10)根据网格调整
等值线锯齿网格太粗检查网格密度网格加密/插值适当的网格密度
数值异常大单位不一致检查Reference Values统一单位系统预设单位模板
流线不光滑步长太大查看Step设置减小Step Factor自适应步长
图片模糊分辨率低检查图片设置提高分辨率预设高质量模板
颜色不当色谱选择查看Colormap选择合适色谱不同变量不同色谱
数据缺失未完全收敛检查残差继续计算设置自动保存
五、后处理最佳实践建议
  1. 系统性流程:先进行结果验证(收敛性、网格独立性),再进行定量分析获取关键数据,最后通过定性分析深入理解流动机理。
  2. 数据管理:建立规范的命名体系,保存所有后处理设置文件(.jou文件),便于重复使用和批处理。
  3. 可视化原则
    • 云图使用合适的色标范围,避免过度饱和
    • 矢量图密度适中,确保清晰可辨
    • 动画帧率选择15-30fps,文件大小与质量平衡
  4. 验证层次
    • 数值验证(残差、守恒性)
    • 物理验证(无量纲数、经验关联式)
    • 实验验证(测试数据对比)
  5. 报告规范:按照"验证-定量-定性"的顺序组织结果,确保科学性和可追溯性。

通过掌握这些后处理技术的操作细节和应用场景,可以全面、准确地分析FLUENT计算结果,为工程设计和科学研究提供可靠依据。

五、完整自检表

5.1 SpaceClaim 几何建模自检表

检查项检查内容检查标准不通过后果
单位设置检查单位系统是否正确单位与实际模型一致,如毫米 (mm) 或米 (m)模型尺寸错误,导致仿真结果完全偏离实际
几何完整性检查模型是否有间隙、重叠或不连续部分无间隙、重叠或不连续,所有面和体完整闭合网格划分失败或仿真结果不准确
模型尺寸检查关键尺寸是否符合设计要求关键尺寸误差在允许范围内(通常 < 1%)仿真结果不符合实际物理问题
开口设置检查进出口开口是否正确设置开口位置和尺寸正确,命名规范(如 inlet、outlet)边界条件设置错误,影响计算结果
布尔操作检查布尔操作(并、交、差)是否正确布尔操作结果符合设计意图,无多余或缺失部分模型结构错误,无法正确模拟实际流动
对称设置检查对称面是否正确设置对称面位置和方向正确,模型对称部分完全一致计算量增加,结果可能不对称
模型简化检查模型简化是否合理简化后的模型保留关键特征,去除非必要细节模型过于复杂导致计算效率低,或简化过度导致结果不准确

5.2 Fluent Meshing 网格划分自检表

检查项检查内容检查标准不通过后果
全局网格尺寸检查全局网格尺寸设置是否合理全局网格尺寸与模型特征尺寸匹配,通常为特征尺寸的 1/10 到 1/20网格过粗导致精度不足,或过密导致计算资源浪费
局部网格加密检查关键区域是否进行了适当加密边界层、狭窄通道、流动分离区域等关键区域网格加密适当无法捕捉流动细节,导致结果不准确
网格类型选择检查网格类型选择是否合适复杂几何使用四面体或混合网格,规则区域使用六面体网格网格质量差,计算不收敛或结果精度低
边界层设置检查边界层网格设置是否合理边界层层数足够(通常 3-5 层),y + 值在合理范围内(湍流模型要求)无法正确模拟壁面附近的流动,影响整体结果
网格质量检查网格质量指标正交性 > 0.2,偏斜度 < 0.9,纵横比 < 5:1(边界层区域 < 10:1)网格质量差导致计算不收敛或结果不准确
网格连通性检查网格是否连通且无孤立区域所有单元正确连接,无孤立节点或单元计算无法进行或结果错误
区域定义检查流体和固体区域定义是否正确流体区域和固体区域正确划分,命名规范物理模型设置错误,影响计算结果
网格统计检查网格统计信息节点数、单元数合理,无负体积或负面积网格存在错误,无法进行计算

5.3 Fluent 求解设置自检表

检查项检查内容检查标准不通过后果
求解器选择检查求解器类型是否正确不可压缩或低马赫数流动选择压力基求解器,高马赫数流动选择密度基求解器计算效率低或结果不准确
时间类型设置检查时间类型(稳态 / 瞬态)是否正确根据流场是否随时间变化正确选择时间类型结果不符合实际物理现象
物理模型选择检查湍流模型、多相流模型等是否正确湍流模型根据雷诺数和流动特性选择,多相流模型根据相态选择无法正确模拟流动特性,结果不准确
材料属性设置检查材料属性是否正确密度、粘度、导热系数等参数设置正确物理特性模拟错误,影响结果
边界条件设置检查边界条件是否正确设置入口、出口、壁面等边界条件符合实际物理情况边界条件错误导致结果偏离实际
初始条件设置检查初始条件是否合理初始条件接近真实解,避免过大的初始扰动收敛困难或结果振荡
求解方法设置检查求解方法参数是否合理压力 - 速度耦合算法、离散格式等选择适当收敛速度慢或结果精度低
松弛因子设置检查松弛因子是否合理主要方程松弛因子在 0.5-1.0 之间,根据收敛情况调整收敛困难或计算发散
残差设置检查残差收敛标准是否合理主要方程残差低于 1e-3,能量方程低于 1e-6结果可能未真正收敛,影响精度
监控设置检查关键物理量监控是否设置设置进出口流量、压力降等关键物理量的监控无法判断计算是否真正收敛
结果验证检查结果是否合理质量守恒、能量守恒误差 < 1%,物理量分布符合预期结果不可靠,需要重新计算

5.4 完整仿真流程检查清单

阶段关键检查点检查标准不通过后果
几何建模单位一致性SpaceClaim 和 Fluent Meshing 单位一致模型尺寸错误
几何完整性无间隙、重叠或不连续部分网格划分失败
开口设置进出口开口位置和尺寸正确边界条件设置错误
网格划分网格质量正交性 > 0.2,偏斜度 < 0.9计算不收敛或结果不准确
边界层设置边界层层数足够,y + 值合理壁面流动模拟错误
网格统计无负体积或负面积计算无法进行
求解设置求解器选择根据流动特性正确选择求解器计算效率低或结果错误
物理模型湍流模型、多相流模型等选择正确流动特性模拟错误
边界条件边界条件符合实际物理情况结果偏离实际
结果分析质量守恒进出口质量流量误差 < 1%结果不可靠
能量守恒能量收支误差 < 1%结果不可靠
物理量合理性压力、速度分布符合物理预期结果可能错误
数据验证与理论或实验数据对比误差在合理范围内结果可信度低

六、常见问题与解决方案

6.1 几何建模常见问题

问题 1:导入几何模型失败

可能原因:

  • 文件路径包含中文字符或特殊字符

  • 文件格式不兼容

  • 文件已损坏

解决方案:

  • 使用英文路径和文件名

  • 尝试其他几何格式(如.stp、.iges)

  • 检查文件完整性,必要时重新生成模型

问题 2:几何模型存在间隙

可能原因:

  • 从 CAD 软件导出时出现错误

  • 模型本身设计存在间隙

  • 简化过程中产生了间隙

解决方案:

  • 使用 SpaceClaim 的 “间隙检测” 功能查找并修复间隙

  • 调整公差设置,自动修复小间隙

  • 手动修复大的间隙

问题 3:开口设置错误

可能原因:

  • 开口位置选择错误

  • 开口面未正确命名

  • 开口尺寸不符合要求

解决方案:

  • 仔细检查开口位置和尺寸

  • 使用 “重命名” 功能正确命名开口面

  • 必要时重新创建开口

6.2 网格划分常见问题

问题 1:网格划分失败

可能原因:

  • 几何模型存在错误

  • 网格尺寸设置不合理

  • 网格类型选择不当

解决方案:

  • 修复几何模型错误

  • 调整网格尺寸参数,增加局部加密

  • 尝试不同的网格类型或划分策略

问题 2:网格质量差

可能原因:

  • 网格尺寸变化过快

  • 曲率或邻近度设置不当

  • 边界层设置不合理

解决方案:

  • 调整 Growth Rate 参数,减小网格尺寸变化率

  • 优化曲率和邻近度设置

  • 调整边界层参数,如层数和增长率

问题 3:计算区域未正确识别

可能原因:

  • 几何模型未正确共享拓扑

  • 流体区域未正确封闭

  • 区域命名错误

解决方案:

  • 使用 “共享拓扑” 功能确保区域正确连接

  • 使用 “封闭流体区域” 功能自动封闭开口

  • 正确命名流体和固体区域

6.3 求解计算常见问题

问题 1:计算不收敛

可能原因:

  • 初始条件不合理

  • 松弛因子设置不当

  • 网格质量差

  • 物理模型选择错误

解决方案:

  • 调整初始条件,使其更接近真实解

  • 减小松弛因子,提高收敛稳定性

  • 改善网格质量

  • 选择更合适的物理模型

问题 2:结果不符合物理预期

可能原因:

  • 边界条件设置错误

  • 物理模型选择错误

  • 网格分辨率不足

  • 求解参数设置不当

解决方案:

  • 仔细检查边界条件设置

  • 确认物理模型选择正确

  • 增加网格密度,特别是关键区域

  • 调整求解参数,如松弛因子和离散格式

问题 3:质量或能量不守恒

可能原因:

  • 边界条件设置错误

  • 源项设置错误

  • 求解控制参数不当

  • 网格质量差

解决方案:

  • 检查进出口边界条件是否正确

  • 确认源项设置无误

  • 调整求解控制参数

  • 改善网格质量,特别是边界层网格

6.4 后处理与分析常见问题

问题 1:结果数据无法正确导出

可能原因:

  • 文件路径包含特殊字符

- 文件权限不足

  • 数据文件过大

解决方案:

  • 使用英文路径和文件名

  • 确保有写入权限

  • 分区域导出数据,减小文件大小

问题 2:云图显示异常

可能原因:

  • 物理量选择错误

  • 显示范围设置不当

  • 网格质量差

解决方案:

  • 确认选择了正确的物理量

  • 调整颜色映射和显示范围

  • 改善网格质量,特别是边界层网格

问题 3:无法正确提取关键数据

可能原因:

  • 监测点或监测面设置错误

  • 数据尚未收敛

  • 物理量命名错误

解决方案:

  • 仔细检查监测点或监测面的位置和类型

  • 确保计算已充分收敛

  • 确认物理量名称正确

七、最佳实践与提高效率建议

7.1 几何建模最佳实践

几何简化策略

对于复杂模型,采用以下简化策略:

  1. 分层次简化:先进行整体简化,再逐步细化关键区域

  2. 去除非关键特征:如小倒角、小凸台等对流动影响不大的特征

  3. 使用 “体积抽取” 功能:快速创建流体区域,简化内部结构

  4. 保留关键几何特征:如狭窄通道、尖锐边缘等对流动有重要影响的特征

命名规范

采用一致的命名规范:

  1. 区域命名:如 fluid、solid、inlet、outlet 等

  2. 开口命名:按功能命名,如 inlet_1、outlet_2 等

  3. 层命名:如 wall、symmetry、interface 等

规范的命名有助于后续的网格划分和边界条件设置,减少错误。

7.2 网格划分最佳实践

网格划分策略

采用以下网格划分策略:

  1. 由粗到精:先创建较粗的网格进行初步计算,根据结果再细化关键区域

  2. 关键区域优先:对边界层、狭窄通道、流动分离区域等关键区域优先加密

  3. 使用 BOI 和 FOI:利用 Body of Influence 和 Face of Influence 控制局部网格密度

  4. 网格质量检查:生成网格后严格检查质量指标,确保符合要求

边界层网格设置

边界层网格设置的最佳实践:

  1. 根据 y + 目标确定层数:对于 k-ω SST 模型,目标 y+≈1,通常需要 3-5 层

  2. 增长率设置:推荐 1.2-1.5,避免增长过快

  3. 总厚度计算:边界层总厚度通常为特征长度的 1-5%

  4. 验证 y + 值:计算完成后检查y + 分布,确保在合理范围内

7.3 求解设置最佳实践

收敛加速技巧

采用以下技巧加速收敛:

  1. 使用混合初始化:比标准初始化更接近真实解,加快收敛

  2. 逐步增加松弛因子:初始计算使用较低的松弛因子,稳定后逐渐增加

  3. 分阶段计算:先使用较粗的网格和简单模型进行初步计算,再细化网格和模型

  4. 监测关键物理量:设置关键物理量的监测,确保其趋于稳定

求解参数优化

优化求解参数以提高计算效率:

  1. 压力 - 速度耦合算法选择:对于复杂几何,使用 SIMPLEC 算法比 SIMPLE 收敛更快

  2. 空间离散格式:优先使用二阶格式,提高计算精度

  3. 库朗数****调整:在密度基求解器中,初始库朗数可设为 5-10,稳定后提高至 50-100

  4. 使用多重网格:启用多重网格加速收敛

7.4 结果验证与分析最佳实践

结果验证方法

采用多种方法验证结果:

  1. 质量和能量守恒检查:确保计算满足基本物理守恒定律

  2. 网格独立性研究:使用至少三种不同密度的网格进行计算,确认结果稳定

  3. 物理合理性****分析:检查结果是否符合基本物理原理和预期

  4. 与实验或理论数据对比:如有可用数据,进行对比验证

后处理技巧

提高后处理效率的技巧:

  1. 创建自定义场函数:使用 Fluent 的表达式功能创建自定义场函数

  2. 自动化脚本:使用 Scheme 或 Python 脚本自动化后处理过程

  3. 结果对比分析:使用 Fluent 的 “Compare” 功能对比不同工况的结果

  4. 数据可视化:使用专业绘图软件(如 CFD-Post、MATLAB 等)进行高级可视化

八、总结与展望

8.1 完整仿真流程回顾

本文详细介绍了使用 SpaceClaim、Fluent Meshing 和 Fluent 进行流体力学仿真的完整流程:

  1. 几何建模:使用 SpaceClaim 创建或处理几何模型,确保模型完整、正确

  2. 网格划分:使用 Fluent Meshing 生成高质量网格,重点关注关键区域的网格质量

  3. 求解设置:在 Fluent 中设置合适的求解器、物理模型和边界条件

  4. 结果分析:对计算结果进行全面分析和验证,确保结果准确可靠

通过遵循本文提供的步骤和自检表,可以系统性地完成流体力学仿真,并提高结果的准确性和可靠性。

8.2 高级技巧与未来发展

高级仿真技术

随着软件和硬件的发展,以下高级技术值得关注:

  1. 大涡模拟 (LES):比 RANS 模型更精确,适用于复杂湍流流动

  2. 分离涡模拟 (DES):结合 RANS 和 LES 的优点,适用于分离流动

  3. 多物理场耦合:流固耦合、热流耦合等复杂多物理场问题

  4. 机器学习辅助仿真:利用机器学习技术加速求解或提高预测精度

未来发展趋势

流体仿真技术的未来发展趋势:

  1. 更高的计算效率:通过并行计算、算法优化等提高计算效率

  2. 更精确的物理模型:发展更精确的湍流模型、多相流模型等

  3. 自动化工作流程:从几何建模到结果分析的全流程自动化

  4. 云仿真平台:基于云平台的仿真服务,提供更便捷的计算资源访问

8.3 持续学习与提高

学习资源推荐

以下资源有助于进一步学习:

  1. ANSYS 官方文档:包括用户手册、理论手册和验证手册

  2. 在线教程:ANSYS Learning Hub 提供丰富的在线教程和案例

  3. 学术论文:阅读流体力学和计算流体力学领域的最新研究论文

  4. 实践项目:参与实际项目,积累经验,提高应用能力

社区与交流

参与社区交流,分享经验:

  1. ANSYS 用户论坛:与其他用户交流经验和问题

  2. 专业学术会议:参加流体力学和计算流体力学相关的学术会议

  3. 行业研讨会:了解最新的行业应用和技术发展

  4. 在线社区:如 CFD Online、ResearchGate 等平台

通过不断学习和实践,您可以逐步提高流体力学仿真的能力和水平,解决更复杂的工程问题。

九、附录:示例仿真项目

9.1 简单管道流动仿真示例

问题描述

模拟不可压缩流体在二维管道中的流动,管道尺寸为 1m×0.1m,入口速度 1m/s,出口压力为大气压。

几何建模步骤

  1. 在 SpaceClaim 中创建矩形:
  • 宽度 1m,高度 0.1m

  • 坐标原点位于左下角

  1. 设置开口:
  • 左侧面为入口 (inlet)

  • 右侧面为出口 (outlet)

  • 上下面为壁面 (wall)

  1. 保存为管道几何模型 (pipe.agdb)。

网格划分步骤

  1. 在 Fluent Meshing 中导入管道几何模型

  2. 设置全局网格尺寸 0.05m

  3. 对壁面附近添加边界层:

  • 3 层,增长率 1.2

  • 目标 y+≈1

  1. 生成非结构化四边形网格

  2. 检查网格质量,确保正交性 > 0.3,偏斜度 < 0.6。

求解设置步骤

  1. 在 Fluent 中导入网格

  2. 选择压力基求解器,稳态计算

  3. 启用 k-ω SST 湍流模型

  4. 设置入口边界条件:

  • 速度入口,速度 1m/s

  • 湍流强度 5%,水力直径 0.1m

  1. 设置出口边界条件:
  • 压力出口,表压 0Pa
  1. 设置壁面边界条件:
  • 无滑移,绝热壁面
  1. 初始化并求解,迭代 500 步。

结果分析

  1. 检查质量守恒,确保进出口流量差 < 1%

  2. 查看速度云图,确认流动充分发展

  3. 查看压力分布,计算压力降

  4. 提取中心线速度分布,与理论解对比。

9.2 带障碍物的管道流动仿真示例

问题描述

模拟不可压缩流体流过管道中的障碍物,管道尺寸为 2m×0.5m,障碍物尺寸 0.2m×0.2m,位于管道中心,入口速度 1m/s,出口压力为大气压。

几何建模步骤

  1. 在 SpaceClaim 中创建管道和障碍物:
  • 管道:2m×0.5m

  • 障碍物:0.2m×0.2m,位于管道中心

  1. 使用布尔差集运算从管道中减去障碍物

  2. 设置开口:

  • 左侧面为入口 (inlet)

  • 右侧面为出口 (outlet)

  • 所有内表面为壁面 (wall)

  1. 保存为带障碍物的管道几何模型 (pipe_obstacle.agdb)。

网格划分步骤

  1. 在 Fluent Meshing 中导入几何模型

  2. 设置全局网格尺寸 0.1m

  3. 对障碍物周围和壁面添加局部加密:

  • 障碍物周围网格尺寸 0.02m

  • 壁面附近添加边界层,3 层,增长率 1.2

  1. 生成混合网格,障碍物周围使用三角形网格,其他区域使用四边形网格

  2. 检查网格质量,确保正交性 > 0.3,偏斜度 < 0.7。

求解设置步骤

  1. 在 Fluent 中导入网格

  2. 选择压力基求解器,稳态计算

  3. 启用 k-ω SST 湍流模型

  4. 设置入口边界条件:

  • 速度入口,速度 1m/s

- 湍流强度 5%,水力直径 0.5m

  1. 设置出口边界条件:
  • 压力出口,表压 0Pa
  1. 设置壁面边界条件:
  • 无滑移,绝热壁面
  1. 初始化并求解,迭代 1000 步。

结果分析

  1. 检查质量守恒

  2. 查看速度矢量图,观察障碍物后的流动分离和回流区域

  3. 查看压力云图,观察压力分布

  4. 提取障碍物表面压力分布,计算阻力系数

  5. 比较不同网格密度下的结果,验证网格独立性。

Logo

AtomGit 是由开放原子开源基金会联合 CSDN 等生态伙伴共同推出的新一代开源与人工智能协作平台。平台坚持“开放、中立、公益”的理念,把代码托管、模型共享、数据集托管、智能体开发体验和算力服务整合在一起,为开发者提供从开发、训练到部署的一站式体验。

更多推荐