四十四、​Fluent 收敛标准-质量和能量守恒

本文介绍如何通过质量守恒和能量守恒来判断CFD软件Fluent的计算是否收敛,包括不同条件下的Fluent Flux通量检查方法。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

1. Fluent质量和能量守恒


Fluent收敛性的判断有时候很让人头疼,我们在三十二、Fluent收敛判断标准及方法中概括的讲述过几种收敛标准。



通常我们判断是否收敛最直接的就是看残差曲线,但是很多时候,残差曲线(尤其连续性曲线)一直是平的,连下降的趋势都没有。这时候最正确的做法应该是监测比较重要的物理量,查看其是否稳定。


如果监测的几个物理量都很稳定,那可以忽略残差曲线而认为计算已经达到收敛。

但是我们计算过程中更多的是遇到不稳定的情况,或者我们根本就不确定计算是否收敛了,那应该怎么办呢??本文提供一点思路来帮助理解计算发散的原因。

任何流动或传热问题的计算都比较质量守恒、能量守恒和动量守恒。当计算发散时,或者我们不确定它是否收敛的时候,就可以查看Fluent是否满足三大守恒,如果不满足,那肯定不收敛。基于此我们可以找到发散的原因。



注:
1. Fluent满足三大守恒也不一定收敛,还要辅助监测物理量来判断收敛性

2. 强推文章四十三、Fluent增强收敛性-伪瞬态计算,能够大大加快收敛速度,尝试了一下,使用SIMPLEC算法10000迭代步可以收敛,使用伪瞬态算法只需要1000步,收敛速度加快了10倍



2. Fluent Flux通量

如何查看Fluent是否满足质量守恒、能量守恒呢??在Fluent后处理中Results-Reports-Fluxes可以查看


但是有几点需要注意:


2.1 Options包含三个选项:

Mass Flow Rate、Total Heat Transfer Rate和Radiation Heat Transfer Rate,其他的以后补充。通过这三个选项可判断计算是否满足质量守恒和能量守恒。


2.2 不涉及源项问题

mass flow rate:所有进出口质量流率之和应该为0,其中进口为正,出口为负。Net Results表示进出口质量流率之和。壁面wall显示为0,计算域内的质量流量显示不为0,但计算Net Results时,内部区域会被跳过,只计算进出口区域。

Net Results:只计算模型进出口的质量流量


Total Heat Transfer Rate:所有进出口及壁面的总传热率之和应该为0,即Net Results=0。计算域内的总传热率显示为0。

Net Results:只计算模型进、出口和壁面的的总传热率;如下图不涉及能源源项,Net Results应该等于0,而图中等于-4.46,说明收敛性不好,应该继续计算或者修改部分参数后计算。


内部边界:如果涉及到内部边界(边界类型 interior),上述三个Fluxes可能会不准确。需要先将内部边界类型更改为wall(用于定位面积向量),不用进行任何操作,然后再更改为interior,计算Fluxes


如果不涉及质量源项,那么模型进出口的mass flow rate应该相等,净值Net Results应该等于0;
如果不涉及能量源项,模型进出口及壁面处的Total Heat Transfer Rate的净值Net Results应该等于0


2.3 涉及源项问题

如果在Cell Zone Conditions---Fluid中设置了质量源项(kg/(m3*s)),则Net results显示的为添加质量源项计算域的积分值。即每个网格的质量源项与计算域网格体积相乘之和,单位恰好为kg/s。


同理如果在Cell Zone Conditions---Fluid中设置了能量源项(W/m3),则Net results显示的为添加能量源项计算域的积分值。



2.4 涉及多相流问题

使用多相流模型,上述三个Fluxes可以对单独相也可对mixture进行计算。如果没有源项,那么每一相的net results=0;


如果对某一相设置了质量或能量源项,那么对这一相的net results就是添加源项计算域的积分值。

如果涉及到相间的传热和传质,那么对于单独相,则mass flow rate和total heat transfer rate的Net Results不为0;

但是对于mixture相,mass flow rate和total heat transfer rate的Net Results=0。比如蒸发冷凝模型、PBM模型。
 


2.5 涉及到颗粒流动

如果使用离散相模型(DPM),则会单独报告粒子注入的质量流量和传热率,并且Net Results结果会将粒子注入的质量流量和传热率也包含进去。

比如上图,Net Results等于Results中各个质量流量之和,再加上DPM Mass Source。Net Results=-1.770709+0.876162 +0.876162+0.01839178=6.367223e-6

也就是说,涉及DPM模型,不涉及源项时,Net result应该也为0。


3. Flux通量总结


单相流且不涉及源项,质量和能量的Net results=0;

多相流不涉及源项,且没有相间传热传质则单独相的Net results=0;
若涉及相间传热传质,则单独相的Net results不等于0,但mixture相的质量和能量的Net results=0;

颗粒流不涉及源项时,Net result也为0。
 

### 查看解释 Fluent 中的残差曲线 在 Fluent 软件中,残差曲线用于评估数值模拟过程中方程式的收敛情况。当观察到残差曲线时,主要关注的是其下降趋势以及最终是否趋于平稳。 #### 残差曲线的查看方式 为了查看 Fluent 的残差曲线,可以在图形界面中的菜单栏选择 `Plot` -> `Residuals...` 或者通过命令面板输入 `/plot/residuals` 命令来激活残差图表窗口[^4]。此图表会动态更新每次迭代后的残差变化状况。 #### 判断收敛标准 判断 Fluent 计算过程是否已经达到了满意的收敛状态可以从以下几个角度入手: -3}\) 至 \(10^{-6}\),具体取决于工程应用的需求)。不过需要注意的是,不同版本 fluent 对于默认显示的比例尺有所不同,因此建议确认实际数值而非仅依赖视觉印象。 - **稳定性波动范围**:即使某些时刻残差不再继续减少而是保持在一个很小范围内上下浮动,只要这种波动幅度非常微弱且整体处于低位,则可初步认定为接近收敛状态[^3]。 - **物理量的一致性检验**:除了单纯依靠数学上的定义外,更重要的是考察所得结果能否合理反映实际情况——即与其他已知条件相符、满足守恒定律等基本原理约束下的预期行为模式一致[^1]。 #### 处理不收敛的情况 面对难以实现理想化快速收敛的情形时,可以尝试调整如下几个方面以改善性能表现: - 提升网格划分的质量并适当加密局部区域; - 审查设定好的边界条件及其合理性; - 尝试改变离散格式的选择顺序由低至高逐步过渡; - 减少松弛因子或控制数(Courant Number),特别是针对稳态问题而言; - 若涉及瞬态仿真则缩短时间步长有助于加快趋向稳定的速度; - 针对复杂的多物理场耦合场景采取渐进加载策略而不是一次性开启全部选项。 ```python # Python 示例代码片段展示如何读取 Fluent 日志文件获取残差信息 (假设日志结构固定) import re def parse_residual_from_log(log_file_path): pattern = r'^(?:\d+\.\d+)\s+(\S+)' with open(log_file_path, 'r') as file: lines = file.readlines() residuals = [] for line in lines: match = re.match(pattern, line.strip()) if match: residual_value = float(match.group(1)) residuals.append(residual_value) return residuals ```
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值