反问题的研究在多个领域中均展现出了重要的价值和广泛的应用前景.然而,地球物理学中的反问题由于其固有的复杂性和不确定性,始终是一个极具挑战性的课题.近年来,量子计算技术的迅速发展和应用为这一难题的解决提供了新的可能性.基于目前流行的两种量子计算范式,本文首先提出了一种利用量子退火的迭代反演算法.该算法结合了量子退火的全局搜索能力和迭代反演的局部优化特性.通过线性和非线性层析成像实验,本文提出的算法可以取得接近于经典方法的反演精度.此外,本文还讨论了变分量子线性系统算法在反问题求解中的潜力.介绍了该算法在反问题求解中的初步尝试,该方法通过构建和优化量子线路来求解反问题,展现出了良好的应用前景.本文研究为量子计算在地球物理学反问题求解中的应用提供了新的思路和方向.随着量子计算技术的不断发展和应用,预计量子计算方法将在未来大尺度地球物理反问题求解方面发挥重要作用.

工业软件作为现代工业体系的核心支撑，正在有力推动制造、能源、交通等关键领域的数字化转型与智能化升级。随着工业4.0和智能制造的快速发展，工业软件在降本增效、优化设计、提质出新等方面发挥着越来越重要的作用。然而，我国在工业软件的核心技术、自主创新和产业化应用方面仍面临诸多挑战。为了促进工业软件领域的学术交流与技术突破，本专辑聚焦于工业软件的创新算法、系统开发及其在工业场景中的应用，展示了该领域的最新研究成果和前沿进展。

本文主要介绍了数值计算的相关发展历程，以及现有的数值计算软件的发展和应用. 此外, 本文介绍了国产计算软件北太天元的发展历程、主要功能介绍、软件开发的困难和重要性，并对未来国产数值软件的发展进行讨论.

格子玻尔兹曼方法是求解玻尔兹曼输运方程的高效算法,在流体流动、声子导热等领域已深入应用. 本文针对声子导热的格子玻尔兹曼方法,基于我们提出的声子弛豫时间有效校正模型,设计了具有不同对称性和离散速度数量的格子模型,通过宏观尺度的热输运过程验证了算法和程序的正确性.进而研究了格子模型的对称性和离散速度数量对纳微尺度方腔热输运模拟的影响.结果表明,仅提高格子模型的对称性无法有效克服模拟过程的非物理温度跃迁,同时增大格子模型的离散速度数量和对称性能够有效削弱该非物理温度跃迁现象.

扫掠曲面生成是计算机辅助几何设计(Computer-Aided Design, CAD)领域中的关键技术, 在工程设计等领域有着广泛应用. 针对此问题, 本文首先给出了与给定B样条曲面几何拼接的单轨和双轨扫掠基础算法, 然后提出了一种拼接误差可控的扫掠曲面生成算法, 随后将这些算法生成的曲面与商业软件Catia和Rhino得到的结果进行了对比分析.在本文拼接误差可控的扫掠曲面生成算法中, 首先对给定曲面边界进行采样, 获取导数和参数等信息, 然后对这些信息进行B样条曲线插值. 通过该B样条曲线进行扫掠操作以获得扫掠曲面. 最后根据预设的扫掠曲面与给定曲面之间的拼接误差阈值, 对采样点进行调整, 以获得拼接误差范围内的扫掠曲面. 实验结果表明, 本文中的算法在生成曲面质量方面优于Rhino, 并且能以较少的控制点数量生成与Catia质量相当的曲面.

航空轮胎是决定飞机起降安全的关键部件,被誉为轮胎工业“皇冠上的明珠”.航空轮胎设计是我国的“卡脖子”技术难题,亟待突破航空轮胎结构数学建模与CAE工业软件开发.航空轮胎结构主要由帘线增强橡胶复合材料构成,需要满足复杂的服役工况. 因此,本文开展了航空轮胎帘线增强橡胶复合材料数学建模方法的研究,针对复合材料变形特点, 提出了帘线增强橡胶复合材料建模方法,建立了帘线与橡胶结构的内埋约束关系.推导了橡胶结构与帘线结构的有限元方程,开发了航空轮胎结构设计工业软件SuperTire, 实现了航空轮胎充气分析,并通过航空轮胎数值算例验证了建模精度与求解效率.

