CN108763683B - 一种三角函数框架下新weno格式构造方法 - Google Patents
一种三角函数框架下新weno格式构造方法 Download PDFInfo
- Publication number
- CN108763683B CN108763683B CN201810472192.4A CN201810472192A CN108763683B CN 108763683 B CN108763683 B CN 108763683B CN 201810472192 A CN201810472192 A CN 201810472192A CN 108763683 B CN108763683 B CN 108763683B
- Authority
- CN
- China
- Prior art keywords
- format
- tweno
- value
- time
- trigonometric function
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 230000004907 flux Effects 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 24
- 238000004088 simulation Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 8
- 230000007704 transition Effects 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000035939 shock Effects 0.000 abstract description 14
- 230000010355 oscillation Effects 0.000 abstract description 9
- 239000002699 waste material Substances 0.000 description 2
- 101000635799 Homo sapiens Run domain Beclin-1-interacting and cysteine-rich domain-containing protein Proteins 0.000 description 1
- 102100030852 Run domain Beclin-1-interacting and cysteine-rich domain-containing protein Human genes 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229910001651 emery Inorganic materials 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种三角函数框架下新WENO格式构造方法,一方面,用三角函数多项式构造的加权基本无振荡格式比用代数多项式构造的经典基本无振荡格式更易于模拟波类或高频振荡类问题,在光滑区域可以获得高阶数值精度,在激波和接触间断处保持基本无振荡的性质;另一方面,新的TWENO格式虽然和经典五阶WENO格式用了同样五个点的信息,但能得到更低的全局L1和L∞范数截断误差。新的TWENO格式采用的线性权不再需要通过繁冗的数值计算得到最优解,可设为满足和为一的任意正数,相比经典的WENO格式更简便,鲁棒性更强,更易于推广到高维空间。该新型TWENO格式有效地数值模拟了几个经典的欧拉问题,充分验证了有效性。
Description
技术领域
本发明属于计算流体力学工程技术领域,具体涉及了一种三角函数框架下新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离散公式将半离散有限差分格式离散成时空全离散有限差分格式;
步骤三、根据时空全离散有限差分格式得到下一时间层上的近似值,依次迭代,得到计算区域内终止时刻流场的数值模拟值。
为优化上述技术方案,采取的具体措施还包括:
所述步骤一中,双曲守恒律方程为:
其半离散格式的形式为:
其中,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以及其周围共五个网格单元组成一个大模板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=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,其计算公式为:
将计算结果代入含有时间导数项的半离散有限差分格式,得到关于时间导数的常微分方程。
所述步骤3中,具体步骤如下:
步骤3.1、在三个模板T1、T2和T3上分别构造三角函数多项式p1(x)、p2(x)和p3(x),使其满足:
步骤3.2、得到每个模板上的三角函数插值多项式p1(x)、p2(x)和p3(x),如下:
所述步骤二中,利用三阶TVD Runge-Kutta离散公式:
得到时空全离散有限差分格式,其中,U(1),U(2)为中间过渡值,Δt为时间步长,上标n表示第n时间层,L(Un),L(U(1)),L(U(2))为-f(U)x的高阶空间离散形式的近似值。
所述步骤三中,时空全离散有限差分格式为关于时间层的迭代公式,初始状态值已知,通过迭代公式求出下一时间层的近似值,依次得到终止时刻计算区域内的数值模拟值。
本发明的有益效果是:相比于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是实施例四中初值条件为(19)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.25时的密度等值线图,分别采用线性权①、②、③。
图7a-7c是实施例四中初值条件为(20)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.25时的密度等值线图,分别采用线性权①、②、③。
图8a-8c是实施例四中初值条件为(21)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.3时的密度等值线图,分别采用线性权①、②、③。
图9a-9c是实施例四中初值条件为(22)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.2时的密度等值线图,分别采用线性权①、②、③。
图10a-10c是实施例四中初值条件为(23)的二维Euler黎曼问题,利用本发明的有限差分TWENO格式得到的t=0.3时的密度等值线图,分别采用线性权①、②、③。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
本发明给出了笛卡尔网格下解可压流场问题的新型五阶有限差分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)处。涡流可看成平均流的速度,温度和熵的扰动,表示为:
其中,τ=r/rc,ε=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给出了初值条件为(18)时二维Euler黎曼问题在t=0.25时刻的密度等值线图。图7a-7c给出了初值条件为(19)时二维Euler黎曼问题在t=0.25时刻的密度等值线图。图8a-8c给出了初值条件为(20)时二维Euler黎曼问题在t=0.3时刻的密度等值线图。图9a-9c给出了初值条件为(21)时二维Euler黎曼问题在t=0.2时刻的密度等值线图。图10a-10c给出了初值条件为(22)时二维Euler黎曼问题在t=0.3时刻的密度等值线图。从图中可以看出基于三角函数多项式空间的有限差分TWENO格式对本发明任意取得的线性权都可以很好的捕捉到黎曼问题的大部分流动特性。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (4)
1.一种三角函数框架下新WENO格式构造方法,在笛卡尔坐标系下,利用TWENO格式对可压流场问题进行数值模拟,其特征在于,包括以下步骤:
步骤一、把双曲守恒律方程离散为空间半离散的有限差分格式,采用TWENO格式重构通量的近似值;
步骤二、对控制方程中的时间导数使用三阶TVD Runge-Kutta离散公式将半离散有限差分格式离散成时空全离散有限差分格式;
步骤三、根据时空全离散有限差分格式得到下一时间层上的近似值,依次迭代,得到计算区域内终止时刻流场的数值模拟值;
所述步骤一中,双曲守恒律方程为:
其半离散格式的形式为:
其中,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以及其周围共五个网格单元组成一个大模板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=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,其计算公式为:
将计算结果代入含有时间导数项的半离散有限差分格式,得到关于时间导数的常微分方程。
4.如权利要求3所述的一种三角函数框架下新WENO格式构造方法,其特征在于:所述步骤三中,时空全离散有限差分格式为关于时间层的迭代公式,初始状态值已知,通过迭代公式求出下一时间层的近似值,依次得到终止时刻计算区域内的数值模拟值。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210373696.7A CN114757070A (zh) | 2018-05-16 | 2018-05-16 | 用于数值模拟的三角函数框架下新weno格式构造方法 |
CN201810472192.4A CN108763683B (zh) | 2018-05-16 | 2018-05-16 | 一种三角函数框架下新weno格式构造方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810472192.4A CN108763683B (zh) | 2018-05-16 | 2018-05-16 | 一种三角函数框架下新weno格式构造方法 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210373696.7A Division CN114757070A (zh) | 2018-05-16 | 2018-05-16 | 用于数值模拟的三角函数框架下新weno格式构造方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108763683A CN108763683A (zh) | 2018-11-06 |
CN108763683B true CN108763683B (zh) | 2022-04-01 |
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 Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210373696.7A Pending CN114757070A (zh) | 2018-05-16 | 2018-05-16 | 用于数值模拟的三角函数框架下新weno格式构造方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN114757070A (zh) |
Families Citing this family (9)
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格式的构造方法 |
CN111859819B (zh) * | 2020-06-16 | 2024-09-06 | 空气动力学国家重点实验室 | 一类高阶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 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于模板构造法的二维边界层网格生成方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091065A (zh) * | 2014-07-03 | 2014-10-08 | 南京信息工程大学 | 一种求解浅水问题模拟间断水流数值的方法 |
CN107220399A (zh) * | 2017-03-23 | 2017-09-29 | 南京航空航天大学 | 基于埃尔米特插值基本加权无振荡格式的全流场模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160306907A1 (en) * | 2012-09-18 | 2016-10-20 | Ming Lu | Numerical method for solving the two-dimensional riemann problem to simulate inviscid subsonic flows |
-
2018
- 2018-05-16 CN CN202210373696.7A patent/CN114757070A/zh active Pending
- 2018-05-16 CN CN201810472192.4A patent/CN108763683B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091065A (zh) * | 2014-07-03 | 2014-10-08 | 南京信息工程大学 | 一种求解浅水问题模拟间断水流数值的方法 |
CN107220399A (zh) * | 2017-03-23 | 2017-09-29 | 南京航空航天大学 | 基于埃尔米特插值基本加权无振荡格式的全流场模拟方法 |
Non-Patent Citations (6)
Title |
---|
"New Finite Volume Weighted Essentially Nonoscillatory Schemes on Triangular Meshes";Jun Zhu et al.;《SIAM Journal on Scientific Computing》;20180131;第1-36页 * |
"A new fifth order finite difference WENO scheme for solving hyperbolic conservation laws";Jun Zhu et al.;《Journal of Computational Physics》;20160507;第110-121页 * |
"Trigonometric WENO Schemes for Hyperbolic Conservation Laws and Highly Oscillatory Problems";Jun Zhu et al.;《Communications in Computational Physics》;20101130;第1-32页 * |
"加权基本无振荡格式研究进展";赵海洋等;《力学季刊》;20050331;第26卷(第1期);第87-95页 * |
"返回舱动态稳定性物理机理分析及被动/主动控制方法研究";赵海洋;《中国博士学位论文全文数据库》;20090715;第1-152页 * |
"高精度WENO格式在三维外流场计算中的应用";刘耀峰等;《弹箭与制导学报》;20051231;第25卷(第3期);第81-83页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114757070A (zh) | 2022-07-15 |
CN108763683A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108763683B (zh) | 一种三角函数框架下新weno格式构造方法 | |
CN109902376B (zh) | 一种基于连续介质力学的流固耦合高精度数值模拟方法 | |
CN107220399A (zh) | 基于埃尔米特插值基本加权无振荡格式的全流场模拟方法 | |
CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
CN108280273B (zh) | 一种基于非等距网格下的有限体积流场数值计算方法 | |
CN110457806A (zh) | 基于交错网格的中心五阶weno格式的全流场模拟方法 | |
Pederson et al. | The Sedov blast wave as a radial piston verification test | |
CN110069854A (zh) | 多重分辨tweno格式对可压流场问题的模拟方法 | |
CN105760602A (zh) | 一种有限体积加权基本无振荡格式的全流场数值模拟方法 | |
Liu et al. | Three-dimensional high-order least square-based finite difference-finite volume method on unstructured grids | |
Daxini et al. | Parametric shape optimization techniques based on Meshless methods: A review | |
CN109740182A (zh) | 一种基于再生核粒子的无网格物理变形仿真方法 | |
CN109657408A (zh) | 一种再生核粒子算法实现结构线性静力学仿真方法 | |
CN111158059A (zh) | 基于三次b样条函数的重力反演方法 | |
CN111177965B (zh) | 基于定常问题求解的多重分辨weno格式定点快速扫描方法 | |
Ojha et al. | Adaptive high-order fluid-structure interaction simulations with reduced mesh-motion errors | |
Zangeneh | Development of a new algorithm for modeling viscous transonic flow on unstructured grids at high Reynolds numbers | |
CN116702571B (zh) | 基于多重光滑度量因子的数值模拟方法及装置 | |
Ranocha et al. | On error-based step size control for discontinuous Galerkin methods for compressible fluid dynamics | |
CN116611369B (zh) | 基于光滑度量量级与候选模板点数的插值方法及装置 | |
CN112307684A (zh) | 一种多重分辨率weno格式结合ilw边界处理的定点快速扫描方法 | |
Nguyen | Model reduction techniques for parametrized nonlinear partial differential equations | |
Deville et al. | Theoretical analysis of brush seals leakage using local computational fluid dynamics estimated permeability laws | |
Qu et al. | Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis | |
Abouri et al. | A stable Fluid-Structure-Interaction algorithm: Application to industrial problems |
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 |