CN112100835A - 一种适用于复杂流动的高效高精度数值模拟方法 - Google Patents
一种适用于复杂流动的高效高精度数值模拟方法 Download PDFInfo
- Publication number
- CN112100835A CN112100835A CN202010925232.3A CN202010925232A CN112100835A CN 112100835 A CN112100835 A CN 112100835A CN 202010925232 A CN202010925232 A CN 202010925232A CN 112100835 A CN112100835 A CN 112100835A
- Authority
- CN
- China
- Prior art keywords
- polynomial
- interface
- template
- reconstruction
- discontinuity
- 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
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
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种适用于复杂流动的高精度数值模拟方法,通过采用二维空间模板插值的方式,完成高阶重构多项式的构造,解决了多维黎曼求解器中所需的重构变量无法由传统适用于结构化网格的高阶格式直接求解的弊端,提高波系结构的分辨率以及计算稳定CFL数;并优选通过采用间断探测技术,有效提高了程序的求解效率。本发明能够在解的光滑区域保持一致的时空高阶精度,基本无震荡地完成对流场间断的捕捉并保证流场解的多维特性保持良好。
Description
技术领域
本发明公开了一种适用于复杂流动的高效高精度数值模拟方法,涉及计算流体力学领域。
背景技术
在航空宇航领域,计算流体力学(CFD)已经成为了一种不可或缺的技术手段,为模拟真实流动现象和降低研究成本提供了有力的支持。相比于传统的风洞实验,CFD技术具有更低的计算成本和更高的飞行大气环境仿真度,从而以较小的计算代价获得较为精确的气动特性,为飞行器设计提供快速、准确的指导。随着飞行器功能需求和任务剖面的日益复杂多样,低耗散性、高分辨率特性以及高精度特性在流动模拟中的需求日益增加,人们对CFD也提出了更高的要求。
CFD最为关键的技术有三种:网格技术,数值离散方法和物理模型。这三个方面中任何一项的进步都能有效的推动CFD的发展以及在工程领域的应用。其中,通量格式和重构格式均属于空间离散格式,对CFD的求解精度和计算效率影响很大。在目前广泛应用于航空航天飞行器设计的较为成熟的CFD软件中,通量格式大多采用一维黎曼求解器,重构格式大多采用二阶格式。其中一维黎曼求解器在计算网格界面的数值通量时只考虑界面法向波系的传播而忽略横向波的影响,无法描述横向波传递至区域边界的流动特征,从而导致多维流动计算中波系结构分辨率降低以及计算允许的CFL数减小。而二阶重构格式虽说基本能够满足对于全机升阻力分析,气动外形设计优化方面的精度需求,但是对于湍流,分离等多尺度流动现象,仍难以给出满意的结果。
Balsara通过推导界面角点处多维黎曼问题的求解公式,提出了一种基于角点框架模式的真正多维黎曼求解器。该求解器通过计算单元界面角点处的数值通量来描述横向波传递的流动特征,从而体现流动的多维效应,提高了波系分辨率以及计算允许的CFL数。该方法具备简单的封闭形式,计算效果良好且算法易于实现。但是,相比于一维黎曼求解器,多维黎曼求解器存在单步求解耗时长、效率低的问题。此外,目前针对多维黎曼求解器高阶格式的研究很少,由于角点处的数值通量求解需要用到角点四周的重构物理量值,传统的适用于结构网格的高阶重构方法无法直接对角点四周的物理量进行重构,从而成为多维求解器向高阶格式推广的主要困难。
发明内容
为进一步提高针对复杂流动的数值模拟精度,本发明提供一种适用于复杂流动的高精度数值模拟方法,以期为更加精准的飞行器设计(例如高性能翼型设计)工作提供一定的技术支撑。本发明主要针对传统CFD求解器计算精度低,稳定CFL数小的弊端,采用多维黎曼求解器以及与之对应的高阶重构方案对流动控制方程进行离散并模拟;同时,针对多维黎曼求解器单步求解耗时长、效率低的问题,采用间断探测技术以提高重构效率、减少计算耗时,最终完成多维复杂流动的高效高精度数值模拟。
为实现上述目的,本发明采用的技术方案包括以下步骤:
步骤1:根据设计任务要求,构建设计对象模型,并建立计算网格,获得需要的网格单元信息;
步骤2:在步骤1建立的网格空间中构造空间离散形式的半离散控制方程:
当采用欧拉方程时,其微分形式的控制方程为
其中,
q是守恒形式的流场变量,f和g表示x和y方向的通量,ρ,u,v,p,E分别表示流体密度、x方向速度,y方向速度,压强和能量;
对通量项进行空间离散得到:
其中i,j是单元节点编号;fi+1/2,j和gi,j+1/2分别为x方向和y方向的界面数值通量;
步骤3、采用HWENO重构方案构造重构多项式
首先对控制方程分别沿x和y进行求导,求导后的控制方程转化为:
选择含有V0~V8在内共9个单元的模板,并将其分为8个子模板:
S1={V0,V1,V2,V8},S2={V0,V2,V3,V4},S3={V0,V4,V5,V6},S4={V0,V6,V7,V8}
S5={V0,V1,V2,V3,V7,V8},S6={V0,V1,V2,V3,V4,V5},S7={V0,V3,V4,V5,V6,V7},S8={V0,V1,V5,V6,V7,V8}
其中V0~V8分别指网格单元Vi,j、Vi+1,j+1、Vi,j+1、Vi1,j+1、Vi-1,j、Vi1,j-1、Vi,j-1、Vi+1,jj、Vi+1,j;
在子模板S1,S2,S3,S4中,关于任意变量q的插值多项式pn需要满足如下约束条件:
其中,
n=1,k=0,1,2,8,kx=8,ky=2;n=2,k=0,2,3,4,kx=4,ky=2;
n=3,k=0,4,5,6,kx=4,ky=6;n=4,k=0,6,7,8,kx=8,ky=6.
在子模板S5,S6,S7,S8中,关于任意变量q的插值多项式pn需要满足如下约束条件:
其中,
n=5,k=0,1,2,3,7,8;n=6,k=0,1,2,3,4,5;
n=7,k=0,3,4,5,6,7;n=8,k=0,1,5,6,7,8.
在每个子模板中,具有三阶精度的插值多项式写为:
pn(x,y)=a0+a1(x-x0)+a2(y-y0)+a3(x-x0)(y-y0)
+a4(x-x0)2+a5(y-y0)2,n=1,2,3.4,5,6,7,8
将插值多项式代入约束条件,从而在每个子模板上得到一组关于多项式系数ak(k=0,1,2,3,4,5)的线性代数方程组,求解该方程组得到各子模板中插值多项式pn的系数ak;
在得到各子模板上的插值多项式后,采用WENO限制器的方法,通过光滑指示因子求得9个多项式的权重,加权组合成最终的重构多项式Pi,j(x,y):
光滑指示因子定义如下:
其中|α|=α1+α2,根据光滑指示因子得到每个多项式的权重如下:
最后通过加权组合的方式得到单元Vi,j上的空间重构多项式:
对于单元Vi,j,首先通过最终的重构多项式求得变量在包含界面中点和角点在内的8个插值点(xi+1,j,yi+1/2,j),(xi,j,yi+1,j),(xi-1,j,yi,j),(xi,j,yi-1,j),(xi+1,j,yi+1/2,j),(xi-1/2,j,yi+1/2,j),(xi-1/2,j,yi-1/2,j),(xi+1/2,j,yi-1/2,j)处的重构值,继而求得界面i+1/2处中点以及上下角点处的多维黎曼求解器所需的状态量 以及其中上标“R”和“L”分别表示界面中点两侧的重构变量值,上标“RU”,“LU”,“LD”和“RD”表示角点四周的重构变量值,下标“i+1/2,j”表示界面中点,下标“i+1/2,j+1/2”和“i+1/2,j-1/2”分别表示界面的上下角点;
步骤5、采用多维黎曼求解器进行界面通量求解:
在步骤2的半离散控制方程中,界面通量的具体求解公式如下:
a是声速,上标“~”表示Roe平均;
步骤6、根据界面通量求解残差,并将半离散有限体积格式转化为时空全离散有限体积格式,全流场进行时间推进求解,得到最终的流场解。
进一步的,步骤3中采用间断探测器判断模板内是否含有间断,在含有间断的模板中采用WENO方法,对各子模板构造得到的插值多项式进行加权平均,得到最终的重构多项式,而在不含间断的模板中,则只求解任一子模板上的插值多项式作为最终的重构多项式。
进一步的,所述间断探测器通过间断探测因子βk判断包含单元Vk及其8个邻域单元在内的模板中是否存在间断;间断探测因子βk定义为
利用以下判断式可以判断单元附近是否存在间断:
其中β0为设定的判断阈值。
进一步的,选取β0=5。
其中上表“c”表示与界面角点相关的物理量,下标“RU”,“LU”,“LD”和“RD”表示角点四周的重构变量值;其中的波速计算采用如下公式:
有益效果
本发明的有益效果是,提供一种适用于复杂流动的高精度数值模拟方法,通过采用二维空间模板插值的方式,完成高阶重构多项式的构造,解决了多维黎曼求解器中所需的重构变量无法由传统适用于结构化网格的高阶格式直接求解的弊端,提高波系结构的分辨率以及计算稳定CFL数;并优选通过采用间断探测技术,有效提高了程序的求解效率。本发明能够在解的光滑区域保持一致的时空高阶精度,基本无震荡地完成对流场间断的捕捉并保证流场解的多维特性保持良好。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明的实现流程图。
图2是HWENO方法的模板示意图。
图3是多维黎曼求解器通量所需的插值点位置示意图。
图4是多维黎曼求解器通量求解的示意图。
图5是实施例中,RAE2822翼型的计算网格示意图。
图6是实施例中,采用本发明方案求得的翼型压力系数等值线图。
图7是实施例中,采用本发明方案求得的翼型表面压力分布和风洞试验结果的对比。
具体实施方式
本发明目的是针对多维黎曼求解器单步求解耗时长、效率低的问题,采用间断探测技术以提高重构效率、减少计算耗时,最终完成对例如高性能翼型设计等应用需求下的多维复杂流动的高效高精度数值模拟。
主要包括以下步骤:
步骤1、根据设计任务要求,构建设计对象模型,如本实施例中RAE2822翼型,模型,并建立复杂流场计算所需的网格,由给定网格得到需要的网格单元信息,如网格尺度,节点坐标等。
步骤2、构造空间离散形式的半离散控制方程。
在步骤1提供的网格空间上构造空间离散形式的半离散控制方程。以欧拉方程为例,其微分形式的控制方程如下所示:
其中,
q是守恒形式的流场变量,f和g表示x和y方向的通量,ρ,u,v,p,E分别表示流体密度、x方向速度,y方向速度,压强和能量。
对通量项进行空间离散可以得到:
其中i,j是单元节点编号;fi+1/2,j和gi,j+1/2分别为x方向和y方向的界面数值通量,由高阶重构格式和多维黎曼求解器求解得到,具体求解过程详见下述步骤。
步骤3、构造重构多项式
重构多项式的构造是高阶重构方案的关键步骤之一。通过构造空间上的二维插值多项式可以获得物理量在网格单元上的分布,从而求得所需位置处的重构变量值。本发明采用HWENO重构方案进行空间二维插值多项式的构造。
为了获得步骤2中界面处的数值通量fi+1/2,j和gi,j+1/2,需要首先获得守恒流场变量在网格单元上的空间分布,进而用于下述步骤5,6中重构变量和界面通量的求解。本发明采用HWENO重构方案进行空间二维插值多项式的构造。
HWENO重构方案采用Hermite插值多项式来构造单元Vi,j上具备三阶精度的空间插值多项式。该方法既需要函数值,也需要函数值的导数,从而加强了单元之间的联系,使得格式更为紧凑和稳定。详细构造过程如下:
首先对控制方程分别沿x和y进行求导,求导后的控制方程转化为:
如图2所示,选择含有V0~V8在内共9个单元的模板,并将其分为8个子模板:
S1={V0,V1,V2,V8},S2={V0,V2,V3,V4},S3={V0,V4,V5,V6},S4={V0,V6,V7,V8}
S5={V0,V1,V2,V3,V7,V8},S6={V0,V1,V2,V3,V4,V5},S7={V0,V3,V4,V5,V6,V7},S8={V0,V1,V5,V6,V7,V8}
其中V0~V8分别指网格单元Vi,j、Vi+1,j+1、Vi,j+1、Vi-1,j+1、Vi-1,j、Vi-1,j-1、Vi,j-1、Vi+1,j-j、Vi+1,j。
在子模板S1,S2,S3,S4中,关于任意变量q的插值多项式pn需要满足如下约束条件:
其中,
n=1,k=0,1,2,8,kx=8,ky=2;n=2,k=0,2,3,4,kx=4,ky=2;
n=3,k=0,4,5,6,kx=4,ky=6;n=4,k=0,6,7,8,kx=8,ky=6.
在子模板S5,S6,S7,S8中,关于任意变量q的插值多项式pn需要满足如下约束条件:
其中,
n=5,k=0,1,2,3,7,8;n=6,k=0,1,2,3,4,5;
n=7,k=0,3,4,5,6,7;n=8,k=0,1,5,6,7,8.
在每个子模板中,具有三阶精度的插值多项式可以写为:
将式(7)代入式(5)和(6),每个子模板上都可得到一组关于多项式系数ak(k=0,1,2,3,4,5)的线性代数方程组,求解该方程组便可得各子模板中插值多项式pn的系数ak。
在得到各子模板上的插值多项式后,类似TWENO重构方案的方法,采用WENO限制器的方法,通过光滑指示因子求得9个多项式的权重,加权组合成最终的重构多项式Pi,j(x,y)。下面采用WENO限制器,根据这五个多项式的光滑程度来赋予其不同的权重。多项式所在的区域越光滑,其所占权重越大,通过该方式能够抑制流场间断附近出现非物理震荡。光滑指示因子定义如下:
其中|α|=α1+α2,根据光滑指示因子得到每个多项式的权重如下:
最后通过加权组合的方式得到单元Vi,j上的空间重构多项式:
步骤4:间断探测器的构造。
步骤3中的WENO格式通过自适应选取候选模板以降低间断所在区域对重构多项式的影响,从而避免流场中出现震荡。但在光滑流场区域,各候选模板的权重相近,继续采用WENO格式会造成不必要的计算消耗。因此,本实施例通过引入间断探测器来进一步重构效率。该方法的核心思想是,通过间断探测技术判断模板内是否含有间断,在含有间断的模板中采用WENO方法,对各子模板构造得到的插值多项式进行加权平均,得到最终的重构多项式P(x,y),而在不含间断的模板中,只需求解任一子模板上的插值多项式(如p0(x,y))作为最终的重构多项式即可,无需再对其它子模板进行多项式构造,从而大幅降低计算耗时。
下面提出间断探测因子βk定义,通过该因子判断包含单元Vk及其8个邻域单元在内的模板中是否存在间断:
显然,利用以下判断式可以判断单元附近是否存在间断:
随着网格尺度h→0,在光滑区域βk→0,而在间断区域βk→∞,因此采用如下判据来监测间断位置:
在本实施例中,选取β0=5。
步骤5、根据重构多项式求解多维黎曼求解器所需的重构状态量;
下面利用步骤3中构造的空间重构多项式对多维黎曼求解器所需的输入重构变量值进行求解,包括如下两个步骤:
步骤5.2如图3所示,令m表示单元内的包含界面中点和角点在内的8个插值点编号,以单元Vi,j为例,m=1,2,…,8分别表示点(xi+1,j,yi+1/2,j),(xi,j,yi+1,j),(xi-1,j,yi,j),(xi,j,yi-1,j),(xi+1,j,yi+1/2,j),(xi-1/2,j,yi+1/2,j),(xi-1/2,j,yi-1/2,j),(xi+1/2,j,yi-1/2,j)。通过可以求得变量在这些插值点处的重构值。
通过求解得到的插值点处的重构值,求得界面i+1/2处中点以及上下角点处的多维黎曼求解器所需的状态量以及其中上标“R”和“L”分别表示界面中点两侧的重构变量值,上标“RU”,“LU”,“LD”和“RD”表示角点四周的重构变量值,下标“i+1/2,j”表示界面中点,下标“i+1/2,j+1/2”和“i+1/2,j-1/2”分别表示界面的上下角点。
例如单元界面两侧的重构变量值为
单元角点C1四周的重构变量值为:
步骤6、采用多维黎曼求解器进行界面通量求解
图4给出了多维黎曼求解器计算通量所需要的初态,除了界面两侧的两个初态“R”和“L”以外,单元的角点处(点C1~C4)还需要四个初态“RU”,“LU”,“LD”和“RD”,这些初态值步骤5得到。下面给出采用多维黎曼求解器进行界面通量求解的详细步骤。
在半离散控制方程(3)中,界面通量的具体求解公式如下:
a是声速,上标“~”表示Roe平均。
其中上表“c”表示与界面角点相关的物理量,下标“RU”,“LU”,“LD”和“RD”表示角点四周的重构变量值,传统适用于结构化网格的高阶重构方法无法求解这些位置的重构变量,但可由本发明步骤3中提出的重构方案求解得到。这里的波速计算采用如下公式:
步骤7、根据界面通量求解残差,并将半离散有限体积格式转化为时空全离散有限体积格式,全流场进行时间推进求解,得到最终的流场解。
步骤6中求得界面通量fi+1/2,j和gi,j+1/2后,由步骤2的式(3)即可求解得到当前第n层时间步的残差L。对时间变量采用三阶Runge-Kutta离散公式将半离散有限体积格式转化为时空全离散有限体积格式:
L表示残差,在得到界面数值通量后,可由式(3)右端项求解得到,如果是NS方程,则额外计算粘性项的影响。上标“n”表示时间步。利用时空全离散有限体积格式求解下一时间步上的流场变量值。重复以上步骤,依次推进,直到得到全流场稳定的数值模拟结果。
本实施例中是以RAE2822翼型跨音速绕流问题进行求解的。图6给出了采用本发明方案求得的压力系数等值线图,能精确捕捉到激波的位置以及激波的强弱。图7给出采用本发明方案求得的翼型表面压力分布和风洞试验结果的对比,可以看出计算结果和试验数据吻合良好。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (5)
1.一种适用于复杂流动的高效高精度数值模拟方法,其特征在于:包括以下步骤:
步骤1:根据设计任务要求,构建设计对象模型,并建立计算网格,获得需要的网格单元信息;
步骤2:在步骤1建立的网格空间中构造空间离散形式的半离散控制方程:
当采用欧拉方程时,其微分形式的控制方程为
其中,
q是守恒形式的流场变量,f和g表示x和y方向的通量,ρ,u,v,p,E分别表示流体密度、x方向速度,y方向速度,压强和能量;
对通量项进行空间离散得到:
其中i,j是单元节点编号;fi+1/2,j和gi,j+1/2分别为x方向和y方向的界面数值通量;
步骤3、采用HWENO重构方案构造重构多项式
首先对控制方程分别沿x和y进行求导,求导后的控制方程转化为:
选择含有V0~V8在内共9个单元的模板,并将其分为8个子模板:
S1={V0,V1,V2,V8},S2={V0,V2,V3,V4},S3={V0,V4,V5,V6},S4={V0,V6,V7V8}
S5={V0,V1,V2,V3,V7,V8},S6={V0,V1,V2,V3,V4,V5},S7={V0,V3,V4,V5,V6,V7},S8={V0,V1,V5,V6,V7,V8}
其中V0~V8分别指网格单元Vi,j、Vi+1,j+1、Vi,j+1、Vi-1,j+1、Vi-1,j、Vi-1,j-1、Vi,j-1、Vi+1,j-1、Vi+1,j;
在子模板S1,S2,S3,S4中,关于任意变量q的插值多项式pn需要满足如下约束条件:
其中,
n=1,k=0,1,2,8,kx=8,ky=2;n=2,k=0,2,3,4,kx=4,ky=2;
n=3,k=0,4,5,6,kx=4,ky=6;n=4,k=0,6,7,8,kx=8,ky=6.
在子模板S5,S6,S7,S8中,关于任意变量q的插值多项式pn需要满足如下约束条件:
其中,
n=5,k=0,1,2,3,7,8;n=6,k=0,1,2,3,4,5;
n=7,k=0,3,4,5,6,7;n=8,k=0,1,5,6,7,8.
在每个子模板中,具有三阶精度的插值多项式写为:
pn(x,y)=a0+a1(x-x0)+a2(y-y0)+a3(x-x0)(y-y0)+a4(x-x0)2+a5(y-y0)2,n=1,2,3,4,5,6,7,8
将插值多项式代入约束条件,从而在每个子模板上得到一组关于多项式系数ak(k=0,1,2,3,4,5)的线性代数方程组,求解该方程组得到各子模板中插值多项式pn的系数ak;
在得到各子模板上的插值多项式后,采用WENO限制器的方法,通过光滑指示因子求得9个多项式的权重,加权组合成最终的重构多项式Pi,j(x,y):
光滑指示因子定义如下:
其中|α|=α1+α2,根据光滑指示因子得到每个多项式的权重如下:
最后通过加权组合的方式得到单元Vi,j上的空间重构多项式:
对于单元Vi,j,首先通过最终的重构多项式求得变量在包含界面中点和角点在内的8个插值点(xi+1,j,yi+1/2,j),(xi,j,yi+1,j),(xi-1,j,yi,j),(xi,j,yi-1,j),(xi+1,j,yi+1/2,j),(xi-1/2,j,yi+1/2,j),(xi-1/2,j,yi-1/2,j),(xi+1/2,j,yi-1/2,j)处的重构值,继而求得界面i+1/2处中点以及上下角点处的多维黎曼求解器所需的状态量 以及其中上标“R”和“L”分别表示界面中点两侧的重构变量值,上标“RU”,“LU”,“LD”和“RD”表示角点四周的重构变量值,下标“i+1/2,j”表示界面中点,下标“i+1/2,j+1/2”和“i+1/2,j-1/2”分别表示界面的上下角点;
步骤5、采用多维黎曼求解器进行界面通量求解:
在步骤2的半离散控制方程中,界面通量的具体求解公式如下:
a是声速,上标“~”表示Roe平均;
步骤6、根据界面通量求解残差,并将半离散有限体积格式转化为时空全离散有限体积格式,全流场进行时间推进求解,得到最终的流场解。
2.根据权利要求1所述一种适用于复杂流动的高效高精度数值模拟方法,其特征在于:步骤3中采用间断探测器判断模板内是否含有间断,在含有间断的模板中采用WENO方法,对各子模板构造得到的插值多项式进行加权平均,得到最终的重构多项式,而在不含间断的模板中,则只求解任一子模板上的插值多项式作为最终的重构多项式。
4.根据权利要求3所述一种适用于复杂流动的高效高精度数值模拟方法,其特征在于:选取β0=5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010925232.3A CN112100835B (zh) | 2020-09-06 | 2020-09-06 | 一种适用于复杂流动的高效高精度翼型绕流数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010925232.3A CN112100835B (zh) | 2020-09-06 | 2020-09-06 | 一种适用于复杂流动的高效高精度翼型绕流数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112100835A true CN112100835A (zh) | 2020-12-18 |
CN112100835B CN112100835B (zh) | 2022-06-14 |
Family
ID=73757800
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010925232.3A Active CN112100835B (zh) | 2020-09-06 | 2020-09-06 | 一种适用于复杂流动的高效高精度翼型绕流数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112100835B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765725A (zh) * | 2020-12-30 | 2021-05-07 | 四川京航天程科技发展有限公司 | 针对多维Euler方程的解析黎曼解算方法 |
CN113408168A (zh) * | 2021-06-18 | 2021-09-17 | 北京理工大学 | 一种基于黎曼问题精确求解的高精度数值模拟方法 |
CN113505443A (zh) * | 2021-09-09 | 2021-10-15 | 南京航空航天大学 | 一种任意外形的三维绕流问题自适应笛卡尔网格生成方法 |
CN114676378A (zh) * | 2021-11-23 | 2022-06-28 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于rad求解器的激波计算方法 |
CN114707254A (zh) * | 2022-06-01 | 2022-07-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于模板构造法的二维边界层网格生成方法及系统 |
CN116611369A (zh) * | 2023-07-17 | 2023-08-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于光滑度量量级与候选模板点数的插值方法及装置 |
CN116702571A (zh) * | 2023-08-07 | 2023-09-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于多重光滑度量因子的数值模拟方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014043846A1 (zh) * | 2012-09-18 | 2014-03-27 | Lu Ming | 求解二维黎曼问题模拟亚音速无粘流的数值方法 |
CN103823916A (zh) * | 2013-10-23 | 2014-05-28 | 沈智军 | 一种基于多维黎曼解的任意拉格朗日欧拉方法 |
CN106682262A (zh) * | 2016-11-21 | 2017-05-17 | 中国航天空气动力技术研究院 | 一种获取飞行器流场的数值模拟方法 |
CN108763683A (zh) * | 2018-05-16 | 2018-11-06 | 南京航空航天大学 | 一种三角函数框架下新weno格式构造方法 |
CN109902376A (zh) * | 2019-02-25 | 2019-06-18 | 北京理工大学 | 一种基于连续介质力学的流固耦合高精度数值模拟方法 |
CN110457806A (zh) * | 2019-08-02 | 2019-11-15 | 南京航空航天大学 | 基于交错网格的中心五阶weno格式的全流场模拟方法 |
-
2020
- 2020-09-06 CN CN202010925232.3A patent/CN112100835B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014043846A1 (zh) * | 2012-09-18 | 2014-03-27 | Lu Ming | 求解二维黎曼问题模拟亚音速无粘流的数值方法 |
CN103823916A (zh) * | 2013-10-23 | 2014-05-28 | 沈智军 | 一种基于多维黎曼解的任意拉格朗日欧拉方法 |
CN106682262A (zh) * | 2016-11-21 | 2017-05-17 | 中国航天空气动力技术研究院 | 一种获取飞行器流场的数值模拟方法 |
CN108763683A (zh) * | 2018-05-16 | 2018-11-06 | 南京航空航天大学 | 一种三角函数框架下新weno格式构造方法 |
CN109902376A (zh) * | 2019-02-25 | 2019-06-18 | 北京理工大学 | 一种基于连续介质力学的流固耦合高精度数值模拟方法 |
CN110457806A (zh) * | 2019-08-02 | 2019-11-15 | 南京航空航天大学 | 基于交错网格的中心五阶weno格式的全流场模拟方法 |
Non-Patent Citations (4)
Title |
---|
FENGQU 等: "A new genuinely two-dimensional Riemann solver for multidimensional Euler and Navier–Stokes Equations", 《COMPUTER PHYSICS COMMUNICATIONS》 * |
FENGQU 等: "A new genuinely two-dimensional Riemann solver for multidimensional Euler and Navier–Stokes Equations", 《COMPUTER PHYSICS COMMUNICATIONS》, 31 October 2019 (2019-10-31), pages 1 - 11 * |
王燕: "微分-差分KdV方程的黎曼v函数周期波解及其渐近性质", 《广东技术师范学院学报》 * |
王燕: "微分-差分KdV方程的黎曼v函数周期波解及其渐近性质", 《广东技术师范学院学报》, no. 05, 15 May 2016 (2016-05-15), pages 14 - 18 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765725A (zh) * | 2020-12-30 | 2021-05-07 | 四川京航天程科技发展有限公司 | 针对多维Euler方程的解析黎曼解算方法 |
CN113408168B (zh) * | 2021-06-18 | 2022-11-08 | 北京理工大学 | 一种基于黎曼问题精确求解的高精度数值模拟方法 |
CN113408168A (zh) * | 2021-06-18 | 2021-09-17 | 北京理工大学 | 一种基于黎曼问题精确求解的高精度数值模拟方法 |
CN113505443A (zh) * | 2021-09-09 | 2021-10-15 | 南京航空航天大学 | 一种任意外形的三维绕流问题自适应笛卡尔网格生成方法 |
CN113505443B (zh) * | 2021-09-09 | 2021-12-14 | 南京航空航天大学 | 一种任意外形的三维绕流问题自适应笛卡尔网格生成方法 |
CN114676378A (zh) * | 2021-11-23 | 2022-06-28 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于rad求解器的激波计算方法 |
CN114676378B (zh) * | 2021-11-23 | 2023-05-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于rad求解器的激波计算方法 |
CN114707254B (zh) * | 2022-06-01 | 2022-08-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于模板构造法的二维边界层网格生成方法及系统 |
CN114707254A (zh) * | 2022-06-01 | 2022-07-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于模板构造法的二维边界层网格生成方法及系统 |
CN116611369A (zh) * | 2023-07-17 | 2023-08-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于光滑度量量级与候选模板点数的插值方法及装置 |
CN116611369B (zh) * | 2023-07-17 | 2023-09-29 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于光滑度量量级与候选模板点数的插值方法及装置 |
CN116702571A (zh) * | 2023-08-07 | 2023-09-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于多重光滑度量因子的数值模拟方法及装置 |
CN116702571B (zh) * | 2023-08-07 | 2023-10-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于多重光滑度量因子的数值模拟方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112100835B (zh) | 2022-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112100835B (zh) | 一种适用于复杂流动的高效高精度翼型绕流数值模拟方法 | |
CN112016167B (zh) | 基于仿真和优化耦合的飞行器气动外形设计方法及系统 | |
Takenaka et al. | Multidisciplinary design exploration for a winglet | |
CN108763683B (zh) | 一种三角函数框架下新weno格式构造方法 | |
CN109726433B (zh) | 基于曲面边界条件的三维无粘低速绕流的数值模拟方法 | |
CN107391891A (zh) | 一种基于模型融合方法的大展弦比机翼优化设计方法 | |
CN113609598B (zh) | 飞行器气动特性模拟的rans/les扰动域更新方法 | |
CN108363843A (zh) | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 | |
CN105678015B (zh) | 一种高超声速三维机翼的非概率可靠性气动结构耦合优化设计方法 | |
CN105718634A (zh) | 一种基于非概率区间分析模型的翼型鲁棒优化设计方法 | |
CN111444661B (zh) | 交互式棱柱网格生成中翘曲现象消除方法 | |
CN109657408A (zh) | 一种再生核粒子算法实现结构线性静力学仿真方法 | |
CN105844025B (zh) | 一种针对高超声速舵面的非概率热气动弹性可靠性设计方法 | |
CN115496006A (zh) | 一种适用于高超声速飞行器的高精度数值模拟方法 | |
Su | Accurate and robust adaptive mesh refinement for aerodynamic simulation with multi‐block structured curvilinear mesh | |
Charvet et al. | Numerical simulation of the flow over sails in real sailing conditions | |
CN109190232B (zh) | 一种飞机平尾区动能损失计算评估方法 | |
CN103530435B (zh) | 一种基于敏感度的船体型线设计方法 | |
Sitaraman et al. | Enhancements to overset methods for improved accuracy and solution convergence | |
CN111625901B (zh) | 一种面向翼型的压力系数曲线智能生成方法 | |
CN108197368A (zh) | 一种适用于飞行器复杂气动外形的几何约束及权函数简捷计算方法 | |
CN114282462B (zh) | 耦合流动及网格分布的重构限制器进行激波捕捉的方法 | |
Sadrehaghighi | Dynamic & Adaptive Meshing | |
Jude et al. | Enhancements to Overset Methods for achieving higher accuracy and solution convergence | |
Biancolini et al. | Validation of structural modeling for realistic wing topologies involved in fsi analyses: ribes test case |
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 |