基于PIV技术的压力场重构算法实现与研究

介绍了有限容积法、直接积分法和Poisson方程法3种基于PIV瞬时速度场重构压力场的基本原理以及相应的计算方法，选取管流突扩流场和偏置方块绕流流场两个不可压缩流场的瞬时速度场数据，采用上述3种压力场重构算法，分别研究了图像噪声、速度场精度、插值算法以及边界条件的类型与精度对重构压力场的影响。最后针对管流突扩过程第20ms的流场，给出了3种重构算法下的压力场云图以及对应的CFD模拟结果。研究表明，有限容积法和直接积分法容易受到噪声的影响而产生剧烈震荡，但是可以在较大的速度场误差范围内保持较高的精度，通过采用双线性插值可以获得更高精度的重构压力场；Poisson方程法不易受到噪声的影响而产生震荡，同时在高精度PIV速度场下的优势较为突出，通过采用双三次差值可以获得更高精度的重构压力场；混合边界条件仅仅测定边界上有限个点的压力值，就获得了接近狄利克雷边界条件下重构压力场的精度，远高于诺依曼边界条件；边界条件的误差严重降低重构压力场的精度，其影响程度比速度场误差还要大。     

基于CFD的直升机旋翼流场及气动力计算

建立了一个基于结构运动嵌套网格的流场求解器,用来精确模拟复杂的旋翼流场,为更好地预测旋翼气动载荷提供一套计算方法.在该求解器中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上采用低数值耗散的Roe格式结合三阶逆风格式(Monotonic upwind scheme for conservation law,MUSCL),湍流模式采用了一方程的Spalart-Allmaras模型.应用所建立的方法,分别针对Caradonna模型旋翼、UH-60A直升机旋翼的悬停流场和7A旋翼、SA349/2直升机旋翼的前飞流场以及旋翼气动载荷进行了数值模拟.计算结果表明,本文的方法在一定程度上提高了对旋翼流场的模拟能力,进而提高了旋翼气动载荷的计算精度.     

弹射座椅系统外高速流场的数值模拟

对绕弹射座椅系统的高速流场进行了数值模拟.控制方程为雷诺平均的Navier-Stokes方程,空间离散为Osher逆风格式和五步龙格库塔格式,湍流模拟采用了基于M-SST两方程紊流模型的DES方法,计算中采用了区域分裂的并行技术.对来流马赫数分别为0.6,1.2的流场进行了计算,计算结果表明该方法可以较好地模拟弹射座椅系统外的高速流场.     

三维流场的计算干涉方法

计算干涉技术通过试验图像和计算图像相比较，可以验证CFD的有效性，是研究复杂流场的有力工具。作者根据M-Z干涉仪的光学成像原理，用计算机仿真实现了将数值计算获得的二维、轴对称和三维的复杂流体的密度场转换为“无限条纹”和“有限条纹”的干涉图像，尤其是在处理三维流场时，采用流场重构的方法，达到了很好的效果。     

用预处理方法及可压Navier-Stokes方程计算低马赫数流动

用可压方法及预处理方法求解了低马赫数紊流流动,其中空间离散格式采用了Roe格式,紊流模型采用了DES模型.以绕NACA0012的紊流流动为例,分析了该方法对提高流动速度和精度的特性,并将其用到多段翼型及绕圆柱的非定常流动中.数值模拟结果表明,该方法能有效求解低速流动及混合流场问题,且具有较高的数值精度.     

边界条件对三维空腔流动振荡的影响

基于全隐式无分裂方法,数值求解三维非定常雷诺平均的N av ier-S tokes方程,紊流模型采用修正的B-L代数模型。运用此方法,本文数值模拟了三维开式空腔的超声速、跨声速、亚声速流动。研究了空腔口处剪切层的非定常运动现象,着重探讨了入口边界条件对三维空腔流自持性振荡的影响。数值计算结果表明,与层流速度型相比,入口条件为紊流速度型和特征边界型得到的结果能够吻合。     

串列叶栅粘性紊流流场有限分析数值研究

本文将近年发展起来的有限分析方法(Finite Analytic Method)用于曲线坐标系上紊流N-S方程的数值计算。计算了来流雷诺数为2.0×10~5,气流攻角分别为0°,10°,-10°三种情况的二维粘性紊流串列叶栅流场。文中用k—ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。数值计算结果与实验结果相比较,吻合程度令人满意。     

