MaximaLabs

Electrowinning / electrochemical metal removal simulation

1 electrolyte -> spent electrolyte (regenerated acid), cathode metal, cell gas

Governing equations

The exact equations the solver works for a electrowinning / electrochemical metal removal — the same math shown in the app's "Theory" panel, not a black box.

n˙M=ηCE IzF\dot n_M = \eta_{CE}\,\frac{I}{zF}
Mz++z e−→M,2 H2O→O2+4 H++4 e−\mathrm{M^{z+}} + z\,e^- \rightarrow \mathrm{M},\qquad 2\,\mathrm{H_2O} \rightarrow \mathrm{O_2} + 4\,\mathrm{H^+} + 4\,e^-
Vcell=Erev+∣ηa∣+∣ηc∣+j δκ+Vhw,Erev=EO2/H2O∘−EMz+/M∘V_{cell} = E_{rev} + |\eta_a| + |\eta_c| + \frac{j\,\delta}{\kappa} + V_{hw},\quad E_{rev} = E^\circ_{O_2/H_2O} - E^\circ_{M^{z+}/M}
P=I Vcell,Esp=zFVcellηCEMMP = I\,V_{cell},\qquad E_{sp} = \frac{z F V_{cell}}{\eta_{CE} M_M}
cout=max⁡ ⁣(cin1+kmA/Q, e(Ec−E∘)zF/RT)(potentiostatic, CSTR)c_{out} = \max\!\left(\frac{c_{in}}{1 + k_m A/Q},\ e^{(E_c - E^\circ) z F/RT}\right)\quad(\text{potentiostatic, CSTR})
n˙M\dot n_M
metal deposited at the cathode [mol/s]
ηCE\eta_{CE}
current efficiency — fraction of the current that plates metal [0..1]
II
cell current [A]
zz
charge of the metal ion
FF
Faraday constant [C/mol]
VcellV_{cell}
cell voltage [V]
ErevE_{rev}
reversible decomposition voltage from the CRC standard potentials [V]
ηa\eta_a
anode overpotential [V]
ηc\eta_c
cathode overpotential [V]
jj
current density I/A [A/m^2]
δ\delta
anode-cathode gap [m]
κ\kappa
electrolyte conductivity [S/m]
VhwV_{hw}
contact and bus-bar voltage [V]
PP
electrical power drawn [W] (reported as duty)
EspE_{sp}
specific energy per kg of metal [J/kg]; also shown in kWh/kg
MMM_M
molar mass of the metal [kg/mol]
cinc_{in}
metal-ion molarity entering [mol/L]
coutc_{out}
metal-ion molarity leaving [mol/L]
kmk_m
mass-transfer coefficient to the cathode [m/s]
AA
cathode area [m^2]
QQ
electrolyte volumetric flow [m^3/s]
EcE_c
cathode potential [V vs SHE]
E∘E^\circ
standard reduction potential of M^z+/M [V vs SHE]
RR
gas constant [J/mol/K]
TT
cell temperature [K]

Parameters

Faraday's-law electrowinning (cathode M^z+ + ze- -> M; anode 2H2O -> O2 + 4H+ + 4e-; the protons leave as regenerated acid); outlets in order spent (electrolyte), cathode (solid metal), gas (O2 + H2); exactly one of: current [A], current_density [A/m^2, with cathode_area], recovery [0..1 of target_metal, current solved], cathode_potential [V vs SHE, with mass_transfer_coefficient and cathode_area: potentiostatic mass-transfer-limited impurity removal; every allowed metal whose Nernst potential is above it plates at its limiting current and never below the Nernst floor]; cathode_area [m^2], mass_transfer_coefficient [m/s], reactor_model ['cstr' or 'plug_flow'], target_metal [metal id]; current_efficiency [0..1 or {metal: share}; REQUIRED, no default], side_reaction ['hydrogen' (default, H2 at the cathode) or 'redox_shuttle' (e.g. Fe3+/Fe2+, no net chemistry)]; cell voltage: cell_voltage [V] or built from CRC standard potentials plus anode_overpotential [V], cathode_overpotential [V], electrolyte_conductivity [S/m], electrode_gap [m], hardware_voltage [V] (none defaulted); optional metals [list of ions allowed to plate; ag cu pb sn ni co (cobalt) cd fe2 (ferrous) zn mn], acid_component [default h2so4 / phosphoric_acid / hno3 if present], protons_per_acid [2 / 1 / 1], temperature [K, isothermal; default feed T]; reports power [W, = duty, electrical work], specific_energy_kwh_per_kg, reversible_voltage, voltage_efficiency, current_efficiency_effective, acid_generated_mol_s, o2_rate, h2_rate

Example flowsheets that use it

Stop fighting legacy software. Build your first flowsheet in 60 seconds.