Particle periodicity in parallel

[Rookie]

Dear OpenLB developers,

After I successfully modified the periodic boundary of the particles, I realized that this boundary can only be used in non-parallel when added.If the periodic boundaries of the particles are added in parallel, an error will be reported to terminate the simulation after the particles have been moving for some time.

terminate called after throwing an instance of ‘std::domain_error’
what(): read only access to data which is not available locally
[banana:883562] *** Process received signal ***
[banana:883562] Signal: Aborted (6)
[banana:883562] Signal code: (-6)
[banana:883562] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x14420)[0x149ca4953420]
[banana:883562] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x149ca432500b]
[banana:883562] [ 2] /lib/x86_64-linux-gnu/libc.so.6(abort+0x12b)[0x149ca4304859]
[banana:883562] [ 3] /lib/x86_64-linux-gnu/libstdc++.so.6(+0x9e8d1)[0x149ca46dc8d1]
[banana:883562] [ 4] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xaa37c)[0x149ca46e837c]
[banana:883562] [ 5] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xaa3e7)[0x149ca46e83e7]
[banana:883562] [ 6] /lib/x86_64-linux-gnu/libstdc++.so.6(+0xaa699)[0x149ca46e8699]
[banana:883562] [ 7] ./channel3d(+0x5535c)[0x55939831535c]
[banana:883562] [ 8] ./channel3d(+0x555b5)[0x5593983155b5]
[banana:883562] [ 9] ./channel3d(+0x53f61)[0x559398313f61]
[banana:883562] [10] ./channel3d(+0x1be36)[0x5593982dbe36]
[banana:883562] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x149ca4306083]
[banana:883562] [12] ./channel3d(+0x1cf5e)[0x5593982dcf5e]
[banana:883562] *** End of error message ***
————————————————————————–
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
————————————————————————–
————————————————————————–
mpirun noticed that process rank 1 with PID 0 on node banana exited on signal 6 (Aborted).
————————————————————————–

I modified this file (olb-1.6r0/src/particles/subgrid3DLegacyFramework/boundaries/periodicBoundary3D.h).My guess is that updating the position after the particle reaches the boundary in parallel causes problems with data exchange. So I realized that this code doesn’t have a setting in parallel, and I thought that I couldn’t get some insight from this by learning how the periodicity of the fluid and the periodicity of the ParticleManager are set up in parallel. My current setup allows me to simulate implementing the periodicity of particles slowly in non-parallel and turning off the periodicity of particles can be simulated quickly in parallel, thanks again for your friendly help.

Regards,
Rookie

[mathias]

Dear Rookie,

if the Spring School doesn’t work for you, then maybe becoming a member of the consortium (https://www.openlb.net/consortium/) is an alternative option, because unfortunately I can’t provide such detailed support here.

Best regards,
Mathias

[Author]

Posts