The objective of the present article is to go further in the simulation of the above extrusion process. Such simulations are challenging notably due to the complex 3D geometry of the TSE as well as due to the choice of a theoretical model used for the simulations. There are indeed several possible choices for the model but in any case, a compromise needs to be found between a refined description (both from the point of view of the number of involved physical actors and on the degrees of freedom of the discretization) and a doable computational cost. The above-mentioned complexity is due in particular to the inherent nature of extrusion paste components, which include constituents that vary in size by several orders of magnitude, so that implies the need for models able to include some sort of multiscale description. Building upon our previous work in , we adopt a DEM method thanks to this already present expertise in the team and the capability to deal with the multiscale nature of this extrusion process. The ultimate goal of this simulation study being naturally to try to optimize the lab experiments in terms of extrusion process parameters, material properties, and the microstructure of extruded filaments. In the present work, we report a more advanced version of our previous model, which is much closer to the original experimental design of Boudeville et al. . We replaced the geometry of the extruder to fit the shape of the physical equipment, applied four material formulations for the particles, and added a recirculation circuit to follow the experimental process, which is also helpful for the DEM model's computational cost. We also explored two injection sequences and analyzed the influence of screw rotation speed on the microstructure of the filaments, which are the outlets of the extrusion process.