2025 AIChE Annual Meeting

(391f) Investigation of the Modelling of Precipitation of Ni-Mn-Co Hydroxides.

The world efforts towards achieving net zero include the development of several technologies aiming at the sustainable production and storage of energy. A share of these efforts is focused on the development of Li-ion batteries, whose utilisation has grown in multiple applications like electric portable devices, EVs and storage due to their ability to store and deliver energy over cycles of use. The performance and cost of the batteries relates to the properties of the cathode materials. Several chemistries have been researched to improve the performance of cathode active materials including the use of lithium-iron phosphate, nickel-manganese-cobalt and nickel-manganese-aluminium. In this work, we look at the synthesis of combined hydroxides of nickel-manganese-cobalt via co-precipitation, particularly the mathematical modelling of this process.

We investigate a semi-batch precipitator applying global sensitivity analysis to determine the relevance of various crystallization processes (e.g. nucleation, growth, agglomeration and breakage) over the particle size distribution and mean particle sizes.

The mathematical model relies on the application of the population balance equation along the mass balances of the dissociated species in solution. The main assumptions of the model are:

  1. The precipitated crystals are spherical.
  2. No seeds in the initial solution or the feed are present.
  3. The precipitator is well mixed.
  4. The precipitation of individual transition metal hydroxides does not occur.
  5. The reaction rate of both the ionic complexes and the hydroxide formation are fast; the controlling mechanisms of the process are the nucleation and growth of the crystals.
  6. The crystalliser operates isothermally.

In addition, the model considered different cases regarding growth and agglomeration/breakage: size independent and size-dependent growth, no agglomeration/breakage, and agglomeration/ breakage occurring. With the previous assumptions and considerations, the model is presented by Equations (1)-(24).

The model equations were implemented and solved by the method of lines for a total of 250 simulations where the values of main operating parameters and the kinetic constants were randomly sampled using the Latin hypercube technique. Figure (2) and (3) are preliminary results of the particle size distribution and the volume-weighted particle size d[4,3] for a few selected cases among considering the model with size-independent growth and no agglomeration/breakage. The effect of the input parameters over the response of the model is visible; in general, for this case, the sensitivity indices (Table 1) obtained from GSA suggested that the model response is primarily sensitive to the growth kinetics, the solution pH and the initial concentration of ammonia. Less relevant is the kinetics of nucleation, whereas the concentration of ammonia feed seems to be irrelevant. Next steps are to analyse the effects of size-dependent growth, and agglomeration/breakage.