基于黎曼解的移动最小二乘粒子动力学数值方法

叙述移动最小二乘粒子动力学(MLSPH)的基本原理,讨论一维MLSPH计算方法,提出两种不同精度的基于无网格的构造黎曼问题初值的方法.通过黎曼解近似粒子间的相互作用,减小接触间断附近的震荡,并给出多介质一维激波管问题的算例,验证此方法的有效性.     

合成展开法解变厚度圆薄板大挠度问题

本文选取圆薄板中心点挠度与中心厚度之比的倒数为第一小参数,设定板厚变化的指数具有一个小参数系数,采用双参数摄动并以 Hencky 薄膜解作为外场解的一级近似,求出外场解的二级近似解,再用内层坐标求得相应的各级内层解,用合成展开法求解变厚度圆薄板在均布载荷作用下的大挠度问题.文中还首次给出了变厚度薄膜解.     

老版本FORTRAN程序的结构化改造

许多程序是在早版本的Fortran77语言或Fortran60环境下设计的，也由于计算机计算存储能力和计算速度限制，程序不是结构化的。这给程序维护与改造造成了很大困难。在充分分析程序特点的基础上，尝试编制软件对老程序进行结构化改造，取得了很好的效果。     

科学计算程序的验证与确认

科学计算程序的验证和确认（Verification＆Validation，简称V＆V）是数值模拟的一个重要组成部分，其目的是为了评估模拟程序和物理模型的可靠性，并量化数值模拟程序计算结果的置信度。这是解决全系统试验无法进行或者代价十分昂贵情况下，确保数值模拟置信度的一种系统方法。     

氧的低温吸附模型与计算

受化学氧碘激光器(COIL)高效率小型化及新型压力恢复技术发展的推动,考虑COIL中气体氧的低温吸附问题。从质量守恒观点出发考察了气体流过沸石多孔床的压力与速度的关系,建立了吸附剂吸收的物质与气流中物质关系的吸附平衡方程。由吸附物较少的特点,得出在确定的压力梯度下速度为常数的近似。由吸附平衡方程与吸附速率方程,构成低温吸附模型方程,其中,亨利等温线方程作为辅助关系。对低温吸附模型方程作数值求解,得到了气流中氧质量浓度及吸收的氧质量浓度随时间变化的规律,以及这些质量浓度分布的长时间渐近特征。     

针对流体大变形问题网格邻域可变技术

结合大量实际问题,研究基于邻域可变的拉氏局部重分新方法,允许网格邻域在计算中局部跟随流体可变,虽然可能会导致任意多边形网格,但比标准的拉氏+重分方法有更好的拉氏特点.网格邻域可变的局部重分策略是基于计算过程中处理网格相交的思想,即在网格相交处改变网格邻域关系,使网格不会出现相交,保持计算正常进行.通过几个数值例子论证方法的特点,并与标准的拉氏+重分方法进行比较,显示其优点.     

辐射流体力学Lagrange方程组一类人为解构造方法

针对一维Lagrange辐射流体力学方程组,基于物理量在Euler空间和Lagrange空间的微分关系,提出了一种人为解构造方法,并构造了一类一维Lagrange辐射流体力学方程组人为解.构造的人为解在整个计算区域光滑可微,质量方程无源项.将构造的人为解应用至二维辐射流体力学Lagrange程序中,从数值误差、收敛阶方面验证了程序的正确性,展示了人为解的可行性和适用性.     

TVD自相容流体力学拉氏格式研究

流体力学数值模拟格式总体上可分为Eulerian(欧氏)、Lagrangian(拉氏)和ALE(Arbitrary LagrangianEulerian),TVD广泛应用于Eulerian格式。本文利用具有TVD保持性质的Runge-Kutta型时间离散方法,构造了流体力学Lagrangian(拉氏)自相容格式,应用von Neumann小扰动技术分析了该格式的稳定性,并进行了相应的数值模拟,较好地抑制了激波波后非物理振荡。     

X射线衍射法(XRD)分析煅烧白云石的物相组成

采用X射线衍射技术(XRD)对白云石原矿和不同温度煅烧的白云石原矿进行了物相分析,白云石原矿中主相为CaMg(CO3)2,相对含量约为98.6%,杂质相为SiO2,相对含量约为1.4%。白云石原矿经过300和500℃,煅烧2h未见主相CaMg(CO3)2分解,在750℃煅烧2h,主相CaMg(CO3)2部分分解为碳酸钙和氧化镁,在1000℃煅烧2hCaMg(CO3)2全部分解为氧化镁和氧化钙。     

任意区域上的粒子均匀分布方法

给出了任意区域内粒子的Voronoi面积均匀划分的方法,对给定的区域边界点,先由边界点得到边界粒子,再利用分布函数结合Delaunay加点算法得到内部粒子,这种方法可以得到一般的区域上粒子的均匀分布。这里给出了一个流体力学中利用SPH方法模拟内爆问题的例子,算例表明,均匀分布可以得到更高的无网格方法数值模拟精度。