The computational workflow starts with the simulation of the slurry phase. To this end, it is known that the rheological properties of the slurry have a great impact on the manufacturing process. The quality and speed of the coating, viscosity and the resulting cathode microstructure all depend on the slurry microstructure and its physiochemical properties. Many studies demonstrate the interdependency between electrode characteristics and the properties of the slurry such as the porosity, conductivity and geometrical tortuosity. The slurry microstructure model is generated based on the experimental weight percentages of 69.0%, 27.6% and 3.4% for AM, SE and CBD, respectively. The diameter of the AM and SE are dependent on the starting materials as well, where the diameter of the AM is 10 µm and the particle size distribution of the SE constitutes of seven different diameters (1.5 µm, 2.5 µm, 3.5 µm, 4.5 µm, 5.5 µm, 6.5 µm and 7.5 µm), the methodology used to obtain these values is explained in the Experimental Section. To avoid overlapping of particles at the initial state, they are randomly placed in a very large simulation box (716 µm x 716 µm x 1000 µm). The initial velocities of the particles were set randomly according to a distribution that corresponds to a temperature of 300 K. For the initial generation of the particles and the slurry equilibration, the number of particles, pressure and temperature are set to be 300 K and 1 atm, respectively (which is known as NPT environment), throughout the CGMD simulation. After that, nonequilibrium CGMD (nE-CGMD) simulations are run to evaluate the shear viscosity (η) of the system against the applied shear rates (γ).