embolism.cpp seems to run without problems.
When I try adding an object (a cylinder), so that I could see the flow around a rigid object, the simulation crashes if I use "forceCoupling3D(...)" in the cpp code. For some reason (probably due to high density or the Peskin stencil method) the calculated forces (fsi[] in ibm3D.hh) diverge to infinity.
Additional info. I created an input file with random particles distributed by the inlet and a cylinder in the middle. The particles are integrated using fix nve, while the cylinder is stationary (fix move linear 0.0 0.0 0.0). If I only use the random distribution of particles without the cylinder the simulation runs fine. Also if I try decreasing the density of the points on the cylinder the simulation runs fine.
embolism.cpp seems to run without problems. When I try adding an object (a cylinder), so that I could see the flow around a rigid object, the simulation crashes if I use "forceCoupling3D(...)" in the cpp code. For some reason (probably due to high density or the Peskin stencil method) the calculated forces (fsi[] in ibm3D.hh) diverge to infinity.
Additional info. I created an input file with random particles distributed by the inlet and a cylinder in the middle. The particles are integrated using fix nve, while the cylinder is stationary (fix move linear 0.0 0.0 0.0). If I only use the random distribution of particles without the cylinder the simulation runs fine. Also if I try decreasing the density of the points on the cylinder the simulation runs fine.
Any ideas?