战术导弹横向喷流数值模拟

为了研究喷流对导弹不同升力面的干扰影响,通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格.通过几种升力面和弹身组合的超声速横向喷流干扰流场的数值模拟,计算的马赫数为2和4.5,计算结果显示具有升力面的喷流干扰比光弹身的喷流干扰大.研究和分析了喷口附近流场的涡系结构和波系结构,并将导弹喷流干扰力放大因子及喷流引起的压心变化的计算结果与目前国外公布的试验值进行比较分析,二者显示了良好的一致性.     

轴对称羽流流场数值模拟和实验验证

针对在轨航天器羽流效应对航天器整体寿命的影响,在"DFH"系列卫星姿控发动机喷管羽流实验研究的基础上,建立喷管的轴对称模型,选取矩形计算区域,设定一定的边界条件和计算条件,采用轴对称单组分直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法对相应实验状态的羽流流场进行数值模拟研究.通过对比数值模拟结果和实验测量结果,表明DSMC数值模拟羽流流场的计算结果符合实际流动规律,实验验证了所编制的轴对称DSMC数值模拟程序的正确性和有效性.      

煤粉燃烧过程的数值模拟

采用计算流体动力学软件对煤粉实际燃烧过程进行数值计算并结合其热态试验数据进行对比分析,已成为验证数学模型和指导工程实践的一种重要研究手段.应用商业软件FLUENT对煤粉燃烧及流场进行了数值模拟分析,采用混合分数/概率密度函数法模拟湍流燃烧,用P-1辐射模型开展辐射传热模拟,利用拉格朗日/欧拉法处理气固两相间的动量、质量和能量交换,对挥发份的析出采用双速率竞争模型,采用动力/扩散反应速率模型模拟煤粉颗粒的表面燃烧.并对模拟结果进行分析与对比,计算结果与实际燃烧过程有较好的一致性.数值模拟计算结果揭示了挥发分释放与焦炭燃烧的过程,为改善和优化煤粉的燃烧提供了重要的参考依据.     

微重力量值对气液两相流相分布及液相湍流统计量影响的数值研究

气液两相流在空间领域具有广阔的应用前景, 深入理解微重力量值(微重力大小)对相分布和液相湍流的影响十分必要. 采用欧拉elax-elax拉格朗日双向耦合模型深入研究了不同微重力量值对相分布和液相湍流的影响. 液相速度场通过直接数值模拟求解, 气泡的运动轨迹由牛顿运动方程跟踪. 研究表明, 气泡分布和液相湍流与微重力量值均具有直接联系. 在低微重力量值下, 气泡近似均匀分布在槽道内, 且对液相湍流统计量几乎没有影响; 然而当微重力量值较高时, 大量气泡聚集在壁面附近, 液相湍流由于气泡的注入受到极大调制.      

两种层板性能对比

针对常规的141型层板,提出了一种增加50％数量扰流柱的设计,即161型层板.采用实验和数值模拟相结合的办法研究了两种类型层板在流阻和换热方面的差别.换热实验部分采用了辐射式加热设备,测量了层板上下表面平均壁温,以及冷气进、出口截面间的焓差,得到了相应进气雷诺数下的层板体积换热努塞尔数.发现161模型流阻降低20％,换热增强5％.对一些实验工况,采用三维计算流体力学程序进行了流-固耦合传热数值模拟,所得结果在趋势上与实验一致,在数值上,流阻和换热与实验结果分别相差5％和30％.利用数值模拟结果分析并比较了两种模型在流场和表面对流换热系数分布方面的细节差别.      

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

歼击机座舱空气流动和传热模拟实验

对某型歼击机座舱的空气流动和传热进行了全尺寸地面模拟实验.模拟了气动热载荷和人体热载荷,并使用ATM2400型36通道多点气流和温度测量系统测量了5个测试面的速度场以及其中两个测试面和邻近服装表面的温度场.观察了供气流量和集中喷嘴开关对舱内空气流动状态和温度分布的影响,评价了该座舱的舒适性.结果表明:舱内速度场 并非相对座舱纵向中心面对称分布,而座舱内壁面温度分布亦极不均匀;集中喷嘴全开和全关时,舱内气流速度场和温度场及邻近服装表面空气温度的分布情况差别显著,但对舱内气流平均速度和平均温度的影响不大;舱内气流平均速度随供气量呈线性变化.实验结果可为数值模拟提供验证依据.     

低速流预处理Roe格式中的数值粘性

使用Roe格式计算低马赫数流动,需采用预处理技术.在预处理技术中,需要选择适合的预处理参数,但有时很难找到,因而得不到合理的计算结果.通过数值实验,发现使用预处理方法时,格式的数值粘性对低速流动的模拟有显著的影响.提出了在预处理Roe格式中加入可调参数,用以控制数值粘性对计算结果的影响.带有可调节数值粘性参数的预处理Roe格式的数值粘性更接近低速流动的真实耗散,反映了低速流动的真实物理现象.对方腔流动和圆柱绕流的低雷诺数定常流动所做的数值实验,得到了比较满意的结果,表明这一改进是有效的.对同一低雷诺数不同低马赫数的定常圆柱绕流所做的数值实验,均得到了满意的结果,表明了方法的适用性.     

气泡雾化喷嘴水平喷射的雾化特性研究

研究了水平喷射、端面注气方式的气泡雾化喷嘴结构及工作参数对雾化及流场特性的影响.试验全部在常温常压下进行,液体采用水,雾化气为压缩空气.用2D PDA测量了液雾的平均直径尺寸分布和速度分布.液体喷射压力变化范围200~600 kPa,气液比变化范围4%~10%左右.试验研究了端面注气方式下,注气孔孔径及数目、混合管长度、液体流动状态及工作参数对雾化特性的影响.结果表明,水平喷射的气泡喷嘴,注气孔的尺寸及数目均对雾化特性产生影响;混合段存在一最佳值范围,在此范围内,喷嘴可获得高的雾化质量;液体旋转流动对雾化特性无显著影响,但可影响两相流动中气泡的分布.     

