求解非线性方程及方程组的模拟退火算法

探讨了用模拟退火算法求解非线性方程及方程组的方法,该方法把方程组求解问题转化为函数优化问题.计算中不需要使用目标函数的导数信息,在Matlab环境下实现了该算法,数值实验结果表明了该算法的有效性.     

液氢贮箱停放过程中的力热分析

为了寻找到运载火箭长时间停放过程中液氢贮箱的最经济液位,用于制定合理的发射流程以及紧急处置方法,采用计算流体力学(CFD)技术,对某型运载火箭停放期间液氢贮箱的力、热情况进行了仿真计算和分析。计算选用了VOF(Volume-of-fluid)两相流模型以及Lee相变模型,为了提高Lee模型在不同压力情况下对相变过程的模拟精度,采用安托因方程修正了该模型。修正后的模型首先由试验数据校验了其精确性,随后开展的液氢贮箱停放过程仿真结果表明：贮箱的竖直方向与径向均存在温度分层的现象,液相内会形成大的漩涡,该漩涡会使得冷热流体不断进行热交换,并导致贮箱内部的液氢出现气化。贮箱停放期间蒸发率最大值超过2 m 3/h,发生在停放4 h左右;而贮箱液位充填至37 m 3以上或17 m 3以下时蒸发率较低,最小值接近1 m 3/h。     

多芯片组件（MCM）用金锡焊带材料电阻率及热导率温度特性研究

对牌号为Au₈₀Sn₂₀的金锡焊带材料在208~423 K的电阻率及热导率与温度的函数关系进行了研究,并对其在多芯片组件（MCM）中的传热效果进行了评估。分别对材料在208~423 K中5个温度点的电阻率及4个温度点的热导率进行了测试,基于理论模型建立电阻率/热导率随温度变化的函数关系,最终采用模拟热扩散数值方法评估材料在高温下的传热能力。结果表明,采用修正函数模型后,金锡焊带材料在208~423 K下热导率与电阻率的关系符合测试结果,随芯片表面温度的边界条件从208 K升高至423 K,采用变温热导率模型得到的热流密度模拟计算结果相比理想化恒定热导率模型的差异性逐渐升高至5.5%。综上,金锡焊带材料热导率与电阻率的关系符合Wiedemann-Franz法则修正后的Smith-Palmer方程,在该材料传热设计时应考虑其热导率温度效应。     

聚氨酯泡沫塑料拉伸本构关系及其失效机理的研究

 通过 3种密度硬质聚氨酯泡沫塑料的拉伸实验,研究了它们的应力 -应变特性及其应变率效应。同时,为确定泡沫塑料拉伸失效机理,对材料进行了扫描电镜下的细观拉伸实验,观察了试件表面胞体的变形与失效过程,再结合宏观拉伸试件断口的扫描电镜分析,进一步讨论了泡沫塑料在拉伸加载下的失效机理。基于泡沫塑料的拉伸应力-应变曲线,用数值方法拟合了泡沫塑料的拉伸本构关系。     

涡喷发动机加速控制系统数学模型及数字仿真

 本文介绍了航空涡轮喷气发动机加速控制系统数学模型建立、数字仿真方法。由于加速过程是一个大偏离动态过程,因而发动机与加速控制器模型均用非线性方程描述,对非线性方程组求解和系统迭代方法本文进行了讨论。算例中分析了在不同飞行条件下发动机的加速规律。     

叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

战斗机弹射救生系统数值模拟研究

对弹射救生系统这种几何外形复杂、不规则非流线型钝体,以Jameson有限体积法为基础发展了一套基于高质量混合网格的Euler方程解算器,在Euler方程解算器基础上,采用S-A湍流模型发展了一套同样基于混合网格的N-S方程解算器,初步建立了适合弹射救生系统这类几何外形复杂物体的数值模拟技术平台,并采用多种加速收敛措施和合适的耗散项模型使得所发展的解算器在软件实现上鲁棒性很强。通过选取合适的计算起始迎角,对弹射救生系统成功地进行了大迎角和大侧滑角绕流计算,获得了与风洞试验比较吻合的计算结果,并且可以清楚地观察流场中马赫数分布、弹射救生系统各部位压力分布等流场细节,对弹射救生系统空气动力特性分析具有重要的指导意义。该计算平台可作为风洞试验的合理补充、加强和完善。     

Hamilton正则方程半解析法的收敛和对称分析

近年来,Hamilton正则方程半解析法在工程问题上的应用越来越广泛,但至今未见有关这种方法收敛性和对称性问题研究的文献。基于Hamilton正则方程的半解析法理论,通过变分原理详细推导了Hamiltonian元素的固支和简支边界公式及对称边界公式。多个实例的数值研究表明：随着网格加密,Hamilton正则方程半解析法的收敛速度快于一般传统位移有限元法,对称解法的效率明显优于整体解法。     

边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.