I am trying to simulate the compaction process of hot mix asphalts and want to do this with mathematica's FEM libraries. How can I define inclusions in meshes (i.e aggregate particles in the binder) while using the FEM and NDSolve capabilities of Mathematica? I am new to FEM analysis so my question may be a bit rudimentary.

  • $\begingroup$ Have a look at the documentation of ToElementMesh and/or the mesh generation tutorial. If you do not find anything there then you might need to clarify you question a bit. $\endgroup$ – user21 Jul 28 '16 at 12:26