气泡雾化喷嘴的试验研究

进行了组不同结构形式的气泡雾化喷嘴的试验研究，初步试验结果表明：这种新型喷嘴对液体介质粘度不敏感，在小气液比条件下获得了很好的雾化效果，雾化性能明显优于从国外引进的某些空气雾化喷嘴，该成果已于１９９３年２月７日被授予中华人民共和国实用型专利。     

后台阶流动的Hybrid RANS/LES模拟

采用雷诺平均的Navier-Stokes和大涡模拟混合的方法对后台阶流动进行了数值模拟,并与实验结果进行了对比;分区尺度采用IDDES(Improved Delayed Detached Eddy Simulation)模式:来流包含背景湍流度时,混合方法转化为WMLES(Wall Modeling in LES)模拟,除了近壁区外的大部分区域都是直接求解;不包含来流湍流度时,转化为DDES(Delayed Detached Eddy Simulation)模式,附面层采用纯RANS模拟,大分离流动区用类脱体涡模拟方法模拟;对流通量的离散采用迎风型低耗散通量分裂格式和中心型四阶偏斜对称格式组合的模式;湍流模型采用Spalart-Allmaras(S-A)模型;亚格子长度尺度取决于网格分布和到壁面的距离.      

气液两相流中气泡速度的电磁互相关测量

气液两相流作为一种工业领域中普遍存在的流态,在测量时需考虑气相和液相间的速度差。相较于流动规律平稳的气相而言,研究不稳定的气泡流动状态具有更高的技术难度。本文针对空气-水为工质的两相流中一项重要参数——气泡速度进行测量,引入基于互相关运算的电磁检测方法,通过提高系统激励频率得到较大的检测信号。实验设计了高频(1 MHz)电磁检测系统,选用垂直上升管道,并在管道的2个平行截面分别安装电磁传感器,每组传感器均包含一个激励线圈及一个接收线圈,对2组接收线圈上的相位信号进行互相关运算,从而求得气泡的速度。实验中,对3种不同速度气泡采集到的信号进行了对比分析,相对误差控制在10%以内,提供了一种完全非接触非侵入的测量气泡速度的方法。此方法可以进行后续补充完善,用于其他工业场合,如金属液中气泡参数的测量等。     

罐车制动时液体晃动的仿真分析

采用Volume of Fluid (VOF) 模型对罐车制动时液体的晃动进行了数值模拟,对单室受力、受力位置及整车轴荷分配进行了计算,并通过与实验结果的对比验证了计算方法的可靠性.计算结果表明,无防波板时,随减速度增加,单室 x,y方向受力峰值增大,整车轴荷比增大;随充液比增加,单室x方向受力峰值先增大后减小,y 方向受力峰值增大,制动初始与结束时充液比越大轴荷比越大,1s左右充液比越大轴荷比越小;单室带防波板时,随防波板面积增加, x,y 方向受力峰值减小,当防波板面积大于横截面的40%时,增加防波板面积能显著改善罐体受力,且防波板面积越大轴荷比峰值越小.     

汽油机进气歧管流通性和均匀性数值模拟

进气歧管的流通性和均匀性对发动机性能产生重要影响.建立了某四缸汽油发动机的GT-power仿真模型,对进气歧管的气体流速进行了仿真计算;利用STAR-CD软件模拟计算了进气歧管不同进气角度时各支管内的气体压力和流速分布,分析了各支管的进气均匀性和充气系数.计算和分析结果表明,通过改变进气口方向和稳压腔容积等几何参数,可以提高进气歧管的气体流通性,改善进气歧管各支管内气体流动质量,使各支管内气体流量趋于均匀.计算结果为该进气歧管的改性设计提供了理论依据.      

失重条件下低温推进剂排放计算

针对CZ-3A系列运载火箭三子级低温贮箱内剩余推进剂排放问题,分别对排放过程中液氢和液氧低温推进剂沿2种不同排放管的流动进行了一维数值计算和分析.计算中认为在排放时沿管路压强降低,从而引起液体推进剂气化,并在管路的某个位置出现气固两相或者全部气化,在计算中,根据不同位置时推进剂的不同状态,分别采用液相、气液两相、气固两相及气相的物理数学模型和控制方程.计算得出了排放的出口参数及排放性能,认为液氢从通道1排放比较好,液氧从通道2排放比较好.为进一步进行排放研究提供了依据.      

基于CCERT与声发射技术的气液固三相流相含率测量

传统的相含率测量方法在测量三相流相含率时只能对某单一相进行检测,针对此问题,基于电容耦合电阻层析成像(CCERT)技术与声发射技术,建立了相含率测量模型,提出了一种三相流各相相含率的非侵入式测量方法。利用偏最小二乘回归法,在静态情况下建立CCERT技术的相含率测量模型,同时在鼓泡床上进行动态实验,采用差压法进行同步测量作为参考值验证模型的有效性,实现了对两相相含率的非接触测量。在此基础上,通过对声发射技术采集到的声音信号进行处理,建立声发射技术的气相相含率预测模型,利用该模型测量出三相体系中的气相相含率,结合CCERT技术测量出三相体系中的不导电相相含率,从而利用非侵入式测量方法得到三相流中的各相相含率。     

固体火箭发动机燃烧室三维流动数值计算

采用k-ε双方程湍流模型,以SIMPLE计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流N-S方程.在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构.计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡.