摘要

基于Lagrange插值多项式,实现了高次边界单元形函数系数的计算机自动生成.在已有闭合单元的基础上,充分利用粒子几何形状的特性构建了高次光滑边界单元.数值算例表明,与传统的二次单元和闭合单元相比,高次光滑边界单元能够较大地提高椭圆和椭球粒子数值模拟的计算精度与效率.