(2) . In order to give a more fluid-like solid behavior for the BC particles, low Young's and high Poisson ratio are defined for this particle type. Experimental measurements show that PP has a very high viscosity (800 Pa.s) during the extrusion (Fig. in SI). In addition, solid carbon nanofibers are part of the BC particle. All this mitigates the inaccuracy of representing this phase as a solid. The cohesion parameter value is chosen to achieve a dense suspension flow, which is assessed by visual inspection of the trajectories. The impact of this value is investigated in the next section. A calibration of these values with experimental data is needed in future work. Rolling friction is also considered by employing the constant directional torque (CDT) model. The Young's Modulus and the Poisson ratio of LFP is found in literature. Using these parameters, the Rayleigh time step can be estimated. The immediate choice of a timestep equal to 20 % of Rayleigh time, as suggested in literature, results in a very long simulation time for filling the extruder. Higher timesteps yielded numerical instability in the simulations. Therefore, a decrease in 3 orders of magnitude of LFP Young's modulus is necessary in order achieve feasible running times and stability. A similar reduction changed the quantitative results for DEM simulations of cohesionless particles but not the trend, which is the essential for our comparative analysis in this work. The same decrease was applied on the Young's modulus of all materials to minimize the effect on the real particle dynamics. The values of the model parameters are given in Table of the SI. Our simulations, involving around 200 000 particles, take between 6-15 days depending mainly on the chosen rotation speed of the extruder. We used 64 cores AMD EPYC 7513 @ 2.60GHz (256 GB of RAM) of the MatriCs platform (Université de Picardie-Jules Verne, France). The following open-source software were used: LIGGGHTS for the DEM simulations, OVITO for visualization and MeshLab for modifying the extruder geometry. The commercial software GeoDict (Math2Market) was employed for post-processing.