一种改进的光滑粒子流体动力学前处理方法

为了便于对任意边界形状的计算域快速地布置均匀粒子,提出了一种改进的光滑粒子流体动力学前处理方法.该方法是在2012年Colagrossi等提出的算法基础上进行改进后得到的.Colagrossi等提出的算法能够计算一些简单外形分布比较均匀的粒子.然而当光滑长度与初始粒子间距的比值较大时该方法在计算过程中会出现较强的数值震荡问题,收敛速度慢;而且在计算过程中可能会遭遇流体粒子穿透固体壁面的问题.本文通过引入未知因素修正的平滑粒子动力学模型来提高计算稳定性,并通过对边界附近的流体粒子施加边界力来避免流体粒子穿透固体壁面.算例验证结果表明,利用改进后的光滑粒子流体动力学前处理方法能够快速地对各种边界形状的计算域分布均匀粒子,并且避免了流体粒子穿透固体壁面的问题.     

基于弱可压与不可压光滑粒子动力学方法的封闭方腔自然对流数值模拟及算法对比

为了对比研究弱可压光滑粒子动力学(WCSPH)方法和不可压光滑粒子动 力学(ISPH)方法在模拟封闭方腔自然对流问题时的特性, 采用粒子位移技术有效地解决了高瑞利数条件下, 拉格朗日型SPH方法模拟封闭方腔自然对流时流体域内的粒子聚集和空穴问题, 将拉格朗日型SPH 方法求解封闭方腔自然对流问题的最高瑞利数提高到了10⁶; 进而通过对比瑞利数分别为10⁴, 10⁵, 10⁶的条件下, 采用拉格朗日型WCSPH、 拉格朗日型ISPH、欧拉型ISPH三种SPH方法模拟得到的封闭方腔速度分布云图、 温度分布云图、壁面努赛尔特数分布曲线和平均努塞尔特数, 分析了三种SPH方法在模拟封闭方腔自然对流时的精度、稳定性和计算效率. 结果表明: 在低瑞利数条件下, 以上三种SPH方法都可以较好地模拟此问题, 在高瑞利数条件下, 欧拉型ISPH方法的模拟结果最为精确; 拉格朗日型WCSPH方法模拟所得结果比拉格朗日型ISPH方法模拟所得结果稍好些. 关键词： 光滑粒子动力学 不可压光滑粒子动力学 粒子位移技术 自然对流     

气泡碰壁受力模型和反弹规律

对Mo=10-8~10-12及Re=5~750范围内的上升气泡与壁面垂直碰撞问题进行了理论求解，研究了不同控制参数下气泡碰壁反弹的规律.气泡上升和碰撞过程的运动方程考虑了浮力、液体阻力、附加质量力和与壁面碰撞时引起的薄膜诱导力.气泡碰壁过程气泡界面与壁面形成的液膜厚度变化规律由Stokes-Reynolds方程计算得到.膜内气泡变形引起的流体压强采用Young-Laplace方程求解.结果表明，基于SRYL方程的薄膜诱导力模型可以很好地预测不同Reynolds数下气泡0到多次的反弹轨迹，计算结果与实验结果吻合良好.气泡在碰壁反弹过程中会形成丰富的薄膜形状，如酒窝状变形，丘疹状变形和涟漪状变形.气泡界面变形会引起膜内压强的变化，压强的分布规律与气泡界面形状有着重要的关系.气泡在与壁面碰撞的过程中，薄膜诱导力会起主导作用，且随着Reynolds数的增加薄膜诱导力最大量级增大.气泡碰撞壁面时，反弹次数与Reynolds数有着直接的联系，不同Morton数下的气泡均在相同Reynolds数附近发生气泡反弹次数的变化.      

兴趣引导型教学模式的改革与实践

