基于OpenFOAM的真实流体模型库及均相求解器开发

为了正确认识高压环境下的流体非理想性特征,及其对于低温射流喷雾过程的影响作用,将多种不同真实流体模型写入到开源CFD程序OpenFOAM仿真平台中,编译形成具有广泛通用性的真实流体热物理模型库;并将高压下真实流体特征的影响引入到压力修正算法中;开发出适用于描述超临界环境下低温流体流动过程的均相求解器。对高压下流体热物理性质的非理想性,以及超临界环境下的低温液氮射流喷注过程进行研究。结果表明:高压下的真实流体效应主要体现在低温区域,本文所建立的模型库可以在任意热力学状态范围内正确预测流体热物性;基于此模型库所开发的均相求解器可正确预测超临界环境下低温射流的详细流动特征和结构,而理想气体模型则会大幅低估射流核心区密度分布,中心线上误差高达68.9%;综合考虑计算精度和效率,PR EoS及相关模型更适合于这类超临界流动问题的求解;低温射流跨临界喷注中局部区域所发生的虚假沸腾现象会导致射流温升速率降低、体积膨胀速率加快;射流核心区压缩因子明显1,在数值仿真时必须考虑高压环境下的真实流体效应。     

基于激波控制的流体推力矢量喷管试验

以二元收扩喷管为对象,开展了基于二次流喷射的流体推力矢量技术研究。基于试验研究,得到了不同喷管落压比、不同的二次流总压比和不同的二次流喷射角度多种工况下的喷管上下壁面中心线压力分布规律以及喷管壁面油流分布图。通过对不同工况下参数变化规律分析,给出了基于二次流喷射的流体推力矢量喷管的主次流气动参数及几何参数对流体推力矢量喷管流场结构和性能影响的关联关系。从试验和分析结果可以看出,喷管落压比、二次流总压比和二次流喷射角度等喷管的主次流气动几何参数对基于流体推力矢量喷管参数变化有明显的影响。     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

流体控制矢量喷管启动过程的瞬态响应

采用三维非定常数值模拟方法对流体矢量喷管启动的瞬态响应过程进行了研究,给出了流场的演化特性及性能参数的响应时间.根据不同时刻的流场和壁面压力分布,分析后得到了注气附近射流扰动波系及漩涡的演化过程,并与上壁面压力分布相互印证.由性能参数(包括流量,推力系数,矢量角)的时间响应曲线发现:流量和推力系数在响应初始阶段存在剧烈波动,而矢量角是平稳变化,三者的响应时间在毫秒量级.      

基于泵式流体回路的小卫星姿控/热控一体化执行机构设计

 主要研究了小卫星姿控/热控一体化执行机构的设计问题.首先,根据流体回路中液体流速变化对小卫星产生力矩实现姿态控制、液体流动吸/散废热实现热控制的原理,提出一种姿控/热控一体化执行机构设计方案.然后针对该设计方案,利用以电机转速为变量的流体回路内压强和电磁力矩方程,推导了一体化执行机构姿控力矩模型;利用散热量随流体回路流速的变化,建立了一体化执行机构热控模型.最后,针对某小卫星设计了基于姿控/热控一体化执行机构的闭环控制系统,并针对该一体化执行机构设计了一种姿控/热控解耦算法,对其姿控/热控能力进行数学仿真验证,仿真结果证明了该一体化执行机构的有效性.     

基于VOF方法分析离心式喷嘴结构参数对性能影响

基于两相界面追踪流体体积(volume of fluid,简称VOF)方法模拟了离心式喷嘴内部的流动过程,在数值方法可靠性验证基础上,得到了15个喷嘴构型的流量系数、液膜厚度以及雾化锥角.定义无量纲影响因子η,比较喷嘴结构参数对喷嘴性能的影响程度.除了喷嘴直径和旋流室直径外,重点研究了出口扩张角、等直段长径比等结构参数的影响,通过与试验相对比,综合以往经验公式,得到如下结论:增大旋流室长度,能够增大雾化锥角;在喷嘴出口增加扩张角,能够显著地减小液膜厚度,减小雾化锥角;增大切向口距喷嘴顶部距离,能够改善喷嘴性能;切向口个数和等直段长度对喷嘴性能影响不大.该工作为下一步优化喷嘴构型打下了基础.      

液体火箭贮箱增压排液过程温度场数值研究

针对某液体火箭贮箱增压排液过程,采用二维数值模拟方法对其温度场进行计算.选用低雷诺数k-ε模型分析流体与固壁间的耦合换热,考虑到气液之间发生热质转移现象,编写了控制相变的用户自定义程序(UDF)并植入Fluent软件.采用文献实验数据对相同工况下的计算结果进行验证,对比结果表明所建立的二维模型能够有效预测气枕温度、壁面温度沿轴向分布规律.数值模拟结果发现:气体扩散器入口方向、入口面积对气枕温度、壁面温度的轴向分布影响较弱,而对靠近增压口附近的温度场影响明显.当增压气体竖直向下进入气枕时,贮箱上封头附近气枕温度较低,有利于保障安全阀的可靠运行.当增压气体水平进入气枕时,扩散器直径变大,贮箱顶端高温区范围相应扩大.      

非交错网格求解传热及流体流动

在非交错网格系统上提出了一种新的以控制容积为基础的有限差分方法，来求解传热及流体流动问题。为了避免压力棋盘问题，压力梯度采用偏迎风压力离散格式进行处理，以保持压力－速度的耦合。利用本算法对数值传热的典型问题进行计算，并将计算结果与采用交错网格的PHOENICS通用程序所得结果相比较，二者吻合甚好。     

准则方程形式的探讨

应用积分方程求解换热问题,不仅数学运算较简单而且计算结果具有很好的直观性,容易看出理论解的物理意义,其计算结果使人们一目了然。下面我们利用积分方程,在第一类边界条件求解Pr=0.001～7500的流体绕流平壁层流换热问题。      

管内振荡流体换热研究

 本文简要总结了管内振荡流体强化换热的重要成果。对分析解,数值结果和量纲分析进行了归纳并与实验数据进行了比较。同时,讨论了强化换热的机理,认为有效导热面积和温度梯度的增加是强化换热的主要原因。这种强化传热技术在宇航工程和电子工程领域有其广阔的远大前景