乙烯燃烧简化化学动力学模型及其验证

采用"准稳态"方法,从乙烯燃烧的详细化学反应动力学模型出发,建立了包含20个组分和16步总包反应的简化化学动力学模型.为检验动力学模型的有效性,在我国同步辐射实验室燃烧站上开展了乙烯/O<,2>/Ar的层流预混火焰组分测量.应用一维数值方法对实验结果进行了模拟,计算的组分分布和实验测量数据进行了对比.结果表明:简化反应动力学模型能有效地再现详细基元反应模型的反应机理,具有较高精度;但采用乙烯详细化学动力学模型计算结果与实验测量还存在一定的差异,需要进一步改进.     

一种基于敏感性分析的RP3替代燃料简化机理

为了研究RP3航空煤油的点火和燃烧特性,基于敏感性分析方法,对一种RP 3数值模拟替代燃 料正葵烷的化学反应机理（62组分344步）进行简化,得到对平衡温度和组分浓度有重 要影响的基元反应,建立一种该替代燃料的36组分62步简化反应机理。利用CHEMKIN 4.1中完 全搅拌反应器对该简化机理进行研究与检验,并在预混燃烧器中对该简化机理进行数值计算 ,与详细机理结果和试验数据进行比较。计算结果表明：在完全搅拌反应器中采用62步机理 计算得到的平衡温度,反应物和主要生成物的摩尔分数与344步详细化学反应机理计算结果 整体变化趋势吻合较好;在预混燃烧器中采用62步机理得到的进口反应物和出口生成物摩尔 分数与详细机理结果和试验数据稍有差异,但基本变化趋势基本一致。因此,文中获得的62 步简化机理能在较大范围内反映正葵烷的燃烧性能,并且大大提高其计算效率。     

激波管在酚醛树脂高温热解动力学研究中的应用

使用激波管作为加热手段,利用其加热速率快的优点,突破了传统方法在加热速率上的限制,研究了酚醛树脂在1400~1700K温度范围内的热解动力学.主要碳氢产物是甲烷、乙烯、乙炔和苯.通过对反应扩散过程的分析,考察了扩散对热解过程的影响.结果表明,实验中酚醛树脂的反应扩散过程迅速达到稳态,扩散影响可以忽略,首次获得了酚醛树脂在芳环开环热解机制下的热解速率常数.     

含化学反应流动的TVD格式及其在燃气射流中的应用

以Harten标准TVD格式为基础，结合固体火箭燃气射流的特点，以数学方法系统地推导了适合于高温、高压和高速流体流动的数值格式，给出了压力偏导数的合理计算公式；利用化学动力学知识，对火箭燃气射流流场中存在的多组分、含有限速率的化学反应问题进行了论述，阐明了化学反应质量源项的求解方法。以12组分9反应方程模型为例，利用编制的计算程序，对某火箭燃气自由射流流场进行了模拟。通过对结果的分析，肯定了数值格式的正确性。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

Roe格式在多组元燃烧流场数值模拟中的应用

对Roe通量分裂格式进行了详细讨论,并将其推广应用到基于非结构网格的高超声速多组元燃烧流场的数值模拟中。控制方程采用三维多组元Euler方程,重点针对Roe格式开展研究;考虑到燃烧流场求解过程中出现的刚性问题,对化学反应源项采用点隐式方法处理。对超声速条件下钝头体激波诱导燃烧流场进行了数值模拟,比较分析了3种经典化学动力学模型对流场结构的影响,获得与数值模拟、试验数据相符的结果。结果表明发展的方法可以用于多组元燃烧流场的数值分析,且化学动力学模型的选取直接影响激波诱导燃烧流场中诱导区厚度。     

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

为研究湍流燃烧数值模拟中化学反应机理计算的加速方法,讨论了动态自适应化学（Dynamic Adaptive Chemistry,DAC）方法和Krylov子空间近似的指数格式的应用情况。在湍流火焰大涡模拟中,使用DAC简化可以加速化学反应计算。然而,在并行燃烧数值模拟中,处理器核心的负载极度不平衡,加速效果有限。而Krylov子空间近似的指数格式的加速效果可以作用于每个处理器核心,更有利于整体计算效率的提高。在同等精度下,相比于隐式格式耦合DAC和MTS加速方法,Krylov子空间近似的指数积分格式对化学反应计算的加速效果更为显著。     