针对力学专业课程设置的特点,进行了兴趣引导型的教学模式的改革与实践. 通过组织概论 型课程与讲座、科学创新实验、力学史教育及实践教育与理论教育相结合等方式,使力学专 业的学生尽早了解本专业的课程设置、专业内容及其分类、专业知识在工程实际中的广泛 应用,从而激发学生学习的兴趣,提高同学学习的积极性与主动性.     

关于Arnold的《经典力学的数学方法》

介绍Arnold的名著《经典力学的数学方法》的主要内容、主要特 色以及它的学术影响, 并给出几点评论.     

辛几何算法在特殊Chaplygin系统中的应用研究

采用辛几何算法对广义力具有广义势的特殊Chaplygin系统的动力学性质进行数值研究,并和传统的Runge-Kutta算法进行比较,得出辛几何算法在这类非完整力学系统的计算中具有优越性. 关键词： 辛几何算法 Chaplygin系统 Helmholtz条件     

无线供能演示装置

根据电磁感应耦合原理,建立电容补偿式电路的相量模型,确定了影响传输效率的主要因素。基于实验得出最优电路参量,制作了简易的无线供电小车演示装置,具有较好的教学演示效果。     

冻土单轴动态加载下的力学性能

研究冻土的动态力学性能对于地下工程人工冻结法施工等具有重要意义。本文应用分离式霍普金森压杆(SHPB),研究了冻土单轴动态加载下的力学性能,涉及-3、-8、-13、-17、-23和-28℃共6个负温的冻土,应变率范围350～1200s-1。获得了相应条件下的冻土应力应变关系。冻土的单轴动态应力应变曲线具有脆性特征。发现冻土具有温度和应变率效应,其强度随温度降低和应变率增大而增大,最终应变随应变率增大而增大。冻土温度越低,应变率敏感性越强;加载应变率越高,冻土的温度效应越显著。文中提出的粘弹性损伤型本构模型能够较好的描述6个温度冻土的应力应变关系。     

滚波演化中聚合过程的数值模拟研究

滚波是一种重力作用下自由液面失稳诱发的水面波动现象, 通常可分为具有相对稳定波形和波速的周期性滚波与波形和波速不断变化的不规则滚波(自然滚波). 不规则滚波的相互作用和发展演化过程十分复杂, 至今对其认识尚不成熟. 本文采用基于雷诺平均Navier-Stokes方程的立面二维数值模型, 对不规则滚波发展过程中的吸收聚合和追赶聚合现象进行了数值模拟研究. 分析了两种聚合模式的演化过程, 给出了滚波聚合过程中完整的波形、波速、速度剖面以及湍流黏性等重要信息. 结果表明滚波的聚合过程是不规则滚波演化和增长的重要机制, 在特定条件下滚波增长由自然增长模式转变为以吸收聚合和追赶聚合为主的增长模式. 滚波聚合过程中, 依次经历后波追赶、爬升、与前波合并、内部流场调制等多个步骤, 最终形成一个具有更大波长和波高的滚波. 本文发现了在3个滚波间距较近的情况下, 会发生二重聚合现象, 即后两个滚波首先聚合, 然后与前波进一步聚合形成一个新的滚波.      

端元提取技术在高光谱图像压缩中的应用

高光谱图像海量数据如何实现大比例有效压缩是限制其应用的主要问题之一,而现有有损压缩方法存在大压缩比与光谱特性信息准确保留的矛盾,即使现有最优有损压缩方法也不能够得到令人满意的结果。文章基于混合像元分解的思想提出基于端元提取技术的数据有损压缩方法来解决该矛盾,首先用顶点成分分析(VCA)方法提取场景中地物的端元光谱,根据各端元与观测像元之间的光谱间余弦角相似性度量方法估计各端元的丰度,接着对端元光谱及丰度数据进行无损压缩,最后利用JPEG2000有损压缩方法对高光谱图像的所有单波段图像进行空间维大比例有损压缩。AVIRIS高光谱图像的仿真结果表明,压缩比得到大幅度提高,光谱信息得到有效恢复。在实现压缩比为50∶1时,大部分像元的光谱角误差在2%左右。