电力能源和交通系统的可用性在很大程度上控制了湍流燃烧过程,如燃气轮机、内燃机、燃烧化石或可再生燃料。进一步优化这样的知名系统只能涉及所有相关过程的一个更好的理解。详细的定量实验需要和非常有用的,但有时不可能的,通常只有少数火焰数量有限。作为补充,详细()直接数值模拟(dns)越来越多地获得理由作为一个可靠的工具的详细调查对基本的认识各种湍流燃烧现象( 1]。数字技术的进展以及计算能力现在允许定量调查的湍流反应流越来越现实的条件下,模拟实际相关的燃料在高雷诺数湍流( 再保险 t ),而同时采用复杂的物理化学模型来描述湍流,分子运输和化学( 2- - - - - - 4]。

实际上只有两种可能性获得统计上显著的结果。系综平均是最直接的方法,基于重复相同的数值实验,”然后平均观测[ 5]。这有时会支持甚至取代了一个空间平均在统计学上均匀的方向,当可用。作为一种替代方法,可以继续仿真很长一段时间,采用样本之间的平均时间间隔足够长的时间间隔。在这个最后的案例中,湍流时间必须保持在恒定水平。因此,时间平均只有空间发展的动荡,通过边界条件通常注射,或通过人为地流在time-decaying模拟,这又与许多理论和实践问题。依靠遍历性,平均方法应该提供相同的结果。考虑数值需求,这两种方法都导致大量增加的计算成本。仿真必须重复(系综平均)或追求在物理时间很长(平均),通常导致一个数量级的增加一个DNS。

更一般的,需要准确的物理模型在实际燃烧计算清楚展示了各种火焰配置(见,例如, 2, 6- - - - - - 8)和引用)。不幸的是,将卷(或大部分)粘度( κ )运输术语仍然是有争议的文献中,因为它已经在几乎所有的传统上被忽视的湍流燃烧的研究到目前为止没有任何令人信服的理由( 9),即使一些出版物已经证明它是错误的忽视它先天的 10, 11]。应该为很多应用程序确实小的影响。另一方面,的价值 κ 为氢,例如,非常高,可能导致明显的火焰的修改。此外,研究的影响 κ 火焰结构而只考虑层流或湍流条件温和可能导致偏见的结论从更现实的,更高的值 再保险 t 很少被评估。

忽视了体积粘性项计算模型似乎逻辑以非常低的马赫数和边界层流动。在实践中,额外的计算成本评估的多组分的混合物也会导致一个更大的程度上排斥。更重要的是,在一定程度上误导的过分自信的建议的斯托克斯( 9可能导致这个系统的逃避。的确,没有体积等稀释单原子气体的粘度影响不( 12)是一个例外,而不是一个规则( 9, 13),因为 κ 出现非常明显的密集的气体和液体。虽然单原子气体显示没有证据表明体积粘性影响,bulk-to-shear粘度比实验测量值( κ / η )相差很大的双原子气体使用现在在学术和实际燃烧实验和数值研究[ 12, 14, 15]。在[ 15),至少说体积粘度是相同的数量级的剪切粘度( η 在300 K),比例 κ / η 对氢高达52 1000 K的温度。我们最近的调查结果显示,湍流预混火焰燃烧hydrogen-containing燃料的不可忽视的影响体积粘性项的湍流火焰结构,即使是平均量( 16]。另一方面,第一个DNS计算火焰燃烧甲烷显示一个小体积粘度的影响考虑情况下全球火焰属性如volume-integrated热释放率和紊流燃烧速度 17]。由此产生的图像是因此非常复杂。一个系统的定量分析的影响 κ 在动荡的甲烷火焰因此有必要完善以前的结果以确定的方式。这是高相关性,特别是燃烧DNS和建模。

因为单个孤立的DNS模拟可能有时会误导,系统的DNS计算已经意识到为了开展ensemble-averaging,卖好统计学意义( 5为观察)。

