适用于气相爆轰数值模拟的WENO-THINC有限顺风格式
A WENO-THINC Finite Downwind Scheme for the Simulations of Gaseous Detonation
钱琛庚 1杨同会 2昝文涛3
作者信息
- 1. 陆军勤务学院 油料系,重庆 401311;北京理工大学 爆炸科学与技术国家重点实验室,北京 100081
- 2. 清华大学 航天航空学院,北京 100084
- 3. 陆军勤务学院 油料系,重庆 401311;中国兵器科学研究院,北京 100089
- 折叠
摘要
针对WENO(Weighted Essentially Non-oscillatory)有限差分格式数值耗散较大,粗网格下模拟爆轰波位置不准确的问题,构造了高精度低耗散WENO-THINC有限顺风格式,分别使用WENO重构和THINC(Tangent of Hyperbola Interface Capturing)重构计算数值通量,选择与有限顺风通量最为接近的数值通量进行空间离散.该格式在接触间断处只需2～3个过渡点,在保证数值格式的稳定性的同时极大地提高了传统WENO格式在接触间断处的分辨率.数值试验表明,即便使用粗网格计算爆轰波传播问题时,该格式也能获得正确的爆轰波位置,且不会出现非物理结构,这说明其数值耗散很小,相比于传统WENO有限差分格式,更适用于工程尺度气相爆轰数值模拟.
Abstract
In response to the problem of high numerical dissipation in the WENO(Weighted Essentially Non-oscillatory)finite differ-ence scheme and inaccurate simulation of detonation wave positions in coarse grids,a high-resolution and low dissipation WENO-THINC finite downwind scheme was constructed.The numerical flux was calculated using WENO reconstruction and THINC(Tan-gent of Hyperbola Interface Capturing)reconstruction,respectively.The numerical flux closest to the finite downwind flux was se-lected for spatial discretization.This scheme only requires 2～3 transition points at contact discontinuities,ensuring the stability of the numerical format while greatly improving the resolution of the traditional WENO scheme at contact discontinuities.Numerical experiments show that the scheme can obtain the correct detonation wave position even when using coarse grid to calculate the deto-nation wave propagation problem,and there is no non-physical structure,which indicates that its numerical dissipation is very small.Compared with the traditional WENO finite difference scheme,it is more suitable for engineering scale numerical simulation of gas-eous detonation.
关键词
WENO有限差分格式/THINC重构/数值模拟/高精度数值格式Key words
WENO finite difference scheme/THINC reconstruction/numerical simulation/high order numerical scheme引用本文复制引用
基金项目
火灾与爆炸安全防护重庆市重点实验室(军用油库安全与设施设备评测实验室)开放基金(LQ21KFJJ03)
国家自然科学基金(12002043)
出版年
2024
国家科技期刊平台
NSTL
万方数据