CN104268322A - Weno差分方法的一种边界处理技术 - Google Patents
Weno差分方法的一种边界处理技术 Download PDFInfo
- Publication number
- CN104268322A CN104268322A CN201410476535.6A CN201410476535A CN104268322A CN 104268322 A CN104268322 A CN 104268322A CN 201410476535 A CN201410476535 A CN 201410476535A CN 104268322 A CN104268322 A CN 104268322A
- Authority
- CN
- China
- Prior art keywords
- weno
- boundary
- difference method
- treatment technology
- carry out
- 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
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种利用WENO差分数值方法进行流场计算时的一种边界处理技术,能够使得进行流场计算时达到一致高精度,并且能够处理复杂边界,结构网格上该发明的主要技术点和步骤为:(1)网格划分,将网格划分为内部区域和边界区域两个区域,边界区域只有一层网格,从而尽量减少边界区域DG的复杂计算。(2)空间离散,对内部区域利用WENO有限差分方法进行空间离散。(3)边界处理,对边界区域利用DG方法进行空间离散。(4)耦合处理,将边界处理部分与内部区域进行耦合,边界区域通过DG多项式逼近提供导数信息。(5)时间离散,通过TVD-Runge-Kutta方法进行显示时间离散。(6)后处理,通过Tecplot或者Paraview等商业软件对流场结果进行可视化模拟。
Description
技术领域
本发明涉及计算流体力学数值方法领域,特别是涉及一种求解双曲守恒律方程的高精度有限差分方法的边界处理技术。
背景技术
飞行器三维复杂流动的数值模拟和相关多目标优化问题是目前计算流体力学中的前沿热点问题,同时也是一个面向工程实际需求的应用问题。然而在目前的计算机规模和求解能力的条件下,目前主流的数值方法在计算效率上还不能满足这一工程实际应用问题的需要,解决这个问题的关键之一是提高流场解算器的效率。
目前已经出现了多种为了充分利用流行的高精度数值方法的优点的耦合方法,例如WENO限制器间断有限元方法,高阶的SV方法,耦合的DG/FV方法,多区域耦合的DG和WENO方法,其中多区域耦合的DG和WENO方法具有高效、高精度且易处理复杂边界的优点,但是目前的所有的高阶方法中效率最高的仍然为有限差分方法,但是该方法一般都需要在结构网格上进行计算,且不能够处理复杂边界,处理边界时也没有紧致性,因而缺少一种能够帮助有限差分方法紧致处理复杂边界的高阶数值方法。
发明内容
为了克服现有技术的不足,我们提出了一种能够处理复杂边界的WNEO有限差分方法作为处理WENO差分方法的一种边界处理技术,该方法发展了多区域耦合的DG和WENO-FD方法,更加充分的利用了有限差分方法的优势和DG的优势,几乎整个区域都利用WENO-FD方法计算,从而效率更高,仅仅在边界的一层网格处使用DG方法,边界处理利用DG方法处理,从而能够紧致处理复杂边界。
本发明的技术方案是:
将所求解的问题进行结构网格剖分,然后将所求解区域分为两种单元区域,一种是内部单元区域,另一种是边界单元区域,而且边界单元区域只有一层计算网格。
利用传统的WENO-FD方法对内部单元区域进行空间离散,利用传统的DG方法对边界单元区域进行空间离散,并且处理边界。
利用多区域耦合DG和WENO方法中的耦合方式实现WENO-FD方法和DG方法的耦合处理。耦合过程中需要构造内部单元区域边界处的数值通量,WENO数值通量的构造方法使用特殊的HWENO构造过程实现。
边界单元处DG方法计算得到该单元上解的n阶多项式逼近结果,通过这个多项式逼近结果可以求出该单元中心点处的值以及不超过n阶的导数值。
在空间离散结束之后,利用高阶的TVD Runge-Kutta方法求解得到下一时刻的数值结果。
利用Tecplot等可视化软件对二维标量双曲守恒律方程和二维Euler方程组等验证算例的结果进行可视化处理。
本发明的有益效果是:
本发明充分利用了高阶WENO有限差分方法高效简单的优点,结合了DG方法能够处理复杂边界并且具有局部性和紧致性的特点,发展了多区域耦合DG和WENO方法,提出了一种基于DG处理边界的WNEO-FD方法,该方法仅仅在边界一层网格处使用DG方法,在处理大规模问题是效率几乎等价于WENO-FD方法,而且也具有处理复杂边界的能力,并且在处理复杂边界时不再需要构造多个虚拟网格点,减小了WENO-FD边界处理时的模板大小从而处理时具有紧致性和局部性。
附图说明
图1为一维算例网格划分示意图
图2为内部单元处传统WENO模板示意图
图3为内部单元边界处HWENO模板示意图
图4为一维标量Burgers方程算例示意图
图5为一维精度验证表格示意图
图6为二维双马赫反射算例结构网格划分示意图
图7为双马赫反射算例描述图
图8为二维双马赫反射算例密度等值线示意图
图9为算法实现过程流程图
具体实施方式
首先我们针对一维标量Burgers方程说明我们的算法主要过程,如图1第一步的区域划分为两个部分,其中左边WENO-FD部分是内部单元区域,右边一层DG网格为一层边界单元区域。
内部区域WENO-FD的具体实施过程如下:
首先我们不考虑通量***,假设f'(u)>0,那么在单元Ii守恒的有限差分格式形式为
其中ui(t)=u(xi,t)是单元Ii中点处的点值,为数值通量可以表示为
如图2除了内部单元边界两个通量外,大部分内部单元通量我们选用图2的模板S=(S1,S2,S3),通量 其中 然后 是由每一个小模板上包含单元的值通过拉格朗日插值得到的,对于5rd的WENO-FD小模板是通过3rd拉格朗日插值得出,最终的插值结果为:
非线性权w1,w2,w3的计算过程中满足wr=dr+O(Δx3):
其中βr被称作“光滑因子”,这里选用光滑因子:
非线性权的计算过程需要计算线性权: 这样的线性系数
如图3,如果我们要计算内部单元边界的两个单元边界通量,我们的方法中使用更小的模板,图中的虚拟网格点我们不再需要,而是在这两个相对小的模板上采用特殊的Hermite插值得到我们要的每一个小模板的插值结果Hermite插值的结果为:
和
然后再通过同样的方式得到我们需要的非线性权和光滑因子,从而可以求解方程得到完整的解。
根据技术方案中的步骤我们还需要通过DG方法得到如图3边界单元I1andIN处单元中点处的值以及高阶导数值,例如IN单元得到的多项式逼近解为:
u(x)=u0p0(x)+u1p1(x)+u2p2(x)
系数u0,u1,u2是由DG方法逼近得到的自由度。高阶导数信息也很容易通过多项式求导得出
u′(x)=u0p′0(x)+u1p′1(x)+u2p′2(x)
u″(x)=u0p″0(x)+u1p″1(x)+u2p″2(x)
根据技术方案中的过程我们需要将两种方法的边界耦合到一起,DG和HWENO耦合过程如下:
其中代表通过通常的WENO重构得到的通量,代表DG方法得到的通量, 是通过上述HWENO重构的方式得到。
为了验证上述算法,如图4我们针对一维标量的Burgers方程计算了强激波算例,并且我们针对线性的一维算例进行了算法的精度验证如图3。
为了验证上述算法,我们针对二维双马赫反射算例进行了验证,双马赫反射问题的网格剖分如图6,其中上面的网格部分代表内部单元网格,最下面一层网格代表边界单元网格。
如图7,假设初始时刻,马赫数Ma=10的正激波向平板运动。激波与平板之间成60°夹角。通过上述描述的间断有限元处理边界的WENO有限差分方法计算流场的演化。
如图3和图7,我们根据我们所述的间断有限元处理边界的WENO有限差分方法的计算结果,通过Tecplot可视化软件对我们的流场计算结果进行了可视化演示,图7中演示了双马赫反射问题中流场密度的变化,计算网格为320*80,初始马赫数为Ma=10,图中演示了从ρ=1.3965到ρ=22.682总共30条密度等值线。
本发明中没有详细说明的内容均属于本技术领域人员熟知的常识。
Claims (7)
1.WENO差分方法的一种边界处理技术,其特征在于所述方法结构网格上的实现步骤如下:
A首先将所求解的问题进行结构网格剖分,然后将所求解区域分为两种单元区域,一种是内部单元区域,另一种是边界单元区域,而且边界单元区域只有一层计算网格。
B利用传统的WENO-FD方法对内部单元区域进行空间离散,利用传统的DG方法对边界单元区域进行空间离散,并且处理边界。
C利用多区域耦合DG和WENO方法中的耦合方式实现WENO-FD方法和DG方法的耦合处理。耦合过程中需要构造内部单元区域边界处的数值通量,WENO数值通量的构造方法使用特殊的HWENO构造过程实现。
D边界单元处DG方法计算得到该单元上解的n阶多项式逼近结果,通过这个多项式逼近结果可以求出该单元中心点处的值以及不超过n阶的导数值。
E在空间离散结束之后,利用高阶的TVD Runge-Kutta方法求解得到下一时刻的数值结果。
F利用Tecplot或者Paraview等可视化软件对利用上述方法求出的流场结果进行可视化实现,包括流场的密度、速度、压力等。
2.根据权利要求1所述的WENO差分方法的一种边界处理技术,其特征在于我们在步骤A中进行网格剖分时,将边界处的一层网格划分为单独的边界单元区域,步骤B中进行空间离散时,边界处理过程完全由DG计算过程承担,并且不再需要构造多个虚拟网格单元,从而保证紧致性。
3.根据权利要求1所述的WENO差分方法的一种边界处理技术,其特征在于我们在步骤C中两种方法进行耦合时,耦合界面处WENO数值通量的构造,我们利用HWENO的思想,利用Hermite插值构造WENO-FD内部单元边界通量fWENO。
4.根据权利要求1中所述的WENO差分方法的一种边界处理技术,其特征在于我们步骤C中进行Hermite插值时需要的信息,一部分由内部单元处WENO-FD方法计算所得的单元中心处的值提供,另一部分由步骤D中DG方法多项式逼近结果得到的边界单元中心处的值以及高阶导数值提供。
5.根据权利要求1中所述的WENO差分方法的一种边界处理技术,其特征在于我们在WENO-FD方法内部单元边界处需要构造边界单元DG方法所需要的DG数值通量fDG(u-,u+),这个数值通量fDG(u-,u+)中的u-则需要通过WENO-FD方法提供,这里我们使用拉格朗日插值的方式提供的u-对应点估计值,其中拉格朗日插值需要的信息完全由内部单元提供。
6.根据权利要求1中所述的WENO差分方法的一种边界处理技术,其特征在于我们在步骤E中使用显示的TVD Runge-Kutta方法进行时间离散,离散过程中为了保持格式是稳定的,选用时间步长有一定的限制,Δx≤CFL*Δx,CFL=min(WENOCFL,DGCFL)。
7.根据权利要求1中所述的WENO差分方法的一种边界处理技术,其特征在于步骤F中我们利用Tecplot等可视化软件对二维标量双曲守恒律方程和二维Euler方程组等验证算例的结果进行可视化处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410476535.6A CN104268322B (zh) | 2014-06-27 | 2014-09-18 | Weno差分方法的一种边界处理技术 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410301536.7 | 2014-06-27 | ||
CN2014103015367 | 2014-06-27 | ||
CN201410301536 | 2014-06-27 | ||
CN201410476535.6A CN104268322B (zh) | 2014-06-27 | 2014-09-18 | Weno差分方法的一种边界处理技术 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104268322A true CN104268322A (zh) | 2015-01-07 |
CN104268322B CN104268322B (zh) | 2018-01-02 |
Family
ID=52159843
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410476535.6A Active CN104268322B (zh) | 2014-06-27 | 2014-09-18 | Weno差分方法的一种边界处理技术 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104268322B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104699909A (zh) * | 2015-03-25 | 2015-06-10 | 厦门大学 | 一种保持强稳定性的变步长多步法时间离散算法 |
CN106156471A (zh) * | 2015-04-16 | 2016-11-23 | 中国科学院计算机网络信息中心 | 一种多层次流场可视化方法 |
CN106201995A (zh) * | 2016-06-24 | 2016-12-07 | 张麟 | 一种图像边界元处理方法 |
CN106528903A (zh) * | 2015-09-15 | 2017-03-22 | 利弗莫尔软件技术公司 | 执行水下***的数字模拟的***和方法 |
CN108052738A (zh) * | 2017-12-13 | 2018-05-18 | 电子科技大学 | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 |
CN108229083A (zh) * | 2018-04-11 | 2018-06-29 | 南京航空航天大学 | 一种基于改进的有限差分格式的水流数值模拟方法 |
CN111563314A (zh) * | 2020-03-23 | 2020-08-21 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102830121A (zh) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | 一种软性磨粒流磨粒群实时检测方法 |
-
2014
- 2014-09-18 CN CN201410476535.6A patent/CN104268322B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102830121A (zh) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | 一种软性磨粒流磨粒群实时检测方法 |
Non-Patent Citations (3)
Title |
---|
JIANXIAN QIU 等: "Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous galerkin method:one dimensional case", 《JOURNAL OF COMPUTATIONAL PHYSICS》 * |
卢耀文: "求解双曲守恒律的非守恒型RKDG和WENO耦合方法", 《科技风》 * |
朱洪强: "简化的WENO重构在h自适应RKDG方法中的应用", 《科学技术与工程》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104699909B (zh) * | 2015-03-25 | 2018-11-02 | 厦门大学 | 一种保持强稳定性的变步长多步法时间离散方法 |
CN104699909A (zh) * | 2015-03-25 | 2015-06-10 | 厦门大学 | 一种保持强稳定性的变步长多步法时间离散算法 |
CN106156471A (zh) * | 2015-04-16 | 2016-11-23 | 中国科学院计算机网络信息中心 | 一种多层次流场可视化方法 |
CN106156471B (zh) * | 2015-04-16 | 2019-05-07 | 中国科学院计算机网络信息中心 | 一种多层次流场可视化方法 |
CN106528903A (zh) * | 2015-09-15 | 2017-03-22 | 利弗莫尔软件技术公司 | 执行水下***的数字模拟的***和方法 |
CN106528903B (zh) * | 2015-09-15 | 2020-10-13 | 利弗莫尔软件技术公司 | 执行水下***的数字模拟的***和方法 |
CN106201995A (zh) * | 2016-06-24 | 2016-12-07 | 张麟 | 一种图像边界元处理方法 |
CN106201995B (zh) * | 2016-06-24 | 2019-10-22 | 张麟 | 一种图像边界元处理方法 |
CN108052738A (zh) * | 2017-12-13 | 2018-05-18 | 电子科技大学 | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 |
CN108052738B (zh) * | 2017-12-13 | 2021-10-15 | 电子科技大学 | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 |
CN108229083A (zh) * | 2018-04-11 | 2018-06-29 | 南京航空航天大学 | 一种基于改进的有限差分格式的水流数值模拟方法 |
CN111563314A (zh) * | 2020-03-23 | 2020-08-21 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
CN111563314B (zh) * | 2020-03-23 | 2023-09-12 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104268322B (zh) | 2018-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104268322A (zh) | Weno差分方法的一种边界处理技术 | |
Li et al. | Wavelet-based numerical analysis: A review and classification | |
Lai et al. | Integrating CAD with Abaqus: a practical isogeometric analysis software platform for industrial applications | |
Gagliardini et al. | Capabilities and performance of Elmer/Ice, a new-generation ice sheet model | |
CN102332046B (zh) | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 | |
CN111859766B (zh) | 可变计算域的拉格朗日积分点有限元数值仿真***及方法 | |
Feng et al. | A compressible lattice Boltzmann finite volume model for high subsonic and transonic flows on regular lattices | |
CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
Gupta et al. | On the enrichment zone size for optimal convergence rate of the Generalized/Extended Finite Element Method | |
CN105426594A (zh) | 一种基于时空场和条件插值的平稳均质风场快速模拟方法 | |
Koric et al. | Evaluation of massively parallel linear sparse solvers on unstructured finite element meshes | |
Berger | Cut cells: Meshes and solvers | |
Xi et al. | A localized collocation scheme with fundamental solutions for long-time anomalous heat conduction analysis in functionally graded materials | |
Yu et al. | A high-order flux reconstruction/correction procedure via reconstruction formulation for unsteady incompressible flow on unstructured moving grids | |
CN104077479B (zh) | 一种基于守恒迎风格式获取参量阵声场空间分布的方法 | |
CN111079326B (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
Cao et al. | A multi-layered point reordering study of GPU-based meshless method for compressible flow simulations | |
CN107423542B (zh) | 一种适用于逐棒计算的非均匀泄漏修正方法 | |
CN107808021B (zh) | 基于cfd的流体设备的阻力计算方法 | |
Ning et al. | A grid generator for 3-D explosion simulations using the staircase boundary approach in Cartesian coordinates based on STL models | |
CN107832482A (zh) | 致密储层多尺度裂缝网络建模及模拟方法 | |
Zimmerman et al. | High-order spectral difference: verification and acceleration using GPU computing | |
Chen et al. | Shape optimization of fluid cooling channel based on Darcy reduced-order isogeometric analysis | |
CN104992046A (zh) | 流体力学计算***及方法 | |
CN103136431A (zh) | 一种自由曲面的保面积参数特性优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |