CN114611421A - 基于模态衰减的人工黏性方法及系统 - Google Patents

基于模态衰减的人工黏性方法及系统 Download PDF

Info

Publication number
CN114611421A
CN114611421A CN202210143246.9A CN202210143246A CN114611421A CN 114611421 A CN114611421 A CN 114611421A CN 202210143246 A CN202210143246 A CN 202210143246A CN 114611421 A CN114611421 A CN 114611421A
Authority
CN
China
Prior art keywords
modal
attenuation
module
variable
grid
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.)
Granted
Application number
CN202210143246.9A
Other languages
English (en)
Other versions
CN114611421B (zh
Inventor
伍彬
于剑
张迪
张青
马万超
秦学宇
李欣益
何威
段旭
贺祥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Electromechanical Engineering
Original Assignee
Shanghai Institute of Electromechanical Engineering
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Electromechanical Engineering filed Critical Shanghai Institute of Electromechanical Engineering
Priority to CN202210143246.9A priority Critical patent/CN114611421B/zh
Publication of CN114611421A publication Critical patent/CN114611421A/zh
Application granted granted Critical
Publication of CN114611421B publication Critical patent/CN114611421B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • 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)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

本发明提供了一种基于模态衰减的人工黏性方法及系统，所述方法包括如下步骤：步骤S1：建立控制方程；步骤S2：使用最小二乘的方法计算衰减率τ；步骤S3：使周围网格守恒变量q组成新的守恒变量qex代替原守恒变量q，并组合得到新的人工黏性模型；步骤S4：通过得到的守恒变量qex和黏性系数μ，得到粘性通量g，基于间断有限元的理论求出随时间t变化的守恒变量q，进而得到整个流场的结果。本发明的人工黏性模型可以很好的捕捉流场中的激波位置，从而有效抑制激波间断处的伪振荡；通过结合当前网格及其相邻网格的变量，增强用于估计衰减率的信息，适用于低阶的计算。

Description

