低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。     

一种基于特征关系式的预处理远场边界条件

高速可压缩流动的计算方法应用于不可压低速问题会因为当地速度和当地声速量级相差较大产生刚性问题并导致数值收敛困难。预处理方法引入预处理矩阵使传统的可压缩方法具备了求解不可压或低速问题的能力。由于预处理方程改变了Navier-Stokes方程的特征值,原Navier-Stokes方程的远场边界条件不再适用。目前被广泛使用的预处理方程的远场边界条件为简化边界条件。本文以预处理方程为基础进行分析,推导了一种基于特征关系式的远场边界条件。通过典型状态算例数值实验,验证了该边界条件及预处理方法的有效性。结果显示,对于定常/非定常低速不可压流动,合理地设置远场边界条件,预处理方法都能够提高计算的收敛性和精度;对于跨声速流动,预处理方法和未加预处理方法两者的计算效率和计算精度相当。本文的预处理远场边界条件与简化预处理远场边界条件相比一是能够进一步提高计算的收敛速度;二是能够有效降低远场边界位置对预处理方法数值计算结果的负面影响。     

计算全马赫数下粘性流动的隐式算法

把预处理方法引入到求解三维粘性流动的Navier-Stokes方程隐式格式算法中,并应用这种方法对平板算例在层流和湍流两种不同流态下及四个不同Ma数下进行了数值模拟,计算结果与理论值及相关公式吻合,同时对比分析了有无预处理情况下的收敛性,证明了本文所采用的预处理隐式算法对于任意马赫数都是适合的,而且可以提高残差精度量级,加快收敛。     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

预处理方法在混合笛卡尔网格中的应用研究

使用基于混合笛卡尔网格的预处理方法,对低马赫数下的定常/非定常流动问题进行了数值模拟研究。网格生成中,在物面附近生成贴体结构网格,其余计算区域使用笛卡尔网格进行填充,两套网格之间通过查找"贡献单元"实现流场数据间的传递。同时,使用基于密度的预处理方法,发展了一套可求解从低速(极低马赫数)到常规马赫数的N-S方程求解器,时间离散使用隐式双时间步LU-SGS格式,空间离散使用具有二阶精度的格心格式有限体积方法。非定常流动问题的数值模拟过程中,使用逆距离插值来进行新现网格单元上流场参数确定。通过对NACA0012翼型在低马赫数下的定常绕流和动态失速问题的数值模拟研究表明:使用预处理方法能够加快低马赫数下数值计算的收敛速度和提高计算的准确性,基于混合笛卡尔网格的N-S方程求解器能够模拟包含运动边界的低速不可压非定常流动问题。耦合预处理技术的混合笛卡尔网格方法给包含低速和常规马赫数的复杂运动边界问题的求解提供了新的思路。     

A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

时间谱方法中的高效GMRES算法

研究了时间谱方法求解周期性非定常流场的计算效率,并对时间谱方法应用于周期性非定常流动的隐式求解方法进行探讨。当采样点数增加或减缩频率增大时,时间谱方法对应的雅可比矩阵对角占优性质迅速恶化,导致很多传统的迭代方法失效。为了解决上述问题,论文采用带预处理的广义极小残差(GMRES)算法来提高雅可比系数矩阵的计算收敛性。使用时间谱方法对NACA0012翼型强迫振荡算例进行计算,并与时域差分方法的计算效率和精度进行对比。研究表明在保证计算精度的同时,时间谱方法普遍可将计算效率提高一个量级左右。对于跨声速周期性流动,广义极小残差算法不论是稳定性还是收敛性都优于对称SGS迭代算法。     

抛物化Navier-Stokes方程在求解高超声速流动稳定性分析 的基本流中的应用