随着工艺的演进和电路集成度的提高, 电路仿真的规模不断扩大, 仿真验证时间也随之延长. 传统硬件架构在性能提升方面的局限性, 给设计工程师带来了巨大的挑战. 近年来, 因人工智能的发展, 通用图形处理单元(GPGPU)得到了快速的发展, 其算力和访存带宽都远超同时代的CPU, 也为电路仿真提供了加速的新机遇. 本文将探讨模拟电路和射频电路仿真所对应的数学模型及工业界常用的求解算法, 分析在异构仿真环境中面临的机遇与挑战. 结合先进的计算平台, 电路仿真可以更高效地予以实现. 尽管如此, 如何有效集成这些新技术以应对复杂电路的仿真需求, 仍是当前研究的重点课题.

投影类算法是一种求解大规模约束优化问题的重要方法,尤其在到约束集合的正交投影容易实现的情况下.在投影方法中,步长的选择对收敛速度有重要影响.然而,理论上的“最优”步长往往在实践场景中表现不佳.已有文献指出,在无约束优化中对最优步长采用一些松弛策略可以对最速下降方法起到加速效果.基于这种见解,本文将这些松弛策略拓展到约束凸二次优化中的投影算法,包括基本格式的投影算法和“预测-校正”型投影算法,得到相应的松弛型算法.数值结果表明,这些松弛策略均提高了上述两种采用最优步长的投影类算法的效果,且在不同问题规模下各有优势.

数值线性代数解法器是科学工程计算与工业软件的核心组件，也是影响这些软件计算效率的主要瓶颈，其高效算法设计与性能优化面临复杂应用特征与超级计算机体系结构特征的双重挑战，一直以来都是学术界和工业界广受关注的问题。自2018年开始，国内相关专家发起并组织了解法器快速算法及应用研讨会（Solver会议），该系列会议至今已成功举办了8届，已成为国内该领域研究人员发布研究成果的交流平台，吸引了工业界和实际应用部门的广泛关注。本系列专辑拟邀请Solver会议组织者担任客座编委，不定期组织活跃于该会议的专家和团队，围绕数值代数解法器的快速算法设计、性能优化、自主软件研发和实际应用撰写文章，展示该领域的最新进展，促进该领域在我国的发展。

本专辑由北京应用物理与计算数学研究所徐小文研究员等人（编委名单见后）负责组织，经过严格同行评审，最终录用了8篇文章，涵盖了自主解法器软件与算法库、应用驱动快速算法、面向国产处理器的性能优化、新型算法等主题，一定程度上反映了我国近年来该领域的研究进展。

基于机器学习的支持向量机(SVM)算法,提出了一种求解病态线性方程组的数值算法．将线性方程组的求解转化为寻求高维空间中的超平面回归问题,分别通过经典SVM和最小二乘SVM回归算法得到病态线性方程组的数值解．由于SVM的求解算法是凸二次优化问题,相较于其他神经网络算法,具有全局唯一解,能够得到误差较小的数值解．数值算例结果表明,该算法获得的数值结果具有较好的精度和稳定性．

随着求解规模不断扩大,直接法求解对称正定的刚度矩阵方程内存需求不断增加,利用充裕的计算机硬盘资源辅助求解,降低内存需求成为一种解决方案.提出3种超节点Cholesky核外分解策略.1)最大速度核外分解,保证每个超节点只被读取和保存一次;2)最小内存核外分解,分解时不参与任务池中任务的超节点一律释放内存.3)限制内存核外分解,根据给定内存大小决定采用策略1)还是策略2).各策略均在共享内存环境下使用OpenMP基于任务池实现异步并行.在求解某4984362阶刚度方程时,MUMPS核外求解失败;FEADS策略1)使用内存18.2GB,相较核内求解所需142.85GB节省87\%,用时2328.07s;PARDISO核外模式使用内存21.3GB,用时4643.18s.

