CN112685978B - 一种适用于五次样条重构格式的自适应人工粘性控制方法 - Google Patents
一种适用于五次样条重构格式的自适应人工粘性控制方法 Download PDFInfo
- Publication number
- CN112685978B CN112685978B CN202110264089.2A CN202110264089A CN112685978B CN 112685978 B CN112685978 B CN 112685978B CN 202110264089 A CN202110264089 A CN 202110264089A CN 112685978 B CN112685978 B CN 112685978B
- Authority
- CN
- China
- Prior art keywords
- interface
- physical quantity
- unit
- state value
- artificial viscosity
- 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 55
- 230000004907 flux Effects 0.000 claims abstract description 25
- 239000006185 dispersion Substances 0.000 claims abstract description 12
- 230000003044 adaptive effect Effects 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000000853 adhesive Substances 0.000 claims description 3
- 230000001070 adhesive effect Effects 0.000 claims description 3
- 239000012530 fluid Substances 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种适用于五次样条重构格式的自适应人工粘性控制方法,包括:步骤1、读取流场数据,求解五次样条重构方程组得到单元界面处物理量的状态值及其一阶导数;步骤2、计算出单元界面处的二阶至五阶导数;步骤3、根据波数识别方法,求出单元界面处流场的等效无量纲波数,从而确定该处流场的人工粘性系数;步骤4、根据单元界面处物理量的状态值计算无粘数值通量,并根据人工粘性系数添加六阶人工粘性项,最后采用相对应的时间离散方法在时间上进行推进。将局部瞬时流场特征与格式的解析能力结合起来,通过判断局部瞬时的流场特征是否在格式的可解析范围内,从而确定格式的人工粘性系数的取值,使格式的人工粘性控制更为合理。
Description
技术领域
本发明涉及计算流体力学中数值计算方法领域,特别涉及一种适用于五次样条重构格式的自适应人工粘性控制方法。
背景技术
在计算流体力学中,数值格式的耗散起着非常重要的作用。数值耗散越大,计算越稳定,但对流场的刻画能力越弱;数值格式的耗散越小,流场的刻画能力就越强,计算越不稳定。因此如何调节格式的数值耗散一直是人们研究的重点。传统调节耗散的方法是调节格式的总体耗散,这种调节耗散的方式并不能完全满足人们的需求。理想的数值格式应当满足以下条件:在流动光滑区域数值耗散尽量小以保证格式精度和分辨率,在流动非光滑区或者梯度较大的区域数值耗散要足够大以保证计算稳定。因此发展根据局部流场信息确定格式耗散大小的自适应耗散方法具有广泛的科研以及应用前景。
发明内容
针对现有技术中存在的问题,提供了一种适用于五次样条重构格式的自适应人工粘性方法;该方法将流场的局部流动特征同格式的色散特性结合起来,由两者共同决定局部流场需要的人工粘性系数的大小,实现人工粘性系数的自适应改变。首先根据波数识别方法,由表征局部流场特征的密度的高阶导数求出该处流场的局部等效无量纲波数,再将等效无量纲波数与格式的修正无量纲波数进行对比。如果等效无量纲波数在格式的完全解析范围内,则人工粘性系数为零或者取一个很小的值;如果等效无量纲波数不在格式的完全解析范围内,则人工粘性系数根据一定的规则逐渐增大,直至最大值。这个判定过程即为自适应人工粘性判据。这种自适应人工粘性方法改变了以往的全流场使用统一的人工粘性系数且人工粘性系数的大小根据经验确定的方式,使得人工粘性系数的取值更为合理。
本发明采用的技术方案如下:一种适用于五次样条重构格式的自适应人工粘性控制方法,包括以下步骤:
步骤1、读取流场数据,求解五次样条重构方程组得到单元界面处物理量的状态值及其一阶导数;
步骤2、根据单元界面处物理量的状态值及一阶导数,计算出单元界面处的二阶至五阶导数;
步骤3、根据波数识别方法,由单元界面处流场密度的高阶导数求出单元界面处流场的等效无量纲波数,代入自适应人工粘性判据确定该处流场的人工粘性系数;
步骤4、根据单元界面处物理量的状态值计算无粘数值通量,并根据人工粘性系数添加六阶人工粘性项,最后采用相对应的时间离散方法在时间上进行推进。
进一步的,在结构网格有限体积方法的基本框架下,不计质量力和源项,在直坐标系下欧拉方程为:
在求解过程中根据单元的平均值重构出单元界面处的物理量的状态值,五次样条重构方程组为:
其中,
下标分别代表着三个方向,和
分别代表单元和单元从单元左侧界面中心至单元右侧界面中心的距离,分别指的五个守恒物理量,和分别指在界面处第个物理量的界面状态值及其一阶导数,和分别指单
元和单元处的第个物理量的单元平均值,通过求解重构方程,可得到单元界面处物理量的状态值及其一阶导数;上述五次样条重构方
程组、界面状态值及一阶导数均为方向,方向和方向需要采用相同的方法求解。
进一步的,波数识别方法为:
进一步的,自适应人工粘性判据为:
进一步的,所述步骤4具体包括:将单元界面处物理量的状态值带入控制方程计算无粘数值通量,无粘通量计算过程为:
与现有技术相比,采用上述技术方案的有益效果为:本发明的方案中人工粘性系数不再是全场统一的值,而是由自适应人工粘性判据根据局部流场的特征来确定,自适应人工粘性判据局部瞬时的流场特征和格式的色散特性共同来决定,此种方法将局部瞬时流场特征与格式的解析能力结合起来,通过判断局部瞬时的流场特征是否在格式的可解析范围内,从而确定格式的人工粘性系数的取值,使格式的人工粘性控制更为合理。
附图说明
图1是本发明的适用于五次样条重构格式的自适应人工粘性方法流程图。
图2是本发明中的等效无量纲波数与修正无量纲波数对比图。
图3是本发明中自适应人工粘性系数随波数变化曲线图。
图4是本发明中与自适应人工粘性系数对应的数值格式耗散曲线。
图5是采用一种固定人工粘性的二维黎曼问题密度等值线图。
图6是采用自适应人工粘性的二维黎曼问题密度等值线图。
图7是采用另一固定人工粘性的二维黎曼问题密度等值线图。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明提出了一种适用于五次样条重构格式的自适应人工粘性方法,包括下述理论依据:
(1)波数识别方法
由式 (1)可得:
将式(2)进行无量纲化,为防止导数为零造成分母为零,等效无量纲波数可写出如下形式:
(2)波数识别在五次样条重构格式中的验证
而根据傅里叶分析得到的五次样条格式的修正无量纲波数为:
图2显示了由高阶导数求出的等效无量纲波数与五次样条格式修正无量纲波数对比图,其中纵坐标表示色散特性,由图中可以看出在修正无量纲波数的可识别范围内,等效无量纲波数与修正无量纲波数符合的很好,说明由波数识别方法得到的等效无量纲波数在一定的波数范围内可以用来代表格式的修正无量纲波数。
(3)自适应人工粘性判据
由于等效无量纲波数在一定波数范围内可以代表格式的修正无量纲波数,因此可
以根据图2进行自适应人工粘性判据的设定:根据修正的无量纲波数的特性设一阈值,
当时,流场等效无量纲波数在五次样条重构格式的精确可识别范围内,人工粘
性系数取零或者很小;当时,此时五次样条重构格式无法精确识别流场波数,
需要增大人工粘性系数来抑制误差。前者根据Lele et al提出的色散误差容限来确定,当时,。后者满足Hu
& Adams提出的色散耗散关系,由此得出的自适应判据如下:
图3为人工粘性系数随波数的变化曲线,可以看出随着波数的增大,人工粘性随之变化。图4为自适应人工粘性系数对应的耗散曲线,其中纵坐标表示格式的耗散特性,可以看出自适应耗散在低波数区域耗散尽量小,不引入耗散误差;在高波数区域尽量大以抑制色散误差带来的波动。
基于上述理论依据,具体方案如下:
如图1所示,一种适用于五次样条重构格式的自适应人工粘性控制方法,包括以下步骤:
步骤1、读取流场数据,求解五次样条重构方程组得到单元界面处物理量的状态值及其一阶导数;
步骤2、根据单元界面处物理量的状态值及一阶导数,计算出单元界面处的二阶至五阶导数;
步骤3、根据波数识别方法,由单元界面处流场密度的高阶导数求出单元界面处流场的等效无量纲波数,代入自适应人工粘性判据确定此处流场的人工粘性系数;
步骤4、根据单元界面处物理量的状态值计算无粘数值通量,并根据人工粘性系数添加六阶人工粘性项,最后采用相对应的时间离散方法在时间上进行推进。
具体的,
步骤(1),
在结构网格有限体积方法的基本框架下,不计质量力和源项,在直坐标系下欧拉方程为:
在求解过程中根据单元的平均值重构出单元界面处的物理量的状态值,五次样条重构方程组为:
其中,
下标分别代表着三个方向,和
分别代表单元和单元从单元左侧界面中心至单元右侧界面中心的距离,分别指的五个守恒物理量,和分别指在界面处第个物理量的界面状态值及其一阶导数,和分别指单
元和单元处的第个物理量的单元平均值,通过求解重构方程,可得到单元界面处物理量的状态值及其一阶导数;上述五次样条重构方
程组、界面状态值及一阶导数均为方向,方向和方向采用相同的方法可进行求解。
步骤(3):式(18)得到单元界面处物理量的二阶至五阶导数后,由式(5)可得到流场界面处的等效无量纲波数,此处求解过程中使用密度的高阶导数。根据自适应人工粘性判据式(7),确定界面处的人工粘性系数。
步骤(4):根据控制方程求解无粘数值通量
将单元界面处物理量的状态值带入控制方程计算无粘数值通量,无粘通量计算过程为:
最后求出通量以后,采用相对应的时间离散方法在时间上进行推进。
在本实施例中,给出以下自适应人工粘性的数值验证,如图5、6、7为采用固定人工
粘性的二维黎曼问题的密度等值线图和使用自适应人工粘性的二维黎曼问题密度等值线
图,由图中可以看出:图5时的耗散太小,无法抑制高频振荡,且计算流场不
对称;图7时的结果耗散太大,流场结构不丰富;图6采用自适应人工粘性
后的流场具有对称且更为丰富的结构。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。如果本领域技术人员,在不脱离本发明的精神所做的非实质性改变或改进,都应该属于本发明权利要求保护的范围。
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
本说明书中公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
Claims (1)
1.一种适用于五次样条重构格式的自适应人工粘性控制方法,其特征在于,包括以下步骤:
步骤1、读取流场数据,求解五次样条重构方程组得到单元界面处物理量的状态值及其一阶导数;
步骤2、根据单元界面处物理量的状态值及一阶导数,计算出单元界面处的二阶至五阶导数;
步骤3、根据波数识别方法,由单元界面处流场密度的高阶导数求出单元界面处流场的等效无量纲波数,代入自适应人工粘性判据确定该处流场的人工粘性系数;
步骤4、根据单元界面处物理量的状态值计算无粘数值通量,并根据人工粘性系数添加六阶人工粘性项,最后采用相对应的时间离散方法在时间上进行推进;
所述步骤1中求解过程具体包括:
在结构网格有限体积方法的基本框架下,不计质量力和源项,在直坐标系下欧拉方程为:
在求解过程中根据单元的平均值重构出单元界面处的物理量的状态值,五次样条重构方程组为:
其中,
下标分别代表着三个方向,和分别代表单元和单元从单元左侧界面中心至单元右侧界面中心的距离,分别指的五个守恒物理量,和分别指在界面处第个物理量的界面状态值及其一阶导数,和分别指单元和单元处的第个物理量的单元平均值,通过求解重构方程,可得到单元界面处物理量的状态值及其一阶导数;上述五次样条重构方程组、界面状态值及一阶导数均为方向;
波数识别方法为:
自适应人工粘性判据为:
所述步骤4具体包括:将单元界面处物理量的状态值带入控制方程计算无粘数值通量,无粘通量计算过程为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110264089.2A CN112685978B (zh) | 2021-03-11 | 2021-03-11 | 一种适用于五次样条重构格式的自适应人工粘性控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110264089.2A CN112685978B (zh) | 2021-03-11 | 2021-03-11 | 一种适用于五次样条重构格式的自适应人工粘性控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112685978A CN112685978A (zh) | 2021-04-20 |
CN112685978B true CN112685978B (zh) | 2021-06-08 |
Family
ID=75458390
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110264089.2A Active CN112685978B (zh) | 2021-03-11 | 2021-03-11 | 一种适用于五次样条重构格式的自适应人工粘性控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112685978B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117966B (zh) * | 2021-12-09 | 2023-04-11 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种物理属性与数据驱动耦合的流声模态分解及预测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102460113A (zh) * | 2009-05-11 | 2012-05-16 | 西安大略大学 | 监测介质的颗粒尺寸分布的超声法方法 |
CN102954910A (zh) * | 2011-08-26 | 2013-03-06 | 航天科工防御技术研究试验中心 | 防热涂层拉伸剪切强度测试模型及其制作方法 |
CN104298797A (zh) * | 2013-07-16 | 2015-01-21 | 中国石油化工股份有限公司 | 一种确定缝洞型油藏高导流通道圈闭剩余油的方法 |
CN108197072A (zh) * | 2017-12-27 | 2018-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于加权守恒变量阶跃的高精度间断Galerkin人工粘性激波捕捉方法 |
CN109032077A (zh) * | 2018-09-05 | 2018-12-18 | 沈阳建筑大学 | 一种基于刀具姿态控制的五轴数控加工指令点插补方法 |
CN109614577A (zh) * | 2018-12-10 | 2019-04-12 | 山东大学苏州研究院 | 一种Burgers方程求解方法及装置 |
US20200109740A1 (en) * | 2018-10-06 | 2020-04-09 | Toyota Jidosha Kabushiki Kaisha | Spline telescopic shaft of vehicular propeller shaft |
WO2020165059A1 (en) * | 2019-02-11 | 2020-08-20 | Miltenyi Biotec B.V. & Co. KG | Generation of human pluripotent stem cell derived artificial tissue structures without three dimensional matrices |
CN112052632A (zh) * | 2020-07-27 | 2020-12-08 | 空气动力学国家重点实验室 | 一种高超声速流向转捩预测方法 |
CN112214869A (zh) * | 2020-09-03 | 2021-01-12 | 空气动力学国家重点实验室 | 一种求解欧拉方程的改进型高阶非线性空间离散方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104317997A (zh) * | 2014-10-17 | 2015-01-28 | 北京航空航天大学 | 一种高负荷风扇/压气机端壁造型优化设计方法 |
CN110209048A (zh) * | 2019-05-20 | 2019-09-06 | 华南理工大学 | 基于动力学模型的机器人时间最优轨迹规划方法、设备 |
-
2021
- 2021-03-11 CN CN202110264089.2A patent/CN112685978B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102460113A (zh) * | 2009-05-11 | 2012-05-16 | 西安大略大学 | 监测介质的颗粒尺寸分布的超声法方法 |
CN102954910A (zh) * | 2011-08-26 | 2013-03-06 | 航天科工防御技术研究试验中心 | 防热涂层拉伸剪切强度测试模型及其制作方法 |
CN104298797A (zh) * | 2013-07-16 | 2015-01-21 | 中国石油化工股份有限公司 | 一种确定缝洞型油藏高导流通道圈闭剩余油的方法 |
CN108197072A (zh) * | 2017-12-27 | 2018-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于加权守恒变量阶跃的高精度间断Galerkin人工粘性激波捕捉方法 |
CN109032077A (zh) * | 2018-09-05 | 2018-12-18 | 沈阳建筑大学 | 一种基于刀具姿态控制的五轴数控加工指令点插补方法 |
US20200109740A1 (en) * | 2018-10-06 | 2020-04-09 | Toyota Jidosha Kabushiki Kaisha | Spline telescopic shaft of vehicular propeller shaft |
CN109614577A (zh) * | 2018-12-10 | 2019-04-12 | 山东大学苏州研究院 | 一种Burgers方程求解方法及装置 |
WO2020165059A1 (en) * | 2019-02-11 | 2020-08-20 | Miltenyi Biotec B.V. & Co. KG | Generation of human pluripotent stem cell derived artificial tissue structures without three dimensional matrices |
CN112052632A (zh) * | 2020-07-27 | 2020-12-08 | 空气动力学国家重点实验室 | 一种高超声速流向转捩预测方法 |
CN112214869A (zh) * | 2020-09-03 | 2021-01-12 | 空气动力学国家重点实验室 | 一种求解欧拉方程的改进型高阶非线性空间离散方法 |
Non-Patent Citations (2)
Title |
---|
Investigation of a Low-Toxicity Energetic Binder for a Solid Propellant: Curing, Microstructures, and Performance;Song Ma等;《ACS Omega》;20201201;第5卷(第47期);30538–30548 * |
三维边界层流动失稳与Bypass转捩模式研究;徐国亮;《中国博士学位论文全文数据库础科学辑》;20121215(第12期);A004-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN112685978A (zh) | 2021-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107272403A (zh) | 一种基于改进粒子群算法的pid控制器参数整定算法 | |
CN110851929B (zh) | 基于自适应网格的二维叶型优化设计方法及装置 | |
CN111079228B (zh) | 一种基于流场预测的气动外形优化方法 | |
CN104778327B (zh) | 基于人工神经网络的飞机翼型优化设计方法 | |
CN106126791B (zh) | 一种考虑几何不确定性的高超声速机翼气动力/热分析方法 | |
CN102540882B (zh) | 一种基于最小参数学习法的飞行器航迹倾角控制方法 | |
CN112685978B (zh) | 一种适用于五次样条重构格式的自适应人工粘性控制方法 | |
CN111006843B (zh) | 一种暂冲式超声速风洞的连续变速压方法 | |
CN113850008B (zh) | 飞行器气动特性预测的自适应网格扰动域更新加速方法 | |
CN112214869B (zh) | 一种求解欧拉方程的改进型高阶非线性空间离散方法 | |
CN107640183B (zh) | 一种基于迭代学习控制的列车运行控制方法 | |
CN113127976B (zh) | 一种边界层分离诱导转捩预测方法、装置、设备及介质 | |
CN113569360B (zh) | 一种风力机叶片抗颤振翼型簇设计方法 | |
CN104834772A (zh) | 基于人工神经网络的飞机翼型/机翼反设计方法 | |
CN112001109A (zh) | 再生核粒子算法实现结构冲击动力学仿真方法 | |
CN110765706A (zh) | 基于ohngbm(1,1)的翼型非定常失速气动系数建模方法 | |
CN106446419A (zh) | 一种火力发电厂燃煤锅炉的建模方法及系统 | |
CN107329131A (zh) | 一种利用粒子滤波的雷达微弱目标检测跟踪方法 | |
CN116373846A (zh) | 一种基于bp神经网络优化的后轮转向车辆稳定性控制方法 | |
CN110852888A (zh) | 一种基于粒子滤波的证券投资组合优化方法 | |
CN114185276B (zh) | 一种基于多维泰勒网的非线性严格系统输出反馈控制方法 | |
CN113977571B (zh) | 一种柔性关节机器人输出力矩控制方法 | |
CN113158339B (zh) | 一种针对sst湍流模型的湍流长度尺度修正方法 | |
Zhuo et al. | Parameter Identification of Tire Model Based on Improved Particle Swarm Optimization Algorithm | |
CN111597741B (zh) | 一种基于改进Hicks-henne算法的仿生蟹滑翔姿态下翼型优化设计方法 |
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 |