高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

基于Deform软件二次开发和BP神经网络的TA15多向锻造微观组织预报

通过热压缩实验构建了TA15合金的热变形真应力-应变曲线,以此为基础分别建立了合金双相区及单相区温度区间的热变形本构方程;基于热压缩试样动态再结晶的统计数据建立了TA15合金的动态再结晶模型。借助Deform提供的二次开发功能实现相关数学模型的程序化,制定正交实验方案,实现了TA15合金多向锻造变形的微观组织仿真。通过正交实验分析得出各项因子的影响对象及强弱差异,提出了双相区及单相区温度区间内的多向锻造最佳因子组合。建立了TA15合金多向锻造变形微观组织的BP （Back Propagation）神经网络预测模型,将预测结果与有限元仿真结果进行比较,结果表明两种方法的预测结果基本一致,但神经网络具备有限元仿真难以实现的良好细节预测能力,能更为细致地实现对微观组织分布状态的划分。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

蜂窝夹层板超高速撞击极限方程分析

为分析撞击极限方程对蜂窝夹层结构的穿透特性的预测能力,调研得到了3类撞击极限方程的表达形式和等效方法,以及131个采用碳纤维复合材料(CFRP)面板的蜂窝夹层板结构的试验数据,并对撞击极限方程的预测能力进行了比较计算。结果发现,MET方程对他源数据的四种(未失效、失效、总体以及安全)预测率均大于80%,进行在轨航天器结构的失效分析时可优先选用;SRL方程对本源数据的安全预测率达到了100%,在他源数据上的安全预测率也很高,适用于航天器防护结构设计。探讨了撞击极限方程中的系数、速度分界值的优化思路,以提高撞击极限方程的预测率。     

三阶幻方问题的代数解法

介绍了中国古代数学家甄鸾和杨辉关于三阶幻方的一个构造方法。用线性代数的方法构造三阶幻方约束方程组时,当自由未知量分别取1,2,3,…,9时,对应的方程组共有56组整数解,为了找出幻方的全部解,需要确定自由未知量取值的限制条件,由此对方程组中的自由未知量的取值进行了分类探讨,进而求出三阶幻方的全部解。利用线性代数的方法去探求三阶幻方的解法种数显然不是最简捷的,但却有完备的理论做保证。此外,提供了一个如何利用线性代数知识来解决实际问题的实例,并简介了三阶幻方所具有的奇妙的特性。     

各向异性Chaboche黏塑性本构方程隐式应力积分算法

针对各向异性Chaboche黏塑性本构方程,利用有效参数控制法推导了该方程的隐式积分算法,该算法把所有未知量都用同一个控制参数表示,并形成控制方程,从而把本构方程应力积分的问题转化成一个非线性控制方程的求解问题.该积分算法程序,与通用有限元软件Marc结合,通过单个单元对单向拉伸、循环、松弛的响应的实验数据进行了数值模拟,通过与显式算法对比,该隐式算法的稳定性和收敛性有所提高.构件的对比分析也表明隐式算法比显式算法具有较高的效率.      

过渡状态倾转旋翼气动力模拟的高效CFD方法

为显著减少倾转旋翼过渡飞行时气动力CFD模拟的计算代价,提出并建立了适合倾转旋翼过渡状态气动特性分析的高效混合CFD方法。首先,提出了适合于过渡状态模拟的嵌套网格系统,并发展了相应的挖洞和贡献单元搜寻方法。在此基础上,结合叶素理论和动量理论建立了旋翼气动力模拟的简化虚拟桨盘模拟方法(Virtual rotor model,VBM)。为了能够分析倾转旋翼气动力的细节特性,多层运动嵌套网格系统和单指令多数据流(Single program multiple data,SPMD)并行技术被引入来建立精确的旋翼模拟方法(Real blade model,RBM)。然后,通过将VBM和RBM方法结合,构建了适合倾转旋翼过渡状态气动特性分析的高效(Hybrid blade model,HBM)方法。最后,通过对有试验值对比的悬停状态典型旋翼和7A旋翼分别验证了VBM和RBM方法的有效性。分别采用3种方法预测了过渡状态不同倾转角下旋翼的气动特性,VBM表现出最优的计算效率,能用于倾转旋翼总体气动性能的分析。HBM方法在保证流场求解精度的基础上,相对于高精度的RBM方法节省了1/3的计算时间。     

超支化聚硅氧烷改性双马来酰亚胺树脂的增韧效果和固化动力学

采用一种超支化聚硅氧烷改性二烯丙基双酚A改性双马来酰亚胺体系,对固化树脂的增韧效果进行了研究,并通过DSC分析了改性树脂体系的固化反应及其动力学。结果表明,当超支化聚硅氧烷的添加量为10%时,浇注体冲击强度达到最大,比未改性体系提高了62.8%;采用非等温DSC研究体系的固化动力学,该体系固化反应分为两步,表观活化能分别为ΔE1=81.14 k J/mol,ΔE2=89.40 k J/mol;固化反应级数分别为n1=0.91,n2=0.91。     

湍流模型对喷流计算精度的影响研究(英文)

采用三维RANS代码MBNS3D对包含了内流、外流及后体干扰的三种简化尾喷流,进行了数值模拟研究。通过壁面压力分布及喷流截面速度型与实验数据的对比,研究了SA和SST两种湍流模型及其修正(SACC,SADES,SSTCC)对喷流计算精度的影响。本文计算表明,SSTCC模拟效果最为精确;利用SA模型进行脱体涡模拟,在同样的网格下优势不明显。湍流模型可压缩修正改善了模拟喷流的膨胀趋势,但对附面层分离位置的预测,在本文的工作中,还不如原始模型。两类模型相比,SST模型的结果好于SA模型。