Solidification Dynamics of Rapidly Solidified Droplets in Battery Powder Synthesis

Atomization

High-velocity gas jets shear liquid metal streams into fine droplets inside the spray chamber. Controlling droplet size during this initial breakup establishes the baseline for thermal extraction rates, phase evolution, and final particle geometry in battery active material precursor powders.

Kinetic Energy and Droplet Disruption

Primary break-up begins where aerodynamic drag overcomes the cohesive surface tension of the liquid filament. High relative velocity between the atomizing gas (typically argon or nitrogen) and the melt stream creates primary capillary waves that amplify along the liquid cylinder until Rayleigh-Taylor instabilities divide the stream into discrete primary ligaments.

Secondary break-up follows immediately as aerodynamic pressure deforms these primary ligaments into thin membranes that burst into micro-droplets. The Weber number quantifies this deformation propensity as the ratio of inertial drag forces to surface tension forces across the droplet diameter:

We = (rho_gas v_rel^2 d_droplet) / sigma_melt

When the local Weber number exceeds a critical value of 12 for low-viscosity precursor melts, bag breakup occurs. Above a Weber number of 80, the mechanism transitions to shear stripping, where boundary-layer shearing pulls smaller droplets from the periphery of larger molten cores.

1. The liquid precursor stream exits the delivery nozzle orifice at controlled volumetric flow rates under pressive pressure.
2. Supersonic inert gas streams impact the liquid column at defined intersection angles between 20 degrees and 45 degrees.
3. Primary stream distortion generates elongated liquid sheets that thin under continuous shear forces.
4. Capillary wave expansion breaks liquid sheets into primary filaments and thread-like structures.
5. Secondary aerodynamic breakup disintegrates filaments into spherical droplets ranging from 5 micrometers to 50 micrometers in diameter.
6. Inert shroud gas conveys airborne droplets through the cooling column to suppress inter-particle collision and agglomeration.

Molten Stream Breakup Mechanics

Instabilities along the jet surface grow exponentially until capillary action pinches off discrete liquid spheres. Disintegration dynamics depend heavily on the Ohnesorge number, which relates viscous forces to inertial and surface tension forces:

Oh = mu_melt / sqrt(rho_melt sigma_melt d_droplet)

Operating at Ohnesorge numbers below 0.1 alongside high Weber numbers places the melt in the atomization regime, triggering immediate breakup at the nozzle tip to generate the fine precursor droplets needed for high-density lithium nickel manganese cobalt oxide cathode powders. Satellite formation occurs when fine secondary droplets form during thread necking, subsequently adhering to larger solidifying particles and degrading powder flowability.

Droplet flight distance before impingement determines whether powder particles weld together or remain discrete units.

High atomization gas pressure always yields finer median particle distributions at the expense of higher satellite attachment rates.

Chill

Rapid thermal extraction from airborne droplets suppresses equilibrium phase transformation during battery precursor synthesis. The heat removal rate dictates whether the freezing droplet solidifies into an amorphous structure, a nanocrystalline matrix, or a micro-segregated crystalline lattice.

Thermal Boundary Resistance at the Gas Interface

Heat removal rates depend primarily on convective heat transfer coefficients through the surrounding inert gas shroud, which dominates thermal extraction during flight. The Nusselt number correlates convective transfer to thermal conductivity of the inert gas shroud:

Nu = 2 + 0.6 (Re^(1/2)) (Pr^(1/3))

Small droplets with diameters below 20 micrometers exhibit Biot numbers well below 0.1, indicating that internal thermal resistance within the droplet is negligible compared to external boundary resistance. The Biot number calculates as:

Bi = (h_conv d_droplet) / (6 k_droplet)

Because internal temperature gradients remain flat under low Biot conditions, Newtonian cooling governs heat extraction, allowing droplets to solidify isothermally from the exterior inward or through uniform volumetric undercooling.

Nucleation Undercooling and Front Velocity

Liquid supercooling drives solid crystal phase transformation across the advancing interface. As the melt drops well below its liquidus temperature, homogeneous nucleation becomes thermodynamically favorable, and the degree of undercooling directly sets interface velocity:

v_front = mu_interface (Delta_T)^n

Where mu_interface represents interface kinetic mobility and Delta_T represents total thermal undercooling. Rapid interface movement during high undercooling conditions traps solute atoms within the advancing crystal matrix, suppressing equilibrium solute rejection. The release of latent heat of fusion during recalescence raises the local temperature behind the solidification front, temporarily slowing front advance until convective cooling dissipates the thermal wave.