通过在边界层内部使用边界层上缘的流向压力梯度代替该处的流向压力梯度,对传统抛物化Navier-Stokes(PNS)方程求解中的流向压力梯度的处理方法进行改进,从而使它可以求解用于稳定性分析的高超声速流动的基本流,并分别计算了平板、零攻角钝锥和小攻角钝锥三种典型算例.对于流向没有分离,且横向流动不强的流动,使用PNS方程计算高超声速流动稳定性分析用的基本流是可靠的.特别是改进后的方法求得的基本流的稳定性分析结果与直接求解Navier-Stokes(N-S)方程的结果非常吻合,但计算的时间消耗和空间占用都减少了一个数量级.     

空化流动隐式求解方法研究

空化流动问题本质上是可压缩流动,应用可压缩方法开展数值模拟研究更符合物理实际。气体和液体压缩性的显著差别使得低速空化流动数值模拟的刚性问题非常突出,通过引入预处理技术解决低速问题中由于特征值量值不一致导致的收敛刚性问题,提高收敛速度。同时,以预处理后的特征值构造Roe格式耗散项,提高低速流动计算精度。鉴于自然空化流动中气液组分转换现象和物质输运现象并存,且气体和液体的密度在常温状态下存在显著差别,预处理后的源项雅克比矩阵的特征值与无粘通量雅克比矩阵的特征值存在量级的差异,这会使得求解过程不稳定或收敛速度极慢,即出现"源项刚性"问题。为此,本文系统推导了预处理框架下的气液两相流隐式求解方法,采用点隐式方法处理源项,通过直接求逆的方式增强算法的稳定性。研究中分别考察了三种不同的算子分裂方式在LU-SGS(Lower-Upper Symmetric-Gauss-Seidel)隐式迭代中的模拟效果,并提出适用于DDADI(Diagonally Dominant Alternating Direction Implicit)的算子分裂方式,在NACA0015水翼空化算例模拟中综合比较了上述四种隐式迭代方式对低速空化问题收敛性的影响。最后,应用本方法对三维尖锥自然空化算例进行了考核模拟,捕捉到自然空泡流场中的主要特征,计算结果与试验测量结果吻合。     

结构非概率可靠性灵敏度分析方法

 结构非概率可靠性的灵敏度分析对工程中非概率问题的分析预测和优化具有重要的指导作用。基于结构非概率可靠性指标理论,提出了结构非概率可靠性的灵敏度分析方法。针对极限状态函数为线性的情况,推导出了非概率可靠性灵敏度的解析解。在此基础上,结合函数线性化理论,提出了求解极限状态函数为非线性情况的非概率可靠性灵敏度的近似解析方法。基于求解非概率可靠性指标的优化方法,并结合差分理论,提出了求解非概率可靠性灵敏度的优化方法。给出了算例分析,结果表明了所提方法的可行性。     

PCI总线驱动程序开发及WDM方法的应用

简单介绍了PCI总线的性能;详细阐述了WDM设备驱动程序的特点并介绍了开发PCI总线驱动程序主要模块的设计方法。     

Theory of formation of the ionosphere

Current knowledge about the solar radiation and absorption and ionization cross sections of atmospheric gases is reviewed. Next the main observed features of ionospheric layers are summarized. Using CIRA 1965 model atmospheres the heights of the peak of the ionization rate are calculated for a number of solar emission lines and it is made clear which of these lines are responsible for the formation of E and F1 layers. The mechanism of electron removal in the F and upper E regions as well as in the lower regions is considered, and the mechanism of formation and some behaviours of each ionospheric layer is discussed. In particular, the equatorial F2 layer is briefly considered. Discrepancies are pointed out between the values of the recombination coefficient and the rate constant for ion-atom interchange reaction obtained from ionospheric observations and from laboratory experiments. Inconsistency of the values of the intensity of solar radiation measured by rocket techniques and inferred from ionospheric considerations is also noted. Some evidence is presented suggesting that corpuscular radiation may be responsible for part of the ionization in the ionosphere even in temperate latitudes.     

图像"三针法"测量锥螺纹中径及测量不确定度评定