基于模态衰减的人工黏性方法及系统
技术领域
本发明涉及计算流体力学的技术领域，具体地，涉及基于模态衰减的人工黏性方法及系统，尤其涉及基于模态衰减的人工黏性模型。
背景技术
在计算流体力学中，高阶格式对于复杂流动的数值模拟具有重要意义。以旋翼流动的模拟为例，对于该类流动，翼梢涡是其主要流场结构之一，采用通常的结构网格计算。而高阶格式对于这类问题则显示出了绝对的优势，采用高阶格式得到的旋翼气动特性和噪声的计算精度都有了显著的提高，并且其计算效率也比在低阶格式的基础上加密网格的效率要高。
高阶格式面临的一大挑战是：在可压缩流动中，以极小的数值耗散来分辨小尺度脉动，在光滑区保持格式的高精度特性不被破坏，同时光滑地捕捉流场中出现的各种间断。混合使用不同的高阶格式来计算，可以解决这一难题。混合格式是将在光滑流场中具有低耗散特性的高阶格式，和可以识别流场间断的激波捕捉格式相组合，在保证低耗散和高精度的同时，也能够准确地捕捉流场中的激波。激波捕捉方法一直是CFD研究的重点和难点，在有限差分和有限体积框架下，研究者们发展了大量经典的激波捕捉方法，包括限制器类方法、重构类方法、人工黏性方法和混合方法。
人工黏性方法是通过添加人工黏性项来抑制非物理振荡，其关键点是构造合适的光滑度量因子和选取适当的人工黏性系数。一种思路是采用某些变量(比如密度、速度等)的导数及其派生变量来作为光滑度量因子。该类方法特点是鲁棒性较强，并且由于针对不同的流动结构设计了不同的度量因子，因此能够灵活的适用于复杂流动。然而，对于更高阶数的导数可以更好的保证光滑区域的计算精度，但这会使得算法复杂度和计算代价大大增加。
而通过模态系数的衰减速度来构造人工黏性，方法构造形式简单，易于理解和实现。其基本思想是光滑流动区域高阶模态系数，应该迅速衰减至较小量级。激波捕捉方法可以在平滑区域保持相应计算格式的精度，对于非平滑的流动，该方法可以很好地抑制激波间断处的伪振荡。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种基于模态衰减的人工黏性方法及系统。
根据本发明提供的一种基于模态衰减的人工黏性方法，所述方法包括如下步骤：
步骤S1：根据所模拟的具体流动，建立控制方程；
步骤S2：由模态衰减的间断指示器，构造模态分量函数
Figure BDA0003507441820000027
然后建立关于衰减率τ的目标函数，并使用最小二乘的方法计算衰减率τ；
步骤S3：使周围网格守恒变量q组成新的守恒变量qex代替原守恒变量q，同时由衰减率τ得到新的人工黏性模型，并计算出黏性系数μ；
步骤S4：通过得到的守恒变量qex和黏性系数μ，得到粘性通量g，基于间断有限元的理论求出随时间t变化的守恒变量q，进而得到整个流场的结果。
优选地，所述步骤S1包括如下步骤：
步骤S1.1：求解一般的双曲方程：
Figure BDA0003507441820000021
其中t为时间变量，q为守恒变量，f为对流通量，g为粘性通量；
步骤S1.2：粘性通量g表示为：
Figure BDA0003507441820000022
其中μ为人工黏性系数；
步骤S1.3：将守恒变量q分解为时间项和空间项：
Figure BDA0003507441820000023
其中Np和r分别表示单元网格中节点个数和坐标向量，m为单元网格中节点的序号，
Figure BDA0003507441820000024
模态分量函数，
Figure BDA0003507441820000025
为一般正交基函数；
步骤S1.4：根据模态衰减的间断指示器，构造模态分量函数为：
Figure BDA0003507441820000026
其中C为已知参数，τ为衰减率。
优选地，所述步骤S2包括如下步骤：
步骤S2.1：构造一个最优的模态衰减系数bm：
Figure BDA0003507441820000031
步骤S2.2：将原模态分量函数
Figure BDA0003507441820000032
与最优模态衰减系数bm组合，构造新的模态系数为
Figure BDA0003507441820000033
Figure BDA0003507441820000034
其中，
Figure BDA0003507441820000035
表示为在网格单元K中取守恒变量q的标准二范数L2；
步骤S2.3：通过如下公式，模态系数
Figure BDA0003507441820000036
依次单调：
Figure BDA0003507441820000037
其中，n表示为网格节点的一组序号；
步骤S2.4：使用最小二乘的方法计算衰减率τ，目标函数选取：
Figure BDA0003507441820000038
优选地，所述步骤S3包括如下步骤：
步骤S3.1：将该网格单元守恒变量q和左右两个网格单元的守恒变量qL、qR构造一个新的守恒变量qex：
qex(r,t)＝γLqL(r,t)+γq(r,t)+γRqR(r,t) (9)
其中，γL，γ，γR表示网格单元守恒变量对应的加权系数；
步骤S3.2：构造的守恒变量值光滑变化，加权系数γL，γ，γR选择如下所示：
Figure BDA0003507441820000039
其中参数r*为：
Figure BDA00035074418200000310
其中，P为计算的精度，P*和Pmax表示计算中所需的精度参数；
步骤S3.3：在计算时如果耗散过度，会破坏高阶格式的精度，计算黏性系数μ的方法为：
Figure BDA0003507441820000041
其中μmax为最大的黏性系数，h为网格尺度，P为计算的精度。
本发明还提供一种基于模态衰减的人工黏性系统，所述系统包括如下模块：
模块M1：根据所模拟的具体流动，建立控制方程；
模块M2：由模态衰减的间断指示器，构造模态分量函数
Figure BDA0003507441820000042
然后建立关于衰减率τ的目标函数，并使用最小二乘的系统计算衰减率τ；
模块M3：使周围网格守恒变量q组成新的守恒变量qex代替原守恒变量q，同时由衰减率τ得到新的人工黏性模型，并计算出黏性系数μ；
模块M4：通过得到的守恒变量qex和黏性系数μ，得到粘性通量g，基于间断有限元的理论求出随时间t变化的守恒变量q，进而得到整个流场的结果。
优选地，所述模块M1包括如下模块：
模块M1.1：求解一般的双曲方程：
Figure BDA0003507441820000043
其中t为时间变量，q为守恒变量，f为对流通量，g为粘性通量；
模块M1.2：粘性通量g表示为：
Figure BDA0003507441820000044
其中μ为人工黏性系数；
模块M1.3：将守恒变量q分解为时间项和空间项：
Figure BDA0003507441820000045
其中Np和r分别表示单元网格中节点个数和坐标向量，m为单元网格中节点的序号，
Figure BDA0003507441820000046
模态分量函数，
Figure BDA0003507441820000047
为一般正交基函数；
模块M1.4：根据模态衰减的间断指示器，构造模态分量函数为：
Figure BDA0003507441820000048
其中C为已知参数，τ为衰减率。
优选地，所述模块M2包括如下模块：
模块M2.1：构造一个最优的模态衰减系数bm：
Figure BDA0003507441820000051
模块M2.2：将原模态分量函数
Figure BDA0003507441820000052
与最优模态衰减系数bm组合，构造新的模态系数为
Figure BDA0003507441820000053
Figure BDA0003507441820000054
其中，
Figure BDA0003507441820000055
表示为在网格单元K中取守恒变量q的标准二范数L2；
模块M2.3：通过如下公式，模态系数
Figure BDA0003507441820000056
依次单调：
Figure BDA0003507441820000057
其中，n表示为网格节点的一组序号；
模块M2.4：使用最小二乘的系统计算衰减率τ，目标函数选取：
Figure BDA0003507441820000058
优选地，所述模块M3包括如下模块：
模块M3.1：将该网格单元守恒变量q和左右两个网格单元的守恒变量qL、qR构造一个新的守恒变量qex：
qex(r,t)＝γLqL(r,t)+γq(r,t)+γRqR(r,t) (9)
其中，γL，γ，γR表示网格单元守恒变量对应的加权系数；
模块M3.2：构造的守恒变量值光滑变化，加权系数γL，γ，γR选择如下所示：
Figure BDA0003507441820000059
其中参数r*为：
Figure BDA00035074418200000510
其中，P为计算的精度，P*和Pmax表示计算中所需的精度参数；
模块M3.3：在计算时如果耗散过度，会破坏高阶格式的精度，计算黏性系数μ的系统为：
Figure BDA0003507441820000061
其中μmax为最大的黏性系数，h为网格尺度，P为计算的精度。
与现有技术相比，本发明具有如下的有益效果：
1、该人工黏性模型可以很好的捕捉流场中的激波位置，从而有效抑制激波间断处的伪振荡，减小计算误差；
2、通过结合当前网格及其相邻网格的变量，在不增加过多计算量下，增强用于估计衰减率的信息，适用于低阶的计算。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明方法完整算法流程图；
图2为本发明方法算例结果图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
本发明的目的是将应用模态衰减模型，包含周围网格单元的信息，扩展到高阶计算精度，构造新的人工黏性。本发明使用了形式简单的模态衰减模型，通过加权函数将周围流场信息组合。本发明不仅能保证准确捕捉流场中的激波位置，也适用与低阶格式的计算。
本发明提出了一种基于模态衰减的人工黏性模型，其流程图如附图1和算例结果如附图2所示，它包含以下四个步骤：
步骤S1：根据所模拟的具体流动，建立控制方程，具体包括如下子步骤：
步骤S1.1：求解一般的双曲方程：
Figure BDA0003507441820000071
其中t为时间变量，q为守恒变量，f为对流通量，g为粘性通量；
步骤S1.2：粘性通量g表示为：
Figure BDA0003507441820000072
其中μ为人工黏性系数；
步骤S1.3：将守恒变量q分解为时间项和空间项：
Figure BDA0003507441820000073
其中Np和r分别表示单元网格中节点个数和坐标向量，m为单元网格中节点的序号，
Figure BDA0003507441820000074
模态分量函数，
Figure BDA0003507441820000075
为一般正交基函数；
步骤S1.4：根据模态衰减的间断指示器，构造模态分量函数为：
Figure BDA0003507441820000076
其中C为已知参数，τ为衰减率。
步骤S2：由模态衰减的间断指示器，构造模态分量函数
Figure BDA0003507441820000077
然后建立关于衰减率τ的目标函数，并使用最小二乘的方法计算衰减率τ，具体包括如下子步骤：
步骤S2.1：构造一个最优的模态衰减系数bm：
Figure BDA0003507441820000078
步骤S2.2：将原模态分量函数
Figure BDA0003507441820000079
与最优模态衰减系数bm组合，构造新的模态系数为
Figure BDA00035074418200000710
Figure BDA00035074418200000711
其中，
Figure BDA00035074418200000712
表示为在网格单元K中取守恒变量q的标准二范数L2；
步骤S2.3：通过如下公式，模态系数
Figure BDA00035074418200000713
依次单调：
Figure BDA00035074418200000714
其中，n表示为网格节点的一组序号；
步骤S2.4：使用最小二乘的方法计算衰减率τ，目标函数选取：
Figure BDA0003507441820000081
步骤S3：使周围网格守恒变量q组成新的守恒变量qex代替原守恒变量q，同时由衰减率τ得到新的人工黏性模型，并计算出黏性系数μ，具体计算步骤如下：
步骤S3.1：将该网格单元守恒变量q和左右两个网格单元的守恒变量qL、qR构造一个新的守恒变量qex：
qex(r,t)＝γLqL(r,t)+γq(r,t)+γRqR(r,t) (9)
其中，γL，γ，γR表示网格单元守恒变量对应的加权系数；
步骤S3.2：构造的守恒变量值光滑变化，加权系数γL，γ，γR选择如下所示：
Figure BDA0003507441820000082
其中参数r*为：
Figure BDA0003507441820000083
其中，P为计算的精度，P*和Pmax表示计算中所需的精度参数；
步骤S3.3：在计算时如果耗散过度，会破坏高阶格式的精度，计算黏性系数μ的方法为：
Figure BDA0003507441820000084
其中μmax为最大的黏性系数，h为网格尺度，P为计算的精度。
步骤S4：将得到的守恒变量qex和黏性系数μ，代入方程(2)可以得到粘性通量g。之后基于间断有限元的理论可以求出随时间t变化的守恒变量q，进而可以得到整个流场的结果。
本发明还提供一种基于模态衰减的人工黏性系统，所述系统包括如下模块：
模块M1：根据所模拟的具体流动，建立控制方程。
模块M1.1：求解一般的双曲方程：
Figure BDA0003507441820000091
其中t为时间变量，q为守恒变量，f为对流通量，g为粘性通量；
模块M1.2：粘性通量g表示为：
Figure BDA0003507441820000092
其中μ为人工黏性系数；
模块M1.3：将守恒变量q分解为时间项和空间项：
Figure BDA0003507441820000093
其中Np和r分别表示单元网格中节点个数和坐标向量，m为单元网格中节点的序号，
Figure BDA0003507441820000094
模态分量函数，
Figure BDA0003507441820000095
为一般正交基函数；
模块M1.4：根据模态衰减的间断指示器，构造模态分量函数为：
Figure BDA0003507441820000096
其中C为已知参数，τ为衰减率。
模块M2：由模态衰减的间断指示器，构造模态分量函数
Figure BDA0003507441820000097
然后建立关于衰减率τ的目标函数，并使用最小二乘的系统计算衰减率τ。
模块M2.1：构造一个最优的模态衰减系数bm：
Figure BDA0003507441820000098
模块M2.2：将原模态分量函数
Figure BDA0003507441820000099
与最优模态衰减系数bm组合，构造新的模态系数为
Figure BDA00035074418200000910
Figure BDA00035074418200000911
其中，
Figure BDA00035074418200000912
表示为在网格单元K中取守恒变量q的标准二范数L2；
模块M2.3：通过如下公式，模态系数
Figure BDA00035074418200000913
依次单调：
Figure BDA00035074418200000914
其中，n表示为网格节点的一组序号；
模块M2.4：使用最小二乘的系统计算衰减率τ，目标函数选取：
Figure BDA0003507441820000101
模块M3：使周围网格守恒变量q组成新的守恒变量qex代替原守恒变量q，同时由衰减率τ得到新的人工黏性模型，并计算出黏性系数μ。
模块M3.1：将该网格单元守恒变量q和左右两个网格单元的守恒变量qL、qR构造一个新的守恒变量qex：
qex(r,t)＝γLqL(r,t)+γq(r,t)+γRqR(r,t) (9)
其中，γL，γ，γR表示网格单元守恒变量对应的加权系数；
模块M3.2：构造的守恒变量值光滑变化，加权系数γL，γ，γR选择如下所示：
Figure BDA0003507441820000102
其中参数r*为：
Figure BDA0003507441820000103
其中，P为计算的精度，P*和Pmax表示计算中所需的精度参数；
模块M3.3：在计算时如果耗散过度，会破坏高阶格式的精度，计算黏性系数μ的系统为：
Figure BDA0003507441820000104
其中μmax为最大的黏性系数，h为网格尺度，P为计算的精度。
模块M4：通过得到的守恒变量qex和黏性系数μ，得到粘性通量g，基于间断有限元的理论求出随时间t变化的守恒变量q，进而得到整个流场的结果。
算例：本算例是检验基于模态衰减构造高阶的激波捕捉方法数值模拟性能的。计算区域为-5＜x＜5，计算网格量为200个，格式的精度阶数P＝3，边界条件为远场边条，计算终止时间为t＝1.8。初始的密度ρ、速度u、压强p为：
Figure BDA0003507441820000105
图2为本实施例模拟结果的密度沿x轴分布曲线。图2中可以看出，本发明清晰地捕捉了流场中的激波(x＜0.5)，说明本发明能够稳定、高精度地捕捉流场中的强间断；对于多尺度湍流结构嵌套区域0.5＜x＜2.5，本发明能够较好地捕捉几个极值点，说明本发明对多尺度湍流结构的分辨能力较好。
本发明提出了一种基于模态衰减的人工黏性模型，而通过模态系数的衰减速度来构造人工黏性，方法构造形式简单，易于理解和实现。包含周围网格单元的信息，扩展到高阶计算精度，构造新的人工黏性。使用了形式简单的模态衰减模型，通过加权函数将周围流场信息组合。不仅能保证准确捕捉流场中的激波位置，也可以提高计算的鲁棒性。构造的人工黏性模型可以在平滑区域保持相应计算格式的精度，对于非平滑的流动，该方法可以很好地抑制激波间断处的伪振荡。
本发明的核心算法是将应用模态衰减模型，包含周围网格单元的信息，扩展到高阶计算精度，构造新的人工黏性。使用了形式简单的模态衰减模型，通过加权函数将周围流场信息组合。在低阶格式的计算时，不会因为模板信息过少而对流场激波捕捉的不够准确，因而构造的人工黏性可以很好的适用于低阶格式的计算。如果对本发明的各种改动和变形不脱离本发明范围，仍然属于本发明的权利要求和同等技术范围之内。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (8)

