一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。     

一种新颖的基于总能量守恒的化学平衡流算法

在推进剂燃烧的建模上,传统的热力计算方法一般基于总焓守恒求解定压绝热燃烧温度和平衡组分,不能考虑壁面传热;在燃气流动的建模上,通常采用的冻结流模型认为本地的组分及热物理性质与燃烧室瞬时一致,忽略了这些参数因来流气体与本网格滞留气体掺混带来的随时间的缓变效应。提出了一种新颖的可以考虑壁面传热的基于总能量守恒的化学平衡流计算方法,运用Fortran2008语言,采用面向对象编程方法建立了化学平衡流燃气发生器管道的模块化仿真模型,并将该模型应用到一个包含42个组件的涡轮试验台气路系统的建模与仿真中。与早期模型仿真结果及试验数据的对比发现,新模型的仿真结果有一定改进,更加接近试验数据。     

地源热泵技术在建筑节能中的应用

本文探讨如何通过地源热泵技术实现建筑节能。研究通过地源热泵来实现一机多用,使其既能用于冬季供暖,又可以用于夏季制冷。同时也讨论了如何与太阳能相配合,来实现地热平衡。文中还给出了通过插值法和曲线拟合法对测量数据进行处理获得变温带和增温带的温度-深度曲线的方法。     

守恒双曲型方程组的一类高效差分格式

本文将文［１］所构造的差分格式推广到方程组情形。数值实验结果表明，格式具有高效分辨激波的能力。     

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

冲压发动机尾喷管实验新方法

阐述冲压发动机尾喷管地面模拟实验的新方法。采用双爆轰技术产生稳定的高焓燃气模拟高马赫数飞行条件下冲压发动机的燃烧气体;利用皮托管测量尾喷管推力,并对测量误差进行了分析;为研究催化复合效应增大推力提供实验基础。     

CE/SE法模拟爆震波点火及传播过程

运用时空守恒元和解元(CE/SE)法模拟一维爆震波的形成及传播过程,分析产生爆震波的不同机理和DDT(爆燃-爆震转捩)过程的影响因素.采取激波点火和热点点火两种点火方式,分别对应SDT(激波-爆震转捩)和DDT过程,热点点火初始场采用线性温度分布和阶梯性温度分布,整个管道均为常压.用隐式梯形方法实现刚性源项积分.结果表明,SDT和DDT对应于不同的爆震波机理,并且DDT受点火区不同参数的影响,为研究爆震波点火和DDT过程提供了理论依据.     

有、无障碍物爆震管起爆过程的CE/SE模拟

用时空守恒元和解元(CE/SE)法模拟有无障碍物情况下爆震管二维爆震波的形成和传播过程,分析爆震波的不同机理,研究在较低能量下,点火区温度、爆震管管径以及障碍物对爆燃-爆震转捩过程(DDT)的影响.控制方程采取无量纲形式,用拓展的隐式梯形法实现刚性源项积分.研究表明,点火能量与点火压力相关,在较低的点火能量下,可以通过提高点火温度、减小爆震管径及加入障碍物加强爆震过程,减小爆燃-爆震转捩过程DDT时间和距离.为研究爆震波点火和DDT过程提供了理论依据.      

一种TVD方法在三维非结构网格中的应用

以三维非结构网格的显式有限体积法为基础,采用了一种TVD方法求解三维Euler方程,使用Roe通量来计算网格单元边界处的守恒量通量.为了验证方法的可行性,用该方法模拟三维爆炸问题,得出的结论是我们的方法可行,稳定且有效.