超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

爆炸波与物体干扰流场的数值模拟

本文构造一类新型的计算效率较高的ＥＮＯ差分格式。利用该格式非定常Ｅｕｌｅｒ方程出发,对爆炸波扫过物体时形成的二维和三维干扰流场进行了数值模拟,计算结果与实验数据符合较好。     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

高速飞行弹丸诱导的三维爆轰波数值模拟

采用三维非定常化学反应控制方程组对驻定在带攻角高速飞行弹丸上的斜爆轰波流场进行了数值模拟。对流项和化学反应项采用时间分裂法,分别用TVD格式和基元反应模型进行处理。计算结果表明三维数值模拟可以更准确的描述弹丸飞行攻角对爆轰波的形成及结构的影响,是进一步研究超驱爆轰推进技术的重要手段。     

一种混合格式的构建及其在多尺度流动中的应用

面向非定常湍流数值模拟方法精度提升问题,构建了实用的矢通量分裂(FVS)形式高阶精度混合格式。该格式平滑区具有五阶精度线性紧致格式的高分辨率特证,而强间断和强梯度区呈二阶耗散型格式特征,进而保证全流场计算的稳定性。基于所建立的混合格式开展有限体积大涡模拟(LES)计算,将其应用于三维高Reynolds数非定常湍流预测中。结果显示,所建立的数值格式能够可靠预测湍流多尺度涡结构特征,且计算量小,具有较高分辨率和鲁棒性。     

高超声速电离绕流的数值模拟

本文从三维N－S方程出发，采用稳式L－U分解算法和7组分15个反应的化学模型，数值模拟高超声速电离空气绕流。首先采用对称TVD格式、AUSMPW＋格式和Van Leer的矢通量分裂格式计算了高超声速球头绕流，并对它们的计算结果做了对比分析。然后用前两种格式，对RAM－C飞行试验模型三个再入高度（81km、71km、61km）的流场进行了数值模拟，计算的流场电子数密度值和试验测量数据符合较好。     

烧蚀外形三维粘性流场数值模拟

本文通过求解层流Navier-Stokes方程,数值模拟了三维奶头型烧蚀外形粘性流场。差分格式为Yee和Harten的迎风隐式TVD格式,在空间方向具有二阶精度。计算得到了高质量的激波和合理的流场结构。物面压力分布与实验结果进行比较,符合较好。     

塞式喷管流场数值模拟

采用ＮＮＤ２Ｍ差分格式,从考虑к－ε湍流模型的理想气体三维薄层近似Ｎ－Ｓ方程出发,对塞式喷管内外流场进行了数值模拟。计算结果与国外文献十分一致。通过流场分析,对塞式喷管流动机理有了清楚的认识。计算表明ＮＮＤ２Ｍ差分格式及相应计算软件适合于模拟复杂喷流流场,可用于塞式喷管发动机的研究。     

飞机三维结冰模型及其数值求解方法

 对飞机结冰外形进行纯三维数值模拟,是目前飞机结冰预测中的一大难点。为了建立三维结冰数值模拟方法,基于Messinger的二维结冰模型,提出了一种考虑液态水溢流效应的三维结冰计算模型,并针对该模型建立了表面单元内溢流水流动的分配方案、发展了相应的迭代求解方法。采用本文方法对MS-317后掠翼结冰进行了计算,并与实验和Lewice3D的计算结果进行了对比。研究结果显示：本文方法具有较好的收敛性,计算的霜冰和明冰外形均与Lewice3D计算的冰形一致;对于霜冰和结冰时间较短的明冰,本文计算的冰形与实验吻合较好;对于结冰时间较长的明冰,本文计算的冰形与实验对比还有一定差异,但冰生长的总体趋势和大致体积与实验一致。同时,对比了MS-317无后掠翼与有后掠翼结明冰的外形,发现机翼后掠导致的三维溢流效应对结冰外形有明显影响,因此,对于三维结冰分析,如果用二维截面的结果来代替三维结果,其合理性还需验证。     

基于锅炉炉膛三维温度场检测算法的研究

鉴于大型锅炉内燃烧火焰变化的多样性,为了实时掌握炉内火焰燃烧状况,提出了一种改进的检测炉内三维温度场的算法.该算法通过对炉膛火焰二维累加辐射能的图像处理,并利用辐射传递方程的逆求解,实现了对炉内火焰三维温度场的确定.仿真研究和实验结果表明此算法是可行的.     

液体火箭发动机推力室再生冷却流动与传热计算研究

为了研究液体火箭发动机推力室再生冷却流动与传热的快速仿真方法,建立了推力室再生冷却的准二维模型,对航天飞机主发动机开展了再生冷却流动与传热计算仿真研究,对比分析了再生冷却准二维模型和三维模型的仿真计算结果。研究表明,两种计算模型均可较好地预测推力室燃气及再生冷却剂的流动和传热。三维模型计算精度高,但计算用时较长。计算得到的航天飞机主发动机的燃气侧壁面最高热流密度为162.2MW/m2,最高壁温为1159.7K,冷却剂温升为244.0K,压降为8.5MPa。准二维模型计算结果精度略有降低,但计算时间较三维模型减小了90%。四个参数与三维模型计算结果的差异分别为0.3%,4.4%,8.6%和4.5%,在可接受范围内。本文的准二维模型计算时间短,适用于液体火箭发动机再生冷却结构的方案筛选和优化设计,三维模型计算精度高,适用于设计完成后的校核计算。     