最后三个主要挑战因此被提出在目前工作:化学和运输的复杂性加大了考虑甲烷火焰粘度和体积效应,同时分别同时增加湍流强度。作者的知识,这样的组合数值,物理,和流动条件从未考虑同时在相同的分析。

在接下来的部分(部分 2),一个轮廓的计算程序,数值方法的描述,问题配置和初始化部分 3。最近一些评论关于模块添加到使用DNS工具对细粒度并行性进行了讨论。中给出的数值结果是最后一节 4结束语(前一节 5)。

DNS方法在于解决尽可能完全所有的物理空间和时间尺度嵌在代表流方程,没有添加任何湍流模型。一个DNS必须提供一个流体动力学和火焰结构的精确解。尽管这种方法需要禁止数值成本实际配置,它提供了一个很好的补充实验为了评估各种物理机制的重要性,对火焰结构获得互补信息,因此提高湍流燃烧模型( 18]。

DNS代码用于这项工作是大规模并行火焰的能手, parcomb( 19, 20.),解决了完整的可压缩无功n - s系统加上详细的化学和多组分传输模型(通过耦合 chemkin( 21), 运输( 22), eglib( 23]库): (1) ∂ ρ ∂ t + ∂ ( ρ u j ) ∂ x j = 0 , ∂ ( ρ u 我 ) ∂ t + ∂ ( ρ u j u 我 ) ∂ x j = - - - - - - ∂ p ∂ x 我 + ∂ τ 我 j ∂ x j , ∂ ( ρ Y k ) ∂ t + ∂ ( ρ u j Y k ) ∂ x j = - - - - - - ∂ ( ρ Y k V k j ) ∂ x j + ω ̇ k ; k = 1、2 , … , N 年代 , ∂ ( ρ e t ) ∂ t + ∂ ( ρ u j e t ) ∂ x j = - - - - - - ∂ 问 j ∂ x 我 + ∂ ( τ 我 j u 我 ) ∂ x j , p ρ = R W ¯ T , 在哪里 ρ 表示混合密度, u j 水动力组件的速度, p 压力, τ 我 j 应力张量, N 年代 物种的总数, V k j 物种的扩散速度的组件 k 的方向, j , ω ̇ k 化工生产的物种 k 和 问 j 的 j 个分量的热通量向量。体积粘度系数 κ 是一个函数的液体本地属性和显式出现在 τ 我 j 动量和能量方程: (2) τ 我 j = - - - - - - p δ 我 j + η ( ∂ u 我 ∂ x j + ∂ u j ∂ x 我 ] + ( κ - - - - - - 2 3 η ) ∂ u k ∂ x k δ 我 j , 在哪里 δ 克罗内克符号和吗 η 是动态(或剪切粘度)。DNS代码已经加上 eglib运输图书馆( 23)体积粘度的实际评价运输项。单一输入数据子程序 EGSKm (),适合fine-grain-distributed架构与大量的处理器,已经被使用。整数 米 连接到子例程的名称指的是各种方法/模型用于评估 κ ( 15, 16]。2和6之间的不同。的价值就越高 米 ,更昂贵的底层算法和更精确的表达式 κ ,因此我们选择 米 = 6 ( EGSK6 ()在目前的研究)的计算。最简单和最便宜的模型( 米 = 2 , EGSK2 ())只考虑交通系统相关矩阵内部的能量,而 米 = 6 ( EGSK6 ())对应不仅最精确的,但也是最昂贵的模型,因为它解决了完整的矩阵系统相关的转化和内部能量碰撞积分和依赖与温度有关的比率,这是评为[中讨论 15, 24, 25]。最后,直接进行反演计算体积粘度对于一个给定的混合物,以温度和物种质量分数。我们感兴趣的读者参考 eglib用户手册( 23为进一步的细节和[] 15)一个广泛讨论的话题。

