转动非平衡玻尔兹曼模型方程统一算法与全流域绕流计算应用

基于过去开展稀薄自由分子流到连续流气体运动论统一算法框架,采用转动惯量描述气体分子自旋运动,确立含转动非平衡效应各流域统一玻尔兹曼模型方程.基于转动能量对分布函数守恒积分,得到计及转动非平衡效应气体分子速度分布函数方程组,使用离散速度坐标法对分布函数方程所依赖速度空间离散降维;应用拓展计算流体力学有限差分方法,构造直接求解分子速度分布函数的气体动理论数值格式;基于物面质量流量通量守恒与能量平衡关系,发展计及转动非平衡气体动理论边界条件数学模型及数值处理方法,提出模拟各流域转动非平衡效应玻尔兹曼模型方程统一算法.通过高、低不同马赫数1:5～25氮气激波结构与自由分子流到连续流全飞行流域不同克努森数（9×10-4～10）Ramp制动器、圆球、尖双锥飞行器、飞船返回舱外形体再入跨流域绕流模拟研究,将计算结果与有关实验数据、稀薄流DSMC模拟值等结果对比分析,验证统一算法模拟自由分子流到连续流再入过程高超声速绕流问题的可靠性与精度.      

GAS-KINETIC UNIFIED ALGORITHM FOR BOLTZMANN MODEL EQUATION IN ROTATIONAL NONEQUILIBRIUM AND ITS APPLICATION TO THE WHOLE RANGE FLOW REGIMES

基于过去开展稀薄自由分子流到连续流气体运动论统一算法框架,采用转动惯量描述气体分子自旋运动,确立含转动非平衡效应各流域统一玻尔兹曼模型方程.基于转动能量对分布函数守恒积分,得到计及转动非平衡效应气体分子速度分布函数方程组,使用离散速度坐标法对分布函数方程所依赖速度空间离散降维;应用拓展计算流体力学有限差分方法,构造直接求解分子速度分布函数的气体动理论数值格式;基于物面质量流量通量守恒与能量平衡关系,发展计及转动非平衡气体动理论边界条件数学模型及数值处理方法,提出模拟各流域转动非平衡效应玻尔兹曼模型方程统一算法.通过高、低不同马赫数1：5～25氮气激波结构与自由分子流到连续流全飞行流域不同克努森数（9×10^-4～10）Ramp制动器、圆球、尖双锥飞行器、飞船返回舱外形体再入跨流域绕流模拟研究,将计算结果与有关实验数据、稀薄流DSMC模拟值等结果对比分析,验证统一算法模拟自由分子流到连续流再入过程高超声速绕流问题的可靠性与精度.     

用格子Boltzmann方法研究Burgers方程

提出了用于Ｂｕｒｇｅｒｓ方程的格子Ｂｏｌｔｚｍａｎｎ模型,应用Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ展开和多重尺度技术,通过选择平衡态分布函数的高阶矩,得出了几种精度的Ｂｕｒｇｅｒｓ方程,模型中的参数通过分析耗散性质和色散给出。     

基于Boltzmann模型方程的气体运动论统一算法研究

模型方程出发,研究确立含流态控制参数可描述不同流域气体流动特征的气体分子速度分布函数方程; 研究发展气体运动论离散速度坐标法, 借助非定常时间分裂数值计算方法和NND差分格式, 结合DSMC方法关于分子运动与碰撞去耦技术, 发展直接求解速度分布函数的气体运动论耦合迭代数值格式; 研制可用于物理空间各点宏观流动取矩的离散速度数值积分方法, 由此提出一套能有效模拟稀薄流到连续流不同流域气体流动问题统一算法. 通过对不同Knudsen数下一维激波内流动、二维圆柱、三维球体绕流数值计算表明, 计算结果与有关实验数据及其它途径研究结果(如DSMC模拟值、N-S数值解)吻合较好, 证实气体运动论统一算法求解各流域气体流动问题的可行性. 尝试将统一算法进行HPF并行化程序设计, 基于对球体绕流及类``神舟'返回舱外形绕流问题进行HPF初步并行试算, 显示出统一算法具有很好的并行可扩展性, 可望建立起新型的能有效模拟各流域飞行器绕流HPF并行算法研究方向. 通过将气体运动论统一算法推广应用于微槽道流动计算研究, 已初步发展起可靠模拟二维短微槽道流动数值算法; 通过对Couette流、Poiseuille流、压力驱动的二维短槽道流数值模拟, 证实该算法对微槽道气体流动问题具有较强的模拟能力, 可望发展起基于Boltzmann模型方程能可靠模拟MEMS微流动问题气体运动论数值计算方法研究途径.      

用格子Boltzmann方程模拟浅水波问题

提出了用格子Ｂｏｌｔｚｍａｎｎ方程（ＬＢＥ）模拟浅水波问题的方法．通过无粘气体运动方程与浅水波方程的比较，确定了ＬＢＥ方法中平衡态的形式，使宏观方程与浅水波方程一致．计算了二维浅水波的一个问题，数值结果与精确解做了比较．     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

一类新的超弹性-循环塑性本构模型

目前,很多经典的超弹性-有限塑性本构模型已被提出,但由于超弹性理论中中间构型的引入使得随动硬化法则相对复杂,故多数文献均采用的是经典的Armstrong-Frederick(A-F)随动硬化法则.本文基于已有的本构理论,利用多机制过程的概念拓展了Lion塑性变形分解理论,明确提出了多重中间构型的概念,并在此基础上,对经典理论中客观性的定义进行了概念上的推广,使其更好地适用于超弹性本构理论分析,同时提出了一类新的超弹性-有限塑性本构模型.这类本构模型满足热动力学法则,且可融合多种小变形循环塑性理论中常用的随动硬化法则(如经典的A-F模型,Chaboche模型,Ohno-Wang(O-W)模型以及Karim-Ohno(K-O)模型等),使得小变形理论中背应力的加法分解性质及其演化的临界面阶跃特性在大变形领域中均有所体现,故本文提出的本构理论可看作是小变形循环塑性模型在大变形理论中的扩展.本文最后以K-O模型为例,对推荐模型进行了详细探讨,并与相应的次弹性模型进行了对比.      

复合材料弹塑性多尺度分析模型与算法

对材料非线性多尺度分析的计算模型与算法进行研究.在构建周期分布单胞分析算法的基础上,发展针对复合材料结构材料非线性多尺度分析的一般有限元方法.方法的特点是将所建立的单胞分析过程作为有限元分析的子程序嵌入到总程序系统当中,完成对应的高斯点应力计算,因而使所发展的方法具有实现方便的特点.给出数值计算结果,验证了方法与所发展的多尺度有限元分析程序的正确与有效性.     

模拟冻-融过程的热-水-应力耦合模型及数值分析

考虑低温条件下出现水-冰相变并有潜热产生的特点,将所开发的热-水-应力耦合模型及其二维有限元程序作了改进,使之可用于岩土体冻融问题的分析.并以Neaupane等的岩石室内冻一融试验为对象进行了数值模拟,将计算、试验和Neaupane的计算中所得的岩石试样上特定点的温度-时间曲线、温度-径向应变曲线作了对比,看到三者的接近程度较好,同时考察了本计算所得的一个岩石试样中不同时刻的温度、位移及应力分布,得出了若干结论.