1.一种基于模态衰减的人工黏性方法，其特征在于，所述方法包括如下步骤：
步骤S1：根据所模拟的具体流动，建立控制方程；
步骤S2：由模态衰减的间断指示器，构造模态分量函数
Figure FDA0003507441810000011
然后建立关于衰减率τ的目标函数，并使用最小二乘的方法计算衰减率τ；
步骤S3：使周围网格守恒变量q组成新的守恒变量qex代替原守恒变量q，同时由衰减率τ得到新的人工黏性模型，并计算出黏性系数μ；
步骤S4：通过得到的守恒变量qex和黏性系数μ，得到粘性通量g，基于间断有限元的理论求出随时间t变化的守恒变量q，进而得到整个流场的结果。
2.根据权利要求1所述的基于模态衰减的人工黏性方法，其特征在于，所述步骤S1包括如下步骤：
步骤S1.1：求解一般的双曲方程：
Figure FDA0003507441810000012
其中t为时间变量，q为守恒变量，f为对流通量，g为粘性通量；
步骤S1.2：粘性通量g表示为：
Figure FDA0003507441810000013
其中μ为人工黏性系数；
步骤S1.3：将守恒变量q分解为时间项和空间项：
Figure FDA0003507441810000014
其中Np和r分别表示单元网格中节点个数和坐标向量，m为单元网格中节点的序号，
Figure FDA0003507441810000015
模态分量函数，
Figure FDA0003507441810000016
为一般正交基函数；
步骤S1.4：根据模态衰减的间断指示器，构造模态分量函数为：
Figure FDA0003507441810000017
其中C为已知参数，τ为衰减率。
3.根据权利要求1所述的基于模态衰减的人工黏性方法，其特征在于，所述步骤S2包括如下步骤：
步骤S2.1：构造一个最优的模态衰减系数bm：
Figure FDA0003507441810000021
步骤S2.2：将原模态分量函数
Figure FDA0003507441810000022
与最优模态衰减系数bm组合，构造新的模态系数为
Figure FDA0003507441810000023
Figure FDA0003507441810000024
其中，
Figure FDA0003507441810000025
表示为在网格单元K中取守恒变量q的标准二范数L2；
步骤S2.3：通过如下公式，模态系数
Figure FDA0003507441810000026
依次单调：
Figure FDA0003507441810000027
其中，n表示为网格节点的一组序号；
步骤S2.4：使用最小二乘的方法计算衰减率τ，目标函数选取：
Figure FDA0003507441810000028
4.根据权利要求1所述的基于模态衰减的人工黏性方法，其特征在于，所述步骤S3包括如下步骤：
步骤S3.1：将该网格单元守恒变量q和左右两个网格单元的守恒变量qL、qR构造一个新的守恒变量qex：
qex(r,t)＝γLqL(r,t)+γq(r,t)+γRqR(r,t) (9)
其中，γL，γ，γR表示网格单元守恒变量对应的加权系数；
步骤S3.2：构造的守恒变量值光滑变化，加权系数γL，γ，γR选择如下所示：
Figure FDA0003507441810000029
其中参数r*为：
Figure FDA00035074418100000210
其中，P为计算的精度，P*和Pmax表示计算中所需的精度参数；
步骤S3.3：在计算时如果耗散过度，会破坏高阶格式的精度，计算黏性系数μ的方法为：
Figure FDA0003507441810000031
其中μmax为最大的黏性系数，h为网格尺度，P为计算的精度。
5.一种基于模态衰减的人工黏性系统，其特征在于，所述系统包括如下模块：
模块M1：根据所模拟的具体流动，建立控制方程；
模块M2：由模态衰减的间断指示器，构造模态分量函数
Figure FDA0003507441810000032
然后建立关于衰减率τ的目标函数，并使用最小二乘的系统计算衰减率τ；
模块M3：使周围网格守恒变量q组成新的守恒变量qex代替原守恒变量q，同时由衰减率τ得到新的人工黏性模型，并计算出黏性系数μ；
模块M4：通过得到的守恒变量qex和黏性系数μ，得到粘性通量g，基于间断有限元的理论求出随时间t变化的守恒变量q，进而得到整个流场的结果。
6.根据权利要求5所述的基于模态衰减的人工黏性系统，其特征在于，所述模块M1包括如下模块：
模块M1.1：求解一般的双曲方程：
Figure FDA0003507441810000033
其中t为时间变量，q为守恒变量，f为对流通量，g为粘性通量；
模块M1.2：粘性通量g表示为：
Figure FDA0003507441810000034
其中μ为人工黏性系数；
模块M1.3：将守恒变量q分解为时间项和空间项：
Figure FDA0003507441810000035
其中Np和r分别表示单元网格中节点个数和坐标向量，m为单元网格中节点的序号，
Figure FDA0003507441810000036
模态分量函数，
Figure FDA0003507441810000037
为一般正交基函数；
模块M1.4：根据模态衰减的间断指示器，构造模态分量函数为：
Figure FDA0003507441810000038
其中C为已知参数，τ为衰减率。
7.根据权利要求5所述的基于模态衰减的人工黏性系统，其特征在于，所述模块M2包括如下模块：
模块M2.1：构造一个最优的模态衰减系数bm：
Figure FDA0003507441810000041
模块M2.2：将原模态分量函数
Figure FDA0003507441810000042
与最优模态衰减系数bm组合，构造新的模态系数为
Figure FDA0003507441810000043
Figure FDA0003507441810000044
其中，
Figure FDA0003507441810000045
表示为在网格单元K中取守恒变量q的标准二范数L2；
模块M2.3：通过如下公式，模态系数
Figure FDA0003507441810000046
依次单调：
Figure FDA0003507441810000047
其中，n表示为网格节点的一组序号；
模块M2.4：使用最小二乘的系统计算衰减率τ，目标函数选取：
Figure FDA0003507441810000048
8.根据权利要求5所述的基于模态衰减的人工黏性系统，其特征在于，所述模块M3包括如下模块：
模块M3.1：将该网格单元守恒变量q和左右两个网格单元的守恒变量qL、qR构造一个新的守恒变量qex：
qex(r,t)＝γLqL(r,t)+γq(r,t)+γRqR(r,t) (9)
其中，γL，γ，γR表示网格单元守恒变量对应的加权系数；
模块M3.2：构造的守恒变量值光滑变化，加权系数γL，γ，γR选择如下所示：
Figure FDA0003507441810000049
其中参数r*为：
Figure FDA0003507441810000051
其中，P为计算的精度，P*和Pmax表示计算中所需的精度参数；
模块M3.3：在计算时如果耗散过度，会破坏高阶格式的精度，计算黏性系数μ的系统为：
Figure FDA0003507441810000052
其中μmax为最大的黏性系数，h为网格尺度，P为计算的精度。
CN202210143246.9A 2022-02-16 2022-02-16 基于模态衰减的人工黏性方法及系统 Active CN114611421B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210143246.9A CN114611421B (zh) 2022-02-16 2022-02-16 基于模态衰减的人工黏性方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210143246.9A CN114611421B (zh) 2022-02-16 2022-02-16 基于模态衰减的人工黏性方法及系统

