盒式机翼布局气动特性研究

采用数值模拟的方法研究了盒式机翼布局的气动特性。在数值模拟过程中,控制方程为雷诺平均N-S方程,采用基于非结构网格的有限体积方法进行求解。空间离散采用了Osher逆风格式,时间离散则采用五步Runge-Kutta格式。紊流模型采用了S-A方程模型。对盒式机翼布局外形在不同状态下的流场进行了数值模拟,给出了每个状态下的气动参数,并进行了分析,结果显示盒式机翼布局具有较好的气动特性。     

跨洞庭湖大桥主梁节段定常气动力研究

采用风洞实验与数值模拟相结合的手段,对跨洞庭湖大桥主梁节段的定常气动力特性进行了研究;流场数值模拟是采用ｋ－ε湍流模型封闭的雷诺时均Ｎ－Ｓ方程,求解方法是建立在交错网格基础上的压力－速度场迭代求解的ＳＩＭＰＬＥ算法;实验是在国防科技大学ＫＤ－０３ 低速风洞进行;研究表明实验与计算结果较为一致     

基于SNGR方法的二维空腔噪声数值模拟

基于SNGR方法,采用有限差分法求解带源项的线化欧拉方程。数值方法采用色散关系保持(DRP)格式。低耗散色散龙格-库塔格式显式时间推进,采用无反射远场边界条件。对二维空腔,长深比为9∶1,在亚声速、跨声速和超声速多个马赫数下进行噪声预测,并与实验结果对比,两者吻合较好。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

底部形状对钝锥体模型气动与电磁散射性能的影响

开展了飞行器气动与隐身综合特性数值研究.分别利用时域有限差分法和数值求解N-S方程的方法对飞行器的电磁散射与气动特性进行了数值模拟,研究了钝锥体模型底部形状对其雷达散射截面(RCS)和零升阻力的影响.由数值计算结果可知:合理地改变钝锥体模型底部形状,可以降低模型的RCS.并且,随着椭球体轴的长度或锥体高度的增加,模型的RCS逐渐减小.当飞行马赫数为5.0,高度为20kin时,底部形状为椭球体或锥体的模型,随着椭球体轴的长度或锥体高度的增加,模型的零升阻力略有下降.     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。     

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

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

高超声速化学非平衡绕流分布式并行算法

讨论了非结构网格高超声速化学非平衡绕流流场的分布式并行数值模拟方法。控制方程采用带化学反应源项的N-S方程,数值离散格式采用Jameson有限体积法。化学反应模型为五组元五反应模型,没有考虑电离效应,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法。通过调用METIS库函数将整体计算区域进行分裂,使用PVM并行机制在PCs-Cluster机群上对二维钝头体高超声速化学非平衡绕流进行了分布式并行计算,得到了结果,并对算法的加速特性进行了对比分析。     

高速列车绕流的数值模拟

给出了用低速N－S方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Train3D，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

高超声速小钝锥尾流转捩规律与光电特性

显式有限差分方法用于求解小钝锥高超声速化学非平衡尾流的轴对称边界层方程。该方法既适用于层流也适用于湍流。应用 Goldburg 转捩准则确定转捩初始位置。详细地计算了小钝锥尺寸、飞行速度、高度对尾流转捩和光电特性的影响。从计算中得到转捩和光电特性变化的一些有用规律。对湍流亚密部分的雷达散射截面进行了规律性分析。计算结果具有实用参考价值。     

多层壁瞬态温度分布计算

风挡玻璃电防冰和除冰装置、飞机机翼及直升机旋翼电除冰装量都是有内加热的多层壁。这些多层壁中温度分布计算不仅在这些装置的设计计算中是必需的，而且在考虑多层壁的热变形及热应力时也是必要的。本文用数值法对有内加热的多层壁温度分布进行了研究，根据导热微分方程及各种边界条件，用克兰克－尼可尔辛法推出了典型结构的有限差表示式，然后用高斯－赛德迭代法求解线性方程组，得到了可用于多层壁瞬态温度分布的计算软件。本文用该软件对实际的多层金属壁及多层玻璃作了计算。计算表明，计算结果与文［３—５］所得结果相符。本文所述的方法可推广到二维及三维瞬态多层壁温度分布计算，也可用于无内加热层的多层壁温度计算。     

一维PSD位置敏感光电器件的应用研究

本文简要叙述了一维ＰＳＤ位置敏感光电器件的工作原理,以及采用光学三角法测量构成的尺寸测量传感器,重点介绍ＰＳＤ、半导体激光器实现铁路的钢轨磨耗自动检测的方法。     

