微可压缩模型(SCM)与可压缩N-S方程数值计算对比研究

本文开展了微可压缩模型与可压缩N-S方程之间的对比研究.通过对层流平板边界层和椭球高雷诺数湍流绕流的数值模拟,比较了微可压缩模型和可压缩N-S方程之间的计算精度和效率.计算结果表明,微可压缩模型和可压缩N-S方程具有相同的计算精度,但是微可压缩模型具有更高的计算效率,内存和CPU时间都节省20%左右.     

一般状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型和高精度多介质Piecewise Parabolic Method(PPM)算法模拟一般形式Mie-Grüneisen状态方程的复杂多流体可压缩流动.根据多流体接触界面无振荡原则设计高精度计算格式,采用双波近似和两层迭代算法求解一般状态方程的Riemann问题,利用Strang分裂方法将一维算法直接推广到高维情况,不需要动能校正.二维多流体接触间断平移问题的模拟结果表明本算法在接触间断附近压力和速度没有振荡,界面数值耗散很小,被控制在两到三个网格之内.本文还模拟了二维Riemann问题和激波同多流体交界面的相互作用,算例结果表明本文方法可以有效地处理接触间断、激波等物理问题,且具有耗散小精度高的特点.     

可压缩流体热线探针校准方法研究

在可压缩流体中利用热线技术进行湍流度测量时,其输出不仅与脉动速度有关,而且也与流体温度、密度紧密相关。因此,需要在与高速可压缩流体特征相似的校准装置中,在使用前对热线探针进行准确校准。本文提出了可压缩流体热线探针校准方法,建立了热线校准(工作)数学模型,并利用自建的热线校准风洞开展了热线校准实验和湍流度测量风洞试验验证,结果表明:热线校准方法可行,校准数学模型合理可靠,热线探针校准精度较高,湍流度测量结果可信,基本可满足高速可压缩流体湍流度测量的工程需要。     

预存气体闭端管路的充填水击研究

为提升补燃循环发动机仿真平台的分布参数计算能力,针对常温液体充填预存气体闭端管路的水击现象,本文使用Modelica语言基于MWorks平台建立了一维有限体积的管路充填模型,并分别使用该模型和集中参数模型对现有文献中的充填水击实验进行了仿真模拟,结果表明：对于预存气体0.1MPa的工况下,一维有限体积模型与集中参数模型对实验的模拟度均较为接近;而在管路接近真空的工况下,一维有限体积模型的计算结果与实验数据吻合度很高,集中参数模型误差则高达76%左右。在不同预存气体压力下的计算结果表明：在管路刚度一定的条件下,预存气体压力0.1MPa以下时,气体的压缩缓冲作用是充填过程中水击压力峰值得到减弱的主要原因;在预存气体压力0.1MPa以上时,气体的缓冲作用与充填流量的减小共同导致了充填水击压力峰值的减弱。对不同弹性管路的计算表明：在闭端管路的充填过程中,管壁弹性对近真空管路的影响相对较大,对预存气体压力较大的管路几乎没有影响。     

可压缩流体恒温热线风速仪校准方法

开展了恒温(CT)热线风速仪校准方法的研究,以满足可压缩流体中湍流度测量需求。将对数函数与多元回归技术相结合,建立了恒温热线风速仪校准的数学模型,并进行了风洞湍流度测量试验验证。在马赫数Ma=0.3~0.55范围内进行了校准试验,根据建立的数学模型对试验数据进行拟合,拟合优度R2=0.999 82,平均绝对速度偏差Δu=0.18m/s,基本满足了恒温热线风速仪校准精度要求。对马赫数Ma=0.55条件下流场进行了湍流度测量,速度测量偏差为0.034%,流场湍流度Tu=0.14%。试验结果证明了恒温热线风速仪应用于可压缩流体速度测量的可行性。     

低温推进剂贮箱压力变化的CFD仿真

为预示低温推进剂贮箱在地面停放阶段的压力变化并研究贮箱内物理过程的相互作用关系,建立了包含液体推进剂和混合气体两相的二维轴对称volume of fluid(VOF)计算流体力学(CFD)模型,并引入了基于热力学平衡假设的推进剂相变模型.对实验液氢贮箱进行仿真得到的压力上升速率与实验结果相差9.1%.通过对地面加压停放阶段下的液氢和液氧贮箱的仿真发现:造成液氢贮箱压力上升的主要因素是壁面漏热对气枕的加热作用,而液氢蒸发影响更小,液氧贮箱在加压停放阶段初期明显受到液氧相变的影响.两个贮箱中液面附近的对流运动在不同的气液传热过程作用下有不同的变化趋势,对流运动会影响推进剂的相变进而影响贮箱的压力变化.      

基于燃气射流冷凝的氧路系统频率特性研究

