Add a new force

Submitted by ecmo on Fri, 03/23/2012 - 20:11

I'm interested in simulating magnetic particles in a media. How can I add this force as a source term in liggghts, which force can be derived analytically, f(x,y,z) or
solved by a FEA solver. For simplicity, I neglect a dipole-dipole effect between particles. Thanks in advance for your guidance.

Forums:

ckloss's picture

ckloss | Wed, 03/28/2012 - 17:30

Hi ecmo,

have a look at fix_addforce.h/cpp. Based on that, you can easily write your own fix that adds custom forces f(x,y,z) onto particles. In case of coupling to FE, the question is if this is highly dynamic or if a one-time FE analysis will do it. In the first case, you will need a fully coupled FE-DEM approach (like we already couple CFD-DEM), in the latter case having the FE solver write a file with the spatial distribution of the forces and having your fix addforce reading it will be sufficient. Good luck :-)

Cheers, Christoph