lammps中怎么添加内部边界条件？

Question

目的：想用dem颗粒模拟气泡的上升和变形。
因为之前有人用mdpd做过，所以想用dem颗粒做一做，用一组颗粒代表周围的液体，另一组颗粒代表气泡，现在的问题是，一开始计算，气泡颗粒和液体颗粒就会混合在一起。个人认为是缺少气液之间的边界条件去阻止互相渗透，lammps中好像没有这种体系内部的边界条件可以用。
有没有人知道怎么在lammps中添加这种边界条件啊？
是不是需要更改源代码，如果需要的话，我应该从哪儿入手呢。因为不熟悉c++，自己直接写肯定很困难，所以想问问有没有可以直接拿过来改的

Answer 1

不知道你这个问题是否已经解决, 如果还没有解决的话:

• 以下回答来自chatgpt:

  对于在LAMMPS中添加用于阻止气液渗透的内部边界条件的问题，我无法给出具体的解决方案。这需要对LAMMPS的源代码和相关领域的知识有深入的了解才能进行修改和添加这种边界条件。若要修改源代码，你需要熟悉C++语言和对LAMMPS的代码结构和算法有一定的了解。

  然而，我可以给出一些推荐的解决方案：

  1.查看LAMMPS的文档和论坛，寻找是否已经有人讨论过类似的问题，并提供了解决方法或代码示例。

  2.与其他研究人员和开发者进行交流，了解是否有现成的插件或模块可以添加这种边界条件。

  3.考虑使用其他分子动力学软件，如ESPResSo，它可能有现成的内部边界条件的实现。

  总之，解决这个问题需要对LAMMPS有深入的了解和编程能力。如果你不熟悉C++和LAMMPS的代码结构，可能需要寻求其他专家的帮助或尝试其他分子动力学软件。

如果你已经解决了该问题, 非常希望你能够分享一下解决方案, 写成博客, 将相关链接放在评论区, 以帮助更多的人 ^-^