旋翼与紊流场干扰噪声计算

本文给出了旋翼与紊流场相互作用产生噪声的理论计算方法。该方法从翼型段在紊流场中作直线运动时的噪声计算分析出发，由理论分析推广到旋翼旋转运动的情况，其中考虑了桨叶上非定常载荷的弦向及展向的非紧致性，叶－叶载荷之间的相关，特别是考虑了由旋翼运动造成的素流场畸变所引起的非各向同性和非均匀性。本文计算了一架模型直升机在＂准悬停＂状态下紊流场的畸变以及相应的干扰噪声谱，并对考虑畸变与非畸变、均匀畸变与非均匀畸变的结果作了比较。结果表明，紊流场的收缩畸变对干扰噪声有显著影响，并对此给出了合理的解释。     

环形燃烧室性能计算

在任意曲线坐标系下对包括扩压器和火焰筒在内的环形燃烧室三维两相反应流进行了数值模拟。燃烧室性能采用多维经验分析法估算所用的数学模型有 κ-ε,紊流模型 ,EBU -Arrhenius紊流燃烧室模型 ,六通量热辐射模型以及颗粒群轨道模型。在非交错网格体系下 ,气相采用 SIMPLE算法和液相采用 PSIC算法求解 ,通过两种工况和三种扩压器进口速度分布计算表明 ,所得的流场各气流参数分布和燃烧室性能较合理 ,本文所建立的程序可为燃烧室优化设计和研制提供有用数据。     

二阶矩-EBU模型用于两相反应流场的数值模拟

采用RNG k-ε紊流模型模拟紊流粘性,二阶矩-EBU紊流燃烧模型计算燃烧速率,六通量热辐射模型考虑辐射传热对两相反应流场的影响,在任意曲线坐标系下对航空发动机环形燃烧室三维两相反应流场进行数值研究。气相流场采用P ISO算法求解,液相采用颗粒轨道模型,并用PS IC考虑气液两相的耦合影响。对两种燃烧室的流场进行计算的结果与实验数据对比,表明该计算方法可靠,可为燃烧室性能的改善提供有用依据。     

一种二元高超声速进气道起动特性的尺度效应研究

对一种二元混压式进气道三维流场进行了数值和实验考察,研究了不同尺度进气道模型自起动性能的变化。结果表明,在相同来流单位雷诺数条件下,随着模型尺度的减小,进气道自起动马赫数有所提高,起动性能有所降低。同时对不同尺度模型进行雷诺数匹配,发现在相同雷诺数下,不同尺度模型的起动性能相近,表明雷诺数是影响不同缩尺模型起动性能不同的主要原因。在可获得的实验结果范围内,数值模拟所得到的自起动结果基本与之相符。此外,对实验中发现在低雷诺数下进气道反而呈现出自起动特征的异常现象进行了初步分析,通过数值模拟比较指出了低雷诺数下来流偏向层流流态,可能会导致进气道呈现一种“起动”状态。     

功率型NiMH动力电池组建模研究

选取合适的动力电池模型用于混合动力汽车控制策略设计和性能仿真.比较分析了基于等效电路的电压-内阻模型和PNGV动态模型,提出了基于测试数据的特性场模型.以功率型NiMH动力电池组为例,进行了模型参数的辨识.仿真结果与实验结果对比表明:特性场模型和PNGV动态模型具有相近的精度,但PNGV动态模型的精度与模型参数的辨识结果密切相关,电压-内阻模型简单适用于稳态仿真.     

电弧加热器超声速湍流平板烧蚀流场变化研究

电弧加热器超声速湍流平板烧蚀试验技术是研究防热材料烧蚀特性的重要手段。为研究超声速湍流平板烧蚀过程中流场变化情况,采用数值求解二维N-S方程的方法进行试验流场模拟。从模拟结果看,未烧蚀模型外形流场模拟得到的模型表面参数结果与试验结果吻合很好。然后对烧蚀过程中的模型外形进行了流场模拟,并与试验流场进行对比,根据模拟结果分析了试验过程中模型表面压力和热流密度分布的变化。根据分析可知,如果平板模型烧蚀量最大的位置在初始高热流区内,可以采用该烧蚀量计算烧蚀速率。     

NUMERICAL SIMULATION ON MULTI-FLUID INTERFACES