随着计算机科学和信息技术的发展,计算机符号计算为复杂计算和数值模拟提供了一种高效的方法和手段.基于计算机符号计算,本文研究一类(2+1)维广义变系数浅水波方程,描述非线性浅水波的动力学特性与非线性作用机制.在Hirota双线性算子的框架下,推导了(2+1)维广义变系数浅水波方程双线性形式,进而将非线性问题转化为双线性问题,实现了问题的可解性.通过运用符号计算软件求解双线性方程,给出了孤子解的表达式、双线性形式Bäcklund变换和团块解.数值模拟展示了孤子解和团块解所描述的非线性浅水波动力学特性,该(2+1)维广义变系数浅水波方程的特殊解存在新颖特性,即团块沿y轴移动,在x轴方向的速度为零.

油藏地质参数一般非均质性强,但探测手段十分有限,导致地质模型通常具有很高的不确定性.无论常规历史拟合还是数据空间反演,均需在考虑地质模型不确定性的前提下预测油藏渗流的未来动态.对此,一般是基于油藏数值模拟进行预测,但对于反映地质模型不确定性的大量实现,油藏数值模拟耗时非常高.为此,本文将动态模式分解(D MD)推广至地质参数空间,即基于一系列反映地层参数不确定性的地质实现进行预测.油藏渗流动态变化随时间往往趋于平缓,因此,构建训练集和测试集,从而通过测试选出动态数据满足局部线性近似的时刻,进而利用DMD替代数值模拟进行接下来的预测,可大幅提升预测效率.利用单相与两相瞬态渗流算例进行了验证与说明.

本文基于求解线性方程组的贪婪随机Kaczmarz算法和贪婪几何概率随机Kaczmarz算法的思想,提出求解矩阵方程的新型贪婪随机算法.而后,利用重要不等式探讨贪婪几何概率随机Kaczmarz算法的收敛性.最后,通过数值实验验证算法的可行性和有效性.数值结果表明:对于大规模矩阵方程,贪婪几何概率随机Kaczmarz算法优于贪婪随机Kaczmarz算法.

在经典Poisson-Nernst-Planck (PNP)模型框架下,研究了外加振荡电场对纳米尺度管道中离子输运性质的影响.利用多尺度平均化方法和数值模拟,我们研究了三种情况:(1)沿管道方向施加时间振荡电场;(2)沿管道方向施加时间振荡电场,并与空间周期性分布的电场相结合;(3)在管道壁存在周期性表面电荷分布的情况下,施加沿管道方向的时间振荡电场.对于高振荡频率的情况,我们推导出了首阶近似的等效模型.结果表明,管道内离子分布和平均输运性质仅与外加振荡电场的时间和空间平均值有关,而与振荡频率无关.

流固耦合(Fluid-Structure Interaction,FSI)的准确模拟对于核反应堆的设计与安全运行至关重要.其中燃料组件的附加质量和阻尼比是重要的工程设计参数.本文利用基于高阶间断有限元(Discontinuous Galerkin,DG)的双向耦合FSI求解器进行数值模拟,通过模拟静水中悬臂梁的振动,计算振动频率和阻尼比,并与理论预测值进行了对比.此外,本文还给出了一个涉及湍流的流致振动数值算例,该算例采用大涡模拟(Large Eddy Simulation,LES)方法解析大尺度涡旋,并利用合成涡方法(Synthetic Eddy Method,SEM)生成入口边界条件.该求解器可以用于计算结构的附加质量和阻尼比,另外通过设置不同的入口条件和湍流模型,该求解器能够给出复杂工况下结构的振动情况,为反应堆的安全运行与设计提供依据.

随机算法近年来发展迅速, 展现出实际应用潜力, 为大规模线性代数方程组的求解提供了新的技术途径. 本文综述稀疏线性代数方程组随机算法的研究现状, 对当前主要的三类随机算法进行了总结, 分析了这些算法的特点、计算复杂度和面临的问题. 在此基础上, 本文对当前随机算法的研究现状进行了评估, 指出了离实际应用需求的差距, 对面向大规模实际应用的随机算法研究进行了展望.

