| 题名 | 复杂流动间断问题的高效高精度算法及应用 |
| 作者 | 李诗尧 |
| 答辩日期 | 2023-05-23 |
| 文献子类 | 博士 |
| 授予单位 | 中国科学院大学 |
| 授予地点 | 北京 |
| 导师 | 申义庆 |
| 关键词 | 复杂流动数值模拟 加权基本无振荡格式 共用权 自适应网格加密方法 凝聚炸药爆轰 |
| 其他题名 | A Study of High-efficiency and High-accuracy Algorithms for Complex Discontinuous Problems |
| 学位专业 | 流体力学 |
| 英文摘要 | 针对流动中存在的间断现象，发展高分辨率高精度的数值求解方法一直是计算流体力学研究的重要内容。加权基本无振荡（Weighted essentially non-oscillatory，简称WENO）格式对光滑解具有较高阶的精度而在间断附近具有基本无振荡的性质，因而在计算流体力学中得到较多的研究与应用，WENO算法也得到不断的发展和完善。本文针对当前WENO格式的一些不足，研究高效、高精度低耗散的WENO格式，并将所发展的方法应用于凝聚炸药爆轰的数值模拟研究。
（1）提出新的权值计算方法，构造了新的WENO-Z型格式（下文简称WENO-ZN格式）,极大提高了临界点的精度。临界点具有较强的数学物理意义，如一阶、二阶、三阶临界点可刻画解的极值、拐点、曲率变化等几何特征，高阶临界点可刻画更小尺度上解的变化。为了提高极值点精度，WENO-Z格式利用整体光滑因子来计算权值。本文通过分析已有的WENO-Z型格式的不足，提出利用整体模板上的最高阶差分的平方作为整体光滑因子，以此使得新格式在尽可能高阶的临界点获得期望的精度；进一步构造了自适应函数代替原权值计算公式的常数1，该函数在光滑区趋于一个远大于1的数，而在含间断的模板其值远小于1，从而获得了一类精度高、耗散低、鲁棒性好的WENO-ZN格式。利用数值算例验证了所发展的系列格式的优越性能。
（2）提出和发展了应用于Euler方程组的共用权WENO（Common-Weights WENO, Co-WENO）格式。当WENO格式应用于方程组时，传统的方法是对方程组中的每个分量分别计算权值及加权组合。本文提出的共用权WENO格式，是指方程组中的每个分量皆利用同一组权值进行加权，因此共用权WENO格式具有两个优点：其一是由于只需计算一组权值，因此该方法具有较高的计算效率；其二是共用权方法保证了同一个模板或者是同一个单元对方程组中的每一个分量具有相同的贡献。本文基于Euler方程组的通量分裂方法，发展了系列的共用权WENO格式。数值算例表明，共用权WENO格式具有效率高、耗散低等性能。
（3）发展了高效的自适应网格加密方法。自适应网格加密方法能够根据需要对流场局部进行加密计算，从而在初始粗网格上通过加密得到与对应均匀细网格几乎一致的计算结果，能显著提升计算效率。在本文的研究工作中，首先提出了一种新的流场加密判定方法，其需要的加密阈值具有较好的普适性；其次，在利用粗网格数值向加密网格的插值过程中，我们提出了共用权WENO插值方法，与通量构造的共用权WENO格式相似，共用权WENO插值极大地提高了网格自适应加密算法的计算效率及稳定性。
（4）凝聚态炸药爆轰研究在国防军事发展和国民经济建设中具有重要的应用价值。由于其复杂的状态方程，针对理想气体状态方程所发展的算法不能直接推广应用。本文发展了适用于凝聚炸药爆轰的Steger-Warming通量分裂方法和半隐式TVD Runge-Kutta时间迭代方法，并结合前文发展的数值格式，模拟了钝感高能炸药的起爆与爆轰波传播过程。数值试验表明本文发展的计算方法在爆轰模拟中也有较高的计算精度和分辨率以及良好的收敛性，能够准确模拟爆轰波衍射过程中的“死区”（Dead Zone）现象和爆轰波碰撞中的正规反射与马赫反射等流动细节。
|
| 语种 | 中文 |
| 内容类型 | 学位论文 |
| 源URL | [http://dspace.imech.ac.cn/handle/311007/92343]  |
| 专题 | 力学研究所_高温气体动力学国家重点实验室 |
| 通讯作者 | 李诗尧 |
| 推荐引用方式 GB/T 7714 | 李诗尧. 复杂流动间断问题的高效高精度算法及应用[D]. 北京. 中国科学院大学. 2023. |
| 个性服务 |
| 查看访问统计 |
| 相关权益政策 |
| 暂无数据 |
| 收藏/分享 |
除非特别说明，本系统中所有内容都受版权保护，并保留所有权利。
修改评论