一种汽车护风罩的气动性能仿真计算方法与流程

本发明是一种汽车护风罩气动性能仿真计算方法，特别涉及一种用于与环形风扇相适配使用的汽车护风罩总成的仿真计算。

背景技术：

汽车发动机运转过程产生大量热量使得发动机内各部件温度上升，降低发动机工作效率，长期处于过热温度增加零件损坏风险，缩短发动机的使用寿命。现代的发动机散热系统在水箱散热的基础上，配以冷却风扇起到更好的散热效果。提升冷却风扇的工作效率能够以更低的能耗带走更多散热量，高效率风扇研究一直是发动机散热领域的一个重点。

在实际应用中，通常有两个措施提升冷却风扇散热效果：一是优化风扇的设计，开发性能更为优异的冷却风扇，这是提升风扇散热效果的主要途径，二是开发与风扇相适配的护风罩，由于风扇周围流场复杂，配备适合的汽车护风罩能充分发挥冷却风扇效果，从而提高散热效率。目前风扇造型研究已经非常充分，仅从风扇造型角度研究提升风扇的性能的效果也非常有限，因此对汽车发动机护风罩的研究有重要的意义。

传统的护风罩设计方法，需要通过大量实验获得数据，研究的周期长，成本的投入也比较高，目前商业化流体仿真软件发展很成熟，仿真结果的精度也有很大提高，因此通过数值计算方法进行汽车护风罩研究具有重要意义。

技术实现要素：

本发明提供一种汽车护风罩模拟计算方法，特别是对环形风扇适配的护风罩的气动性能计算，建立仿真分析模型，并与实验结果对标，在保证模型准确性基础上，对不同结构参数的护风罩仿真，优化护风罩设计。

本发明的技术方案是：

一种汽车护风罩的气动性能仿真计算方法，包括以下步骤：

步骤一、使用三维建模软件catia建立冷却风扇与护风罩的几何模型，取z轴为风扇的旋转轴，风扇旋转方向为正转方向，z轴负方向为来流方向；

步骤二、计算区域模型建立与网格划分，计算模型的尺寸与管道实体尺寸相一致，为提高fluent仿真计算精度，风扇叶片与护风圈面网格采用三角形网格，最大尺寸控制在10mm以内，风扇前、后缘尺寸较小，网格尺寸控制在1mm以内；整个计算域模型体网格分为5个部分，顺着流线方向依次为入口区，入口过渡区，旋转区，出口过渡区，出口区域。风扇面网格为三角形网格，过渡区流场复杂，在旋转区和过渡区采用四面体网格可以提升网格质量和计算精度；入口区与出口区流量较为平稳，在入口区域出口区网格采用六面体网格可以减小网格数量，保证计算精度；

步骤三、在gambit中设立面边界条件，包括入口边界条件，出口边界条件，监测面，交界面，风扇表面边界条件；

步骤四、采用基于有限体积法的商业软件fluent计算给定风扇转速下流体的纳维-斯托克斯方程，设置求解参数、边界条件、初值条件、控制方程求解器、离散方法、参考系类型、压力梯度的解耦方法；

步骤五、根据数值计算的结果与实验数据对标，误差在可接受范围内时，可认定模型正确，若误差超出接受范围，则返回步骤二，重新调整网格再次计算；基于涡动力学和风扇周边流场特性，建立具有不同筋条几何尺寸参数护风圈模型，仿真计算得到最佳筋条布置的护风圈设计方案。

进一步地，所述步骤二中还包括进行网格无关性验证的步骤，若两套网格计算结果的偏差在1％之内，则认为计算结果与网格数量无关，无需进一步增加网格数量。

进一步地，所述步骤三中的入口边界条件采用质量入口边界条件，出口边界条件采用压力出口边界条件。

进一步地，所述步骤四中，采用基于压力求解器，先求解动量方程获得速度场，继而进行压力修正以满足连续性方程，从而保证流场同时满足动力方程与连续性方程；采用rngk-epsilon湍流模型；采用基于格林高斯节点的速度梯度插值方案；采用simple算法进行压力-速度耦合方程求解，采用标准格式离散求解连续性方程，采用一阶迎风格式离散求解动量方程，采用二阶迎风格式离散求解湍流动能与湍流耗散率；采用多参考系方法，旋转区采用动参考系，其余区域采用静坐标系。

进一步地，所述步骤五中，通过编辑journal文件给定的不同质量流量入口边界条件，设定fluent在每个流量点自动保存计算的结果，并设定自动读取下个流量点的入口边界条件，避免大量的手动操作。

进一步地，所述步骤五，若仿真计算的结果与实验数据偏差≤6％，可认为模型具有足够准确性。

相比现有技术，本发明有以下积极效果：

(1)网格无关性验证验证仿真计算模型网格划分尺寸的合理性，减小计算量，提高优化设计效率；