混合动力复合翼应急迫降在线航迹规划与制导

针对混合动力复合翼飞行器巡航模式下空中停车后无动力应急迫降（VTOL）问题,提出在线航迹规划与制导方法。根据复合翼空中停车时初始位置/航向不确定散布,发展一套满足动力学约束、终端约束的三维航迹在线规划方法：利用几何规划方法快速生成扩展Dubins二维航迹,再根据下滑性能约束进行三维扩展。针对低速无动力下滑航迹跟踪更易受风干扰以及三维航线分段连接处曲率不连续的特性,发展一种基于非线性模型预测控制的三维制导算法。将纵横解耦的制导律嵌入到预测模型框架内,跟踪误差、外界风扰动、航迹曲率不连续等非线性因素则通过系统输出建立目标约束,其后利用滚动优化实时求解制导指令。最后对航迹在线规划方法与三维制导律的适用性进行仿真分析与验证,结果表明所提出的航迹规划方法适用于不确定初始位置/航向散布的应急迫降在线规划,所提的制导算法具备抵抗风扰、提高三维制导精度的能力。     

跨声速风扇的流固耦合的颤振分析(英文)

发展了时间推进的叶片颤振的数值方法,采用了时间推进求解流体和固体相互作用的过程。气动模型是基于求解三维雷诺平均的N-S方程,采用了多块结构化网格的有限体积格式,对流和耗散通量的计算使用了二阶迎风格式和中心格式。在变形的动网格上流体运动的守恒型方程的求解采用双时间步,隐式格式和多重网格方法。叶片振动采用了振型叠加的线性气动弹性模型。该方法在气动弹性标准算例4进行了验证,并用于求解跨音速风扇的颤振问题。     

固体火箭发动机体空间缺陷最大直径矢量特征测量方法

针对固体火箭发动机工业CT(computed tomography)三维扫描数据,从计算时间和测量精度两方面考虑,结合固体火箭发动机内部缺陷体空间数据场的特征,通过改进传统的空间最大距离求解法——擂台法,提出了基于分类种子点法的体空间缺陷最大直径矢量特征测量方法.设计了预置缺陷的固体火箭发动机,经实验验证,相比传统擂台法,该方法能够提高测量精度和缩短计算时间,最大直径及其与轴向锐角夹角的最大测量误差在10%以下,为固体火箭发动机三维可视化故障诊断奠定了基础.      

Range-Doppler Imaging of Rotating Objects

During the integration time required to obtain fine Dopplerfrequency resolution in a range-Doppler imaging radar, a point on a rotating object may move through several range and Doppler resolution cells and produce a smeared image. This motion can be compensated by storing the appropriately processed return pulse, and the angular coordinates are determined by the angular coordinates of the radar antenna. The resulting stored data represents the three-dimensional Fourier transform of the object reflectivity density, and hence can be processed by an inverse Fourier transformation. Also included is an analysis of the three-dimensional radar/object geometry with separate source and receiver locations. The effects of various system aberrations are investigated and experimental results from a microwave test range which demonstrate the image improvement are presented.     

不定常差分方法中的二阶三层格式

本文介绍一种求解三维理想可压流问题的二阶三层格式并给出它的稳定性证明。通过大量二维和三维问题的实际计算,表明此方法的优点是精度好,时间省,公式简单,程序量少,便于实际应用。文中对应用此法如何求解间断流场的问题,针对具体算例做了较详细的讨论。对不同外形在不同流动条件下,如何设计差分网格也做了分析。最后给出部分复杂外形的二维和三维计算结果。     

Three-dimensional guidance law based on adaptive integral sliding mode control

For the terminal guidance problem of missiles intercepting maneuvering targets in the three-dimensional space, the design of guidance laws for non-decoupling three-dimensional engage-ment geometry is studied. Firstly, by introducing a finite time integral sliding mode manifold, a novel guidance law based on the integral sliding mode control is presented with the target acceler-ation as a known bounded external disturbance. Then, an improved adaptive guidance law based on the integral sliding mode control without the information of the upper bound on the target accel-eration is developed, where the upper bound of the target acceleration is estimated online by a designed adaptive law. The both presented guidance laws can make sure that the elevation angular rate of the line-of-sight and the azimuth angular rate of the line-of-sight converge to zero in finite time. In the end, the results of the guidance performance for the proposed guidance laws are pre-sented by numerical simulations. Although the designed guidance laws are developed for the con-stant speed missiles, the simulation results for the time-varying speed missiles are also shown to further confirm the designed guidance laws.     

三维流场的干涉层析技术及进展状况

本文概要地介绍了三维流场干涉层析技术及其进展状况。对三维流场的多方向干涉投影数据的获取方法做了简要的评述；着重介绍了三维图像建方法。对莫尔偏折法在三维流场中的应用也做了一些探讨。