湍流燃烧模拟中化学反应的加速算法研究进展

刘再刚, 孔文俊

刘再刚, 孔文俊. 湍流燃烧模拟中化学反应的加速算法研究进展[J]. 实验流体力学, 2019, 33(4): 1-10. DOI: 10.11729/syltlx20180123
引用本文: 刘再刚, 孔文俊. 湍流燃烧模拟中化学反应的加速算法研究进展[J]. 实验流体力学, 2019, 33(4): 1-10. DOI: 10.11729/syltlx20180123
Liu Zaigang, Kong Wenjun. Progress on acceleration algorithm of the computation for chemical reactions in turbulent combustion simulation[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(4): 1-10. DOI: 10.11729/syltlx20180123
Citation: Liu Zaigang, Kong Wenjun. Progress on acceleration algorithm of the computation for chemical reactions in turbulent combustion simulation[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(4): 1-10. DOI: 10.11729/syltlx20180123

湍流燃烧模拟中化学反应的加速算法研究进展

基金项目: 

国家自然科学基金项目 91441131

国家自然科学基金项目 U1738113

详细信息
    作者简介:

    刘再刚(1991-), 男, 河北石家庄人, 博士研究生。研究方向:燃烧化学反应机理简化加速计算。通信地址:北京市北四环西路11号中国科学院工程热物理研究所(100190)。E-mail:zaigangliu@foxmail.com

    通讯作者:

    孔文俊, E-mail: wjkong@buaa.edu.cn

  • 中图分类号: TK16

Progress on acceleration algorithm of the computation for chemical reactions in turbulent combustion simulation

  • 摘要: 为研究湍流燃烧数值模拟中化学反应机理计算的加速方法,讨论了动态自适应化学(Dynamic Adaptive Chemistry,DAC)方法和Krylov子空间近似的指数格式的应用情况。在湍流火焰大涡模拟中,使用DAC简化可以加速化学反应计算。然而,在并行燃烧数值模拟中,处理器核心的负载极度不平衡,加速效果有限。而Krylov子空间近似的指数格式的加速效果可以作用于每个处理器核心,更有利于整体计算效率的提高。在同等精度下,相比于隐式格式耦合DAC和MTS加速方法,Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著。
    Abstract: To investigate the methods of accelerating the computation of chemical kinetics in turbulent combustion, the application of Dynamic Adaptive Chemistry (DAC) and exponential integrator with Krylov subspace approximation is discussed. In the large eddy simulation of a turbulent flame, using DAC can accelerate the computation of chemical kinetics. However, in the context of parallel combustion simulation, the loads on the different processors are extremely imbalanced, which limits its performance. For the exponential integrator with Krylov subspace approximation, the acceleration effect acts on each processor, which is beneficial to improve the global computational efficiency. Compared to that of the implicit scheme coupled with DAC and MTS, the accelerating performance of the exponential integrator with Krylov subspace approxi-mation under the same level of accuracy is more obvious.
  • 燃烧是发动机中燃料能量转换的主要形式,准确地描述燃烧过程,从而得到精确的燃烧数值模拟结果, 可以高效指导发动机的设计与实验。实际应用中,发动机中的燃烧大多是湍流的,燃烧与湍流之间的相互作用机理复杂,使得湍流燃烧数值模拟成为了重要的研究课题。

    随着计算机能力的提升,大涡模拟(Large Eddy Simulation, LES)逐渐成为湍流燃烧数值模拟的重要研究手段。在湍流燃烧LES中需要使用燃烧模型来封闭标量方程中的化学反应源项。目前精度较高的燃烧模型,如小火焰模型[1-5]、线性涡模型[6-7]、基于概率密度函数的模型[8-12]、条件矩模型[13-14]、加厚火焰面模型[5, 15]、涡耗散模型[16]和部分预混反应器模型[17]等,均需要借助详细的化学反应机理来描述燃烧时的复杂反应过程,如燃烧场中的局部熄火和再燃、化学反应速率随压力的非线性变化、低温化学反应的影响、污染物排放等。

    化学反应机理是指描述化学过程的基元反应的集合。在燃烧时,燃料分解产生大量化学分子和自由基,称为组分,组分间最基本的化学反应为基元反应。要准确模拟燃烧过程,所需的化学反应机理的规模是巨大的,可以包含数量成百上千的组分和基元反应[18]。而且化学反应组分常微分方程组(Ordinary Differential Equations,ODEs)通常是刚性的,如果使用常规的显式格式求解,为了计算稳定必须使用极小的时间步长。常见的解决方法是使用基于隐式格式的刚性求解器,如DVODE[19]和DASSL[20]等。因此,在湍流燃烧LES中,求解化学反应方程的时间常占据主要部分[21]

    化学反应机理规模的增加将导致化学反应源项计算量的迅速增长。应用小火焰类燃烧模型,如FGM(Flamelet Generated Manifolds)或FPV(Flamelet/Progress-Variable)等,可以通过预先建表的方法回避在线求解化学反应问题,类似的还有PRISM(Piecewise Reusable Implementation of Solution Mapping)[22]、RS-HDMR(Random Sampling-High Dimensional Model Representation)[23-24]等映射方法。相比于直接使用总包反应机理,此类方法虽然在一定程度上考虑了详细化学反应机理的作用,但化学反应表中的变量数量严重受限于计算机内存空间,当燃烧的工况不处于预先存储结果的工况范围内时,就会缺乏灵活性和适应性,与直接在线求解化学反应相比准确性较低。为了降低直接求解化学反应的计算时间,虽然可以使用在线的建表方法,如ISAT(In Situ Adaptive Tabulation)[25-29],但是,ISAT方法对内存需求仍较大,并且在线建表的速度快慢也依赖于直接积分的求解速度。而提高直接求解的计算效率,可以通过使用化学反应机理简化或采用高效的积分求解器来实现。

    对于不同的燃烧问题,反应机理的组分和基元反应并不都是活跃的,可能仅仅其中的某一部分是重要的,而其他部分的反应速率则可以忽略不计。因此,研究人员发展了许多化学反应机理简化方法,如敏感性分析法[30-32]、准稳态假设法[33]、局部平衡法[33]、计算奇异摄动法[34]、ILDM(Intrinsic Low Dimensional Manifold)法[35]和DRG(Direct Relation Graph)法[36]及其衍生算法,如DRGEP(Directed Relation Graph Error Propagation)方法[37]、PFA(Path Flux Analysis)[38-40]、EFA(Element Flux Analysis)[41-43]和GPS(Global Pathway Selection)方法[44]等。

    上述反应机理简化方法需要在流场计算前完成,并在计算过程中直接使用预先简化好的反应机理。为了保证简化机理的准确性,这种预先简化的机理常包含了较多的组分。然而,对于同一燃烧场中的不同时刻或当地,反应机理可以进一步简化。因此研究人员发展了一类动态自适应化学法(Dynamic Adaptive Chemistry,DAC)[39, 45-48],可以在计算过程中动态简化化学反应机理。动态简化主要基于DRG、PFA或GPS算法,因为这些方法能以较小的额外时间花费完成机理简化。DAC方法已应用于多种简单模型的燃烧计算,如HCCI(Homogeneous Charge Compression Igni-tion)模型[47, 49]、部分预混搅拌反应器模型[27, 50-52]、自着火问题和一维火焰传播[48, 53-54]等。

    DAC方法也在不断地改进,并与其他简化方法耦合应用。为了进一步降低机理简化本身所占用的时间,Sun等[55]提出了自相关动态自适应化学(Correlated Dynamic Adaptive Chemistry,CoDAC)方法。该方法将具有相近状态的不同空间位置和时间节点网格点进行归类后统一处理,这与分区模型和CA(Cell Agglomeration)方法[56-57]的思想相同,自相关技术有效地降低了化学反应机理简化所需要的计算时间。Oluwole等[54]发展了一种DSRR(Decoupled Species and Reaction Reduction)简化方法,相比于DRG简化,具有同等的简化计算速度并提供了改进的误差控制方法。DSRR的显著优点是不需要指定初始组分,但简化结果可能会导致元素不守恒[54]。Xie等[58]提出的TSRA(Time-Scale and Jacobian-aided Rate Analysis)将组分分为活跃的(Active)、直接耦合的(Directly coupled)和不重要的(Inconsequential)3类,通过迭代计算来保证元素守恒。此外,DAC方法和ISAT方法结合,构成TDAC(Tabulated Dynamic Adaptive Chemistry)[27, 51, 59-60]方法,可以进一步提高化学反应计算速度。除了将DAC和ISAT二者结合,Xie等[60]提出的DAAM(Dynamic Adaptive Acceleration Method)可以根据DAC和ISAT的计算耗时自动选择二者之中较快的方法进行计算。最后,分区方法或CA方法也可以和DAC及ISAT结合[52, 61-62],实现更高效的计算。

    通过提高求解器计算效率也可以加速化学反应求解,如使用外推法[63],点隐式求解器[64]、带预处理的Krylov求解器[65]、刚性去除法[66]和基于多时间尺度(MTS,Multi-Timescale Method)分析的混合多时间尺度(HMTS,Hybrid Multi-Timescale Method)法[67],HMM(Heterogeneous Multiscale Method)方法[68]和AHI(Adaptive Method for Hybrid Integra-tion)[69]方法等。

    求解化学反应常微分方程组,除了使用广泛应用的显式或隐式格式外,还可以使用指数格式积分方法(Exponential Integrator)[70]。这种方法构造含有Jacobian矩阵指数函数形式的积分格式,其优势在于可以使用显式大步长时间推进得到稳定的计算结果,适用于刚性常微分方程问题和振荡系统问题。因此,指数格式被广泛应用于金融与统计学[71]、正则化理论[72]、电磁场模拟[73]及浅水方程模拟[74]等领域。在指数积分方法中,需要计算Jacobian矩阵的指数函数。在多种计算矩阵指数函数的方法[75]中通用性较好、计算效率较高的方法是收缩乘方和Padé近似[76]方法。然而,这种方法所需的计算量与矩阵维数的3次方成正比,因此适用于较小规模的矩阵,对于大型矩阵,尤其是大型稀疏矩阵,其计算效率很低。一种广泛采用的解决方法是利用Krylov子空间近似[77],将原问题投影到其Krylov子空间,从而实现降维,减少问题规模[74, 78-81]。Krylov子空间近似也常被应用于GMRES (Generalized Minimum Residual)算法和BiCGSTAB (Biconjugate Gradient Stabilized)算法等线性方程组求解方法中。

    在燃烧数值模拟中,刚性化学反应问题也适合用指数格式求解。Bisetti[82]和Curtis等[83]的研究表明, 四阶指数积分格式和Krylov子空间近似方法可应用于求解自着火和HCCI的化学反应机理。然而,四阶指数积分格式每步需要计算3次矩阵指数函数,计算量较大。Niesen和Wright[79]发展了一种自适应时间步长和Krylov子空间大小的算法:Phipm算法。该求解器已被应用于自着火问题数值模拟,与DVODE求解器相比,显著提高了化学反应求解效率[84]

    为了研究湍流燃烧数值模拟中化学反应机理计算的加速方法,本文阐述了DAC和指数格式在湍流燃烧数值模拟中的应用情况。首先,探讨了DAC在湍流燃烧数值模拟中的加速效果。然后,考察了Krylov子空间近似的指数格式求解湍流-化学反应系统的加速效果。

    Liu等[21, 85]使用53组分/325反应的GRI-Mech 3.0反应机理和基于PFA的CoDAC简化方法对Sandia Flame D进行了LES数值模拟。Sandia Flame D是美国Sandia国家实验室发展的经典湍流射流火焰,由主射流、值班射流和空气伴流组成,其中主射流由25%甲烷和75%空气组成,值班射流为甲烷-空气燃烧后的平衡产物。具体参数可见Barlow等[86]的研究。数值模拟结果如图 1所示,通过对比发现,使用CoDAC简化的LES计算结果与未使用简化得到的结果符合得很好,表明使用CoDAC可以准确模拟非预混湍流火焰特性。使用CoDAC后,LES中化学反应计算时间总体下降了29%,化学反应速率的求解仍然占据总求解时间的主要部分。Yang等[64, 87-89]同样使用基于PFA的CoDAC方法简化化学反应计算,并结合CoTran方法降低组分输运系数的计算。对湍流预混火焰的DNS模拟[64]结果表明, 简化方法的加速因子可以达到2.7,在湍流非预混火焰的DNS模拟[89]中可以达到4.0。CoDAC简化方法也被用于模拟实验室尺度湍流火焰Sandia Flame D[87]和Flame E[88]的LES,使用的反应机理为20组分/84反应的简化机理,该反应机理是通过对GRI-Mech3.0机理使用GPS简化后获得的。该研究对比了有限速率化学模型和小火焰/进度变量方法,结果表明,2种模型都较好地反映了实验结果的变化趋势,但使用有限速率化学模型和CoDAC简化相比于FPV模型,可以得到更准确的统计结果。

    图  1  Sandia Flame D火焰在x/d=3、7.5、15、30、45和60截面处时均温度的周向分布[21]
    Fig.  1  Radial profiles of time-average of temperature for Sandia Flame D at x/d = 3, 7.5, 15, 30, 45 and 60[21]

    DAC还被应用于超声速燃烧数值模拟中。Wu等[90]对比了DAC、TDAC方法和骨架机理在超声速燃烧室数值模拟中的表现。结果表明DAC和TDAC可以准确模拟燃烧室的特性和火焰结构,而骨架机理不能准确预测火焰的稳定特性。在计算效率方面,在57组分/269反应的UCSD反应机理的基础上使用DAC可以获得与使用30组分/143反应的骨架机理同等的加速效果,而在DAC的基础上使用TDAC可以将计算效率再提高1倍。

    在实际燃烧设备的数值模拟中,DAC方法也得到了应用,Li等[59]对比研究了TDAC方法和5种不同的化学反应机理简化方法:DRG、DRGEP、Contino等的DAC方法[91]、PFA和EFA。研究使用非定常雷诺平均方法(URANS)数值模拟了DJHC(Delft Jet-in-Hot-Coflow)燃烧器的燃烧特性,结果表明DRGEP、DAC和EFA方法比DRG和PFA方法准确。对于小规模反应机理,TDAC方法中的制表方法贡献了主要的计算加速,而对于大规模反应机理,DAC简化起了主要作用。Zhou和Lu等[92-93]使用URANS和LES方法模拟了Spray H正庚烷湍流喷雾火焰,结果表明使用DAC的加速因子可达到2.0。在化工应用方面,Huang等[94]使用基于DRG的DAC方法模拟了石英玻璃反应炉的SiCl4/H2/O2燃烧系统,加速因子为2.2。

    DAC在内燃机数值模拟中受到了广泛关注,DAC方法及其他多种化学反应机理简化方法被应用于加速计算。Shi等[95]使用基于DRGEP的DAC进行了DI(Direct-Injection)发动机的二维数值模拟,并改进了确定初始组分的方法。计算速度的对比表明,在多维发动机数值模拟中的DAC简化加速效果不如在单区绝热发动机模型模拟中显著,但仍有效加速了化学反应计算。对于34组分和61组分的正庚烷反应机理,计算时间分别节省了30%和50%。Contino等[91]使用TDAC模拟了HCCI发动机的燃烧过程,加速因子可高达500.0,而对于柴油发动机,加速因子在9.0左右。Xie等[60]将DAAM方法用于PCCI(Premixed Charge Compression Ignition)发动机的数值模拟,模拟结果显示DAAM相比于ISAT方法使计算加速了50%,其计算效率也高于ISAT-DAC方法,且数值模拟的准确性处于同一水平。此外,DAC和ISAT方法也可以和CA方法联合[61-62],在分区的基础上再进行DAC简化或ISAT建表,可以进步一提高计算速度。Zhou等[61]的研究表明,对于约40个组分的反应机理,CA联合DAC的方法可以使计算加速一倍,而对于约140个组分的反应机理,加速因子可达到9.0左右。

    DAC在湍流燃烧数值模拟中的加速效果不尽相同,其影响因素是多方面的:

    (1) DAC简化的加速效果与数值模拟结果误差是相关的,用户所能容忍的误差越大,DAC简化能实现的加速就越大。误差的调节方式与DAC中具体应用的简化方法有关,如DRG、DRGEP、PFA等通过简化阈值调节,合理阈值的大小还根据所使用简化机理的不同而不同,最终的误差需要通过结果来验证;而DSRR则通过给定容许误差来调节。

    (2) 加速效果与燃烧特性有关。Liu等[21]的研究表明,对于燃烧反应不活跃的情况,如自着火过程发生前的慢速反应区域和自着火完成后的平衡区域,以及一维预混火焰面上游的未燃区和火焰面下游的已燃区,燃烧过程中活跃的组分数较小,使用DAC简化可以显著降低局部简化机理中的组分,从而大大提高计算速度。而对于反应较为剧烈的区域(如自着火时刻或火焰锋面处),化学反应过程复杂,涉及组分多,因此DAC简化产生的局部简化机理规模较大,甚至与输入的反应机理规模相当,此时DAC的加速效果不明显。

    (3) 加速效果与输入反应机理有关,应用规模越大的反应机理通常能获得越明显的加速效果。例如Wu等[90]的研究表明,同等网格数量的情况下,对于57组分、75组分和111组分的反应机理,DAC的加速因子分别为1.9、2.2和2.5。

    在实际应用中,DAC如果与其他加速手段耦合,加速效果也会受到影响。例如DAC方法与ISAT方法耦合时[51],DAC的加速效果随ISAT误差阈值的增大而略有降低。在湍流燃烧的并行计算中,化学反应计算的并行方式也对DAC的加速效果有较大影响。Liu等[21]的研究表明,由于并行LES中化学反应的计算效率受火焰面附近处理器核心的计算效率控制,在并行计算中,化学反应求解完成后,核与核之间需要进行一次同步并通信,然后再进行下一步计算,因此总体的化学反应计算时间由各个核心中最慢的决定。例如,图 2所示是不同核心上的化学反应计算时间tODE, p,可以看到不同计算区域中的tODE, p差别很大,火焰反应区附近的tODE, p显著高于其他反应活性较低的区域。放大图 2(a), 显示出最高的tODE, p出现在火焰中r=d, 10d < x < 25d的区域。这正是湍流-化学相互作用最强烈的区域,湍流火焰中的局部熄火和再燃频繁发生。CoDAC和DAC的主要区别在于是否进行自相关算法运算。图 2(d)所示是当地各个网格点上的简化机理中通过PFA算法得到而不是通过自相关算法得到的占比ϕPFA, p,分布结果同样表明在火焰反应区附近自相关算法带来的计算加速效果不明显。在自着火的模拟结果中(如图 3所示),同样可以发现在反应活跃的点火时刻,计算时间也出现了显著升高的现象。因此,要进一步提高计算效率,一方面应平衡处理器之间的计算效率,如使用x2f_mpi算法[96]等;另一方面,可以发展更高效的化学反应积分方法,从而提高整体计算效率。

    图  2  在各处理器核心上的(a) ODE时间(tODE, p)以及(b) ODE时间、(c)温度T、(d) PFA占比(ϕPFA, p)和(e)简化机理中组分数(NS)的局部放大视图[21]
    Fig.  2  Contour plot of (a) ODE time on each processor (tODE, p) together with zoom-in view of (b) ODE time on each processor, (c) temperature T, (d) percentage of cells performing PFA (ϕPFA, p) and (e) number of selected species in locally reduced mechanism (NS) [21]
    图  3  使用DAC和CoDAC在初始条件为ϕ =1.0, T0=1300K和简化阈值εr=0.1时相对计算时间和温度随时间的变化[21]
    Fig.  3  Evolution of relative computational time tr and temperature for auto-ignition simulation with DAC and CoDAC at ϕ=1.0, T0=1300K and reduction threshold εr=0.1[21]

    对于如下的燃烧化学反应常微分方程问题:

    (1)

    二阶的指数积分格式可以表示为[70]

    (2)

    式中ϕ表示包括温度和组分质量分数在内的微分方程变量,ωk表示各个变量的化学反应速率,h为时间步长,当时间步长较短时可认为Jacobian矩阵Jk=∂ωk/∂ϕk是近似恒定的,矩阵的φ函数定义为:

    (3)

    此时,常微分方程的求解变成了对一般形式为φ1(τA)v的向量的求解,其中A为在n×n维实空间Rn×n中的矩阵,v为在n维实空间Rn中的向量。

    n很大时,φ函数的计算需要巨大的计算量,利用Krylov子空间的性质,组合φ1(τA)v可以方便地在更小的m维Krylov子空间中得到近似。Krylov子空间是由向量序列v, Av, A2v, …, Am-1v张成的空间,构造Krylov子空间的算法称为Arnoldi算法[77]。原向量φ1(τA)v在Krylov子空间中的近似可以表示为[77]

    (4)

    式中β=‖v‖,Vm=(v1, v2, …, vm)为Krylov子空间的基向量,Hm为上Hessenberg矩阵,是原矩阵A在Krylov子空间中的投影,e1em分别表示第1维和第m维单位向量。当矩阵A在Krylov子空间中投影为Hm时,Hm的特征值称为Ritz值,随着m的增加,收敛于原矩阵A模最大的特征值。这一重要性质使Krylov子空间广泛应用于近似求解矩阵特征值问题[97]。通过近似式(4),原n维问题降为m维问题,使指数格式求解所需的计算量显著降低。

    Bisetti[82]使用四阶Rosenbrock指数格式exp4[98]计算了火焰的自着火过程,结果表明, 尽管使用大时间步长时非线性系统刚性很大,但指数格式方法仍是稳定的。Curtis等[83]对比了CPU和GPU计算时exp4和expbr43指数积分格式的计算效率,结果表明2种指数积分格式在CPU和GPU上都慢于隐式求解器CVODE和Radau-IIA。然而,四阶指数积分格式每推进1个时间步长需要计算3次矩阵指数函数。而矩阵指数函数的计算十分耗时,为了实现加速计算,应尽量减少计算次数。因此,刘再刚等[84]使用二阶指数格式,如式(2),并引入Schur分解方法有效控制了算法的舍入误差[99],发展了用于加速燃烧化学反应计算的指数积分求解器。

    指数格式算法的计算效率可以使用简单的燃烧模型来估计。带配对混合(Pairwise mixing)的局部搅拌反应器(Partially-Stirred Reactor, PaSR)是一个假设的反应器[100]。这种反应器可以用于模拟湍流燃烧数值模拟中的小尺度混合和反应过程,常用于检验燃烧化学反应计算方法和ISAT方法的准确性和计算效率。模拟使用的燃料和参数与Liu等[21]的研究相同,PaSR入流由3部分组成,分别为甲烷、甲烷-空气燃烧后的平衡产物以及空气,并分别对应Sandia Flame系列火焰的主射流、值班射流和空气伴流。在Liu等[21]的研究中,CoDAC在PaSR模拟中得到了与LES模拟相近的加速结果,表明了PaSR能合理估计LES的计算效率。对于燃烧化学反应问题,还可以同时使用一些加速算法。CoDAC简化应用于指数格式可以减少算法中矩阵Av的维数,SpeedCHEM(SC)库[101]中的稀疏矩阵算法和解析Jacobian矩阵算法可以降低矩阵存储所需的内存空间,并加速Jacobian矩阵和Arnoldi算法中矩阵-向量乘法Av的计算。MTS方法[69]可以分离刚性问题中的刚性部分和非刚性部分,从而降低刚性方程组的数量。在指数格式中,作用与CoDAC相同,可以减少算法中矩阵Av的维数。Liu等[99]将指数积分格式与SC库、CoDAC简化和MTS方法结合,并与DVODE求解器的计算结果作了详细对比,不同的组合算例列于表 1

    表  1  使用不同加速方法的算例设置
    Table  1  Simulated cases for different acceleration methods
    加速方法\求解器 DVODE 指数格式
    Non R -
    SC As Bs
    SC + CoDAC Asc Bsc
    SC+MTS Asm Bsm
    SC + CoDAC + MTS Ascm Bscm
    下载: 导出CSV 
    | 显示表格

    图 4所示是分别使用GRI-Mech 3.0[102]和USC Mech Ⅱ[103](111组分/784反应)计算得到的结果。图中的加速因子和误差均以表 1中的参考工况R为基准,加速因子和误差的详细定义可参考文献[99]。从图中可以看出,在同等精度并同时采用了SC的情况下,使用指数格式的计算速度均显著快于使用隐式格式求解器DVODE,并且加速效果显著优于使用DVODE结合3种加速算法。但指数格式结合CoDAC和MTS加速算法时,不仅计算精度下降,计算效率也有所下降。从图 5可以看出,效率下降主要是由CoDAC和MTS带来的额外开销导致的。因此只使用含SC的指数格式可以得到最佳的加速效果。

    图  4  使用反应机理(a)GRI-Mech 3.0[102]和(b)USC Mech Ⅱ[103]得到的表 1中算例的加速因子γ与最大平均相对误差ε的关系。Xsc(εr),Xsm(α)和Xscm(εr, α)表示表 1所列工况算例,括号中的数值εrα分别表示CoDAC简化阈值和MTS安全因子[99]
    Fig.  4  The speedup factor γ as functions of maximum averaged relative error ε for the cases in Table 1 with (a) GRI-Mech 3.0[102] and (b) USC Mech Ⅱ [103]. Xsc (εr), Xsm(α) and Xscm(εr, α) represent the cases in Table 1 and numbers εr and α in the parentheses represent the reduction threshold of CoDAC and the safety factor of MTS, respectively[99]
    图  5  使用USC Mech Ⅱ[103]时,表 1中的算例中用于计算快方程、慢方程和CoDAC简化的壁面时间[99]
    Fig.  5  Wall time spent on fast equation integration, slow equation integration and CoDAC reduction for the cases in Table 1 with USC Mech Ⅱ [99, 103]

    指数格式的加速效果一方面来源于大步长显式时间推进的优势,另一方面是由于与Krylov子空间近似的降维效果。从图 6中可以看到,Krylov子空间的大小增长速度远低于刚性方程组规模的增速。由此可知,在化学反应计算中应用指数格式在保证计算精度的同时可以显著提高计算速度。因此,其在燃烧LES中有很好的应用前景。指数格式求解器直接求解化学反应方程组,如果与查表类加速方法,如ISAT相结合,可以进一步降低化学反应计算时间,从而实现快速高精度的湍流燃烧数值模拟,从而使得详细化学反应机理用于实际发动机燃烧室湍流燃烧数值模拟成为可能,并以此指导湍流燃烧研究和燃烧器设计与实验。

    图  6  使用GRI-Mech 3.0[102]和USC Mech Ⅱ [103]机理时表 1中算例Ascm和Bscm的Krylov子空间维数(mKrylov)随快方程数量(nfast)的变化[99]
    Fig.  6  The dimension of the Krylov subspace (mKrylov as a function of the number of fast equations (nfast) for cases Ascm and Bscm shown in Table 1 with the mechanisms of GRI-Mech 3.0 [102] and USC Mech Ⅱ [99, 103]

    在湍流燃烧数值模拟中使用详细化学反应机理可以显著提高燃烧反应计算精度,但是由此导致的计算量激增是阻碍大规模高精度燃烧数值模拟的主要因素之一。在计算过程中动态、自适应地进行化学反应机理简化可以获得一定的加速效果。然而,在并行燃烧数值模拟中,由于燃烧场中火焰反应区域的分布不均匀,处理器核心的负载极度不平衡,且核心之间需要进行必要的同步通信,加速效果在很大程度上不能充分发挥。相比之下,使用刚性求解器得到的加速将直接作用于每个处理器核心,更有利于整体计算效率提高。相比于常见的隐式格式以及CoDAC和MTS加速方法,在同等精度下,Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著,加速主要来源于指数格式的显式大时间步推进能力和Krylov子空间近似的降维效果。

    为了进一步提高化学反应的求解速度,首先,可行的方法是以高效积分求解器为基础,紧密结合查表类方法,互相弥补2种方法的不足:用查表代替重复性计算,用高效积分补充表中的元素。其次,使用自适应网格以及自适应平衡处理器核心负载的算法,可以有效避免计算资源的空闲,提高整体计算效率。最后,燃烧化学反应的求解涉及的浮点数运算密集,较适合使用图形处理器架构进行运算,这也是燃烧数值模拟加速的一个发展方向。

    致谢: 感谢国家自然科学基金项目(91441131和U1738113)的支持。
  • 图  1   Sandia Flame D火焰在x/d=3、7.5、15、30、45和60截面处时均温度的周向分布[21]

    Fig.  1   Radial profiles of time-average of temperature for Sandia Flame D at x/d = 3, 7.5, 15, 30, 45 and 60[21]

    图  2   在各处理器核心上的(a) ODE时间(tODE, p)以及(b) ODE时间、(c)温度T、(d) PFA占比(ϕPFA, p)和(e)简化机理中组分数(NS)的局部放大视图[21]

    Fig.  2   Contour plot of (a) ODE time on each processor (tODE, p) together with zoom-in view of (b) ODE time on each processor, (c) temperature T, (d) percentage of cells performing PFA (ϕPFA, p) and (e) number of selected species in locally reduced mechanism (NS) [21]

    图  3   使用DAC和CoDAC在初始条件为ϕ =1.0, T0=1300K和简化阈值εr=0.1时相对计算时间和温度随时间的变化[21]

    Fig.  3   Evolution of relative computational time tr and temperature for auto-ignition simulation with DAC and CoDAC at ϕ=1.0, T0=1300K and reduction threshold εr=0.1[21]

    图  4   使用反应机理(a)GRI-Mech 3.0[102]和(b)USC Mech Ⅱ[103]得到的表 1中算例的加速因子γ与最大平均相对误差ε的关系。Xsc(εr),Xsm(α)和Xscm(εr, α)表示表 1所列工况算例,括号中的数值εrα分别表示CoDAC简化阈值和MTS安全因子[99]

    Fig.  4   The speedup factor γ as functions of maximum averaged relative error ε for the cases in Table 1 with (a) GRI-Mech 3.0[102] and (b) USC Mech Ⅱ [103]. Xsc (εr), Xsm(α) and Xscm(εr, α) represent the cases in Table 1 and numbers εr and α in the parentheses represent the reduction threshold of CoDAC and the safety factor of MTS, respectively[99]

    图  5   使用USC Mech Ⅱ[103]时,表 1中的算例中用于计算快方程、慢方程和CoDAC简化的壁面时间[99]

    Fig.  5   Wall time spent on fast equation integration, slow equation integration and CoDAC reduction for the cases in Table 1 with USC Mech Ⅱ [99, 103]

    图  6   使用GRI-Mech 3.0[102]和USC Mech Ⅱ [103]机理时表 1中算例Ascm和Bscm的Krylov子空间维数(mKrylov)随快方程数量(nfast)的变化[99]

    Fig.  6   The dimension of the Krylov subspace (mKrylov as a function of the number of fast equations (nfast) for cases Ascm and Bscm shown in Table 1 with the mechanisms of GRI-Mech 3.0 [102] and USC Mech Ⅱ [99, 103]

    表  1   使用不同加速方法的算例设置

    Table  1   Simulated cases for different acceleration methods

    加速方法\求解器 DVODE 指数格式
    Non R -
    SC As Bs
    SC + CoDAC Asc Bsc
    SC+MTS Asm Bsm
    SC + CoDAC + MTS Ascm Bscm
    下载: 导出CSV
  • [1]

    Ihme M, Pitsch H. Prediction of extinction and reignition in nonpremixed turbulent flames using a flamelet/progress variable model:2. Application in LES of Sandia flames D and E[J]. Combustion and Flame, 2008, 155(1-2):90-107. DOI: 10.1016/j.combustflame.2008.04.015

    [2]

    Benim A C, Iqbal S, Meier W, et al. Numerical investigation of turbulent swirling flames with validation in a gas turbine model combustor[J]. Applied Thermal Engineering, 2017, 110:202-212. DOI: 10.1016/j.applthermaleng.2016.08.143

    [3]

    Olbricht C, Stein O T, Janicka J, et al. LES of lifted flames in a gas turbine model combustor using top-hat filtered PFGM chemistry[J]. Fuel, 2012, 96(1):100-107. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1376d56355d324a1e57c710827c8c249

    [4]

    Donini A, Bastiaans R J M, van Oijen J A, et al. A 5-D implementation of FGM for the large eddy simulation of a stratified swirled flame with heat loss in a gas turbine combustor[J]. Flow Turbulence and Combustion, 2017, 98(3):887-922. DOI: 10.1007/s10494-016-9777-7

    [5]

    Zhang H, Ye T, Wang G, et al. Large eddy simulation of turbulent premixed swirling flames using dynamic thickened flame with tabulated detailed chemistry[J]. Flow Turbulence and Combustion, 2017, 98(3):841-885. DOI: 10.1007/s10494-016-9791-9

    [6]

    Zhou X J, Jiang X, Martinez D M. The effects of chemical kinetic mechanisms on large eddy simulation (LES) of a nonpremixed hydrogen jet flame[J]. International Journal of Hydrogen Energy, 2016, 41(26):11427-11440. DOI: 10.1016/j.ijhydene.2016.04.079

    [7]

    Srinivasan S, Menon S. Linear eddy mixing model studies of high karlovitz number turbulent premixed flames[J]. Flow Turbulence and Combustion, 2014, 93(2):189-219. DOI: 10.1007/s10494-014-9542-8

    [8]

    Jangi M, Zhao X Y, Haworth D C, et al. Stabilization and liftoff length of a non-premixed methane/air jet flame dischar-ging into a high-temperature environment:An accelerated transported PDF method[J]. Combustion and Flame, 2015, 162(2):408-419. https://www.sciencedirect.com/science/article/abs/pii/S0010218014002375

    [9]

    Elbahloul S, Rigopoulos S. Rate-Controlled Constrained Equilibrium (RCCE) simulations of turbulent partially premixed flames (Sandia D/E/F) and comparison with detailed chemistry[J]. Combustion and Flame, 2015, 162(5):2256-2271. DOI: 10.1016/j.combustflame.2015.01.023

    [10]

    Raman V, Pitsch H. A consistent LES/filtered-density function formulation for the simulation of turbulent flames with detailed chemistry[J]. Proceedings of the Combustion Institute, 2007, 31:1711-1719. DOI: 10.1016/j.proci.2006.07.152

    [11]

    Popov P P, Pope S B. Large eddy simulation/probability density function simulations of bluff body stabilized flames[J]. Combustion and Flame, 2014, 161(12):3100-3133. DOI: 10.1016/j.combustflame.2014.05.018

    [12]

    Hiremath V, Lantz S R, Wang H F, et al. Computationally-efficient and scalable parallel implementation of chemistry in simulations of turbulent combustion[J]. Combustion and Flame, 2012, 159(10):3096-3109. DOI: 10.1016/j.combustflame.2012.04.013

    [13]

    Zhang H W, Garmory A, Cavaliere D E, et al. Large Eddy Simulation/Conditional Moment Closure modeling of swirl-stabilized non-premixed flames with local extinction[J]. Proceedings of the Combustion Institute, 2015, 35:1167-1174. DOI: 10.1016/j.proci.2014.05.052

    [14]

    Zhang H, Mastorakos E. Modelling local extinction in Sydney swirling non-premixed flames with LES/CMC[J]. Proceedings of the Combustion Institute, 2017, 36(2):1669-1676. https://www.sciencedirect.com/science/article/pii/S1540748916303091

    [15]

    Rittler A, Proch F, Kempf A M. LES of the Sydney piloted spray flame series with the PFGM/ATF approach and different sub-filter models[J]. Combustion and Flame, 2015, 162(4):1575-1598. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=6dea79083482522efe8ae0414b4698ca

    [16]

    Minotti A, Bruno C, Cozzi F. A LES simulation of a CH4/air microcombustor with detailed chemistry[J]. Combustion Science and Technology, 2011, 183(6):554-574. DOI: 10.1080/00102202.2010.523031

    [17]

    Fooladgar E, Chan C K, Nogenmyr K J. An accelerated computation of combustion with finite-rate chemistry using LES and an open source library for In-Situ-Adaptive Tabulation[J]. Computers & Fluids, 2017, 146:42-50. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=d7c627c35e600c9655ed9526e19782e9

    [18]

    Lu T, Law C K. Toward accommodating realistic fuel chemistry in large-scale computations[J]. Progress in Energy and Combustion Science, 2009, 35(2):192-215. DOI: 10.1016/j.pecs.2008.10.002

    [19]

    Brown P N, Byrne G D, Hindmarsh A C. VODE, a variable-coefficient ODE solver[J]. Siam Journal on Scientific and Statistical Computing, 1989, 10(5):1038-1051. DOI: 10.1137/0910062

    [20]

    Petzold L R. A description of DASSL: A differential/algebraic system sover[R]. Sandia National Laboratories, SAND82-8637, 1982.

    [21]

    Liu Z, Han W, Kong W, et al. LES modelling of turbulent non-premixed jet flames with correlated dynamic adaptive chemistry[J]. Combustion Theory and Modelling, 2018, 22(4):694-713. DOI: 10.1080/13647830.2018.1447148

    [22]

    Tonse S R, Moriarty N W, Frenklach M, et al. Computational economy improvements in PRISM[J]. International Journal of Chemical Kinetics, 2003, 35(9):438-452. DOI: 10.1002/kin.10140

    [23]

    Li G Y, Hu J S, Wang S W, et al. Random sampling-high dimensional model representation (RS-HDMR) and orthogona-lity of its different order component functions[J]. Journal of Physical Chemistry A, 2006, 110(7):2474-2485. DOI: 10.1021/jp054148m

    [24]

    Li G, Rabitz H, Hu J, et al. Regularized random-sampling high dimensional model representation (RS-HDMR)[J]. Journal of Mathematical Chemistry, 2008, 43(3):1207-1232. DOI: 10.1007/s10910-007-9250-x

    [25]

    Pope S B. Computationally efficient implementation of combustion chemistry using in situ adaptive tabulation[J]. Combustion Theory and Modelling, 1997, 1(1):41-63. DOI: 10.1080/713665229

    [26]

    Yang B, Pope S B. Treating chemistry in combustion with detailed mechanisms-In situ adaptive tabulation in principal directions premixed combustion[J]. Combustion and Flame, 1998, 112(1-2):85-112. DOI: 10.1016/S0010-2180(97)81759-2

    [27]

    Contino F, Jeanmart H, Lucchini T, et al. Coupling of in situ adaptive tabulation and dynamic adaptive chemistry:An effective method for solving combustion in engine simulations[J]. Proceedings of the Combustion Institute, 2011, 33(2):3057-3064. https://www.sciencedirect.com/science/article/pii/S1540748910003627

    [28]

    Sun W, Won S H, Ju Y. In situ plasma activated low temperature chemistry and the S-curve transition in DME/oxygen/helium mixture[J]. Combustion and Flame, 2014, 161(8):2054-2063. DOI: 10.1016/j.combustflame.2014.01.028

    [29]

    Lefkowitz J K, Uddi M, Windom B C, et al. In situ species diagnostics and kinetic study of plasma activated ethylene dissociation and oxidation in a low temperature flow reactor[J]. Proceedings of the Combustion Institute, 2015, 35(3):3505-3512. https://www.sciencedirect.com/science/article/pii/S154074891400385X

    [30]

    Peters N, Kee R J. The computation of stretched laminar methane air diffusion flames using a reduced 4-step mechanism[J]. Combustion and Flame, 1987, 68(1):17-29. https://www.sciencedirect.com/science/article/abs/pii/0010218087900629

    [31]

    Chen J Y. A general procedure for constructing reduced reaction mechanisms with given independent relations[J]. Combustion Science and Technology, 1988, 57(1-3):89-94. DOI: 10.1080/00102208808923945

    [32]

    Ju Y, Niioka T. Reduced kinetic mechanism of ignition for nonpremixed hydrogen/air in a supersonic mixing layer[J]. Combustion and Flame, 1994, 99(2):240-246. https://www.sciencedirect.com/science/article/abs/pii/0010218094901279?via=ihub&cc=y

    [33]

    Law C K. Combustion physics[M]. New York:Cambridge University Press, 2006.

    [34]

    Lam S H, Goussis D A. The CSP method for simplifying kinetics[J]. International Journal of Chemical Kinetics, 1994, 26(4):461-486. DOI: 10.1002/kin.550260408

    [35]

    Maas U, Pope S B. Simplifying chemical kinetics:Intrinsic low-dimensional manifolds in composition space[J]. Combustion and Flame, 1992, 88(3):239-264. DOI: 10.1016-0010-2180(92)90034-M/

    [36]

    Lu T F, Law C K. A directed relation graph method for mechanism reduction[J]. Proceedings of the Combustion Institute, 2005, 30:1333-1341. DOI: 10.1016/j.proci.2004.08.145

    [37]

    Lu T F, Law C K. On the applicability of directed relation graphs to the reduction of reaction mechanisms[J]. Combustion and Flame, 2006, 146(3):472-483. DOI: 10.1016/j.combustflame.2006.04.017

    [38]

    Sun W T, Chen Z, Gou X L, et al. A path flux analysis method for the reduction of detailed chemical kinetic mechanisms[J]. Combustion and Flame, 2010, 157(7):1298-1307. DOI: 10.1016/j.combustflame.2010.03.006

    [39]

    Gou X L, Chen Z, Sun W T, et al. A dynamic adaptive chemistry scheme with error control for combustion modeling with a large detailed mechanism[J]. Combustion and Flame, 2013, 160(2):225-231. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f3ec5c8ad533b33eeafa960c21541644

    [40]

    Wang W, Gou X L. An improved path flux analysis with multi generations method for mechanism reduction[J]. Combustion Theory and Modelling, 2016, 20(2):203-220. DOI: 10.1080/13647830.2015.1117660

    [41]

    He K, Androulakis I P, Ierapetritou M G. On-the-fly reduction of kinetic mechanisms using element flux analysis[J]. Chemical Engineering Science, 2010, 65(3):1173-1184. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=66f189777c234022e451d3fc1d5607c1

    [42]

    Androulakis I P, Grenda J M, Bozzelli J W. Time-integrated pointers for enabling the analysis of detailed reaction mechanisms[J]. Aiche Journal, 2004, 50(11):2956-2970. DOI: 10.1002/aic.10263

    [43]

    Zhang S, Broadbelt L J, Androulakis I P, et al. Comparison of biodiesel performance based on HCCI engine simulation using detailed mechanism with on-the-fly reduction[J]. Energy & Fuels, 2012, 26(2):976-983. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ded211921909b7aa7276d8bcdb445647

    [44]

    Gao X, Yang S, Sun W. A global pathway selection algorithm for the reduction of detailed chemical kinetic mechanisms[J]. Combustion and Flame, 2016, 167:238-247. DOI: 10.1016/j.combustflame.2016.02.007

    [45]

    Schwer D A, Lu P, Green W H. An adaptive chemistry approach to modeling complex kinetics in reacting flows[J]. Combustion and Flame, 2003, 133(4):451-465. DOI: 10.1016/S0010-2180(03)00045-2

    [46]

    Liang L, Stevens J G, Farrell J T. A dynamic multi-zone partitioning scheme for solving detailed chemical kinetics in reactive flow computations[J]. Combustion Science and Technology, 2009, 181(11):1345-1371. DOI: 10.1080/00102200903190836

    [47]

    Liang L, Stevens J G, Raman S, et al. The use of dynamic adaptive chemistry in combustion simulation of gasoline surrogate fuels[J]. Combustion and Flame, 2009, 156(7):1493-1502. DOI: 10.1016/j.combustflame.2009.02.008

    [48]

    Ren Z, Xu C, Lu T, et al. Dynamic adaptive chemistry with operator splitting schemes for reactive flow simulations[J]. Journal of Computational Physics, 2014, 263:19-36. DOI: 10.1016/j.jcp.2014.01.016

    [49]

    Liang L, Stevens J G, Farrell J T. A dynamic adaptive chemistry scheme for reactive flow computations[J]. Proceedings of the Combustion Institute, 2009, 32(1):527-534. https://www.sciencedirect.com/science/article/pii/S1540748908000941

    [50]

    Yang H T, Ren Z Y, Lu T F, et al. Dynamic adaptive chemistry for turbulent flame simulations[J]. Combustion Theory and Modelling, 2013, 17(1):167-183. DOI: 10.1080/13647830.2012.733825

    [51]

    Ren Z Y, Liu Y F, Lu T F, et al. The use of dynamic adaptive chemistry and tabulation in reactive flow simulations[J]. Combustion and Flame, 2014, 161(1):127-137. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1026da2ed0980362210dd0176c0f6768

    [52]

    Liang Y W, Pope S B, Pepiot P. A pre-partitioned adaptive chemistry methodology for the efficient implementation of combustion chemistry in particle PDF methods[J]. Combustion and Flame, 2015, 162(9):3236-3253. DOI: 10.1016/j.combustflame.2015.05.012

    [53] 刘再刚, 杨帆, 孔文俊.基于CO-DAC的合成气反应极限分析[J].燃烧科学与技术, 2017, 23(1):61-67. http://d.old.wanfangdata.com.cn/Periodical/rskxyjs201701010

    Liu Z G, Yang F, Kong W J. Reaction limit analysis of syngas based on CO-DAC method[J]. Journal of Combustion Science and Technology, 2017, 23(1):61-67. http://d.old.wanfangdata.com.cn/Periodical/rskxyjs201701010

    [54]

    Oluwole O O, Ren Z, Petre C, et al. Decoupled species and reaction reduction:An error-controlled method for Dynamic Adaptive Chemistry simulations[J]. Combustion and Flame, 2015, 162(5):1934-1943. DOI: 10.1016/j.combustflame.2014.12.012

    [55]

    Sun W, Gou X, El-Asrag H A, et al. Multi-timescale and correlated dynamic adaptive chemistry modeling of ignition and flame propagation using a real jet fuel surrogate model[J]. Combustion and Flame, 2015, 162(4):1530-1539. DOI: 10.1016/j.combustflame.2014.11.017

    [56]

    Aceves S M, Martinez-Frias J, Flowers D L, et al. A decoupled model of detailed fluid mechanics followed by detailed chemical kinetics for prediction of iso-octane HCCI combustion[R]. SAE Technical Paper, 2001-01-3612, 2001.

    [57]

    Goldin G M, Ren Z, Zahirovic S. A cell agglomeration algorithm for accelerating detailed chemistry in CFD[J]. Combustion Theory and Modelling, 2009, 13(4):721-739. DOI: 10.1080/13647830903154542

    [58]

    Xie W, Lu Z, Ren Z, et al. Dynamic adaptive chemistry via species time-scale and Jacobian-aided rate analysis[J]. Proceedings of the Combustion Institute, 2017, 36(1):645-653. https://www.sciencedirect.com/science/article/pii/S1540748916303637

    [59]

    Li Z, Lewandowski M T, Contino F, et al. Assessment of on-the-fly chemistry reduction and tabulation approaches for the simulation of moderate or intense low-oxygen dilution combustion[J]. Energy & Fuels, 2018, 32(10):10121-10131.

    [60]

    Xie W, Lu Z, Ren Z, et al. Dynamic adaptive acceleration of chemical kinetics with consistent error control[J]. Combustion and Flame, 2018, 197:389-399. DOI: 10.1016/j.combustflame.2018.08.018

    [61]

    Zhou D, Tay K L, Li H, et al. Computational acceleration of multi-dimensional reactive flow modelling using diesel/biodiesel/jet-fuel surrogate mechanisms via a clustered dynamic adaptive chemistry method[J]. Combustion and Flame, 2018, 196:197-209. DOI: 10.1016/j.combustflame.2018.06.008

    [62]

    Zhou L, Zhong L, Qin W, et al. Application of cell agglomeration algorithm coupled with dynamic adaptive chemistry for transient engine simulation of diesel fuel[J]. Fuel, 2018, 234:1313-1321. DOI: 10.1016/j.fuel.2018.07.155

    [63]

    Imren A, Haworth D C. On the merits of extrapolation-based stiff ODE solvers for combustion CFD[J]. Combustion and Flame, 2016, 174:1-15. DOI: 10.1016/j.combustflame.2016.09.018

    [64]

    Yang S, Ranjan R, Yang V, et al. Parallel on-the-fly adaptive kinetics in direct numerical simulation of turbulent premixed flame[J]. Proceedings of the Combustion Institute, 2017, 36(2):2025-2032. DOI: 10.1016/j.proci.2016.07.021

    [65]

    Perini F, Galligani E, Reitz R D. A study of direct and Krylov iterative sparse solver techniques to approach linear scaling of the integration of chemical kinetics with detailed combustion mechanisms[J]. Combustion and Flame, 2014, 161(5):1180-1195. DOI: 10.1016/j.combustflame.2013.11.017

    [66]

    Lu T F, Law C K, Yoo C S, et al. Dynamic stiffness removal for direct numerical simulations[J]. Combustion and Flame, 2009, 156(8):1542-1551. DOI: 10.1016/j.combustflame.2009.02.013

    [67]

    Gou X L, Sun W T, Chen Z, et al. A dynamic multi-timescale method for combustion modeling with detailed and reduced chemical kinetic mechanisms[J]. Combustion and Flame, 2010, 157(6):1111-1121. DOI: 10.1016/j.combustflame.2010.02.020

    [68]

    Zhou D, Yang W. A heterogeneous multiscale method for stiff combustion chemistry integration in reactive flows[J]. Combustion and Flame, 2018, 188:428-439. DOI: 10.1016/j.combustflame.2017.09.039

    [69]

    Gao Y, Liu Y, Ren Z, et al. A dynamic adaptive method for hybrid integration of stiff chemistry[J]. Combustion and Flame, 2015, 162(2):287-295. DOI: 10.1016/j.combustflame.2014.07.023

    [70]

    Hochbruck M, Ostermann A. Exponential integrators[J]. Acta Numerica, 2010, 19:209-286. DOI: 10.1017/S0962492910000048

    [71]

    Tangman D Y, Gopaul A, Bhuruth M. Exponential time integration and Chebychev discretisation schemes for fast pricing of options[J]. Applied Numerical Mathematics, 2008, 58(9):1309-1319. DOI: 10.1016/j.apnum.2007.07.005

    [72]

    Hochbruck M, Honig M, Ostermann A. Regularization of nonlinear Ⅲ-posed problems by exponential integrators[J]. ESAIM-Mathematical Modelling and Numerical Analysis, 2009, 43(4):709-720. DOI: 10.1051/m2an/2009021

    [73]

    Tokman M, Loffeld J. Efficient design of exponential-Krylov integrators for large scale computing[J]. Procedia Computer Science, 2012, 1(1):229-237. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Open J-Gate000004239135

    [74]

    Gaudreault S, Pudykiewicz J A. An efficient exponential time integration method for the numerical solution of the shallow water equations on the sphere[J]. Journal of Computational Physics, 2016, 322:827-848. DOI: 10.1016/j.jcp.2016.07.012

    [75]

    Moler C, Van Loan C. Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later[J]. Siam Review, 2003, 45(1):3-49. DOI: 10.1137/S00361445024180

    [76]

    Higham N J. The scaling and squaring method for the matrix exponential revisited[J]. Siam Journal on Matrix Analysis and Applications, 2005, 26(4):1179-1193. DOI: 10.1137/04061101X

    [77]

    Saad Y. Analysis of some Krylov subspace approximations to the matrix exponential operator[J]. Siam Journal on Numerical Analysis, 1992, 29(1):209-228.

    [78]

    Sidje R B. Expokit:A software package for computing matrix exponentials[J]. Acm Transactions on Mathematical Software, 1998, 24(1):130-156. DOI: 10.1145/285861.285868

    [79]

    Niesen J, Wright W M. Algorithm 919:A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators[J]. Acm Transactions on Mathematical Software, 2012, 38(3):1-19. http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_0907.4631

    [80]

    Lord G J, Stone D. New efficient substepping methods for exponential timestepping[J]. Applied Mathematics and Computation, 2017, 307:342-365. DOI: 10.1016/j.amc.2017.02.052

    [81]

    Vo H D, Sidje R B. Implementation of variable parameters in the Krylov-based finite state projection for solving the chemical master equation[J]. Applied Mathematics and Computation, 2017, 293:334-344. DOI: 10.1016/j.amc.2016.08.013

    [82]

    Bisetti F. Integration of large chemical kinetic mechanisms via exponential methods with Krylov approximations to Jacobian matrix functions[J]. Combustion Theory and Modelling, 2012, 16(3):387-418. DOI: 10.1080/13647830.2011.631032

    [83]

    Curtis N J, Niemeyer K E, Sung C-J. An investigation of GPU-based stiff chemical kinetics integration methods[J]. Combustion and Flame, 2017, 179:312-324. DOI: 10.1016/j.combustflame.2017.02.005

    [84] 刘再刚, 孔文俊.求解刚性燃烧化学反应系统的Krylov子空间中的指数积分法[J].工程热物理学报, 2018, 39(9):2062-2071. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gcrwlxb201809027

    Liu Z G, Kong W J. An exponential integrator in Krylov subspace to solve stiff combustion chemical reaction system[J]. Journal of Engineering Thermophysics, 2018, 39(9):2062-2071. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gcrwlxb201809027

    [85] 刘再刚, 隋春杰, 韩文虎, 等.基于CO-DAC反应机理简化的非预混射流火焰的大涡模拟[J].燃烧科学与技术, 2018, 24(1):39-45. http://d.old.wanfangdata.com.cn/Periodical/rskxyjs201801007

    Liu Z G, Sui C J, Han W H, et al. LES for turbulent non-premixed jet flame with CO-DAC reduction[J]. Journal of Combustion Science and Technology, 2018, 24(1):39-45. http://d.old.wanfangdata.com.cn/Periodical/rskxyjs201801007

    [86]

    Barlow R S, Frank J H, Karpetis A N, et al. Piloted methane/air jet flames:Transport effects and aspects of scalar structure[J]. Combustion and Flame, 2005, 143(4):433-449. DOI: 10.1016/j.combustflame.2005.08.017

    [87]

    Yang S, Wang X J, Yang V, et al. Comparison of flamelet/progress-variable and finite-rate chemistry LES models in a preconditioning scheme[R]. AIAA-2017-0605, 2017.

    [88]

    Yang S, Wang X J, Yang V, et al. Comparison of finite-rate chemistry and flamelet/progress-variable models Ⅱ: Sandia Flame E[R]. AIAA-2018-1427, 2018.

    [89]

    Yang S, Ranjan R, Yang V, et al. Sensitivity of predictions to chemical kinetics models in a temporally evolving turbulent non-premixed flame[J]. Combustion and Flame, 2017, 183:224-241. DOI: 10.1016/j.combustflame.2017.05.016

    [90]

    Wu K, Contino F, Yao W, et al. On the application of tabulated dynamic adaptive chemistry in ethylene-fueled supersonic combustion[J]. Combustion and Flame, 2018, 197:265-275. DOI: 10.1016/j.combustflame.2018.08.012

    [91]

    Contino F, Lucchini T, D'errico G, et al. Simulations of advanced combustion modes using detailed chemistry combined with tabulation and mechanism reduction techniques[J]. SAE International Journal of Engines, 2012, 5:185-196. DOI: 10.4271/2012-01-0145

    [92]

    Zhou L, Lu Z, Ren Z Y, et al. Numerical analysis of ignition and flame stabilization in an n-heptane spray flame[J]. International Journal of Heat and Mass Transfer, 2015, 88:565-571. DOI: 10.1016/j.ijheatmasstransfer.2015.05.003

    [93]

    Lu Z, Zhou L, Ren Z Y, et al. Effects of spray and turbulence modelling on the mixing and combustion characteristics of an n-heptane spray flame simulated with dynamic adaptive chemistry[J]. Flow Turbulence and Combustion, 2016, 97(2):609-629. DOI: 10.1007/s10494-015-9702-5

    [94]

    Huang Y, Lu Z, Zheng L. Study of SiCl4/H2/O2 chemical kinetics and its application to fused silica glass synthesis[J]. Combustion Science and Technology, 2018, 190(10):1861-1885. DOI: 10.1080/00102202.2018.1476349

    [95]

    Shi Y, Liang L, Ge H-W, et al. Acceleration of the chemistry solver for modeling DI engine combustion using dynamic adaptive chemistry (DAC) schemes[J]. Combustion Theory and Modelling, 2010, 14(1):69-89. DOI: 10.1080/13647830903548834

    [96]

    Lu L, Lantz S R, Ren Z, et al. Computationally efficient implementation of combustion chemistry in parallel PDF calculations[J]. Journal of Computational Physics, 2009, 228(15):5490-5525. DOI: 10.1016/j.jcp.2009.04.037

    [97]

    Stewart G W. A Krylov-Schur algorithm for large eigenproblems[J]. Siam Journal on Matrix Analysis and Applications, 2002, 23(3):601-614. DOI: 10.1137/S0895479800371529

    [98]

    Hochbruck M, Lubich C, Selhofer H. Exponential integrators for large systems of differential equations[J]. Siam Journal on Scientific Computing, 1998, 19(5):1552-1574. DOI: 10.1137/S1064827595295337

    [99]

    Liu Z G, Consalvi J L, Kong W J. An exponential integrator with Schur-Krylov approximation to accelerate combustion chemistry computation[J]. Combustion and Flame, 2019, 203:180-189. DOI: 10.1016/j.combustflame.2019.01.031

    [100]

    Lu L Y, Pope S B. An improved algorithm for in situ adaptive tabulation[J]. Journal of Computational Physics, 2009, 228(2):361-386. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=331d60857be5dc1d0c837e9e87da1b43

    [101]

    Perini F, Galligani E, Reitz R D. An analytical Jacobian approach to sparse reaction kinetics for computationally efficient combustion modeling with large reaction mechanisms[J]. Energy & Fuels, 2012, 26(8):4804-4822. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=fd64e2068cd2c5edca40b7a696028b74

    [102]

    Smith G P, Golden D M, Frenklach M, et al. GRI-Mech 3.0.[EB/OL]. (1999-07-30)[2018-08-27]. http://www.me.berkeley.edu/gri_mech.

    [103]

    Wang H, You X, Joshi A V, et al. USC Mech Version Ⅱ: High-Temperature Combustion Reaction Model of H2/CO/C1-C4 Compounds.[EB/OL]. (2007-05)[2018-08-27]. http://ignis.usc.edu/USC_Mech_II.htm

  • 期刊类型引用(5)

    1. 王珏,王誉超,季辰. 超声速风洞带舵机状态全尺寸舵颤振亚临界试验. 空天防御. 2023(02): 77-83 . 百度学术
    2. 郭力,刘杰平,吕计男,王昕江,季辰,刘子强. 前缘钝度对高速翼型颤振边界的影响. 宇航总体技术. 2020(01): 47-55 . 百度学术
    3. 杨超,赵黄达,吴志刚. 吸气式高超声速飞行器热气动弹性研究进展. 北京航空航天大学学报. 2019(10): 1911-1923 . 百度学术
    4. 周波,高川,杨洋. 2m超声速风洞流场变速压控制方法研究. 实验流体力学. 2019(06): 72-77 . 本站查看
    5. 季辰,刘子强,李锋. 钝前缘梯形翼高超声速风洞颤振试验. 气体物理. 2018(01): 54-63 . 百度学术

    其他类型引用(3)

图(6)  /  表(1)
计量
  • 文章访问数:  242
  • HTML全文浏览量:  115
  • PDF下载量:  25
  • 被引次数: 8
出版历程
  • 收稿日期:  2018-08-29
  • 修回日期:  2019-03-13
  • 刊出日期:  2019-08-24

目录

/

返回文章
返回
x 关闭 永久关闭