一种涡特征驱动的旋转湍流PANS模型制造技术

本发明专利技术公开了一种涡特征驱动的旋转湍流PANS模型。其实现过程主要包括：引入伪涡度数以区分局部流体的运动形式；计算涡特征指数以实现流场特性的自适应追踪；计算旋转湍流中基于涡特征指数的湍动能滤波因子及耗散率滤波因子；计算近壁区涡粘阻尼以保证旋转湍流一致性条件；以经典的RNG k‑ε模型为基底，通过修正输运方程和经验系数确定新PANS模型。本发明专利技术引入了涡特征驱动参数，仅通过联系流场涡特征即可实现对旋转湍流的高效解析，其湍动能和耗散率滤波因子的动态调节无需联系网格空间尺度，能有效降低数值模拟对网格的要求，能更好地兼顾计算精度与计算效率，可为流体机械内旋转湍流的高效计算提供新思路。

A pan model of rotational turbulence driven by vortex characteristics

The invention discloses a pan model of rotational turbulence driven by vortex characteristics. The implementation process mainly includes: introducing pseudo vorticity to distinguish the motion form of local fluid; calculating the vorticity characteristic index to realize the adaptive tracking of flow field characteristics; calculating the turbulent kinetic energy filter factor and dissipation rate filter factor based on the vorticity characteristic index in rotating turbulence; calculating the eddy viscosity damping near the wall to ensure the consistency of rotating turbulence; taking the classical RNG K \u2011 \u03b5 model as the base, The new pans model is determined by modifying the transport equation and empirical coefficient. The invention introduces the vortex characteristic driving parameter, which can realize the efficient analysis of the rotating turbulence only by connecting the vortex characteristic of the flow field. The dynamic adjustment of the turbulent kinetic energy and the dissipation rate filtering factor does not need to contact the space scale of the grid, which can effectively reduce the requirements of the numerical simulation on the grid, give better consideration to the calculation accuracy and efficiency, and improve the efficient calculation of the rotating turbulence in the fluid machinery For new ideas.

