蚯蚓蠕动过程中非光滑波纹形体表的力学分析

蚯蚓体表的形状是蚯蚓适应其生存环境和自然进化的结果．用力学分析方法，解释蚯蚓的多体节非光滑形态在其运动过程中所表现出的省力、节省能量的特性．计算结果表明具有波纹形表面的蚯蚓在蠕动过程中可以节省1．12%的纵向肌肉的收缩力，在运动过程中所消耗的能量也可以节省1％左右．     

用格子Boltzmann方法数值模拟混合层中旋涡的合并过程

用格子Boltzmann方法计算混合层中的流动问题。在流场的入口处加不同频率、振幅和相位的小扰动,观察混合层中旋涡的演进机理,模拟二维混合层中旋涡合并现象。在基本扰动波的基础上,又加入频率为基本波频率一半的亚谐波,得到了两个涡合并的计算结果,当加入的亚谐波频率为基本波频率的三分之一时,得到了三个涡合并的计算结果。这些计算结果与已有文献的结果基本一致,显示用格子Boltzmann方法模拟混合层问题是可行的。     

用格子Boltzmann方法模拟非线性热传导方程

对具有非线性源项和非线性扩散项的热传导方程建立格子Boltzmann求解模型.在演化方程中增加了两个关于源项分布函数的微分算子,对演化方程实施Chapman-Enskog展开.通过对演化方程的进一步改进,恢复出具有高阶截断误差的宏观方程.对不同参数选取下的非线性热传导方程进行了数值模拟,数值解与精确解吻合得很好.该模型也可以用于同类型的其他偏微分方程的数值计算中.     

用离散速度方法计算浅水长波方程

用离散速度法计算浅水波方程，将空气动力学方程和浅水波方程作了比较，用Nadiga提出的近平衡流动方法模拟浅水波方程的连续和间断解。计算了一维的溃坝波问题和Thacker提出的连续解问题，结果与精确解作了比较，并且计算了水流跃过障碍物的问题。     

蚯蚓爬行的生物力学分析与数值计算

蚯蚓是典型的多体节动物。本文从生物力学的角度,分析了蚯蚓在静止和运动状态下的体壁应力,又分析了蚯蚓爬行过程中的变形以及与地面之间的摩擦力,给出了蚯蚓在地表爬行的动力学模型。进一步,又数值计算了一定质量范围的蚯蚓在静止和运动状态下体壁承受的最大应力和动力学模型方程,得到蚯蚓爬行的平均速度与质量的关系。计算结果与实验观测值基本一致。      

实数型格子气模型及方柱绕流的模拟

在4Bits的HPP格子气模型基础上提出了用实数表示节点每一速度方向的粒产密度的实数型HPP模型,建立了粒子碰撞的运算关系式(代替碰撞规则),用此方法对方柱绕流问题进行了计算并介绍了边界条件的处理方法,计算结果较好地反映了真实流动特性。     

动力学的流矢量分裂法模拟浅水波方程

Based on the analogy to gas dynamics, the kinetic flux vector splitting (KFVS) method is used to stimulate the shallow water wave equations. The flux vectors of the equations are split on the basis of the local equilibrium Maxwell-Boltzmann distribution. One dimensional examples including a dam breaking wave and flows over a ridge are calculated. The solutions exhibit second-order accuracy with no spurious oscillation.     

慢周期性流动中毕托管和热线测速的对比

本文采用微机控制调速系统控制低速风洞,在风洞内获得流速按周期性变化的气流.用毕托管和热线风速仪同时测量风洞内按周期性变化的气流速度.毕托管所测结果与热线测得结果相比有较大的相位滞后和振幅衰减.本文分析了毕托管的动态响应特性,采用一阶传递函数对毕托管测速系统的动态特性进行描述,并确定了测速系统的惯性时间和纯滞后时间.     

徐州矿区采空塌陷区公路工程地质灾害与防治对策分析

采用格子Boltzmann方法模拟可变形膜与周围流体的相互作用. 分析了格子Boltzmann 方法中的边界处理方法和边界受力的计算方法，并且用此方法计算流场中可变形膜的受力. 可将离散化后的膜看作一系列的质点，从而得到膜的动力学方程. 将可变形膜在流场中受到 的力引入方程中，可以计算膜的变形. 求解了几种不同情况下，膜的形状随时间的变化. 发现，如果可变形膜非常软或者非常硬，经过足够长的时间后，膜的形状会接近一 条直线，即回到初始状态. 模拟过程是二阶精度的.     

骰子中灌铅对掷骰子结果的影响

骰子是一种常用的赌博工具，段子灌铅后，重心可能偏离其几何中心，将骰子视为刚体，将掷骰子过程看作刚体自由下落然后与桌面碰撞的刚体动力学过程．本计算出重心偏离的骰子下落后，各个面向下的概率。