物质点法中如何实现周期性边界？ 10

路过的相关方向的大佬们，求求给点建议，可以有偿。目前想用物质点法模拟一个无限长斜面上的稳定颗粒流算例。如下左图： 斜面倾角为θ，长度无限长，斜面上覆盖有一定厚度的颗粒体，在自重条件下加速直到最终达到稳定速度（类似于μ(I)模型）。选取了一个代表单元建模，底部采用固定边界，两侧需要设置为周期边界，对粒子施加一个斜重力对应于倾斜的平面，如下中间图： 边界设置如右图：有颜色的是物质点，代表物质点x方向的速度，没有颜色的是背景网格。目前采用的求解格式为MUSL，底部蓝色的固定边界实现方式为：该边界上网格的速度和加速度设为0；红色的周期边界实现方式为：在粒子将质量，力，动量投射到网格后，依次找到periodic1和periodic2上除外同一高度的两个对应节点，两两节点的动量，力和质量相加。这样每组网格都有同样的初速度，加速度。然后求解动量方程，再将网格上的信息映射回物质点，更新物质点的坐标和速度，如果物质点x坐标越过了periodic2，则将其坐标减掉两个周期边界的x距离，使物质点从另一侧进入。 我通过这样的处理方式一直无法得到一个正确稳定的速度剖面。所以想请教各位路过的大佬！帮忙看看是哪里出了问题。感激不尽！！！！