Cooling rates exceeding 100,000 K per second suppress dendritic solute partitioning in nickel-rich precursor droplets.

• Inter-particle Sintering occurs when droplets strike chamber walls before completing bulk solidification, forming dense metallic crusts that destroy powder sizing distribution.
• Gas Entrapment Voids form when atomizing gas dissolves into the molten droplet under pressure and becomes trapped as internal pores during ultrafast planar front advance.
• Asymmetric Grain Growth develops when directional convective gas flow induces uneven temperature distribution across the droplet surface during nucleation.
• Amorphous Phase Decomposition happens when local cooling rates fall below critical threshold levels, triggering unintended phase separation during subsequent calcination steps.

Insufficient cooling rates allow secondary phase precipitation during freezing, which degrades discharge capacity and accelerates capacity decay across battery cycle life.

Partition

Rapid dendritic growth generates steep concentration gradients in multi-element transition metal droplets. In cathode materials such as LiNi0.8Mn0.1Co0.1O2 or high-entropy olivine phosphates, non-uniform distribution of nickel, manganese, and cobalt across individual powder grains induces micro-cracking during electrochemical cycling.

Solute Redistribution in Multicomponent Systems

Transition metal ions migrate across the advancing solid-liquid boundary according to non-equilibrium distribution coefficients. The non-equilibrium solute partition coefficient k depends directly on interface velocity v_front:

k(v) = (k_0 + (v_front / v_diff)) / (1 + (v_front / v_diff))

Where k_0 represents equilibrium partition coefficient and v_diff represents diffusive speed of solute species across the interface boundary. As the solidification front approaches the solute diffusive speed, k approaches unity, suppressing compositional divergence between the solid core and interdendritic liquid.

Will Ultrasonic Atomization Reduce Dendritic Segregation?

High-frequency acoustic agitation disrupts growing dendrite arms by inducing severe localized micro-cavitation within liquid spheres. Secondary dendrite arm spacing lambda_2 shrinks as cooling rate increases according to a power-law empirical relation:

lambda_2 = A (dT / dt)^(-n)

Where A represents material constant and exponent n ranges between 0.3 and 0.45 for transition metal oxide precursors. Refining secondary dendrite arm spacing below 100 nanometers shortens the diffusion distance required for chemical homogenization during subsequent solid-state calcination, allowing shorter thermal dwelling times that avoid lithium evaporation while achieving fully ordered layered structures.

Supplier compliance with particle sizing standard ISO 13320 shifts the financial burden of out-of-spec powder fractions to the chemical refiner.

1. Audit Raw Solute Ratios by measuring raw feed liquid stoichiometry via inductively coupled plasma optical emission spectroscopy prior to atomization stream injection.
2. Monitor Solidification Velocity using high-speed pyrometry to verify that liquid cooling rates exceed critical solute trapping thresholds.
3. Inspect Dendritic Arm Spacing through cross-sectional scanning electron microscopy on representative powder samples from each production lot.
4. Calculate Homogenization Index by assessing elemental variance using energy-dispersive X-ray spectroscopy line scans across grain cross sections.
5. Reject Out-of-Spec Batches whenever local transition metal segregation ratios exceed 3 percent relative variance across primary grain cores.

Minor transition metal segregation does not homogenize completely during downstream high-temperature calcination without extending thermal dwell times.

Sphericity

Surface tension forces contract deformed droplets into spheres before solid phase conversion locks their geometry in place. Highly spherical powder particles maximize packing density and improve slurry viscosity.

Rheological Dynamics and Packing Density

Spherical particles reduce inter-particle friction in slurries during roll-to-roll electrode casting. Irregular geometries require higher solvent volumes, depressing the achievable solids loading in cathode slurries, which in turn raises drying energy demands and worsens binder migration during solvent evaporation.

Sphericity Psi calculates as the ratio of surface area of an equivalent volume sphere to actual surface area of the particle:

Psi = (pi^(1/3) (6 V_particle)^(2/3)) / A_particle

When particle sphericity drops below 0.85, powder tap density decreases significantly. Low tap density degrades volumetric energy density in finalized battery cells by reducing active material mass per unit volume of coated electrode foil.

Surface Tension Driven Spheroidization Window