生物组织中光场分布的有限元分析

基于对脑组织多层结构特殊性的认识,将有限元方法应用到脑组织的光场分布模拟研究中,即通过区域划分和分片插值确定基函数,并根据变分原理推导了辐射传输方程在边界条件下的泛函,进而得到每个单元的特征式和整个场的联立方程式,可求解联立方程式得到光场分布的模拟结果.最后,进一步讨论了有限元方法的合理性选择问题,给出有限元方法的适用条件是吸收系数远小于散射系数的生物组织.将此方法和蒙特卡罗方法模拟的光场分布进行比较,结果表明,有限元方法对脑组织内的光场分布模拟相当精确,同时模拟所用的时间大大缩短.     

一种快速计算电磁各向异性介质目标RCS的新方法

从电磁理论中的体等效原理出发,推导了任意电磁各向异性介质目标电磁散射的体积分方程。在此基础上,由SWG基函数和伽略金法建立其矩阵方程。在迭代求解过程中,应用了快速计算阻抗元素的等效偶极子法和加速矩阵向量积的快速偶极子法两种快速算法,有效地降低了传统矩量法计算电磁各向异性介质目标RCS的时间。同时,在快速偶极子法中,远场组间互阻抗元素无需显式计算,并且形式简单的聚集、转移和发散函数可现用现算,从而大量节省了计算机内存。数值仿真结果表明:本文使用的方法在分析电磁各向异性介质目标的电磁散射特性中是非常有效的,不仅有高的计算效率和低的内存需求,而且有良好的数值精度。     

非结构网格高超声速绕流数值模拟及涵道构型减阻特性分析

在非结构网格上对二维高超声速化学非平衡粘性绕流进行了数值模拟，并应用此方法对涵道构型高超声速减阻特性进行了数值分析。空间离散采用VanLeer逆风通量分裂格式，时间推进采用显式的Runge—Kutta格式。化学非平衡动力学模型为五组元五反应模型，对化学反应源项进行了点隐式格式处理，温度场的计算采用牛顿迭代法。分别对二维类弹头体涵道构型、二维半圆柱模型的高超声速绕流流场进行了数值模拟，得到了数值结果，并与完全气体绕流计算结果进行了对比分析。结果表明本文方法可以应用于高超声速绕流计算，并且涵道构型具有高超声速减阻特性。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

重根特征向量导数的统一迭代法

在模态法之外,作者提出过好几种计算重特征向量导数的方法,如直接扰动法、广义逆法、动柔度法和混合移频动柔度法等。这些方法是不同的,其中有两种方法仅用于许多特征向量导数才是有利的．还有一种方法不能用于密集根情况。本文企图将这些方法统一成同一迭代格式。这就大大简化了这些方法的编程和执引过程。这种统一迭代格式不仅可用于密集根状态,还可以经济地用于许多特征向量导数的计算。数值结果表明,所有演变、改造过的方法都是成功的。     

CALCULATION OF THREE DIMENSIONAL BOUNDARY LAYER EQUATIONS USING DIFFERENTIAL QUADRATURE METHOD

ＣＡＬＣＵＬＡＴＩＯＮＯＦＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＢＯＵＮＤＡＲＹＬＡＹＥＲＥＱＵＡＴＩＯＮＳＵＳＩＮＧＤＩＦＦＥＲＥＮＴＩＡＬＱＵＡＤＲＡＴＵＲＥ　ＭＥＴＨＯＤＴａｎｇＤｅｎｇｂｉｎ；ＣｈｅｎｇＧｕｏｗｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏ...     

设计灵敏度分析的迭代模态法

在结构动态有限元数学模型修正、结构动力学修改以及结构动态设计中，灵敏度分析是十分重要的一环，往往也是主要计算工作量所在。如何在保证精确度前提下减少灵敏度分析工作量和计算时间，具有重要意义．本文提出的迭代模态法，可以大大提高特征向量对设计参数导数（模态灵敏度）的计算效率。本文还从理论上证明经典模态法和修改模态法只是迭代模态法的特殊情况。文中还给出了ＦＯＲＫ，３Ｄ－ＦＲＡＭＥ两种结构为实例的计算机仿真，并与已有两种模态方法进行比较。理论分析与计算机仿真表明，本方法不仅能充分保证特征向量灵敏度分析的精度，而且大大提高了计算速度．