基于求解偏微分方程边值问题的传统谱配点法,建立了用矢量形式表示的分数阶Chebyshev多项式的整数和分数阶积分及微分算子,构造了考虑固结方程及其边界和初始条件的解函数,将其代入分数阶黏弹性固结方程后,通过谱配点将其转化为一个代数方程组,建立了求解饱和黏土一维分数阶黏弹性固结方程谱配点法的数值方法.该算法是一种在时间域内求解分数阶黏弹性固结方程的直接解法,孔压、有效应力及沉降量等均为有限项级数组成的显式函数,可方便地计算随时间变化的荷载、分级加载等情况下的一维分数阶黏弹性固结问题.算例分析表明,该方法具有较高的计算精度和有效性.

本文针对均匀网格上的非线性退化抛物方程构造一种新型六阶不等距模板混合WENO格式（HUS-WENO）.该格式仅利用三个不等距空间模板上的信息,在光滑区域达到六阶精度,在间断处保持本质无振荡特性.此外,HUS-WENO格式的线性权可以任取和为1的正数.为降低计算成本,基于不等距模板WENO格式中六点模板的重构多项式设计一种混合策略,能自动、准确且高效地识别问题单元,并且不包含任何人工参数.通过数值算例,验证HUS-WENO格式在计算效率和高分辨率等方面的优越性.

本文给出求解一类垂直线性互补问题的模系矩阵分裂迭代法.通过把垂直线性互补问题重新转化为一个等价的非线性方程组,创建一类新的基于模的矩阵分裂迭代法,并在一定条件下证明了算法的收敛性.最后提供两个数值算例证明所提算法的有效性.

大规模代数特征值问题广泛应用于材料科学与工程结构分析等领域. 在多尺度复杂系统中, 随着问题规模的不断扩大以及尺度差异引发的误差累积, 传统特征值算法在精度、稳定性以及计算效率方面逐渐显现不足, 在部分场景下甚至难以实现有效求解. 本文围绕特征值解法器GCGE的研究进展, 系统开展了高效并行算法的设计与实现, 致力于提升其在复杂系统中的稳定性与适应性. 针对传统解法器在处理复杂小规模问题及大规模病态问题时存在的效率与稳定性瓶颈, 本文对GCG算法进行了优化改进, 引入对角归一化预处理策略以增强收敛性和数值稳定性, 并提出面向质量矩阵不正定问题的$L^2$-正交与刚度矩阵正交策略以提升数值性能. 此外, GCGE已完成与SLEPc软件环境的集成, 支持Hermitian特征值问题和最大特征值问题的求解, 并在北太天元高性能计算平台实现模块化部署, 展现出良好的工程应用前景.

整体代数多重网格(BAMG)法是一种求解偏微分方程组离散系统的高效预条件方法, 可用于求解多孔介质中的多相渗流模型. 本文针对BAMG方法, 面向多相多组分油气藏问题研究了块矩阵粗化中不同范数的选取及基于对角块的经典插值策略对收敛性的影响. 针对启动阶段的并行瓶颈, 引入了"延迟更新"和"结构保持更新"两种策略, 并基于可容忍迭代增长阈值的判据设计了一种自适应启动BAMG预条件方法. 进一步, 对求解阶段中的磨光算子与残量计算等热点进行了深度优化与向量化. 数值实验表明, 所提出方法在收敛性、高效性及并行可扩展性方面具有优势; 例如, 在亿级网格规模的测试中, 当使用8192个CPU核时, 提出方法的计算时间相比于传统方法减少了57.6%, 并行效率提升了38.1%.

