CN114757070A - 用于数值模拟的三角函数框架下新weno格式构造方法 - Google Patents
用于数值模拟的三角函数框架下新weno格式构造方法 Download PDFInfo
- Publication number
- CN114757070A CN114757070A CN202210373696.7A CN202210373696A CN114757070A CN 114757070 A CN114757070 A CN 114757070A CN 202210373696 A CN202210373696 A CN 202210373696A CN 114757070 A CN114757070 A CN 114757070A
- Authority
- CN
- China
- Prior art keywords
- format
- tweno
- trigonometric function
- value
- flux
- 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.)
- Pending
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/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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
本发明涉及一种用于数值模拟的三角函数框架下新WENO格式构造方法,用三角函数多项式构造的加权基本无振荡格式比用代数多项式构造的经典基本无振荡格式更易于模拟波类或高频振荡类问题,在光滑区域可以获得高阶数值精度,在激波和接触间断处保持基本无振荡的性质;新的TWENO格式虽然和经典五阶WENO格式用了同样五个点的信息,但能得到更低的全局L1和L∞范数截断误差。新的TWENO格式采用的线性权不再需要通过繁冗的数值计算得到最优解,可设为满足和为一的任意正数,相比经典的WENO格式更简便,鲁棒性更强,更易于推广到高维空间。该新型TWENO格式有效地数值模拟了几个经典的欧拉问题,充分验证了有效性。
Description
本申请是名为《一种三角函数框架下新WENO格式构造方法》的专利申请的分案申请,原申请的申请日为2018年05月16日,申请号为201810472192.4。
技术领域
本发明属于计算流体力学工程技术领域,具体涉及一种用于数值模拟的三角函数框架下新WENO格式构造方法。
背景技术
在工程应用中,流场问题常常出现,比如气体动力系统和浅水建模等。因此,制定解决这类问题的鲁棒的、精确的、高效的数值模拟方法至关重要,也吸引了很多研究者的兴趣。1959年,Godunov为解流场问题提出了一阶精度的数值模拟格式。一阶精度的数值模拟方法在捕捉激波时不会出现非物理的数值振荡但会过度抹平强间断,而往往强间断对问题的后续研究有着重要意义,因此需引进高精度数值计算格式模拟强间断类问题。
为了提高格式的精度,模拟解的结构以及准确的捕捉激波位置,Harten于1983年首次提出了TVD(Total Variation Diminishing)格式,并在此基础上与Osher于1987年提出了ENO(Essentially Non-Oscillatory)高精度格式。ENO格式的主要思想是在逐次扩展的模板中选用最光滑的模板构造多项式求出单元边界处的值,进而在光滑区域达到高阶精度、高分辨率,同时在间断附近实现基本无振荡的效果。但是,在方法的实现过程中,ENO格式最终只选用所有候选模板中最优的模板,造成计算结果的浪费,且构造的精度越高浪费的越多,导致计算效率不高。因此,Liu,Osher和Chan等于1994年提出了WENO(WeightedEssentially Non-oscillatory)格式,提高了计算结果的利用率并且使得r阶精度的ENO格式提高到r+1阶精度。1996年,Jiang和Shu进一步改善了WENO格式,使得数值精度能够提高到2r-1阶,并设计出新光滑因子和非线性权的构造框架。WENO格式的主要思想是通过低阶重构通量的线性凸组合获得高阶近似。但该经典WENO格式的实现过程中,线性权依赖于母模板,且其求解过程相当复杂,因此,2016年Zhu和Qiu改善了该WENO格式,在维持精度不减的情况下,随机选取大于零且总和为一的线性权。这些格式已被成功地用到很多应用领域,特别是包含激波和复杂解结构的问题,比如模拟可压缩湍流系统和空气声学系统等。
波类和高频振荡类问题在工程应用中常常出现。然而,对模拟此类问题的更适合的三角函数多项式插值WENO格式的研究较少。虽然Baron于1976年研究三角函数插值呈现了Neville类方法,Muhlbach提出了牛顿三角函数插值,但是这些成果不能直接应用于ENO类插值格式。为此,Christofi于1996年提出了能直接用于ENO格式中的三角函数重构方法。Zhu和Qiu于2010年提出了用三角函数多项式重构WENO格式的方法,但计算复杂,不易实现。
发明内容
本发明针对现有技术中的不足,提供一种用于数值模拟的三角函数框架下新WENO格式构造方法,能针对各种可压流场问题,进行高精度数值模拟。
为实现上述目的,本发明提供了如下方案:
一种用于数值模拟的三角函数框架下新WENO格式构造方法,在笛卡尔坐标系下,利用TWENO格式对可压流场问题进行数值模拟,包括以下步骤:
步骤一、把双曲守恒律方程离散为空间半离散的有限差分格式,采用TWENO格式重构通量的近似值;
步骤二、对控制方程中的时间导数使用三阶TVD Runge-Kutta离散公式将空间半离散的有限差分格式离散成时空全离散有限差分格式;
步骤三、根据时空全离散有限差分格式得到下一时间层上的近似值,依次迭代,得到终止时刻计算区域内流场的数值模拟值;
其中,采用TWENO格式重构通量的近似值包括:
其中,L(U)表示-f(U)x的空间离散形式,U=(ρ,ρu,E)T表示守恒变量,f(U)=(ρu,ρu2+p,u(E+p))T表示通量,f(U)x表示f(U)对x求导,x表示空间变量,ρ、u、p、E分别表示流体密度、速度、压强、能量,T表示转置,Ui(t)表示U在目标网格单元Ii内点xi处的值U(xi,t),和分别表示通量f(U)在目标网格单元Ii的边界和处的五阶近似的数值通量;
步骤2、将目标网格单元Ii以及其周围共五个网格单元组成一个大模板T1=[Ii-2,Ii-1,Ii,Ii+1,Ii+2],从大模板中选择两个包含两个网格单元的小模板T2=[Ii-1,Ii]和T3=[Ii,Ii+1];
步骤3、在T1、T2、T3每个模板上分别重构三角函数多项式p1(x)、p2(x)和p3(x),使得:
p2(x),p3(x)∈span{1,sin(x-xi)};
步骤4、任意取三组线性权,每组线性权的和均为1;
步骤5、根据三角函数多项式计算光滑指示器,用于衡量三角函数多项式pl(x)在目标网格单元上的光滑度;l=1,2,3表示对应模板序号;
步骤6、通过线性权和光滑指示器计算非线性权;
将计算结果代入含有时间导数项的空间半离散的有限差分格式,得到关于时间导数的常微分方程。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
相比于WENO格式,该TWENO格式通过把三角函数多项式而不是代数多项式作为有限差分TWENO格式的构建模块,模拟了波类和高频振荡类可压流场问题,同时在光滑区域能够达到高阶精度;相比于已有三角函数多项式重构格式,该TWENO格式得到的全局L1截断误差与L∞截断误差更小,同时也避免了在强激波和接触间断处产生非物理振荡,该新五阶TWENO格式中的相关线性权不再需要通过复杂的计算得到而是被设为和为一的任意正数,因此该新TWENO格式具有更简便更易拓展到高维空间的优势。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1a-1c是实施例一中的台阶问题,利用本发明的有限差分TWENO格式得到的密度等值线图,分别采用线性权①、②、③。
图2a-2c是实施例二中的双马赫问题,利用本发明的有限差分TWENO格式得到的密度等值线图,分别采用线性权①、②、③。
图3a-3c是实施例三中的激波和涡流相互干扰问题,利用本发明的有限差分TWENO格式得到的t=0.35时的压强等值线图,分别采用线性权①、②、③。
图4a-4c是实施例三中的激波和涡流相互干扰问题,利用本发明的有限差分TWENO格式得到的t=0.6时的压强等值线图,分别采用线性权①、②、③。
图5a-5c是实施例三中的激波和涡流相互干扰问题,利用本发明的有限差分TWENO格式得到的t=0.8时的压强等值线图,分别采用线性权①、②、③。
图6a-6c是实施例四中初值条件为(15)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.25时的密度等值线图,分别采用线性权①、②、③。
图7a-7c是实施例四中初值条件为(16)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.25时的密度等值线图,分别采用线性权①、②、③。
图8a-8c是实施例四中初值条件为(17)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.3时的密度等值线图,分别采用线性权①、②、③。
图9a-9c是实施例四中初值条件为(18)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.2时的密度等值线图,分别采用线性权①、②、③。
图10a-10c是实施例四中初值条件为(19)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.3时的密度等值线图,分别采用线性权①、②、③。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种用于数值模拟的三角函数框架下新WENO格式构造方法,能针对各种可压流场问题,进行高精度数值模拟。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明给出了笛卡尔网格下解可压流场问题的新型五阶有限差分TWENO高精度数值计算格式的构造过程,相比于经典WENO格式,该TWENO格式通过把重构的三角函数多项式而不是重构的代数多项式作为有限差分WENO格式的构建模块,解决了波类和高频振荡类的可压流场问题的数值模拟,且在光滑区域能够达到高阶精度近似,捕捉到尖锐和无振荡激波的转换。新的TWENO格式采用的线性权不再需要通过繁冗的数值计算得到,可设为满足和为一的任意正数,此格式方法简单、精度高,易于推广到多维空间中。该方法在笛卡尔坐标系下,利用TWENO格式对可压流场问题进行数值模拟,具体步骤如下:
步骤一、把双曲守恒律方程离散为空间半离散的有限差分格式,采用TWENO格式重构通量的近似值。
考虑一维双曲守恒律方程:
其半离散格式的形式为:
其中,U=(ρ,ρu,E)T表示守恒变量,f(U)=(ρu,ρu2+p,u(E+p))T表示通量,Ut表示U对t求导,f(U)x表示f(U)对x求导,t表示时间变量,x表示空间变量,ρ、u、p、E分别表示流体密度、速度、压强、能量,T表示转置,U0表示初始状态值,L(U)表示-f(U)x的空间离散形式。
步骤2、将目标网格单元Ii以及其周围共5个网格单元组成一个大模板T1=[Ii-2,Ii-1,Ii,Ii+1,Ii+2],从大模板中选择两个包含两个网格单元的小模板T2=[Ii-1,Ii]和T3=[Ii,Ii+1],其中Ii为对应序号的网格单元。
步骤3、在每个模板上分别重构三角函数多项式p1(x)、p2(x)和p3(x),使得:
p2(x),p3(x)∈span{1,sin(x-xi)}。
其具体过程如下:
步骤3.1、在三个模板T1、T2和T3上分别构造三角函数多项式p1(x)、p2(x)和p3(x),使其满足:
步骤3.2、得到每个模板上的三角函数插值多项式p1(x)、p2(x)和p3(x),如下:
步骤4、任意取三组线性权:
①γ1=0.98,γ2=0.01,γ3=0.01;
②γ1=1/3,γ2=1/3,γ3=1/3;
③γ1=0.01,γ2=0.495,γ3=0.495。
步骤5、计算光滑指示器βl,用于衡量重构多项式pl(x)在目标网格单元上的光滑度,计算公式为:
步骤6、通过线性权γl和光滑指示器βl计算非线性权ωl,其计算公式为:
其次,将计算结果代入含有时间导数项的半离散有限差分格式,得到关于时间导数的常微分方程。
步骤二、对控制方程中的时间导数使用三阶TVD Runge-Kutta离散公式将半离散有限差分格式离散成时空全离散有限差分格式。
利用三阶TVD Runge-Kutta离散公式:
得到时空全离散有限差分格式,其中,U(1),U(2)为中间过渡值,Δt为时间步长,上标n表示第n时间层,L(Un),L(U(1)),L(U(2))为-f(U)x的高阶空间离散形式的近似值。
步骤三、根据时空全离散有限差分格式得到下一时间层上的近似值,依次迭代,得到终止时刻计算区域内流场的数值模拟值。
时空全离散有限差分格式为关于时间层的迭代公式,初始状态值已知,通过迭代公式求出下一时间层的近似值,依次得到终止时刻计算区域内流场的数值模拟值。对于二维问题,逐维用上面的重构过程。
下面给出几个算例作为本发明所公开方法的具体实施例。
实施例一、台阶问题。该问题是Emery于1968年提出的一个用于检验非线性双曲型守恒律格式的经典算例。初始数据为水平来流马赫数为3,密度为1.4,水平速度为3,竖直速度为0,压强为1,管道区域为[0,3]×[0,1],在距离左边界0.6处有一高度为0.2的台阶,且台阶延伸到管道的尽头。上下边界为反射边界,左边界为来流边界,右边界为出流边界。图1a-1c给出了t=4时的密度等值线图。
实施例二、双马赫反射问题。该问题描述了一个与x轴成60°角的强激波射到反射墙上发生的变化,来流是马赫数为10的强激波。计算区域为[0,4]×[0,1],区域底部从y=0处开始为反射边界条件,其它的底部边界(从x=0到那部分)为波前条件。图2a-2c给出了t=0.2时在[0,3]×[0,1]区域的密度等值线图。
实施例三、激波和涡流相互干扰问题。马赫数为1.1的激波位于x=0.5处且垂直于x轴,激波初始状态是小涡流位于该激波的左边且其中心位于(xc,yc)=(0.25,0.5)处,涡流可看成平均流的速度,温度和熵的扰动,表示为:
其中,ε=0.3,rc=0.05,α=0.204,γ=1.4,计算区域为[0,2]×[0,1]。图3a-3c给出了t=0.35时在[0,1]×[0,1]区域的压强等值线图。图4a-4c给出了t=0.6时在[0.4,1.45]×[0,1]区域的压强等值线图。图5a-5c给出了t=0.8时在[0,2]×[0,1]区域的压强等值线图。
实施例四、二维Euler黎曼问题。计算区域为[0,1]×[0,1],初值条件分别设为:
图6a-6c给出了初值条件为(15)时二维Euler黎曼问题在t=0.25时刻的密度等值线图。图7a-7c给出了初值条件为(16)时二维Euler黎曼问题在t=0.25时刻的密度等值线图。图8a-8c给出了初值条件为(17)时二维Euler黎曼问题在t=0.3时刻的密度等值线图。图9a-9c给出了初值条件为(18)时二维Euler黎曼问题在t=0.2时刻的密度等值线图。图10a-10c给出了初值条件为(19)时二维Euler黎曼问题在t=0.3时刻的密度等值线图。从图中可以看出基于三角函数多项式空间的有限差分TWENO格式对本发明任意取得的线性权都可以很好的捕捉到黎曼问题的大部分流动特性。
本说明书中每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种用于数值模拟的三角函数框架下新WENO格式构造方法,在笛卡尔坐标系下,利用TWENO格式对可压流场问题进行数值模拟,其特征在于,包括以下步骤:
步骤一、把双曲守恒律方程离散为空间半离散的有限差分格式,采用TWENO格式重构通量的近似值;
步骤二、对控制方程中的时间导数使用三阶TVD Runge-Kutta离散公式将空间半离散的有限差分格式离散成时空全离散有限差分格式;
步骤三、根据时空全离散有限差分格式得到下一时间层上的近似值,依次迭代,得到终止时刻计算区域内流场的数值模拟值;
其中,采用TWENO格式重构通量的近似值包括:
其中,L(U)表示-f(U)x的空间离散形式,U=(ρ,ρu,E)T表示守恒变量,f(U)=(ρu,ρu2+p,u(E+p))T表示通量,f(U)x表示f(U)对x求导,x表示空间变量,ρ、u、p、E分别表示流体密度、速度、压强、能量,T表示转置,Ui(t)表示U在目标网格单元Ii内点xi处的值U(xi,t),和分别表示通量f(U)在目标网格单元Ii的边界和处的五阶近似的数值通量;
步骤2、将目标网格单元Ii以及其周围共五个网格单元组成一个大模板T1=[Ii-2,Ii-1,Ii,Ii+1,Ii+2],从大模板中选择两个包含两个网格单元的小模板T2=[Ii-1,Ii]和T3=[Ii,Ii+1];
步骤3、在T1、T2、T3每个模板上分别重构三角函数多项式p1(x)、p2(x)和p3(x),使得:
p2(x),p3(x)∈span{1,sin(x-xi)};
步骤4、任意取三组线性权,每组线性权的和均为1;
步骤5、根据三角函数多项式计算光滑指示器,用于衡量三角函数多项式pl(x)在目标网格单元上的光滑度;l=1,2,3表示对应模板序号;
步骤6、通过线性权和光滑指示器计算非线性权;
将计算结果代入含有时间导数项的空间半离散的有限差分格式,得到关于时间导数的常微分方程。
5.如权利要求1所述的方法,其特征在于:所述步骤4中,三组线性权为:
γ1=0.98,γ2=0.01,γ3=0.01;
γ1=1/3,γ2=1/3,γ3=1/3;
γ1=0.01,γ2=0.495,γ3=0.495。
10.如权利要求1所述的方法,其特征在于:所述步骤三中,时空全离散有限差分格式为关于时间层的迭代公式,初始状态值已知,通过迭代公式求出下一时间层的近似值,依次得到终止时刻计算区域内流场的数值模拟值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210373696.7A CN114757070A (zh) | 2018-05-16 | 2018-05-16 | 用于数值模拟的三角函数框架下新weno格式构造方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810472192.4A CN108763683B (zh) | 2018-05-16 | 2018-05-16 | 一种三角函数框架下新weno格式构造方法 |
CN202210373696.7A CN114757070A (zh) | 2018-05-16 | 2018-05-16 | 用于数值模拟的三角函数框架下新weno格式构造方法 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810472192.4A Division CN108763683B (zh) | 2018-05-16 | 2018-05-16 | 一种三角函数框架下新weno格式构造方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114757070A true CN114757070A (zh) | 2022-07-15 |
Family
ID=64008294
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210373696.7A Pending CN114757070A (zh) | 2018-05-16 | 2018-05-16 | 用于数值模拟的三角函数框架下新weno格式构造方法 |
CN201810472192.4A Active CN108763683B (zh) | 2018-05-16 | 2018-05-16 | 一种三角函数框架下新weno格式构造方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810472192.4A Active CN108763683B (zh) | 2018-05-16 | 2018-05-16 | 一种三角函数框架下新weno格式构造方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN114757070A (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110069854A (zh) * | 2019-04-22 | 2019-07-30 | 南京航空航天大学 | 多重分辨tweno格式对可压流场问题的模拟方法 |
CN110781626A (zh) * | 2019-10-31 | 2020-02-11 | 南京航空航天大学 | 有限差分多重分辨三角函数weno格式的模拟方法 |
CN111159956B (zh) * | 2019-12-10 | 2021-10-26 | 北京航空航天大学 | 一种基于特征的流场间断捕捉方法 |
CN111177965B (zh) * | 2019-12-25 | 2022-06-17 | 南京航空航天大学 | 基于定常问题求解的多重分辨weno格式定点快速扫描方法 |
CN111563314B (zh) * | 2020-03-23 | 2023-09-12 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
CN112163312B (zh) * | 2020-08-17 | 2022-08-02 | 空气动力学国家重点实验室 | 一种通过高阶weno格式降阶对可压缩流动问题进行数值模拟的方法 |
CN112100835B (zh) * | 2020-09-06 | 2022-06-14 | 西北工业大学 | 一种适用于复杂流动的高效高精度翼型绕流数值模拟方法 |
CN114707254B (zh) * | 2022-06-01 | 2022-08-26 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于模板构造法的二维边界层网格生成方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014043846A1 (zh) * | 2012-09-18 | 2014-03-27 | Lu Ming | 求解二维黎曼问题模拟亚音速无粘流的数值方法 |
CN104091065A (zh) * | 2014-07-03 | 2014-10-08 | 南京信息工程大学 | 一种求解浅水问题模拟间断水流数值的方法 |
CN107220399A (zh) * | 2017-03-23 | 2017-09-29 | 南京航空航天大学 | 基于埃尔米特插值基本加权无振荡格式的全流场模拟方法 |
-
2018
- 2018-05-16 CN CN202210373696.7A patent/CN114757070A/zh active Pending
- 2018-05-16 CN CN201810472192.4A patent/CN108763683B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108763683B (zh) | 2022-04-01 |
CN108763683A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108763683B (zh) | 一种三角函数框架下新weno格式构造方法 | |
Renson et al. | Numerical computation of nonlinear normal modes in mechanical engineering | |
Loseille et al. | Achievement of global second order mesh convergence for discontinuous flows with adapted unstructured meshes | |
CN105718634B (zh) | 一种基于非概率区间分析模型的翼型鲁棒优化设计方法 | |
Chen et al. | A sample-efficient deep learning method for multivariate uncertainty qualification of acoustic–vibration interaction problems | |
CN102436550B (zh) | 复杂边界及实际地形上溃坝洪水的自适应模拟方法 | |
Pederson et al. | The Sedov blast wave as a radial piston verification test | |
CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
Liu et al. | Three-dimensional high-order least square-based finite difference-finite volume method on unstructured grids | |
CN111177965B (zh) | 基于定常问题求解的多重分辨weno格式定点快速扫描方法 | |
Jouhaud et al. | A kriging approach for CFD/wind-tunnel data comparison | |
Chen et al. | Towards high-accuracy deep learning inference of compressible turbulent flows over aerofoils | |
Zangeneh | Development of a new algorithm for modeling viscous transonic flow on unstructured grids at high Reynolds numbers | |
Zhu et al. | An isogeometric approach to Biot-Cosserat continuum for simulating dynamic strain localization in saturated soils | |
CN111159956B (zh) | 一种基于特征的流场间断捕捉方法 | |
Song et al. | A Surrogate Model for Shallow Water Equations Solvers with Deep Learning | |
Ranocha et al. | On error-based step size control for discontinuous Galerkin methods for compressible fluid dynamics | |
Šmolík et al. | 3D vector field approximation and critical points reduction using radial basis functions | |
CN103425835B (zh) | 一种用于研究水-气界面演化过程的数值方法 | |
Yang et al. | An immersed boundary method based on parallel adaptive Cartesian grids for high Reynolds number turbulent flow | |
Breton et al. | On the prediction of tip vortices in the near wake of the MEXICO rotor using the actuator surface method | |
Zhiping et al. | Study on effects of thickness on airfoil-stall at low Reynolds numbers by cusp-catastrophic model based on GA (W)-1 airfoil | |
Kasturi Rangan et al. | A face-based immersed boundary method for compressible flows using a uniform interpolation stencil | |
CN116702571B (zh) | 基于多重光滑度量因子的数值模拟方法及装置 | |
CN105488265A (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 |