基于基元反应的总包机理建模及算法优化

以煤油/氧的化学动力学体系为例,介绍了基于基元反应的总包机理建模、反应速率计算等国际先进方法,归纳了针对化学动力学过程由建模到求解的一套完整流程。首先选用正十二烷作为煤油的替代燃料,通过直接关系图法(DRG)结合计算奇异值摄动法(CSP)的机理简化方法将正十二烷203组分、738步的详细反应机理化简为了32组分、36步的总包反应机理。为了计算总包反应速率,采用线性准稳态假设(LQSSA)计算准稳态组分浓度,解耦了准稳态组分间的非线性耦合,在保证计算精度的同时提高算法稳定性和效率。最终,通过VODE刚性积分求解器求解化学动力学方程组,并形成了一套可与CFD主程序耦合使用的化学动力学计算子程序。     

基于划分检测模型的终端区异常轨迹检测方法

针对异常轨迹对轨迹聚类效果的影响以及给进离场程序管制适用性的量化分析提供分析数据的考虑，在结合数据挖掘中异常检测理论的基础上，利用划分检测框架及将距离和密度结合起来，设计了一种有效检测算法，运用到终端区飞行轨迹的异常检测当中，检测出飞行轨迹中的异常轨迹，从而改善轨迹的聚类效果。实例仿真结果表明：方法可以准确地检测出异常轨迹。     

A Parallel Adaptive 3D MHD Scheme for Modeling Coronal and Solar Wind Plasma Flows

A parallel adaptive mesh refinement (AMR) scheme is described for solving the governing equations of ideal magnetohydrodynamics (MHD) in three space dimensions. This solution algorithm makes use of modern finite-volume numerical methodology to provide a combination of high solution accuracy and computational robustness. Efficient and scalable implementations of the method have been developed for massively parallel computer architectures and high performance achieved. Numerical results are discussed for a simplified model of the initiation and evolution of coronal mass ejections (CMEs) in the inner heliosphere. The results demonstrate the potential of this numerical tool for enhancing our understanding of coronal and solar wind plasma processes. This revised version was published online in June 2006 with corrections to the Cover Date.     

针栓式喷注单元膜束撞击雾化混合过程数值模拟

为了全面认识针栓式喷注器喷雾场结构，基于自适应网格加密技术和分三相计算的PLIC VOF（Piecewise Linear Interface Calculation Volume of Fluid）方法对针栓式喷注单元膜束撞击雾化混合过程进行了仿真分析，通过对两路推进剂分别进行界面追踪，获得了膜束撞击雾化混合过程的详细结构特征，与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。以此为基础对膜束撞击的喷雾场结构、撞击变形过程、流场涡结构、雾化破碎典型特征及破碎后的雾化混合分布特征进行了识别分析，结果表明：膜束撞击形成了液束未穿透液膜和液束穿透液膜2种不同的喷雾扇结构。膜束撞击形成的喷雾扇呈"Ω"形，膜束同时发生弯曲变形和横截面变形。另外，膜束撞击同时受到正压和剪切应力作用，导致了一系列复杂涡流现象，使得相互作用增强，雾化混合均增强，这也是膜束撞击喷注构型优于膜膜撞击的本质原因。最后，还发现膜束撞击喷雾场液滴分布呈现分区结构特征，分别是液束控制主导的上雾化区、液膜控制主导的下雾化区及夹在中间的混合区，实际中应兼顾雾化特性和混合特性，选取中等动量比膜束撞击，这可为针栓式喷注器的理论研究和工程设计提供重要参考。     

目标函数误差估算及网格自适应处理

发展了一种基于目标函数误差估算的网格自适应准则,进而通过网格自适应处理提高目标函数计算的效率和准确性。首先描述了目标函数的误差估算及修正的方法,该方法通过伴随方程将原方程的残值误差与目标函数联系起来。然后,基于误差估算建立网格自适应准则,以减少目标函数修正后的剩余误差,提高目标函数计算结果的准确性,并将方法进一步发展至多目标问题。最后将该准则应用于Euler方程控制的NACA0012翼型无粘可压流动的网格自适应模拟。数值计算成功地捕获了与升力、阻力和力矩等目标函数相关的特征流动区域,计算出符合指定精度要求的目标函数值,验证了本文所发展的方法。     

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

Transient simulation of regression rate on thrust regulation process in hybrid rocket motor

The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB(hydroxyl-terminated polybutadiene) hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.     

PIV技术在复杂二相流场中的应用

光学元器件随飞行器在大气中飞行时,其工作性能越来越多地受到大气悬浮汇聚微粒的影响。大气微粒在复杂流场中呈现何种运动汇聚效应,对于合理准确评估机载光学元器件的工作效能具有十分重要的工程意义,而复杂气动流场中微粒分布状态的预估一直是飞行器外界环境研究中的一个难点。气动问题的复杂性、大气中微粒的多样性一直是制约各种试验手段展开、数值模型建立的主要因素。利用先进的激光粒子图像技术,在风洞中对舵面旋涡主导的复杂流场中的微粒速度及分布特性进行了实验研究。在测量舵面翼梢脱落旋涡特性的基础上,通过激光片光扫描流场全域,同时高帧频CCD相机同步曝光,利用PIV 拍摄到的流场中涡流截面内微粒分布的瞬态图像。结合图像后处理技术,对原始粒子图像进行互相关、二值化处理,通过对图像区域内的灰度值计算,统计相对流场截面内的粒子浓度系数,得到在复杂旋涡结构流场内瞬态粒子的分布特性规律。研究结果表明,利用大气中微粒在激光片光下的米氏散射原理,可以有效地拍摄到复杂流场结构下粒子光学散射及分布的特性图像,解决了传统环境测试设备无法对复杂条件下流场内粒子分布进行实时测量的缺陷;在旋涡为主导的流场中,大气中的微粒由向心力牵引,在涡核周围达到平衡运动状态,微粒环绕涡核形成一条环状带,这一区域中的粒子浓度系数要远大于自由流场中的微粒,涡核中心粒子呈“空洞”状态。     

轻质隔热纳米孔结构耐烧蚀酚醛材料

在一定温度下,通过改性剂组分与热固性钡酚醛树脂之间热共聚反应,制备了改性酚醛树脂。采
用模压工艺制备出具有目标尺寸和形状的改性酚醛树脂固化物,对该样品分别进行了400 和700℃的炭化处
理,得到了相应的酚醛炭化样品,并测试其密度、失重率、收缩率等特征参数。研究改性剂用量对碳化物特征参
数和性能的影响,如密度、失重率形貌、比表面积、隔热性能、压缩强度等。SEM 研究表明,700℃炭化后样品形
成了“珊瑚状”碳结构,形成了30 ~150 nm 纳米孔。比表面积测试(BET)结果表明其比表面积可达464 m2 / g。
隔热性能与力学性能的测试结果表明,相对于未改性酚醛树脂碳化物,改性酚醛碳化物的热导率降低了
45． 57%,并且保持了较好的压缩强度。     

三维非对称斜劈高速流动计算方法研究

交叉激波是高超声速进气道发展的关键气动问题之一。为实现对此类高速流动的准确数值模拟,基于RANS方程的自研多块结构网格计算软件,对经典三维非对称斜劈构型的斜激波交叉干扰流动特征开展了计算研究。通过对比分析,系统评估了计算网格、湍流模型及高精度格式对空间流场、底板摩擦力流线、压强分布、绝热壁温分布等的影响,获得了适用于三维非对称高速流动问题的网格生成准则和计算策略,同时验证了自研软件的模拟能力。