(2)通过分析旋转区大小对结果的影响，保证计算结果的具有高精度；

(3)利用fluent的journal文件，自动读取不同入口流量大小，保存结果数据，减小的重复性的工作；

(4)本发明根据计算结果，给出了一种新型护风罩的设计方案，能够起到更好的护风作用。

附图说明

图1是流场区域网格。

图2是改进前护风罩的正视图。

图3是改进后护风罩的正视图。

图4是改进前护风罩的左视图半剖示意图。

图5是改进后护风罩的左视图半剖示意图。

具体实施方式

下面通过具体实施例对本发明的目的作进一步详细地描述，实施例不能在此一一赘述，但本发明的实施方式并不因此限定于以下实施例。

一种汽车护风罩的气动性能仿真计算方法，为适配的环形风扇相护风罩的优化设计提供一气动性能数值计算方法，按以下步骤实现：

步骤一：风扇护风罩实物模型建立；

使用三维建模软件catia建立冷却风扇与护风罩的几何模型，取z轴为风扇的旋转轴，风扇旋转方向为正转方向，z轴负方向为来流方向。本实施例建模采用7叶片、直径为286mm，护风罩内径为292mm。

步骤二：仿真模型网格划分；

计算区域模型建立和网格划分，计算模型尺寸与风管实体尺寸相一致。风管入口直径为550mm，风管长6000mm，自由流场区域为圆柱形，直径为4000mm，风扇叶片与护风圈面网格采用三角形网格，最大尺寸控制在10mm以内，风扇前、后缘尺寸较小，面网格尺寸控制在1mm以内(见图1)。

整个计算域模型体网格分为5个部分，顺着流线方向依次为入口区，入口过渡区，旋转区，出口过渡区，出口区域。风扇面网格为三角形网格，过渡区流场复杂，在旋转区和过渡区采用四面体网格，要求四面体网格的volskew低于0.75；入口区与出口区流量较为平稳，在入口区域出口区网格采用六面体网格可以减小迭代误差，要求六面体网格jacobian高于0.6。

网格无关性验证：

影响仿真计算结果准确性的因素主要有几个方面：模型误差、不同湍流模型、离散误差、迭代误差等几个方面。离散误差随着网格总数增加而降低，但同时迭代误差随网格总数增加而增加，网格数量过多可能导致模型收敛性变差。因此在满足计算准确性前提下，减小网格数量，加快收敛速度可以缩短计算周期。通常两套疏密不同的网格模型计算结果偏差在1％以内可认为计算结果与网格数量总量无关。网格无关性验证结果如表1，可见当网格总数大于128万后，两套网格偏差为0.26％，因此之后计算网格数量控制在128万以上可以保证有较高的计算精度。

表1网格无关性验证结果

旋转区大小确定：

采用多坐标系方法进行计算时，在动坐标系中求解动量方程需要考虑科氏力和旋转离心力作用，而静坐标系无需考虑这两个力作用。通常认为旋转区划分时不宜太小，旋转区尺寸必须要将风扇周边复杂的流线区包含在内，否则计算结果会有明显偏差。定义lu为前旋转区到风扇轴向距离，ld为后旋转区到风扇轴向距离，r为风扇半径，旋转区大小分析结果如表2所示，可见当lu/r＝0.03时，大小流量点的静压偏差都比较大，但lu/r取值为0.1和0.2时，静压的偏差都在6％以下，且当lu/r值取为0.1时，大流量点的误差为-4.2％，小流量点误差为5.8％，该方案在大流量点的误差更小，当lu/r为0.2时，这种旋转区在小流量点的准确性更高，考虑到该风扇的设计流量点更接近0.514(m³/s)，在之后计算中选取lu/r＝0.1左右，可认为仿真结果的准确性满足计算要求。

表2不同旋转区大小误差分析

步骤三：边界条件确定；

在gambit中设立面边界条件，包括入口边界条件，出口边界条件，监测面，交界面，风扇表面边界条件。相比于压力入口，质量入口有更快收敛速度，故之后计算中边界条件选用质量入口条件，出口边界条件选用压力出口；在四面体网格与六面体网格交接面设为interface面。

步骤四：fluent仿真计算；

采用基于有限体积法的商业软件fluent计算给定风扇转速下流体的纳维-斯托克斯方程，湍流模型选用rngk-epsilon，利用simple算法求解，离散格式为一阶精度，求解器具体设置如下：

表3求解器参数设置

步骤五：利用fluent的journal脚本文件实现仿真的自动计算，只要启动fluent后读取journal，按照设定入口边界条件以及自动保存计算结果，计算结果通过文本文件读取。根据fluent数值计算结果，基于涡动力学理论和叶片周围流场特性优化护风圈结构。

经过上述计算后的到的护风罩的结构如图3和图5所示，优化前的护风罩结构如图2和图4所示。

本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。