由火焰聚心和激波聚焦诱导的爆震波研究

通过对几种不同结构形式的爆震发生器进行数值模拟,研究了激波聚焦和火焰聚心现象、气动特性及其机理.数值计算采用多组分理想气体详细的化学反应机理、二维轴对称非定常流动Navier-Stokes方程来模拟流体动力学和化学动力学过程.数值计算表明用较低的点火能量对射流火焰燃烧器中可燃混合物点火,层流火焰在狭窄管壁作用下加速,射流火焰在轴线上汇聚过程有利于激波的加强,强激波加速火焰,在多重激波与火焰反复作用下,激波和火焰面之间出现热点,热点迅速放大并形成压力很高的过驱爆震波,而后衰减为稳定的爆震波,不同的激波聚焦腔爆震波的形成过程不同.通过对数值计算的结果进行分析,得到了起爆距离和稳定爆震距离,为进一步试验提供参考.     

基于混合网格上的高温非平衡流的数值模拟

新一代再入飞行器及空间传输系统需要了解其飞行时的气动问题及热力学问题。数值模拟高空、低密度、高熵及非平衡流动是一具有挑战性的问题。本文针对带座舱飞船高超声速再入大气过程中存在的严重气动加热现象,利用混合网格及Osher逆风格式,数值求解了三维化学非平衡Navier-Stokes方程,其中化学模型是5组分1 7个化学反应的空气化学模型,对带座舱飞船再入高度为40 km和马赫数为2 0 ,1 0的化学非平衡流场进行了数值模拟,给出了迎角为0°和2 0°情况下的各个组分的密度分布、压力等参数,并与量完全气体的计算结果进行了比较。     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

城市污水处理厂污泥热解机理

本研究采用热重分析仪对南昌市青山湖污水污泥进行热重分析实验,分别获得了20 K/min和40 K/min加热速率下污泥的TG-DTG曲线,得出了污水处理厂污泥的热解特性。根据实验结果对干燥污泥热解过程进行了分析且利用积分法确定了热解机理,求出反应动力学参数-频率因子A、活化能E。     

几何非线性固体壳单元新列式的研究

通过定义广义应力，提出了一个改进的刚度矩阵，以克服固体壳元的厚度自锁问题；由一个基于广义应力的新的非线性变分泛函，推导了一个用于几何非线性分析的十八节点固体壳单元，精心选择低、高阶应力插值模式，保证二者正交，且对应于低阶项的刚度阵与减缩积分单元相当，对应高阶项的刚度阵用于克服单元的零能模式且可推得显式，从而显著提高了计算效率，此外该单元还拥有优秀的收敛性能，即使在非常大的载荷步下也能取得好的收敛结果。     

截面形式对复合材料杆件固化变形的影响

热固性复合材料构件应用广泛,但其在加热固化成型过程中易产生变形。通过优化铺层设计,可以减小碳纤维增强环氧树脂细长杆件固化过程的变形,但是忽略了截面形式对固化变形的影响程度。本文提出在复合材料大长径比构件设计时应采用有限元分析(Finite element analysis, FEA)方法进行热变形分析,充分考虑截面对称性对产品固化变形的影响,有效降低该类产品固化过程中的变形。     

基于可靠性的弹性折叠机构设计优化方法研究

描述了基于可靠性的弹性折叠机构设计优化方法。采用特征造型方法建立了弹性折叠机构几何模型,CAD二次开发技术实现了折叠机构模型的参数化自动建模;考虑折叠机构几何参数不确定性,建立了折叠机构可靠性动态分析模型;基于功能函数的极小变换法将机构动态可靠性分析问题简化为静态问题;采用一次二阶矩法分析了机构折叠时间、最大应力和机构寿命的可靠度;以折叠机构几何参数为设计变量,采用修正的可行方向法驱动优化使折叠机构质量最轻且寿命更长。数值分析结果显示功能函数的极小变换法将机构动态可靠性分析简化为静态可靠性分析,可以在保证分析精度的前提下,提高可靠性分析效率;同时与传统优化方法比较表明,可靠性优化设计方法可以通过牺牲部分优化目标性能换取优化设计结果可靠性的提高。     