Publications (2)

Publication Number Publication Date
CN114611421A true CN114611421A (zh) 2022-06-10
CN114611421B CN114611421B (zh) 2023-07-07

Family

ID=81859233

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210143246.9A Active CN114611421B (zh) 2022-02-16 2022-02-16 基于模态衰减的人工黏性方法及系统

Country Status (1)

Country Link
CN (1) CN114611421B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120065950A1 (en) * 2010-09-09 2012-03-15 Ming Lu Numerical method for simulating subsonic flows based on euler equations in lagrangian formulation
CN108153984A (zh) * 2017-12-27 2018-06-12 中国空气动力研究与发展中心计算空气动力研究所 一种基于流场密度阶跃的高精度间断迦辽金人工粘性激波捕捉方法
CN108197072A (zh) * 2017-12-27 2018-06-22 中国空气动力研究与发展中心计算空气动力研究所 一种基于加权守恒变量阶跃的高精度间断Galerkin人工粘性激波捕捉方法
CN108197367A (zh) * 2017-12-27 2018-06-22 中国空气动力研究与发展中心计算空气动力研究所 一种基于流场通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120065950A1 (en) * 2010-09-09 2012-03-15 Ming Lu Numerical method for simulating subsonic flows based on euler equations in lagrangian formulation
CN108153984A (zh) * 2017-12-27 2018-06-12 中国空气动力研究与发展中心计算空气动力研究所 一种基于流场密度阶跃的高精度间断迦辽金人工粘性激波捕捉方法
CN108197072A (zh) * 2017-12-27 2018-06-22 中国空气动力研究与发展中心计算空气动力研究所 一种基于加权守恒变量阶跃的高精度间断Galerkin人工粘性激波捕捉方法
CN108197367A (zh) * 2017-12-27 2018-06-22 中国空气动力研究与发展中心计算空气动力研究所 一种基于流场通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法