【技术实现步骤摘要】
一种涡特征驱动的旋转湍流PANS模型
本专利技术涉及工程流体力学计算领域，特别是涉及一种涡特征驱动的旋转湍流PANS模型。
技术介绍
旋转湍流是水轮机、水泵、汽轮机和压气机等诸多流体机械中的典型流动现象。受强旋转、大曲率、多壁面等因素的影响，流体机械中的旋转湍流表现出非线性、强瞬态、高脉动、宽尺度等特点，在各种流体动力学研究中，此种流动被认为是最复杂的流动之一。因此，开发可兼顾计算精度与计算效率的湍流模型一直是工程计算流体力学领域的研究热点。近年来，PANS模型作为一种可实现从RANS到DNS光滑过渡的混合方法，在流体机械中的应用越来越多。该方法能在控制计算效率的前提下有效提高旋转湍流涡结构的解析尺度与精度，对流体机械内复杂流动的研究具有重要意义。然而，目前PANS模型的动态计算基本都依赖于网格空间尺度和泰勒尺度的相对大小，即在联系流场特征的同时还与计算域的网格条件相关，这将对网格数量和质量提出一定的要求，进而限制数值计算的效率和鲁棒性。为此，开发仅依赖于流场特征的PANS模型是十分必要的，它能自然地降低对网格的要求，进而更好地实现流体机械内旋转湍流的高效计算。
技术实现思路
(一)要解决的技术问题本专利技术的目的是提供一种涡特征驱动的旋转湍流PANS模型，旨在降低PANS方法对网格的要求，在保证计算精度的同时进一步提高旋转湍流的计算效率，为流体工程的高效推进奠定科学基础。(二)技术方案为了解决上述技术问题，本专利技术提供一种涡特征驱动的旋转湍流PANS模型，包括如下步骤：步骤一、引入伪涡度数Vq以区分局部流体的主导运动形式；步骤二、计算涡特征指数EC以实现流场特性的自适应追踪；步骤三、计算旋转湍流中基于所述涡特征指数EC的湍动能滤波因子fk及耗散率滤波因子fε；步骤四、计算近壁区涡粘阻尼fw以保证旋转湍流一致性条件；步骤五、以经典的RNGk-ε模型为基底，推导得出新PANS模型所需的模化湍动能输运方程、模化耗散率输运方程和经验系数，并与雷诺方程结合得到所述新PANS模型；所述伪涡度数Vq由下式得出：式中，||TR||F为旋转率张量TR的F-范数，||TD||F为应变率张量TD的F-范数，为速度梯度张量，γ为一正小量以保证分母非零，可直接取10-10；所述涡特征指数EC由下式得出：式中，Vq为所述伪涡度数；所述湍动能滤波因子fk由下式得出：fk＝lg[0.0852·cos(4.1·EC)+2.7644·sin(4.1·EC)+3.8959]所述耗散率滤波因子fε由下式得出：fε＝lg[1.9822·cos(1.5·EC)+1.4255·sin(1.5·EC)+7.3504]所述近壁区涡粘阻尼fw由下式得出：fw＝1-exp[1301.51·cos(0.001·y*)-0.1788·sin(0.001·y*)-1301.51]式中，y*为近壁距比率，其中，所述近壁距比率由下式得出：y*＝30·δ/Y1式中，δ为近壁计算节点至壁面的实际距离，Y1为第一层计算节点至壁面的实际距离；Y1需保证在对数律层内缘，δ需保证为Y1的4倍；所述新PANS模型所需的模化湍动能输运方程、模化耗散率输运方程和新经验系数由下式关系推导得出：式中，k为所述经典的RNGk-ε模型中的原湍动能，ε为所述经典的RNGk-ε模型中的原耗散率，Cμ为所述经典的RNGk-ε模型中的原经验系数，km为所述新PANS模型中所需的模化湍动能，εm为所述新PANS模型中所需的模化耗散率，Cμm为所述新PANS模型中所需的新经验系数。(三)有益效果本专利技术提供了一种涡特征驱动的旋转湍流PANS模型，本专利技术引入了涡特征驱动参数，仅通过联系流场涡特征即可实现对旋转湍流的高效解析，其湍动能和耗散率滤波因子的动态调节无需联系网格空间尺度，能有效降低数值模拟对网格的要求，能更好地兼顾计算精度与计算效率，可为流体机械内旋转湍流的高效计算提供新思路。附图说明图1为旋流突扩管算例计算域示意图；图2为监测线上轴向速度平均值对比图；图3为监测线上轴向速度脉动值对比图；图4为SSTk-ω模型计算所得湍流涡结构图；图5为本专利技术实施例提供的PANS模型计算所得湍流涡结构图；图6为所述涡特征驱动的旋转湍流PANS模型实现方法流程图。具体实施方式下面结合附图和实施例，对本专利技术的具体实施方式作进一步详细描述。以下实例用于说明本专利技术，但不用来限制本专利技术的范围。在本专利技术的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本专利技术中的具体含义。本专利技术提供一种涡特征驱动的旋转湍流PANS模型，包括如下步骤：步骤一、引入伪涡度数以区分局部流体的运动形式；在这一步骤中，按下式引入伪涡度数Vq以区分局部流体的主导运动形式，其中，伪涡度数由下式得出：式中，||TR||F为旋转率张量TR的F-范数，||TD||F为应变率张量TD的F-范数，为速度梯度张量，γ为一正小量以保证分母非零，可直接取10-10；步骤二、计算涡特征指数以实现流场特性的自适应追踪；在这一步骤中，按下式引入涡特征指数EC以计算得出实现流场特性的自适应追踪；式中，Vq为伪涡度数；步骤三、计算旋转湍流中基于涡特征指数的湍动能滤波因子及耗散率滤波因子；在这一步骤中，湍动能滤波因子由下式得出：fk＝lg[0.0852·cos(4.1·EC)+2.7644·sin(4.1·EC)+3.8959]式中，fk为湍动能滤波因子，lg[]为常用对数函数，cos()为余弦函数，sin()为正弦函数；耗散率滤波因子由下式得出：fε＝lg[1.9822·cos(1.5·EC)+1.4255·sin(1.5·EC)+7.3504]式中，fε为耗散率滤波因子，lg[]为常用对数函数，cos()为余弦函数，sin()为正弦函数；步骤四、计算近壁区涡粘阻尼以保证旋转湍流一致性条件；在这一步骤中，按下式确定近壁区涡粘阻尼fw以保证旋转湍流一致性条件：fw＝1-exp[1301.51·cos(0.001·y*)-0.1788·sin(0.001·y*)-1301.51]exp[]为自然指数函数，cos()为余弦函数，sin()为正弦函数，其中，y*为近壁距比率，其计算式为：y*＝30·δ/Y1式中，δ为近壁计算节点至壁面的实际距离，Y1为第一层计算节点至壁面的实本文档来自技高网...

【技术保护点】
1.一种涡特征驱动的旋转湍流PANS模型，其特征在于，包括如下步骤：/n步骤一、引入伪涡度数V

【技术特征摘要】
1.一种涡特征驱动的旋转湍流PANS模型，其特征在于，包括如下步骤：
步骤一、引入伪涡度数Vq以区分局部流体的主导运动形式；
步骤二、计算涡特征指数EC以实现流场特性的自适应追踪；
步骤三、计算旋转湍流中基于所述涡特征指数EC的湍动能滤波因子fk及耗散率滤波因子fε；
步骤四、计算近壁区涡粘阻尼fw以保证旋转湍流一致性条件；
步骤五、以经典的RNGk-ε模型为基底，推导得出新PANS模型所需的模化湍动能输运方程、模化耗散率输运方程和经验系数，并与雷诺方程结合得到所述新PANS模型；
所述伪涡度数Vq由下式得出：



式中，||TR||F为旋转率张量TR的F-范数，||TD||F为应变率张量TD的F-范数，▽V为速度梯度张量，γ为一正小量以保证分母非零，可直接取10-10；
所述涡特征指数EC由下式得出：



式中，Vq为所述伪涡度数；
所述湍动能滤波因子fk由下式得出：
fk＝lg[0.0852·cos(4.1·EC)+2.7644·sin(4.1·EC)...

【专利技术属性】
技术研发人员：王福军，王超越，陶然，汤远，赵浩儒，王本宏，闫婷婷，何庆儒，
申请(专利权)人：中国农业大学，
类型：发明
国别省市：北京;11