采用CARS试验技术与UFPV数值方法研究航空发动机燃烧室

在自主开发的软件平台上,采用基于URANS的方法计算航空发动机燃烧室的三维两相燃烧流动,考虑了液态燃油从液膜-液滴-燃气-燃烧的完整物理化学过程。其中,颗粒相采用LISA一次破碎模型,KH-RT二次破碎模型和标准的蒸发模型,湍流燃烧模型采用可以考虑非稳态燃烧特性的非稳态火焰面/反应进度变量方法,得到了航空发动机燃烧室中温度、组分浓度和燃油液滴的颗粒直径分布规律。同时,采用CARS光学手段测量燃烧室主燃区的温度分布,并将数值计算结果与光学试验测量值进行比较,数值计算结果和试验值吻合较好,数值计算误差小于7.3%。说明了本文的数值计算方法和UFPV方法在计算航空发动机燃烧室的两相燃烧流动时具有较高的精度。     

煤粉热解研究的激波管技术

为开拓激波管技术在煤粉热解研究中的应用,对一座激波管实体的系统、运行原理、煤粉撒播及热解气体的回收和分析技术等作了介绍。在激波管装置上共作了八种典型煤的热解实验。文中报导了一种典型国产煤在1700K 热解温度下,逸出的十种碳氢化合物和其它气体的产量数据,与澳大利亚类似煤种的结果基本吻合。     

桨尖下反对旋翼BVI噪声特性的抑制分析

基于计算流体力学(Computational fluid dynamics,CFD)方法,分析研究了新型下反桨尖对前飞状态旋翼桨/涡干扰(Blade-vortex interaction,BVI)噪声特性的影响。首先,通过求解Navier-Stokes方程获得准确的噪声声源信息,湍流模型采用一方程S-A模型。采用双时间法进行时间推进,为了提高流场的收敛速度,在伪时间方向上使用高效的隐式LU-SGS格式推进,并采用并行算法进行加速。然后,在获得的可靠声源信息基础上,采用基于FW-H方程的Farassat 1A公式求解噪声。对有试验结果的算例进行了对比计算分析,验证了本文噪声预测方法的有效性。在此基础上,针对具有典型BVI噪声特征的前飞斜下降状态,开展了不同下反角度桨尖新型旋翼噪声辐射特性的计算分析,通过噪声辐射球的对比结果表明,选择适当的下反桨尖,可以有效地降低前飞斜下降状态下旋翼BVI噪声,从而达到较好的旋翼噪声抑制效果。     

基于CFD方法的主动襟翼控制旋翼翼型涡特性研究

基于计算流体力学(Computational fluid dynamics,CFD)技术,建立了一个适用于旋翼二维主动襟翼控制(Active flap control,AFC)数值模拟的方法。在满足对翼型参数化分析的前提下,使用Euler方程求解以提高计算速度,并采用嵌套网格方法对AFC旋翼后缘襟翼进行运动控制。应用所建立的方法,首先进行了算例验证计算,然后着重对AFC旋翼翼型进行了数值模拟。在此基础上,进一步开展了AFC旋翼翼型主要参数对后缘涡影响的计算分析。结果表明:提高桨尖马赫数、增加后缘小翼摆动频率能加快涡产生速度;而提高桨尖马赫数、增大后缘小翼摆动幅度和后缘小翼长度能增大涡的强度;但增大后缘小翼与主桨叶缝隙间距仅在一定范围内能够增加涡强度。     

航天飞行器脉动压力数值计算方法综述

航天飞行器以高马赫数飞行时,由于边界层的复杂流动将产生强烈的脉动压力环境。由此会引发仪器设备强烈的振动与噪声环境,从而大大降低了系统的可靠性,甚至导致飞行失败。长期以来,国内外在飞行器脉动压力预示方面开展了大量的理论分析、数值计算与试验研究,取得了一定的成果。随着计算流体力学和脉动压力计算方法的日趋成熟,数值方法正日益成为脉动压力环境预测的主要工具。针对飞行器脉动压力环境的基本特征,对目前常用的脉动压力数值计算方法和现有研究成果进行了较为全面的介绍,对新一代航天飞行器的研制提供了一定的参考依据。