线性方程组解法器库X-Solver旨在实现并优化Krylov子空间求解方法与预条件算法，用于在配备GPU加速卡的分布式内存机群上高效求解大规模稀疏线性方程组. 结合实际应用需求与硬件架构趋势，我们在以下三个方面实现了突破: (1) 提供用户友好、易于集成的分布式计算支持; (2) 深度适配异构架构特征，充分挖掘硬件计算潜力; (3) 保持跨多种硬件平台的高性能可移植性. 数值实验表明，在目标硬件上，X-Solver相较于其它先进解法器库在多个典型应用中展现出优势.

代数多重网格(AMG)是一种高效的线性方程组求解预条件方法.半结构 AMG 利用结构化信息高效计算, 且支持存在非结构信息, 因此可以同时达到高性能和高灵活性, 从而广泛应用于科学与工程计算的各个场景中.然而, 目前主流半结构 AMG 求解器在绝对速度和可扩展性上仍然具有明显缺陷, 为此我们研发了 Semi-StructMG 求解器. 一方面, 它利用多维粗化, 降低了复杂度, 提高了单步运行速度和可扩展性; 另一方面, 它在光滑器和插值算子中考虑块间连边, 改善了在各种复杂问题中的收敛性.我们在基准测试和多个真实应用中对 Semi-StructMG 进行了测试, 相比 hypre 中的 SSAMG, Split 和 BoomerAMG 达到了 5.97x, 15.2x 和 3.85x 的加速比.

三温辐射扩散方程组能够准确地描述辐射能量在介质中传播、散射、吸收和发射等物理过程.文[Yue X, He J, Xu X, Shu S, Wang L. Commun. Comput. Phys., 2022, 32: 829-849]基于一种物理量重叠分解策略提出了精确限制型加性和乘性Schwarz预条件算法和迭代算法,但未给出相关收敛性支撑和数值验证.本文将为它们所含的Schur补矩阵中子块逆矩阵提出合理的逼近假设,据此构造相应的非精确限制型加性和乘性Schwarz预条件算法,并受范数等价和数值域等价的数学启发,在保对称有限体元离散格式(或其系数矩阵)的稳定性条件下,开展相应预处理广义极小残差迭代方法的最优收敛性分析,且通过源自激光驱动球形内爆减速阶段中流体动力学不稳定实际模拟实验证实理论结果的有效性.

辐射扩散问题广泛出现在天体物理学和惯性约束聚变等多物理耦合领域, 基于问题的物理与代数特征的 AMG 法已成为当今多重网格法研究领域的热点.本文重点针对三温辐射扩散方程组的线性化离散系统, 首先给出了一种常见的 UA-AMG 预条件算法及相应的 PGMRES 解法器 T2T2-ILU(0)-V-FGMRES. 进一步, 为改善该解法器的计算性能, 对不同离散系统凝练了若干物理和代数特征, 设计了基于这些特征的自适应 UA-AMG 预条件算法, 并研制了相应的 PGMRES 解法器 Adapt-UA-AMG-FGMRES. 数值实验表明: 新解法器具有更好的稳健性和计算效率, 与解法器 T2T2-ILU(0)-V-FGMRES 和 HMIS-V-FGMRES(在基于几种常见非聚集型粗化算法的 AMG 预条件子中计算性能最好)相比 CPU 时间分别减少了约 49.1% 和 25.3%.上述算法设计思想容易推广到多群辐射扩散方程组等更一般的模型问题中.

在大规模线性系统求解中, 传统稀疏直接法解法器往往采用单一精度计算方式, 难以灵活平衡计算效率与数值精度. 为解决该问题, 基于分布式稀疏直接法解法器PanguLU提出了一种面向异构众核处理器MT-3000的混合精度优化算法. 该算法依据矩阵块的空间位置与数值敏感性, 动态选择块存储精度, 从而在数值分解阶段实现混合精度计算. 同时, 针对解法器中的通用矩阵乘子任务, 设计了一种计算与存储精度分离的流水线机制. 实验结果表明, 所提出的方法在数值分解阶段实现了1.04倍至1.19倍的性能提升, 同时将相对残差较单精度方案降低了1.97倍至4.15倍, 在提升求解速度的同时有效控制了精度损失.