基于热力学不平衡两流体六方程模型模拟了液氧煤油补燃循环发动机燃气射流在氧路系统泵间管路液氧流体中的冷凝过程，获得了燃气射流冷凝特征参数沿管路流向的分布规律。根据Rayleigh-Plesset方程和燃气射流特性建立了泵间管路燃气冷凝过程的传递函数模型，并与氧路管路、泵等组件模型联立求解，分析了发动机氧路系统频率特性。研究了在不同氧路入口压力和液氧温度边界条件下，泵间管路燃气射流冷凝过程对氧路系统频率特性的影响，仿真结果表明高入口压力和过冷液氧改变了射流气体的惯性和柔度，使得氧路系统特征频率增大。不同边界条件下发动机热试车结果表明，提高氧路入口压力或降低液氧温度使得氧路系统频率从8.3 Hz提高至11 Hz，与数值仿真结果一致，验证了泵间管燃气射流冷凝过程是影响氧路系统频率特性的重要环节。     

低温液体推进剂充填管路的数值模拟

研究了低温液体推进剂供应管路预冷充填过程的计算方法, 利用一维均相平衡态流体动力学模型和涵盖预冷过程中主要传热工况的传热模型, 考虑了低温液体推进剂的可压缩性, 用有限容积法求解管流方程, 用有限差分法求解管壁内的一维非稳态导热方程.计算了某型低温液体推进剂火箭发动机实验台系统供应管路的预冷充填过程, 分析比较了仿真与实验的结果, 为发动机和实验台系统的改进及新系统的设计提供了依据, 仿真结果及分析结论已应用于现有发动机实验台系统的改造和长距离液氢输送管道的设计中.      

V锥流量计火箭发动机液氢液氧推进剂测量性能

为探究火箭发动机液氢液氧低温推进剂流量测量新方法,通过数值模拟研究了V锥流量计低温流体的测量性能。湍流模型采用Realizableκ-ε模型,空化模型为Schnerr-Sauer模型,并通过自行编写UDF程序,在能量方程中考虑汽化潜热等热力学效应的影响。获得了V锥流量计流出系数和压力损失系数的变化规律,并分析了V锥流量计的测量误差。研究结果表明,存在一个雷诺数"稳定区",该区域内流出系数和压力损失系数基本为常数,液氢液氧和常温水对应的平均流出系数基本相等,且稳定区雷诺数下限值也基本相同;不同流体稳定区的平均流出系数对应的雷诺数范围差别较大,低温流体尤其是液氢的雷诺数上限值明显高于常温水。此外,空化较轻时,对流出系数和压力损失系数影响较小,当空化区域对锥尾低压口附近的压力分布产生较大影响时,则会导致流出系数迅速下降和压力损失系数增大。在稳定区对应的雷诺数范围内,液氢、液氧和水的质量流量均具有较高的测量精度,其相对误差在±0.5%之内,尤其对于液氢和液氧,其在很宽的测量范围内也可以保持较高的测量精度,空化的产生亦对V锥流量计测量精度影响较小。     

液氧煤油补燃发动机泵间管路高温富氧燃气掺混冷凝特性数值研究

针对液氧煤油补燃发动机液氧预压泵和主泵间管路富氧燃气掺混冷凝现象,建立了大过热度下富氧燃气和液氧两相流动掺混冷凝特性的全三维数值仿真方法,并以常温制冷剂R123为工质,通过气液掺混冷凝实验验证了数值仿真模型对管内两相流型和气液再液化性能的精确预测能力。仿真结果表明：弯管段气液两相在离心力作用下发生横向相对流动,强化了相间热质交换;在较低的液体流速（1m/s）下,气体水平注入管路后形成一个与气孔相连接的局部气腔,注气速率低于80m/s时,气腔一侧贴在管路内壁上,注气速率超过100m/s后气腔脱离管路内壁面。气相在气腔下端被撕裂成离散的气泡,随液体向下游流动并逐渐冷凝。在实际工况下管路的富氧燃气没有全部完成再液化过程,此时流体状态会对液氧主泵造成气蚀影响。     

可压缩流动SIMPLE算法压力修正方程边界条件处理探析

可压缩粘性流体Nivier-Stokes方程SIMPLE算法求解中压力修正方程起着关键的作用.提出从连续性方程本质出发确定压力修正方程边界条件的方法.通过对双圆弧(凸包)通道内流场的计算验证了方法的正确性.     

准一维可压缩瞬变管流的有限体积模型(Ⅰ)流场的有限体积模型

针对现有流体系统一维动力学模型缺乏统一的理论推导过程和能够涵盖众多模型的基本方程组的现象,从适用于连续介质控制体的欧拉型积分形式的流动守恒方程出发,经过合理模化推导得到了可压缩变截面管流准一维瞬变流动积分形式和微分形式的守恒方程.通过保留能量方程中的变体积项考虑了由于控制体体积的改变所输出的膨胀功.通过对守恒方程在空间交错网格上的有限体积离散建立了一种计算准一维可压缩变截面管流瞬变流场的有限体积模型,该模型可考虑变体积、变物性、轴向热传导、重力场影响,一方面,在深化和阐明了有限元状态变量模型体系的理论基础的同时,拓宽了其对瞬变流动的适用范围,另一方面,结合阀芯节流模型,从此模型出发可推导出管路系统常见元件的流场模型.      

