I am working on coupling the fiber model with CFD and have hit a snag with a simulation.
I have a long fiber fixed at one end (using fix freeze) while all other atoms are free to interact.
I have 0.5 m/s air traveling up from the bottom of the domain and all other boundaries are zero gradient.
After a few moments an explosion occurs at the spheres and I do not know where to start to look to find the problem.
I have tried IB solvers and the piso solver but both cause this type of problem.
video of what is happening in the simulation
I thought maybe it was by grid size but I do not think this is the case as the grid in the video is fairly fine (far right side of video is the grid).
I don't think it has to do with the the bonding by its self. I know that care has to be taken with multisphere particles but there are no overlaps
in the spheres that are bonded.
Any ideas on where I should start looking would be a great help.
Attached are the input files. If you wish to run this, you will need my fiber model code off of github.