针对一种基于CCD图像的锥螺纹参数测量系统的开发,介绍了图像"三针法"的测量原理,并对图像"三针法"测量锥螺纹中径进行了测量不确定度的评定.     

空中加油问题的最优化研究

通过对空中加油在作战使用中的问题进行分析,针对受油机往返各加１次油的情况,给出了空中加油的数学模型,并采用跗算法对问题进行了优化计算,实际计算表明,所建模型具有良好的适应性,能够适合不同的作战环境和不同作战要求,采用跗优化算法研究空中加油问题,优化效率高、具有全局性,并易于计算机实现,所得的加油方案是合理的。     

整体壁板结构弹塑性弯曲中性层位置分析

中性层位置是壁板结构设计和分析及其弯曲成型工艺的重要参数,工程上将截面形心作为弹塑性弯曲中性层位置的方法存在很大误差。通过弹塑性分析,考虑了材料塑性强化效应,建立了一种便捷的计算整体壁板结构弹塑性弯曲中性层位置的方法,并和数值仿真结果进行了对比,分析结果表明:弹性弯曲时,中性层位于截面形心;而进入塑性弯曲阶段,中性层位置随着弯曲半径的减小向蒙皮方向显著移动;在整体壁板常用的相对弯曲半径下(ρ/t≥10),壁板的变薄可以忽略不计,应力中性层和应变中性层几乎重合。     

一种飞船返回舱落点选优算法的设计与实现

设计一种选优算法。算法依据概率论,以寻找落点最密集区域为思路,从已知众多飞船返回舱落点数据中计算产生一个“优秀”落点,即尽可能接近返回舱的实际着陆点。     

伞舱系统跟踪引导数据迭代算法

针对在返回舱乘伞下降阶段,测控链中弹道测量设备丢失目标,以及跟踪结束后无法外推生成可靠引导数据,光学设备难以快速捕获目标等问题,以风场修正模型为基础,通过风场提前修正与采用测量数据实际修正相结合的方法,建立伞舱系统跟踪引导数据迭代算法.利用多次返回段任务数据对该算法进行实验验证,并与实际任务中目前使用的方法进行比较,结果表明：通过该算法推算较长时间后的引导数据,其平均误差为现有方法的6％;推算较短时间后的引导数据,其平均误差为现有方法的20％.因此,利用该算法计算的引导数据精度较高,有助于光学设备及其他测控设备快速捕获或重捕目标.     

基于CAM/CNC集成的航空大型薄壁件数控加工在机刀轨调整方法

航空大型薄壁零件刚性弱,加工过程中易产生超差变形,需要对零件进行在机检测,并将检测数据反馈给CAM系统进行拟合并调整刀轨。以上过程需要在不同部门和系统间进行切换,严重影响了生产效率。针对以上问题,提出了一种基于CAM/CNC集成的航空大型薄壁件数控加工在机刀轨调整方法。建立了CAM/CNC集成框架,并构建了用于CAM/CNC间信息传递的动态特征信息模型,基于CAM系统离线生成特征理论中间状态,基于CNC系统建立特征实际中间状态检测点同理论中间状态的映射关系。通过分析理论厚度和实际厚度,确定特征中间状态的变形情况。针对变形导致特征中间状态加工余量无法包络符合公差要求的最终加工状态的情况,进行在机刀轨调整。基于所提出的方法开发的航空大型薄壁件数控加工在机刀轨调整系统在某大型航空企业得到了应用验证,可有效避免信息的反复传递,减少了工艺人员的重复劳动,提高了生产效率。     

飞行模拟训练器中发动机起动过程模型的建立和仿真

在飞行模拟训练器的研制中 ,发动机的起动过程是需要重点研究和模拟的过程之一。本文根据发动机起动工作过程的具体特点 ,采用三段积分法建立了仿真模型并进行了仿真研究。所得结论基本达到了训练器的要求。