粒子穿过可变形界面的问题

[凤雏岂能消逝]

大家好，我在课题中遇到大量颗粒穿过一个2D可变形界面的case。我需要计算穿过界面的法向速度，来计算流通量。界面通过immersed boundary method 来表达，我的想法是：穿过界面上相邻lagrangian点做矢量为基准求每个颗粒的相应的法向速度。但总的颗粒量500万左右，这个方法凭空想想，需要的耗费很多资源，时间上不允许。
可否有大神能提供一些思路，或者建议，万分感谢！

[bestucan]

是否可以通过一个时间步后粒子的坐标相对于界面的位置来判断粒子是否穿过界面

https://www.baidu.com/s?wd=判断点在线的某侧

给粒子在界面某侧声明个变量和一个flag，每次修改变量的时候给flag扶植1，然后每次计算完统计flag为1的粒子数量，同时给flag赋值0.

然后看每个时间步统计的粒子数量差就知道通量了

[凤雏岂能消逝]

@bestucan 真的非常感谢您的建议。对于二维的case尝试了一下，我设置一个不变形的界面，进入filter区域的粒子数量大约是等同于流通量。但我手上的计算资源和代码的优化方面并不乐观，耗时久。如果界面变成可变形的情况，计算时间应该会加倍，时间上等不起。