上述系统的控制方程是解决 parcomb在传统的笛卡尔网格高阶准确和nondissipative数值方案。空间sixth-order中央计划,逐步减少到四阶边界附近使用。在升级的过程中从原来的2 d版的代码当前3 d状态( 4, 20.),称为斜对称的一种改进方案制定( 26)已经实现了对流条件以减少进一步数值耗散和增加稳定性。时间集成在一个显式执行方式与四阶龙格-库塔方案。扩展的n - s特征边界条件(NSCBCs [ 27, 28]),应用nonreflecting边界和压力放松沿着流出的面孔。

parcomb并行使用一个三维域分解和MPI消息传递协议,提供一个良好的峰值性能和近乎完美的并行扩展为一个完整的大规模的三维运行4096计算核心的 IBM BlueGene / P。更多细节关于代码结构、优化应用程序,和最近的升级/添加模块可以在[ 4, 29日]。

为了获得更高的价值 再保险 t 在精密并联系统中,两个湍流生成技术基于数字滤波器( 30.和随机噪声扩散 31日)杂化,执行,并行大规模并行计算机上( 16]。通过这个简单、灵活、准确的方法,限制问题的大小由先前基于逆FFT实现发电机( 32)已经缓解,铺平了道路模拟大域相当高雷诺数湍流。对结果提出后,生成的初始湍流场已经系统地使用这种混合动力技术。

第二个主要问题需要注意对细粒度并行性是有效的,完全并行数据输入/输出(I / O)。传统的顺序I / O方法已经被完全并行I / O(通过MPI-I / O),多个进程并行程序的访问数据从一个共同的(读/写),共享文件。这提供了更高的性能(加速在写作/阅读所有文件所需时间至少3倍)和单(重启/解决方案)数据文件。

化学计量球形预拌沼气空气火焰被认为是在所有的计算,在一个立方计算域的两侧 l = 3所示。0 4.0厘米,一个统一的网格间距的26岁 μ米,允许足够的分辨率的空间尺度上对所有配置。的点火和随后的扩张/发展这样一个预混合火焰内核的影响下的湍流流场是一个非常有趣的配置允许研究湍流火焰远离外部约束的影响,如墙壁和人工边界条件。它有另外一个直接相关的工业案例包括火花点燃式内燃机和燃气轮机再起以及安全问题。

甲烷氧化由25-step骨架建模方案由16种(CH₄阿,₂H₂H₂O, CH₂O,有限公司有限公司₂,何₂,哦,H、O, CH₃、HCO H₂O₂,CH₃O, N₂)[ 33]。这里的反应机理是保留由于其简单性和稳定性,并提供足够准确的结果为精益化学计量的条件。它已经成功地用于大规模多维nonpremixed甲烷喷射火焰的直接计算( 34),最近,高度湍流预混火焰( 4, 17]。然而,它将对甲烷火焰不足由于缺乏C₂和更高的碳链反应的原因 Φ ≤ 1。0 目前的研究。