Also Published As

Publication number Publication date
CN114611421B (zh) 2023-07-07

Similar Documents

Publication Publication Date Title
US8775140B2 (en) Time and space scaled S-model for turbulent fluid flow simulations
CN108153984B (zh) 一种基于流场密度阶跃的高精度间断迦辽金人工粘性激波捕捉方法
CN111859529B (zh) 一种飞行器绕流数值模拟的多重网格扰动域更新加速方法
CN113723615A (zh) 基于超参优化的深度强化学习模型的训练方法、装置
Richmond et al. Evaluation of an offshore wind farm computational fluid dynamics model against operational site data
CN106650086A (zh) 基于fluent软件的间接空冷系统数值模拟平台
CN105678002A (zh) 等离子体粒子-场自洽系统长期大规模高保真模拟方法
Chronis et al. Generative fluid dynamics: integration of fast fluid dynamics and genetic algorithms for wind loading optimization of a free form surface
CN113864112B (zh) 风力发电机组的尾流流场的确定方法、装置及系统
CN115438598A (zh) 基于一般时间根方尺度的雷诺应力湍流模型的数值方法
CN117217062B (zh) 基于刚度矩阵的流体仿真方法及装置
CN115879335A (zh) 一种基于图生成神经网络的流体多物理场参数预测方法
Quagliarella et al. An open-source aerodynamic framework for benchmarking multi-fidelity methods
CN111159956B (zh) 一种基于特征的流场间断捕捉方法
CN114611421B (zh) 基于模态衰减的人工黏性方法及系统
Xie et al. A new solution to the spherical particle surface concentration of lithium-ion battery electrodes
CN109960841B (zh) 一种流体表面张力的仿真方法、终端设备及存储介质
Sonar Classical finite volume methods
CN110909511B (zh) 一种无曲面体积分的无粘低速绕流数值模拟方法
CN114239437A (zh) 一种精细化瞬态中子输运计算的时间步加速方法
CN104679953B (zh) 户内变电站主变室流场温度场数值模拟快速计算方法、系统
CN112989739B (zh) 一种Trap-Gear时间离散格式的时间步长设定方法
CN117172161B (zh) 一种流场模拟方法、装置、计算机设备及存储介质
Safta et al. Toward using surrogates to accelerate solution of stochastic electricity grid operations problems
CN116306038B (zh) 一种气象模式和同化系统的垂直层匹配方法和装置

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
CB02 Change of applicant information

Address after: No. 1333-1 Zhongchun Road, Minhang District, Shanghai, 201109

Applicant after: SHANGHAI INSTITUTE OF ELECTROMECHANICAL ENGINEERING

Address before: No. 3888, Yuanjiang Road, Minhang District, Shanghai, 201100

Applicant before: SHANGHAI INSTITUTE OF ELECTROMECHANICAL ENGINEERING

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant