连续爆轰发动机多波面现象的三维并行数值模拟

来源 :北京大学 | 被引量 : 0次 | 上传用户:xpzcz1987
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
针对连续爆轰发动机实验研究中出现的多波面现象,本文通过三维数值模拟对连续爆轰发动机不同进气方式条件下燃烧室内流场进行比较,观测波面产生的过程,并对其生成的原因作出解释,为实验研究提供参考。同时为了今后进行燃料、氧化剂和燃烧产物混合程度对连续爆轰发动机燃烧室内爆轰波传播的影响,尝试了使用基元化学反应对二维连续爆轰发动机进行了模拟,并对结果进行了分析和评价。此外,由于三维和二维数值模拟计算量巨大,本文编写了基于MPI+CUDA的并行程序,以达到快速高效计算的目的。  本文主要内容如下:  一、综合叙述了文章所使用的计算格式和技术。详细介绍了数值模拟所采用广义坐标系下控制方程的形式,叙述了WENO差分格式的使用理由,列出对应的矢通量分裂和3阶Rung-Kutta时间步迭代的表达式,说明了文中使用的各种边界条件。此章还对一步化学反应和基元化学反应进行了说明,全面给出了本文使用的化学反应模型的表达式和参数,以供参考。最后还对基于GPU的CUDATM架构进行了介绍,说明了由于硬件特点而在CUDA编程中所需注意的问题。  二、初步探索了利用基元化学反应模型对连续爆轰发动机进行数值模拟。本文首先通过一维基元化学反应模型模拟验证了模型和程序的正确性,然后利用二维矩形流场进行了燃烧室的近似模拟。模拟结果表明,若直接使用真实发动机参数进行模拟,模拟结果与事实相差很远,即基元化学反应模型对具有较高压强的混合燃料模拟会出现严重失真的现象。讨论认为基元化学反应对连续爆轰发动机的模拟存在不足,若采用基元化学反应模型分析流程需要做折中的处理。  三、介绍了本文数值模拟的CUDA并行程序的编写。程序根据函数的数据结构特点对不同的函数进行了不同的优化。本文形象刻画了优化算法,并对GPU单卡、多卡在加速性能上的体现进行了比较。研究证明GPU利用CUDA编程可以大幅提升程序运行速度。  四、对传统进气模式下多处点火产生的多波面现象进行研究。利用一步化学反应模型进行连续爆轰发动机燃烧室的模拟,通过设定均匀分布的多个初始点火点,人工建立多波面起爆的初始条件,并分析此后流场的变化。影响波面个数上限和下限的因素通过流场各个时间点上的发展变化得以阐述,并体现了不同数量波面达到稳定时间上的差异,得出对于给定尺寸和进气条件的燃烧室,存在一个最易稳定的波面个数的结论。  五、研究不同进气方式对燃烧室内流场的影响。文章在传统进气方式的基础上建立了四种新的进气方式。通过对各个进气方式下流程的分析,得到不同进气条件对流场截然不同的影响,并初步解释了燃料混合气体与燃烧产物相互作用对流场的影响。在四种进气方式中,其中的放射间隔进气方式得到了实验中多波面自发形成的现象。在对四种模拟结构和传统进气方式数值模拟的比较之后,初步得到了燃烧室内可以存在稳定爆轰波传播和不稳定爆轰波传播的现象。  六、对多波面自发形成现象进行深入讨论。本文首先对一个形成稳定多波面传播的算例进行了详细分析,根据流场的变化将点火后波面形成的过程分阶段进行讨论。数值研究结果展现了点火中期阶段新波面反复生成和消失的过程,解释了新波面产生的成因。其次将运算结果与他人实验结果进行比较,解释了实验结果中的一些现象,并证明了数值模拟的高逼真度。最后通过改变算例的条件,探讨燃烧室条件对新波面形成和波面稳定传播的影响关系。讨论认为波面的形成对各种条件都很敏感,波面个数、传播方向和稳定性与燃烧室各个条件没有定量的关系,同时再次验证了不稳定爆轰波传播属于爆轰发动机正常的工作状态。  
其他文献
该文主要研究了陶瓷衬管在热镶嵌及射击过程中,产生端部效应的原因及相应的2解决方法,对AlO及SiC陶瓷衬管进行了热镶嵌试验,并通过射击实验考核,采用X-光应力分析方法进行了
由于大量不规则结构的出现,使得结构抗震性能不能满足规范要求。为了解决这一问题,便提出了消能减震技术,该技术能够使不规则结构满足结构抗震性的要求。本文首先介绍了介绍
该论文根据采场实际情况,利用渗流理论和管路漏流流动规律、线性弥散定律和质量守衡定律分别建立了采场风流流动和采场瓦斯浓度分布数学模型.通过伽辽金有限元法离散了所建立
产权风险,主要指由于企业的产权安排不合理或不能适应企业发展,而给企业带来的治理风险。主要表现有:“一股独大”造成的强权制度、股权过于分散造成的高管人失控、关联交易
利用代数图论、稳定性理论和线性矩阵不等式(LMI)方法,研究了连续多智能体系统的一致性问题和包围控制问题。由于时滞现象在实际环境中广泛存在,对具有可变时滞多智能体系统的
既有民用建筑结构的加层改造,便捷且经济实用,其中以既有RC框架结构钢结构加层传力明确、构造简单、经济实用而较为常见。然而,加层后,既有结构与加层体系的协同作用机理(地震响应
随着电力企业精益化管理的深入和用户对电力供应可靠性要求的提高,电网运行已经由原来“有电可供”的粗放式管理向目前“提高经营效益、电能质量”的集约化方式转变。电力系统无功电压是电网重要的研究指标之一,关系到电网稳定、电能质量及运行经济性。因此,开展无功电压优化研究是电网运行管理单位共同面临的重要课题之一。建立以系统有功损耗最小为目标的模型,采用原-对偶内点法进行求解,对于计算机求解则按照给定的误差和步
本文以内蒙古电力公司鄂尔多斯电业局于2010年上线运行的用电信息采集与监控系统为写作背景。通过对鄂尔多斯电业局用电信息采控系统的监视运行,为今后进一步加大开发力度全面推进智能化电网建设奠定基础。用电信息采集与监控系统的建设是采集客户端实时用电数据的一个平台,是不断提升营销管理水平,堵漏增收、降损增效的重要举措。本文通过对用电信息采集与监控系统设计原则和思想的分析,对系统软件和硬件部分的深入研究,最
写作水平的高低,作文成绩的好坏,直接影响学生的语文素质和语文成绩.因此,写作教学成为中学语文教学中最重要的环节.而当前在语文教学中教师往往只是强调文本阅读,忽略了培养
为了将弹塑性断裂力学从确定性分析推广至随机分析以解决工程中延性金属结构可靠度计算问题,该文借助前人的"弹塑性断裂力学工程方法"和工程随机分析理论,建立了弹塑性断裂随