混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

基于近似技术的高亚声速运输机机翼气动/结构优化设计

探索基于近似技术的高亚声速运输机机翼气动/结构多学科设计优化方法,建立了基于近似技术的多学科设计优化框架。气动学科采用全速势方程加黏性修正进行翼身组合体跨声速流动的气动计算,结构学科采用有限元分析方法进行应力与变形计算。采用均匀设计法给出若干样本点,分别采用二次响应面、Kriging模型和径向基神经网络等多种近似技术,构造气动学科和结构学科的近似分析模型,并对几种近似模型精度进行了分析和比较。研究发现,Kriging模型和二次响应面具有几乎等同的较高的近似精度,神经网络的近似精度则较差,由于二次响应面计算量更小,故最终选定为机翼设计优化的近似方法。以升阻比和结构重量为目标,考虑升力、机翼面积以及应力和应变约束条件,对运输机机翼4个外形参数和4个结构参数进行多目标、多约束优化设计。优化后的机翼具有较好的气动/结构综合性能,表明本文方法是可行的。     

可靠性灵敏度分析的一种新方法

基于极限状态函数矩估计的失效概率计算,提出一种新的可靠性灵敏度分析方法。推导极限状态函数的矩对基本变量分布参数的偏导数,并进而利用失效概率与极限状态方程矩的关系,推导失效概率对基本变量分布参数的偏导数,从而得到可靠性灵敏度。与改进一次二阶矩可靠性灵敏度分析方法相比,所提方法不用求极限状态方程的设计点,因而不需用到极限状态函数对基本变量的梯度函数,适用于隐式极限状态方程的可靠性灵敏度分析,算例结果也充分显示所提方法的合理性和精度。     

一种新的有限长等截面航空发动机短舱声学模型及数值结果

基于航空发动机短舱声学设计的需要,发展一种与以往工作不同的有限长管道声学模型声场的计算方法——直接边界元和传递单元法相结合。本方法对管道内部几何条件和声衬的敷设情况有更少的限制,可以很好地用于分段声衬的优化设计,并且提供求解有限长变截面管道声学模型的基础。     

直升机振动响应计算方法研究

在直升机研制阶段,预测全机振动响应是控制直升机振动水平的关键技术。本文建立了利用频响传递函数实测数据,对全机振动分析有限元模型进行修正的方法,机体的振动响应预测结果与实测结果具有很好的一致性,达到了较高的精度。通过研究,探索出了一套直升机结构振动响应模型修改和分析技术,为直升机振动响应预测提供了一种有效的分析手段,也为深入开展直升机振动响应预测研究奠定了基础。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

Preparation,Growth Mechanisms and Characterizations of ZnSe Films via the Solvothermal Method

With diethylamine as a solvent, ZnSe films were formed on the Si substrate directly from zinc and selenium through the modified solvothermal method. The effects of holding temperature, deposition time and substrate surface treatment on the quality and morpholo-gies of the ZnSe films were investigated. The growth mechanism of ZnSe films was proved to be a layer-nucleation growth process, which was tied in with the Stranski-Krastanov (SK) model. ZnSe films were identified by the X-ray diffraction pattern (XRD), the scan-ning electron microscope (SEM), the X-ray photoelectron spectroscope (XPS) and the photoluminescence (PL) techniques. The results indicate that the modified solvothermal method with diethylamine as a solvent is suitable to prepare high quality ZnSe films.     

安全关键软件可靠性验证测试方法研究

为了在不降低安全关键软件可靠性验证测试结果可信性的前提下减少测试用例量,在分析经典统计假设测试和无先验贝叶斯统计方法的基础上,提出了一种先验知识动态整合的贝叶斯推断统计测试方法;并提供了软件失效概率的概率密度函数先验分布参数的详细求解办法。实验表明,所提供的安全关键软件可靠性验证测试方法可以用较少的测试用例获得同样的结果可信性。     

蜂窝填充薄壁盒式梁的约束扭转

在刚性剖面假设的基础上分析具有蜂窝芯材薄壁盒式梁约束扭转问题。应用分离变量法,建立并求解了两个常微分支配方程。边界条件可以精确满足,同时,解答是用本征函数展开式表达的。扭矩图是由三角级数表示的,并用于决定上述本征函数展开式的待定系数。计算结果给出了蒙皮正应力与剪应力以及蜂窝芯材剪应力沿翼弦和翼展的分布规律。     

高速铁路隧道初始压缩波一维流动模型的数值分析方法

基于隧道内空气流通截面是时间和距离的二元函数条件与一维流动理论,提出了高速列车驶入设有缓冲结构隧道端口时产生初始压缩波的数值分析方法。一维流动理论采用考虑摩擦、传热以及空气质交换等因素的一维可压缩非定常不等熵流动模型。数值分析方法采用广义黎曼变量特征线法。与国外有关未设缓冲结构、开孔型与未开孔型缓冲结构的隧道端口三种情况的现车试验结果进行了比较,表明本文所建方法正确、合理,能够模拟计算多种类型缓冲结构条件下的压缩波波动规律。