某型涡轴发动机等压差活门建模分析

考虑某型涡轴发动机燃油流量调节器通道截面的局部能量损失,推导出流经计量油针的燃油流量方程,在此基础上详细分析等压差活门的功用、结构组成及其工作原理,以连续方程和力平衡方程为基础,采用线性化处理方法,建立等压差活门的数学模型。通过分析等压差活门系统结构,得到等压差活门系统的传递函数,进而对其稳定性进行分析,计算单位阶跃输入的稳态误差。结果表明:以流量方程和力平衡方程为基础建立等压差活门数学模型的方法可行,分析出等压差活门稳定工作条件和稳态误差的主要影响因素,并且为燃油调节器数学模型的建立以及仿真计算奠定了基础。     

高超声速湍流高效模拟算法

采用了直接过滤的Navier-Stokes(N-S)方程组对高速可压缩湍流进行研究。针对高超声速湍流的非线性流动特性,对N-S方程直接过滤推导了大尺度湍流流场的控制方程,更精确地反映高速湍流的可压缩性,建立了可压缩湍流的大涡模拟TDM模型。使用传统的Smagorinsky模型的非线性推广,采用基于非Favre过滤的超声速可压缩湍流的大涡模拟模型,应用Caylay-Hamilton定理,建立可压缩湍流大涡模拟的非线性亚格子模型,并发展为动力学模式,模型中的两个常数通过当地流场动态的确定,消除了可调经验常数的影响。针对构造的高超声速湍流大涡模型开发相应的高效并行算法。     

超声速燃烧室内氢气燃烧的三维数值研究

通过隐式耦合求解可压缩N-S方程及组分方程,对氢气的超声速燃烧进行三维数值模拟,对燃烧室内的冷热态流场,以及氢气喷射压力对燃烧的影响进行了系统分析和研究.研究表明,燃烧放热以及氢气喷射压力提高都会导致燃烧室内波系结构变化,沿周向的流向涡增强,压缩波和膨胀波与氢气核心区的相互作用增强,促进氢气和空气的掺混,更利于氢气的燃烧.      

高空试车台供气压缩机组动态仿真建模方法研究

针对国内高空试车台研究中缺少对供气压缩机组动态仿真研究方法的现状，提出了一种适用于该热力系统的动态仿真建模方法。基于守恒方程，获得了压缩机、管网、换热器、调节阀的动态仿真计算模型；将Greitzer压缩机模型应用到复杂热力系统的动态建模中；在容腔模型的基础上增加流动损失计算模块，提出了“容腔-损失”管网模型。在Matlab/Simulink平台上搭建了某试车台供气压缩机组动态仿真模型。将仿真结果与实验结果对比，验证压缩机、管网等部件动态仿真模型的可靠性。其中，压力、温度的仿真结果的平均误差均小于1%；与传统的容腔模型相比，“容腔-损失”管网模型误差明显减小。证明了该建模方法的精确度与可靠性。     

可压缩管流中真实气体效应的理论分析与数值模拟

为研究超燃冲压发动机冷却通道内复杂碳氢燃料的可压缩流动特性及相关影响因素,本文在典型工况条件下讨论了超声速膨胀流中真实气体效应的影响。基于平衡热力学理论,明确了气体动力学基础导数对可压缩流动的关键影响,并从气体动力学基础导数的角度出发考察了多种碳氢燃料在临界点附近的特性。同时,在数值模拟程序中耦合了立方型状态方程以准确描述超临界碳氢燃料的热物性参数。最后,对超临界条件下的二维超声速膨胀流进行了数值模拟,捕捉到了超声速扩张管道流动中马赫数下降的现象并揭示了其相应的物理机制。此外,参考发动机冷却通道内的典型工况条件,系统考察了不同入口马赫数、压力与温度等不同来流条件下真实气体效应的影响规律。     

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

本文采用直接过滤的Navier-Stokes方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似,但仅保留级数的一阶导数项。这样产生的误差相当于丢失了模型的部分耗散作用,本文用一种动力学非线性亚格子模型来补偿丢失的耗散影响。本文根据Caylay-Hamilton定理导出了一个非线性亚格子模型,模型中包含的系数由动力学模式确定。与传统的动力学Smagorinsky模式相比,这种动力学非线性模型(DNM)稳定性更好,且不需要在统计均匀方向进行统计平均来计算模型系数,因此减少了计算开销。本文用这种非线性亚格子模型对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验值、计算值及理论值进行了对比。结果表明,本文模型可以有效地模拟可压缩湍流流场。