最初的系统是一个固定热( T b ≈ 2 200年 K)完美的球形层流火焰内核的初始半径 r o = 5.0 mm,位于中心的计算框和新的大气包围预混合的甲烷和空气的混合物在一个未燃的温度 T u = 300年 K。最初的原始轮廓度变量 Φ 根据规定 (3) Φ = Φ o + Δ Φ 2 ( 1 - - - - - - 双曲正切 ( 年代 ⋅ ( r - - - - - - r o r o ) ) ] , 在哪里 Δ Φ 之间的区别是初始值( Φ o 在新鲜和燃烧气体的混合物。常数 年代 是新鲜的刚度的测量/燃烧气体接口。初始质量分数值 Y C H 4 = 0.055 和 Y O 2 = 0.220 在 T u 外的内核 Y C O 2 = 0.120 和 Y H 2 O = 0.137 在 T b 内部。一个适当的补充氮添加无处不在。

对湍流雷诺数的影响进行调查 再保险 t 基于积分尺度火焰结构,同样的计算重复了相同的初始pseudoturbulent结构和初始成分,但先后高紊流速度波动。均方根速度波动 u ′ 湍流雷诺数 再保险 t 艾迪的人事变动 τ 用于量化火焰/湍流交互time-decaying动荡( 35),Karlovitz号码 卡 ≈ ( ( u ′ / 年代 l ) 3 ( l t / δ t h ) - - - - - - 1 ] 1 / 2 给出了在表 1。混合物的特性粘度, ν = 1.56 × 1 0 - - - - - - 5 米²/ s,层流火焰速度 年代 l = 0.504 m / s,层流火焰热厚度 δ t h = ( T b - - - - - - T u ) / 马克斯 | ∇ T | = 0.36 毫米,整体规模 l t / δ t h = 8.89 这里给出所有模拟都是常数。

基于这些湍流特性、火焰认为这里属于薄反应区(TRZ)政权根据修改后的政权彼得斯(图 36, 37),如图 1。它是政权相信小湍流漩涡能够穿透和扰乱预热区,但达不到这样的反应区,因为一个数量级差异在这些特征层的厚度。

的一般视图配置如图 2,iso-surface温度与H₂O₂片和有限公司₂片和流行(薄白线)显示,揭示了严重皱后火焰前锋与三维交互time-decaying均匀各向同性合成湍流速度场为0.8毫秒。

总数计算核心介于512和096是用来解决这个问题 IBM POWER5集群(HPCx)系统EPCC(提供总峰值性能6 Tflop / s维持)。通常,所需的计算时间是10天到达 t / τ = 1 。最长的计算需要大约两个月。下面展示的结果中,无量纲时间 t = 1。2 τ ,需要得到适当的湍流和化学之间的平衡条件( 35),是系统保留的分析。

获得DNS数据集使用内部后处理库处理 AnaFlame( 38, 39]。为 例3和 4,每个实验都进行两次完全相同的初始和边界条件包括湍流特性,除了一个事实,即体积粘度方面在另一但激活释放。此外,系综平均的使用是为了改善结果的统计学意义通过重复每一对( κ ≠ 0 和 κ = 0 ) 案例3不同湍流字段初始化随机种子不同但拥有相同的初始强度( u ′ ),然后平均观测。因为一个随机数发生器是用来确定速度阶段 40),每一对DNS与湍流的相同的全局属性(频谱、相关性、波动、雷诺数等)但对应不同的初始空间特性,从而不同的进化(即。实现)的火焰。这些进一步的实现将为这里的结果提供统计学意义( 5]。接下来,所有资料显示 案例3因此6实现的平均值。

瞬时解的条件分析了统计量与建模相关如温度、热释放,选择质量分数说明turbulence-impaired火焰结构。体积粘度的影响是通过比较上述评估概要文件为每个的双重计算 例3和 4。

数字技术的进展以及计算能力现在允许定量调查的湍流反应流越来越现实的条件包括详细的物理化学模型。化学计量预混甲烷火焰的直接数值模拟参数研究中被认为是包括与六对一个案例的计算 再保险 t = 615年 2 460(因子4)通过访问高端并行计算机在欧洲水平。后与I / O相关的解决问题和初始湍流生成,详细分析使用专用的基于matlab的图书馆已经完成 AnaFlame( 39]。

火焰速度随气流搅拌增加而迅速增加,火焰明显增厚。同时,放热峰值条件概况和主要和次要物种质量分数有系统地降低湍流强度增加。在所有考虑情况下,预混合的甲烷火焰体积粘度的影响可以忽略不计。任何差异在层流计算,证实了理论结果低马赫数条件。以前的观测对湍流氢火焰,揭示一个明确的体积粘度的影响,不适用目前甲烷火焰雷诺数在研究中考虑。节省计算资源,包含多组分体积粘度影响的多维计算湍流火焰燃烧甲烷和可能更高的碳氢化合物以及因此气馁,只要低马赫数。

这个研究表明也重复的重要性DNS实现为了获得统计上显著的数据。单一的实现可能会导致虚假的差异,迅速平滑时平均几个结果,观察到这里当比较的结果 案例3(多个DNS实现) 例4(单DNS实现)。