CN117574029B - 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 - Google Patents
高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 Download PDFInfo
- Publication number
- CN117574029B CN117574029B CN202410078856.4A CN202410078856A CN117574029B CN 117574029 B CN117574029 B CN 117574029B CN 202410078856 A CN202410078856 A CN 202410078856A CN 117574029 B CN117574029 B CN 117574029B
- Authority
- CN
- China
- Prior art keywords
- representing
- viscosity
- resolution
- vortex
- rans
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 238000000926 separation method Methods 0.000 claims abstract description 16
- 230000035939 shock Effects 0.000 claims abstract description 13
- 238000004088 simulation Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 2
- 230000002452 interceptive effect Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 8
- 238000005457 optimization Methods 0.000 abstract description 7
- 230000000452 restraining effect Effects 0.000 abstract description 3
- 208000012322 Raynaud phenomenon Diseases 0.000 abstract 3
- 238000010586 diagram Methods 0.000 description 5
- 230000006978 adaptation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及数据驱动的雷诺平均湍流场优化技术领域,尤其涉及高分辨率雷诺应力与雷诺平均Navier‑Stokes方程求解器的兼容方法,包括:S1:获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据;S2:耦合SST湍流模型的RANS求解器获得低分辨率的时均场数据,S3:计算高分辨率雷诺应力各向异性张量,S4:利用最小二乘法反算每个网格单元的涡粘,S5:再通过加权计算更新湍流粘性,并约束其上、下界,S6:运行RANS求解器,更新速度场和平均应变率张量;S7:循环S3‑S6,直至RANS速度场和DNS速度场之间的均方误差收敛。实现了高分辨率雷诺应力数据和基于SST模型低分辨率RANS流场的相容性处理,确保优化过程中迭代计算的稳定性,实现力系数的再次收敛,提高了RANS流场的准确度。
Description
技术领域
本发明涉及数据驱动的雷诺平均湍流场优化技术领域,尤其涉及高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法。
背景技术
目前,雷诺平均Navier-Stokes(RANS)方法仍然是航空航天工程领域计算湍流的主流方法。但是由于湍流模型的局限性,该方法对激波/边界层干扰等复杂分离流的预测偏差较大,难以满足工程实践的需求,基于高分辨率计算方法获取的雷诺应力数据在一定程度上可以改善RANS结果的准确度。现有技术基于“冻结”假设和最小二乘法实现了低速流动下的RANS流场的优化,但并不适用于超声速等高速流动,在高分辨率雷诺应力和RANS时均场迭代过程中容易出现误差积累造成求解器发散问题。
有鉴于此,如何基于高分辨率雷诺应力数据实现超声速激波/边界层干扰复杂分离流的RANS流场改善是本领域亟待解决的问题。
发明内容
有鉴于此,本发明提供高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,旨在实现高分辨率雷诺应力和RANS时场均迭代过程的稳定以及流场的收敛,进而实现改善RANS流场准确度的目的。
为解决以上技术问题,本发明的技术方案为提供高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,包括:
S1:获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据;
S2:将所述高分辨率雷诺应力数据耦合SST湍流模型的RANS求解器获得低分辨率的时均场数据,并将所述时均场数据设定为初始值,所述时均场数据包括第一涡粘和第一平均应变率张量;
S3:计算高分辨率雷诺应力各向异性张量并将数据读取至RANS求解器;
S4:利用最小二乘法反算每个网格单元的涡粘,得到涡粘反算值;
S5:将涡粘反算值与所述第一涡粘进行加权计算并约束其上、下界,得到约束后的第二涡粘;
S6:将所述第二涡粘赋值到所述RANS求解器中,并进行质量方程、动量方程、能量方程的迭代,更新速度场和平均应变率张量;
S7:循环S3-S6,直至RANS速度场和DNS速度场之间的均方误差收敛。
作为一种实施方式,所述获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据的方式包括:
通过直接数值模拟获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据,或,
通过大涡模拟获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据。
作为一种实施方式,所述SST湍流模型的控制方程被配置为通过方程1和方程2计算得到SST涡粘,其中,
方程1:
,
方程2:
,
SST涡粘:
,
其中,表示偏导运算,/>表示密度,/>表示湍功能,/>表示时间,/>表示速度,/>表示张量形式空间坐标轴,下标/>表示空间维度索引,/>表示动量粘性,/>表示t时刻的运动粘性,/>表示Bradshaw常数,/>表示湍流动量粘性,/>表示动能生成项,/>表示比耗散率生成项,/>表示涡量大小,/>,/>,/>,/>均为因变量,由相应下标1和2的常数参数通过混合函数得到,/>、/>为常数参数,/>和/>为经验函数,/>表示比耗散率。
作为一种实施方式,所述高分辨率雷诺应力各向异性张量的计算方程被配置为:,其中,/>表示雷诺应力张量,/>表示雷诺应力各向异性张量,/>表示湍功能,/>表示Kronecker符号,下标/>和/>表示空间维度索引。
作为一种实施方式,所述利用最小二乘法反算每个网格单元的涡粘,包括:
采用公式反算每个网格单元的涡粘,其中,/>表示第k步迭代的平均应变率张量,/>表示涡粘反算值,高分辨率雷诺应力各向异性张量/>,/>表示雷诺应力张量,/>表示湍功能,/>表示Kronecker符号,下标/>和/>表示空间维度索引。
作为一种实施方式,所述将涡粘反算值与所述第一涡粘进行加权计算并约束其上、下界之前,还包括:
设定加权系数的数值,在需要优化的流场区域设为1,在无需优化的流场区域设为0,交界处通过tanh函数予以光顺连接。
作为一种实施方式,所述第二涡粘被约束的上、下界为,其中表示约束后的第二涡粘,/>表示约束前的第二涡粘。
作为一种实施方式,所述均方误差的计算公式被配置为:,其中,上标DNS表示DNS速度场相关参数,RANS表示RANS速度场相关参数,/>、/>分别表示水平和垂直速度,/>表示优化网格单元的总数,下标/>表示优化网格单元索引。
本发明的首要改进之处:能够实现高分辨率雷诺应力数据和基于SST模型低分辨率RANS流场的相容性处理,确保优化过程中迭代计算的稳定性,实现力系数的再次收敛,提高了RANS流场的准确度,最终达到提高RANS方法对分离泡、力系数等关键信息预测准确度的目的。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图做简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1是本发明一实施例提供的高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法的步骤示意图;
图2是本发明一实施例提供的计算网格示意图;
图3是本发明一实施例提供的阻力系数演化效果图;
图4是本发明一实施例提供的优化开始后速度场均方误差变化图;
图5是本发明一实施例提供的壁面压力系数分布对比图;
图6是本发明一实施例提供的壁面摩擦阻力系数分布对比图。
具体实施方式
为了使本领域技术人员更好的理解本发明实施例,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1,是本发明一实施例提供的高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法的步骤示意图,具体包括以下步骤:
S1、获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据。
进一步的,通过直接数值模拟获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据,或,通过大涡模拟获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据。
S2、将高分辨率雷诺应力数据耦合SST湍流模型的RANS求解器获得低分辨率的时均场数据,并将时均场数据设定为初始值,时均场数据包括第一涡粘和第一平均应变率张量。
进一步的,SST湍流模型的控制方程被配置为通过方程1和方程2计算得到SST涡粘,其中,
方程1:
,
方程2:
,
SST涡粘:
,
在方程1、方程2和SST涡粘中,动能生成项和比耗散率生成项/>:
,/>,
参数和/>:
,/>,
参数、/>、/>和/>由相应下标1和2的常数参数通过混合函数计算得到,/>表示参数/>、/>、/>和/>,/>和/>表示相应下标1和2的常数参数,混合函数:
,
关于参数相应下标1和2的常数参数/>和/>:
,/>,
关于参数、/>和/>相应下标1和2的常数参数:
=0.85,/>=1.0,/>=0.075,/>=0.0828,/>=0.5,/>=0.856,
涡量大小:
,
在涡量大小中,旋转率张量/>:
,
在参数和参数/>中,参数/>和参数/>:
,/>,
在参数和参数/>中,双曲正切函数/>和交叉扩散项/>:
,/>,
上述方程和公式中,表示偏导运算,/>表示密度,/>表示湍功能,/>表示时间,/>表示速度,/>表示张量形式空间坐标轴,/>表示比耗散率,下标/>和/>表示空间维度索引,动量粘性/>,/>表示Kronecker符号,/>表示t时刻的运动粘性,/>表示物面距离,表示自然底数,常数参数/>=0.09,/>=0.41,Bradshaw常数/>,湍流动量粘性。
需要说明的是,通过上述计算获得的低分辨率的时均场数据,包括第一涡粘和第一平均应变率张量/>,将其作为初始值,有/>。
S3、计算高分辨率雷诺应力各向异性张量并将数据读取至RANS求解器。
进一步的,高分辨率雷诺应力各向异性张量的计算方程被配置为:,其中,/>表示雷诺应力张量,/>表示雷诺应力各向异性张量,/>表示湍功能,/>表示Kronecker符号,下标/>和/>表示空间维度索引。
S4、利用最小二乘法反算每个网格单元的涡粘,得到涡粘反算值。
进一步的,采用公式反算每个网格单元的涡粘,其中,/>表示第k步迭代的平均应变率张量,/>表示涡粘反算值,高分辨率雷诺应力各向异性张量,,/>表示雷诺应力张量,/>表示雷诺应力各向异性张量,/>表示湍功能,/>表示Kronecker符号,下标/>和/>表示空间维度索引。需要说明的是,本实施例采用的网格是如图2所示的24°压缩拐角计算网格。
S5、将涡粘反算值与所述第一涡粘进行加权计算并约束其上、下界,得到约束后的第二涡粘。
进一步的,将涡粘反算值与所述第一涡粘进行加权计算并约束其上、下界之前,需要设定加权系数的数值,在需要优化的流场区域设为1,在无需优化的流场区域设为0,交界处通过tanh函数予以光顺连接。
进一步的,第二涡粘被约束的上、下界为,其中/>表示约束后的第二涡粘,/>表示约束前的第二涡粘。由此得到约束后的第二涡粘为。
S6、将第二涡粘赋值到所述RANS求解器中,并进行质量方程、动量方程、能量方程的迭代,更新速度场和平均应变率张量。
将第二涡粘赋值到所述RANS求解器中,运行Navier-Stokes方程求解器,更新速度场和应变率张量。
S7、循环S3-S6,直至RANS速度场和DNS速度场之间的均方误差收敛。
进一步的,循环步骤S3-S6,直至RANS速度场和DNS速度场之间的均方误差收敛,均方误差的计算公式被配置为:,其中,/>、/>分别表示水平和垂直速度,/>表示优化网格单元的总数,下标/>表示优化网格单元索引。
需要说明的是,结合图3-图6所示的效果图和比对图可以很明显的发现,采用本发明提供的兼容方法可以进一步对速度场的均方误差实现二次收敛,提高了RANS流场的准确度。
本发明提供的高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,通过以下步骤S1:获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据;S2:将高分辨率雷诺应力数据耦合SST湍流模型的RANS求解器获得低分辨率的时均场数据,并将时均场数据设定为初始值,时均场数据包括第一涡粘和第一平均应变率张量;S3:计算高分辨率雷诺应力各向异性张量并将数据读取至RANS求解器;S4:利用最小二乘法反算每个网格单元的涡粘,得到涡粘反算值;S5:将涡粘反算值与所述第一涡粘进行加权计算并约束其上、下界,得到约束后的第二涡粘;S6:将第二涡粘赋值到RANS求解器中,并进行质量方程、动量方程、能量方程的迭代,更新速度场和平均应变率张量;S7:循环S3-S6,直至RANS速度场和DNS速度场之间的均方误差收敛。能够实现高分辨率雷诺应力数据和基于SST模型低分辨率RANS流场的相容性处理,确保优化过程中迭代计算的稳定性,实现力系数的再次收敛,提高了RANS流场的准确度,最终达到提高RANS方法对分离泡、力系数等关键信息预测准确度的目的。
需要说明的是,本发明中,“第一方面”、“第二方面”、“第三方面”、“第四方面”等中,术语“第一”、“第二”、“第三”、“第四”等仅用于描述目的,不能理解为指示或暗示相对重要性或数量,也不能理解为隐含指明所指示的技术特征的重要性或数量。而且“第一”、“第二”、“第三”、“第四”等仅起到非穷举式的列举描述目的,应当理解并不构成对数量的封闭式限定,亦或者是对特定名词的封闭式理解。
此外,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
以上仅是本发明的优选实施方式,应当指出的是,上述优选实施方式不应视为对本发明的限制,本发明的保护范围应当以权利要求所限定的范围为准。对于本技术领域的普通技术人员来说,在不脱离本发明的精神和范围内,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,其特征在于,包括:
S1:获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据;
S2:将所述高分辨率雷诺应力数据耦合SST湍流模型的RANS求解器获得低分辨率的时均场数据,并将所述时均场数据设定为初始值,所述时均场数据包括第一涡粘和第一平均应变率张量,所述SST湍流模型的控制方程被配置为通过方程1和方程2计算得到SST涡粘,其中,
方程1:
,
方程2:
,
SST涡粘:
,
其中,表示偏导运算,/>表示密度,/>表示湍功能,/>表示时间,/>表示速度,/>表示张量形式空间坐标轴,下标/>表示空间维度索引,/>表示动量粘性,/>表示t时刻的运动粘性,/>表示Bradshaw常数,/>表示湍流动量粘性,/>表示动能生成项,/>表示比耗散率生成项,/>表示涡量大小,/>,/>,/>,/>均为因变量,由相应下标1和2的常数参数通过混合函数得到,/>、/>为常数参数,/>和/>为经验函数,/>表示比耗散率;
S3:计算高分辨率雷诺应力各向异性张量并将数据读取至RANS求解器,所述高分辨率雷诺应力各向异性张量的计算方程被配置为:,其中,/>表示雷诺应力张量,/>表示雷诺应力各向异性张量,/>表示湍功能,/>表示Kronecker符号,下标/>和/>表示空间维度索引;
S4:利用最小二乘法反算每个网格单元的涡粘,得到涡粘反算值,其中,采用公式反算每个网格单元的涡粘,/>表示第k步迭代的平均应变率张量,/>表示涡粘反算值;
S5:将涡粘反算值与所述第一涡粘进行加权计算并约束其上、下界,得到约束后的第二涡粘;
S6:将所述第二涡粘赋值到所述RANS求解器中,并进行质量方程、动量方程、能量方程的迭代,更新速度场和平均应变率张量;
S7:循环S3-S6,直至RANS速度场和DNS速度场之间的均方误差收敛。
2.根据权利要求1所述的高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,其特征在于,所述获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据的方式包括:
通过直接数值模拟或大涡模拟获取高速激波/边界层干扰复杂分离流的高分辨率雷诺应力数据。
3.根据权利要求1所述的高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,其特征在于,所述将涡粘反算值与所述第一涡粘进行加权计算并约束其上、下界之前,还包括:
设定加权系数的数值,在需要优化的流场区域设为1,在无需优化的流场区域设为0,交界处通过tanh函数予以光顺连接。
4.根据权利要求3所述的高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,其特征在于,所述第二涡粘被约束的上、下界为,其中/>表示约束后的第二涡粘,/>表示约束前的第二涡粘。
5.根据权利要求1所述的高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法,其特征在于,所述均方误差的计算公式被配置为:,其中,上标DNS表示DNS速度场相关参数,RANS表示RANS速度场相关参数,/>、/>分别表示水平和垂直速度,/>表示优化网格单元的总数,下标/>表示优化网格单元索引。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410078856.4A CN117574029B (zh) | 2024-01-19 | 2024-01-19 | 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410078856.4A CN117574029B (zh) | 2024-01-19 | 2024-01-19 | 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117574029A CN117574029A (zh) | 2024-02-20 |
CN117574029B true CN117574029B (zh) | 2024-04-26 |
Family
ID=89890462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410078856.4A Active CN117574029B (zh) | 2024-01-19 | 2024-01-19 | 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117574029B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010182287A (ja) * | 2008-07-17 | 2010-08-19 | Steven C Kays | 適応型インテリジェント・デザイン |
CN107273593A (zh) * | 2017-06-01 | 2017-10-20 | 北京航空航天大学 | 一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 |
CN109299569A (zh) * | 2018-10-24 | 2019-02-01 | 广州市香港科大霍英东研究院 | 一种基于相干结构的不可压缩黏性流体的大涡模拟方法 |
CN111324993A (zh) * | 2020-02-21 | 2020-06-23 | 苏州浪潮智能科技有限公司 | 一种湍流场更新方法、装置及其相关设备 |
AU2020102748A4 (en) * | 2020-10-16 | 2020-12-10 | Tianjin University | Intelligent feedback real-time control system and method for hydrodynamic circulation under external interference |
CN114595643A (zh) * | 2022-03-04 | 2022-06-07 | 重庆大学 | 耦合多点测风仪和微尺度稳态模拟的高分辨率山地风场测量方法 |
CN115186608A (zh) * | 2022-07-12 | 2022-10-14 | 北京航空航天大学 | 一种基于湍流能谱耦合rsm模型的网格自适应湍流模拟方法 |
CN115438598A (zh) * | 2022-09-02 | 2022-12-06 | 中国人民解放军国防科技大学 | 基于一般时间根方尺度的雷诺应力湍流模型的数值方法 |
CN115438597A (zh) * | 2022-09-02 | 2022-12-06 | 中国人民解放军国防科技大学 | 基于一般时间根方尺度的两方程涡黏湍流模型的数值方法 |
CN115859781A (zh) * | 2022-11-15 | 2023-03-28 | 重庆大学 | 基于注意力和卷积神经网络编解码器的流场预测方法 |
CN115859760A (zh) * | 2022-04-29 | 2023-03-28 | 中广核风电有限公司 | 流场湍流的模拟方法和装置 |
CN115935546A (zh) * | 2022-12-16 | 2023-04-07 | 扬州大学 | 一种基于分岔法的旋转曲率修正的自适应ddes模型构建方法 |
CN116309919A (zh) * | 2023-04-03 | 2023-06-23 | 长沙理工大学 | 多切片组多被试fMRI张量的约束CCPD算法 |
-
2024
- 2024-01-19 CN CN202410078856.4A patent/CN117574029B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010182287A (ja) * | 2008-07-17 | 2010-08-19 | Steven C Kays | 適応型インテリジェント・デザイン |
CN107273593A (zh) * | 2017-06-01 | 2017-10-20 | 北京航空航天大学 | 一种用于高马赫数强激波流场气动热预测的湍流模型及其建立方法 |
CN109299569A (zh) * | 2018-10-24 | 2019-02-01 | 广州市香港科大霍英东研究院 | 一种基于相干结构的不可压缩黏性流体的大涡模拟方法 |
CN111324993A (zh) * | 2020-02-21 | 2020-06-23 | 苏州浪潮智能科技有限公司 | 一种湍流场更新方法、装置及其相关设备 |
AU2020102748A4 (en) * | 2020-10-16 | 2020-12-10 | Tianjin University | Intelligent feedback real-time control system and method for hydrodynamic circulation under external interference |
CN114595643A (zh) * | 2022-03-04 | 2022-06-07 | 重庆大学 | 耦合多点测风仪和微尺度稳态模拟的高分辨率山地风场测量方法 |
CN115859760A (zh) * | 2022-04-29 | 2023-03-28 | 中广核风电有限公司 | 流场湍流的模拟方法和装置 |
CN115186608A (zh) * | 2022-07-12 | 2022-10-14 | 北京航空航天大学 | 一种基于湍流能谱耦合rsm模型的网格自适应湍流模拟方法 |
CN115438597A (zh) * | 2022-09-02 | 2022-12-06 | 中国人民解放军国防科技大学 | 基于一般时间根方尺度的两方程涡黏湍流模型的数值方法 |
CN115438598A (zh) * | 2022-09-02 | 2022-12-06 | 中国人民解放军国防科技大学 | 基于一般时间根方尺度的雷诺应力湍流模型的数值方法 |
CN115859781A (zh) * | 2022-11-15 | 2023-03-28 | 重庆大学 | 基于注意力和卷积神经网络编解码器的流场预测方法 |
CN115935546A (zh) * | 2022-12-16 | 2023-04-07 | 扬州大学 | 一种基于分岔法的旋转曲率修正的自适应ddes模型构建方法 |
CN116309919A (zh) * | 2023-04-03 | 2023-06-23 | 长沙理工大学 | 多切片组多被试fMRI张量的约束CCPD算法 |
Non-Patent Citations (4)
Title |
---|
Considerations of Stress Limiter for the SST Turbulence Model in Dual Throat Nozzle Predictions;El-Samanoudy R. Tharwat1 M. 等;《Ninth International Conference on Computational Fluid Dynamics (ICCFD9) 》;20160715;1-15 * |
Quantification of Reynolds-averaged-Navier–Stokes model-form uncertainty in transitional boundary layer and airfoil flows;Chu Minghan 等;《Physics of Fluids》;20221003;第34卷(第10期);1-23 * |
旋转流体的并行计算与分析;赵运业;《中国优秀硕士学位论文全文数据库基础科学辑》;20230115(第01期);A004-66 * |
智能空气动力学若干研究进展及展望;唐志共 等;《空气动力学学报》;20230715;第41卷(第07期);1-35 * |
Also Published As
Publication number | Publication date |
---|---|
CN117574029A (zh) | 2024-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marsh et al. | The influence of turbulence model and two and three-dimensional domain selection on the simulated performance characteristics of vertical axis tidal turbines | |
Alkishriwi et al. | A large-eddy simulation method for low Mach number flows using preconditioning and multigrid | |
Mannion et al. | A two and three-dimensional CFD investigation into performance prediction and wake characterisation of a vertical axis turbine | |
Breuer et al. | Flow around a surface mounted cubical obstacle: Comparison of LES and RANS-results | |
Chen | Applications of lattice Boltzmann method to turbulent flow around two-dimensional airfoil | |
Coppel et al. | Simulating the fluid dynamic behaviour of oar blades in competition rowing | |
Zha et al. | An improved higher-order moving particle semi-implicit method for simulations of two-dimensional hydroelastic slamming | |
Borsboom et al. | A generating‐absorbing boundary condition for dispersive waves | |
Zhang et al. | Research on data assimilation strategy of turbulent separated flow over airfoil | |
CN117574029B (zh) | 高分辨率雷诺应力与雷诺平均Navier-Stokes方程求解器的兼容方法 | |
Duc et al. | Predicting aerodynamic performance of savonius wind turbine: An application of generalized k-ω turbulence model | |
Wang et al. | Simulation of turbulent flow around a surface-mounted finite square cylinder | |
Myong | Numerical investigation of fully developed turbulent fluid flow and heat transfer in a square duct | |
Ferrandis et al. | Influence of viscosity and non-linearities in predicting motions of a wind energy offshore platform in regular waves | |
CN113158339B (zh) | 一种针对sst湍流模型的湍流长度尺度修正方法 | |
Zhang et al. | A finite particle method based on a Riemann solver for modeling incompressible flows | |
JiSeok et al. | Flow around a flexible plate in a free stream | |
Yang et al. | Calculation of propeller’s load noise using LES and BEM numerical acoustics coupling methods | |
Amaral et al. | Comprehensive verification and validation of a CFD analysis | |
Boschitsch et al. | Aeroelastic Analysis Using Deforming Cartesian Grids | |
Löwe et al. | hp-adaptation for a discontinuous Galerkin method in an unstructured flow solver | |
Lin et al. | Prediction of aerodynamic performance of NREL offshore 5‐MW baseline wind turbine considering power loss at varying wind speeds | |
Tong et al. | A Cartesian grid method with improvement of resolving the boundary layer structure for two‐dimensional incompressible flows | |
Karbon et al. | Turbulence Modeling Using ZF, RSM, LES and WMLES for Flow Analysis in Z-Shape Ducts | |
Tupkar et al. | CFD Analysis of Water Replenishment Holes in an Offshore Wind Turbine Foundation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |