求助：lammps中如何实现在体系内部添加边界以阻止不通原子之间互相渗透？

李海波

目的：想用dem颗粒模拟气泡的上升和变形。
因为之前有人用mdpd做过，所以老师想让我用dem颗粒做一做（当然能不能成功还需要验证），
用一组颗粒代表周围的液体，另一组颗粒代表气泡，不做任何处理，运行一段时间后，气泡颗粒和液体颗粒就会混合在一起，如图（周围代表液体的颗粒设置成透明了）。
应该是缺少气液之间的边界条件（蓝色圆圈处）去阻止互相渗透，lammps中好像没有这种体系内部的边界条件可以用。考虑过fix wall/region，但是region不能跟随气泡（中间的颗粒团）发生变化

我应该怎么在lammps中添加这种边界条件啊？
是不是需要更改源代码，如果需要的话，我应该从哪儿入手呢。因为不熟悉c++，自己直接写很困难，所以想问问有没有可以直接拿过来改的