Abstract We present an updated version of the open-source Hypersonics Task-based Research (HTR) solver for hypersonic aerothermodynamics. The solver, whose first version was presented in Di Renzo et al. (Comput. Phys. Commun. 255, 2020), is designed for direct numerical simulation (DNS) of canonical hypersonic flows at high Reynolds numbers in which thermo-chemical effects induced by high temperatures are relevant. The solver relies on high-order spatial discretization on structured meshes and efficient time integrators for stiff systems within the Regent/Legion software stack, which makes the code highly portable and scalable in CPU and GPU-based supercomputers. The new version herein presented includes several optimizations and new tools for data analysis, along with novel user option for hybrid skew-symmetric/targeted essentially non-oscillatory numerics, to offer higher computational efficiency and lower numerical dissipation at moderate Mach numbers, inclusion of a new combustion mechanism for methane and oxygen, and new recycling-rescaling inlet boundary conditions targeted to the simulation of fully developed turbulent boundary layers. Program summary Program Title: Hypersonics Task-based Research solver CPC Library link to program files: Developer’s respository link: https://github.com/stanfordhpccenter/HTR-solver.git Licensing provisions: BSD 2-clause Programming language: Regent, C++, and CUDA Journal Reference of previous version: Di Renzo, M., Fu, L., & Urzay, J. (2020). HTR solver: An open-source exascale-oriented task-based multi-GPU high-order code for hypersonic aerothermodynamics. Computer Physics Communications, 107262. Does the new version supersede the previous version?: Yes Reasons for the new version: Release of new features Summary of revisions: • New optional sixth-order hybrid scheme has been implemented (activated by the flag “hybridScheme” in the input file). The scheme combines the energy-preserving properties of a sixth-order skew-symmetric central difference scheme [1] in smooth flow regions with the shock-capturing properties of a sixth-order targeted essentially non-oscillatory (TENO) scheme at points where shocks are involved. The switch between the two schemes is controlled by a TENO sensor whose cutoff value is adapted based on the maximum value of a modified Ducros sensor [2] across the reconstruction stencil. If the flag “hybridScheme” is set to false, the numerical scheme will revert to the TENO6-A scheme released in the previous version of the solver [3]; • New recycling-rescaling inflow boundary conditions [4,5] for the simulation of turbulent compressible boundary layers are now available to the user; • Support for Legion tracing, which significantly improves the strong scalability of the solver, has been implemented; • A diagnostic tool to monitor the time evolution of the flow variables in a subvolume of the computational domain is now available; • A single-step chemistry mechanism for methane/oxygen combustion has been added to the mixtures handled by the HTR solver; • Sample scripts for strong scaling have been added to the “testcases” directory; • Unit test and regression test suites have been added to the repository; • The input file scheme has been modified in order to reduce verbosity and increase flexibility in specifying the boundary conditions and type of gas mixture; • Hyperbolic sine stretching functions have been made available to users during the grid generation process; • Computationally intensive tasks have been ported to C++ and CUDA in order to achieve higher efficiency on all hardware; • Several optimizations of the tasks body and mapper have been implemented in order to increase the computational efficiency and reduce the memory footprint. Nature of problem: This code solves the Navier–Stokes equations at hypersonic Mach numbers including finite-rate chemistry for dissociating air and multicomponent transport. The solver is designed for direct numerical simulations (DNS) of transitional and turbulent hypersonic turbulent flows under high-enthalpy conditions, and it accounts for thermochemical effects (vibrational excitation and chemical dissociation). Solution method: This code uses a low-dissipation sixth-order schemes for the spatial discretization of the conservation equations on Cartesian stretched meshes. Time advancement is carried out by either an explicit method if chemistry is slow, hence not introducing additional stiffness, or by an operator-splitting algorithm whereby chemical production rates are handled implicitly. Additional comments including restrictions and unusual features: The HTR solver builds on the runtime Legion [6,7] and is written in the programming language Regent [8,9] developed at Stanford University. Instructions for installation of the components are provided in the README file enclosed with the HTR solver and in the Legion repository [6]. References [1] S. Pirozzoli, Journal of Computational Physics 229 (2010) 7180–7190. https://doi.org/10.1016/j.jcp.2010.06.006 . F. Ducros, F. Laporte, T. Souleres, V. Guinot, P. Moinat, B. Caruelle, Journal of Computational Physics 161 (2000) 114–139. https://doi.org/10.1006/jcph.2000.6492 . M. Di Renzo, L. Fu, J. Urzay, Computer Physics Communications 255 (2020) 107262. https://doi.org/10.1016/j.cpc.2020.107262 . T. S. Lund, X. Wu, K. D. Squires, Journal of Computational Physics 140 (1996) 233–258. https://doi.org/10.1006/jcph.1998.5882 . S. Pirozzoli, M. Bernardini, F. Grasso, Journal of Fluid Mechanics 657 (2010) 361–393. https://doi.org/10.1017/S0022112010001710 . Legion web page, 2020. URL: https://legion.stanford.edu . M. Bauer, S. Treichler, E. Slaughter, A. Aiken, Legion: Expressing locality and independence with logical regions, International Conference for High Performance Computing, Networking, Storage and Analysis, SC (2012), IEEE. Regent web page, 2020. URL: http://regent-lang.org . E. Slaughter, W. Lee, S. Treichler, M. Bauer, A. Aiken, SC ’15: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis (2015) 1–12,. https://doi.org/10.1145/2807591.2807629 .

中文翻译：

---

HTR-1.2 求解器：基于高超声速任务的研究求解器 1.2 版

