Issue with Read_Data command during coupled simulation only

Submitted by RussellHolt on Wed, 12/08/2021 - 12:03

Hi all

I am attempting to simulate an unresolved coupled CFDDEM channel flow of particles using a fixed particle boundary for the walls instead of a meshed boundary. Currently I am running into an issue when using the READ_DATA command to bring in the wall particles, I have set the wall particles as type 2, and my flowing particles as type 1.

I have used two DEM input scripts, first an initialization (in.Test_init) where the wall particles and flowing particles are read in, and then the running simulation (in.Test). When running both of these scripts in LIGGGHTs individually there seems to be no problem, likewise running the CFD code separately also works fine. Additionally if I simply remove the READ_DATA command and omit the 2nd atom type the simulation runs fine.

I think I must missing something as I cannot spot the error here, can anyone please help? I have attached the log file with the error message (MPI message that confuses me!), by printing out a simple "here" in the input script I can confirm the error occurs when the script hits the 'run 100' line.

I'm hoping I have simply made a mistake in the data file.
Note: I have moved the generated restart file into the appropriate folder to be read before running, I have run the _init file prior to running the coupled simulation which only using the in.Test file.

RussellHolt | Wed, 12/08/2021 - 13:46

Issue was that the atoms in the data file extended beyond the periodic domain, this apparently was not an issue for LIGGGHTS which surprises me, but once coupled it caused problems.

By this I mean my simulation box (in the x direction) went from 0 to 50, and I had particles at both 0 and 50.

One to look out for