Capillary stress drives fluid motion from high-curvature lobes toward depressed regions across the liquid surface. The characteristic time required for a distorted liquid filament to relax into a spherical droplet depends on viscosity, surface tension, and droplet diameter:

tau_relaxation = (5 mu_melt d_droplet) / (8 sigma_melt)

Complete spheroidization occurs only when relaxation time remains shorter than the total solidification time required for liquid phase freezing:

tau_solidification = (rho_melt L_fusion d_droplet) / (6 h_conv (T_melt – T_ambient))

Consider a worked construction evaluating high-nickel NMC precursor droplets. Take a droplet batch with average diameter d_droplet = 25 micrometers, melt viscosity mu_melt = 0.005 Pa s, surface tension sigma_melt = 1.2 N/m, latent heat L_fusion = 280,000 J/kg, density rho_melt = 4,500 kg/m^3, convective heat transfer coefficient h_conv = 4,500 W/m^2 K, and temperature difference Delta_T = 1,100 K.

Calculating relaxation time yields:

tau_relaxation = (5 0.005 25 10^-6) / (8 1.2) = 6.51 10^-8 seconds

Calculating solidification time yields:

tau_solidification = (4,500 280,000 25 10^-6) / (6 4,500 1,100) = 1.06 10^-3 seconds

Because solidification time exceeds relaxation time by four orders of magnitude, surface tension forces fully round the liquid droplets before phase transformation locks in particle shape. Assume now an atomization scenario producing larger 120-micrometer agglomerated droplets under reduced gas pressure where h_conv drops to 800 W/m^2 K. Solidification time extends to 2.1 10^-2 seconds, but internal thermal gradients induce rapid skin formation before surface tension flattens high-aspect lobes, resulting in low-sphericity particles with Psi = 0.76.

High particle sphericity lowers slurry viscosity and permits higher active material solids loading during electrode coating.

• Particle Size Distribution Curve certified by dynamic light scattering according to ISO 13320 standard specifications.
• Morphology Characterization Dossier including automated static image analysis parameters defined under ISO 9276-6 guidelines.
• Specific Surface Area Certification measured via multi-point Brunauer-Emmett-Teller gas adsorption analysis per ISO 9277 rules.
• Tap Density Verification Certificate established through standard mechanical tapping procedures outlined in ASTM B527 testing protocols.
• Internal Porosity Quantification verified by mercury intrusion porosimetry or cross-sectional X-ray computed tomography mapping.

Incorporating standard compliance terms under ISO 9276-6 fixes the minimum acceptable sphericity index at 0.88, transferring the commercial risk of poor slurry rheology back to the powder vendor.

Audit

Chemical and physical characterization across powder sizing batches establishes baseline performance stability for commercial battery manufacturing. Purchasing specifications must enforce strict metrics on particle size distribution metrics, compositional tolerances, and tap density values across incoming shipments.

Specification Tolerances for Active Material Precursors

Stoichiometric variations across particle size fractions alter secondary phase formation during high-temperature calcination. Finer powder fractions under 5 micrometers experience higher relative cooling rates during atomization, driving higher solute retention and lower surface segregation relative to coarse fractions exceeding 30 micrometers. Blending asymmetric fractions leads to inconsistent lithium absorption during calcination, leaving unreacted lithium carbonate residues on particle surfaces.

Commercially delivered lots require tight particle size distribution spans defined as:

Span = (D_90 – D_10) / D_50

Spans exceeding 1.4 indicate wide size distributions that introduce packing inconsistencies during electrode calendering. High calendering pressure applied to non-uniform powder layers breaks coarse particles while under-compacting fine regions, causing localized impedance variations across the coated electrode sheet.

Verification Protocols for Powder Batches

Statistical sampling plans enforce strict chemical limits before precursor powder lots enter calcination kilns. Incoming inspection protocols isolate samples across top, middle, and bottom sections of shipping containers to detect size segregation induced by vibration during transit. Sieve analysis coupled with inductively coupled plasma testing identifies compositional drift across particle size cuts.

Batch acceptance criteria define strict upper and lower control limits for critical quality attributes:

C_p = (USL – LSL) / (6 sigma_process)

Process capability indices C_p exceeding 1.33 ensure that powder generation processes run within operational limits without generating off-spec material. Powder lots failing process capability benchmarks face immediate quarantine or price-discounted rework processing.

Whether sub-micron precursor fines generated during rapid solidification can be fully re-melted without introducing chemical contamination remains an open operational question for high-volume powder refiners.