摘要 我们提出了用于高超声速空气热力学的开源高超声速基于任务的研究 (HTR) 求解器的更新版本。求解器，其第一个版本在 Di Renzo 等人中提出。(Comput. Phys. Commun. 255, 2020)，专为高雷诺数下规范高超声速流动的直接数值模拟 (DNS) 而设计，其中高温引起的热化学效应是相关的。求解器依赖于结构化网格上的高阶空间离散化和 Regent/Legion 软件堆栈中刚性系统的高效时间积分器，这使得代码在基于 CPU 和 GPU 的超级计算机中具有高度的可移植性和可扩展性。此处介绍的新版本包括多项优化和用于数据分析的新工具，以及混合偏斜对称/目标基本非振荡数值的新颖用户选项，以在中等马赫数下提供更高的计算效率和更低的数值耗散，包括新的甲烷和氧气燃烧机制，以及新的循环重新缩放入口边界用于模拟完全发展的湍流边界层的条件。程序概要 程序名称：Hypersonics Task-based Research solver CPC Library 程序文件链接：Developer's respository 链接：https://github.com/stanfordhpccenter/HTR-solver.git 许可条款：BSD 2-clause 编程语言：Regent、C++ , 和 CUDA 之前版本的期刊参考：Di Renzo, M., Fu, L., & Urzay, J. (2020)。HTR 求解器：用于高超音速空气热力学的开源的面向百亿亿次任务的多 GPU 高阶代码。Computer Physics Communications, 107262。新版本是否取代了以前的版本？：是 新版本的原因：新功能的发布 修订摘要： • 新的可选六阶混合方案已实施（由标志“hybridScheme”激活）在输入文件中）。该方案将平滑流动区域中六阶偏对称中心差分方案 [1] 的能量保存特性与六阶定向基本非振荡 (TENO) 方案在冲击点处的冲击捕获特性相结合。涉及。两种方案之间的切换由一个 TENO 传感器控制，该传感器的截止值根据重建模板上修改后的 Ducros 传感器 [2] 的最大值进行调整。如果标志“hybridScheme”设置为false，数值方案将恢复到以前版本求解器中发布的TENO6-A方案[3]；• 用户现在可以使用用于模拟湍流可压缩边界层的新的再循环-重新缩放流入边界条件 [4,5]。• 支持 Legion 追踪，显着提高了求解器的强大可扩展性；• 用于监测计算域子体积中流动变量时间演变的诊断工具现已可用；• 甲烷/氧气燃烧的单步化学机制已添加到 HTR 求解器处理的混合物中；• 强扩展的示例脚本已添加到“testcase”目录中；• 单元测试和回归测试套件已添加到存储库中；• 修改了输入文件方案，以便在指定边界条件和气体混合物类型时减少冗长并增加灵活性。• 在网格生成过程中，用户可以使用双曲正弦拉伸函数；• 计算密集型任务已移植到 C++ 和 CUDA，以便在所有硬件上实现更高的效率；• 对任务主体和映射器进行了多项优化，以提高计算效率并减少内存占用。问题性质：此代码以高超音速马赫数求解 Navier-Stokes 方程，包括用于分离空气和多组分传输的有限速率化学。该求解器专为高焓条件下过渡和湍流高超声速湍流的直接数值模拟 (DNS) 而设计，并考虑了热化学效应（振动激发和化学解离）。求解方法：此代码使用低耗散六阶方案对笛卡尔拉伸网格上的守恒方程进行空间离散化。如果化学反应缓慢，则通过显式方法执行时间推进，因此不会引入额外的刚度，或者通过算子分裂算法隐式处理化学生产率。其他评论，包括限制和不寻常的功能： HTR 求解器建立在运行时 Legion [6,7] 上，并用斯坦福大学开发的编程语言 Regent [8,9] 编写。HTR 求解器附带的 README 文件和 Legion 存储库 [6] 中提供了安装组件的说明。参考文献 [1] S. Pirozzoli，Journal of Computational Physics 229 (2010) 7180–7190。https://doi.org/10.1016/j.jcp.2010.06.006。F. Ducros、F. Laporte、T. Souleres、V. Guinot、P. Moinat、B. Caruelle，计算物理学杂志 161 (2000) 114–139​​。https://doi.org/10.1006/jcph.2000.6492。M. Di Renzo、L. Fu、J. Urzay，计算机物理通信 255 (2020) 107262。https://doi.org/10.1016/j.cpc.2020.107262。TS Lund, X. Wu, KD Squires, Journal of Computational Physics 140 (1996) 233–258。https://doi.org/10.1006/jcph.1998.5882。S. Pirozzoli、M. Bernardini、F. Grasso，流体力学杂志 657 (2010) 361–393。https://doi.org/10.1017/S0022112010001710 。军团网页，2020。网址：https：//legion.stanford.edu 。M. Bauer、S. Treichler、E. Slaughter、A. Aiken，Legion：用逻辑区域表达局部性和独立性，高性能计算、网络、存储和分析国际会议，SC（2012 年），IEEE。摄政网页，2020 年。网址：http://regent-lang.org。E. Slaughter、W. Lee、S. Treichler、M. Bauer、A. Aiken，SC '15：高性能计算、网络、存储和分析国际会议论文集（2015 年）1-12，。https://doi.org/10.1145/2807591.2807629。M. Bauer, A. Aiken, SC '15: 高性能计算、网络、存储和分析国际会议论文集 (2015) 1–12,. https://doi.org/10.1145/2807591.2807629 。M. Bauer, A. Aiken, SC '15: 高性能计算、网络、存储和分析国际会议论文集 (2015) 1–12,. https://doi.org/10.1145/2807591.2807629。