INTRODUCTIONThe study on phenomenon of multi- fluid in-terface driven by shock waves include:( 1 ) totrack the motion of multimaterial interfaces;( 2 )to simulate the chemical reaction on the inter-faces. As the pioneering work,chemical reactionis neglected in this paper. The objective is totrack the interfaces. Of course,the fluids aretreated as compressible medium under such ahigh impact.Although Eulerian schemes workwell for most gas flows,they have been shownto be too dissipative nea…     

超声速流动现象的水法流模拟与仿真

根据流体流动控制方程组,研究发现明渠中的浅水波和气体中的扰动波具有相同的传播特征,因而可用直观形象的水流实验来模拟复杂的气体动力学现象.运用收缩扩张形实验明渠得到的水流实验结果与气体流场数值仿真结果具有良好的一致性,验证了明渠流模拟气体动力学现象的有效性.研究结果表明:拉瓦尔喷管中的马赫数M与同尺寸明渠中的弗劳德数Fr具有相同的变化规律;明渠中的超临界流通过楔形体和圆柱体时的流动现象可以比较理想地模拟出超声速气流在相应条件下的膨胀波和激波现象;通过改变出口外部水位的方式在明渠内形成的水跃波,可以非常形象地模拟出拉瓦尔喷管在相应条件下形成的管内激波.     

从一维正投影光流场恢复二维刚体运动和结构

研究从一维正投影光流场恢复二维刚体运动和结构。证明利用一维正投影光流场中最少3个光流点可以完全确定刚体的二维旋转运动参数，一般情况下有2个解，尽管不能确定刚体的二维平移运动参数，但利用一维正投影光流场可以完全确定二维刚体的结构。文中给出了具体算例。     

NUMERICAL SIMULATION OF FLOW PATTERNS IN PLANAR JETS

INTRODUCTIONTurbulent jet flow is the most widely usedflow type in engineering applications especiallyfor combustion engines which provide most ener-gy sources for industrial production and trans-portation,and a large proportion of energy trans-portation and conversion processes is dominatedby transient evolution of large- scale structures,or so called coherent structures,in turbulentjets.As the flow control and optimization pro-cesses are closely related to the understanding ofdetailed in…     

从壁面脉冲到湍流相干结构演化的数值模拟

发展了高效求解基于Navier-Stokes方程的扰动方程高阶紧致差分方法，研究了壁面强脉冲到湍流相干结构的复杂演化过程。该方法采用高精度、高分辨率的三维耦合中心差分格式，以及和高阶时间分裂法相结合因而更为有效，比一般谱方法更适用于较为广泛的流动区域和边界条件的湍流相干结构研究。中模拟了槽道流动中壁面强脉冲诱导的湍流相干结构的产生、发展和演化过程，显示了湍流相干结构的基本特征和发展规律．计算结果表明，所用壁面强脉冲模型比共振三波模型更有效．     

短试验段内格栅湍流场风洞试验研究

为在轴向距离较短的风洞中调制出局部高湍流场,利用风洞试验的方法,对5种格栅（2种方型格栅、3种竖条格栅;3 cm × 3 cm截面铝制型材组装）后的近流场区域进行热线风速仪测量,得到格栅后近场区域的湍流参数分布及各向同性特性,依据相关湍流参数的变化规律,引入无量纲量拟合近场区域湍流强度变化规律经验公式,拟合优度为0.96。分别利用经典和现代谱估计对湍流功率谱密度进行分析,对比发现不同形式的谱估计均能准确预测近场区域湍流功率谱密度,仅在低频处存在偏差。通过改变格栅与测点距离、来流速度以及栅条结构的形式,可以改变格栅湍流场中的能量结构。     

非烧蚀防热材料地面模拟试验有效性分析

流场化学非平衡度与材料表面催化反应的耦合控制着服役于化学非平衡流场中非烧蚀防热材料表面的气动热载荷,若在该类防热材料性能模拟研究中忽略这一耦合效应,地面模拟试验将无法获得材料的有效使用性能。为此,本文依据钝头体高超声速飞行器边界层驻点热流关系式,分析了影响驻点热流的主要流场参数、地面高焓模拟设备所提供的高焓超声速流场特点以及与飞行热环境之间的主要差异,采用CFD分析了"三参数"模拟方法的有效性。针对化学非平衡边界层驻点传热分析,提出"四参数"模拟方法并分析了"四参数"模拟方法中离解焓无法模拟时的防热材料性能并提出初步解决方法。