Fick's first law states that the diffusive flux J is proportional to the concentration gradient: J = -D (dC/dx). In the shrinking core model used here, the gradient across the boundary layer of thickness delta is approximated as (Cs - Cbulk) / delta, giving a flux J = D * (Cs - Cbulk) / delta. The total molar dissolution rate is then J multiplied by the crystal surface area (2 pi r for the 2-D circular model shown). As the radius r shrinks, the surface area decreases and the rate falls even if all other parameters stay constant.

Pharmaceutical tablets must dissolve at a controlled rate to release active ingredients in the bloodstream at the correct time and dose — this is governed entirely by Fick diffusion through the dissolution medium. Geologists use dissolution kinetics to model how limestone karst landscapes form over thousands of years. In construction, cement clinker particles dissolve in water to form calcium silicate hydrate, the binding phase that gives concrete its strength. Sugar and salt production both rely on optimized dissolution and crystallization cycles.

Stirring speeds up dissolution only when diffusion through the boundary layer is the rate-limiting step, which is the most common case for sparingly soluble compounds. When dissolution is surface-reaction limited rather than diffusion limited, stirring has little effect because the bottleneck is the detachment of ions from the crystal lattice, not their transport away from the surface. The simulation models the diffusion-controlled regime, where stirring compresses the boundary layer thickness and directly increases the diffusion flux.

The shrinking core model was formalized in chemical engineering by Octave Levenspiel in his landmark 1962 textbook "Chemical Reaction Engineering." The model describes a reacting or dissolving particle that shrinks from the outside inward as the unreacted or undissolved core gets smaller. Earlier theoretical groundwork was laid by Adolf Fick, who published his diffusion laws in 1855 based on experiments with salt diffusing through water, drawing an analogy with Fourier's law of heat conduction.

Crystal dissolution is the reverse of crystallization, and both are governed by the same thermodynamic driving force: the difference between actual and equilibrium concentration. Related phenomena include Ostwald ripening, where small crystals dissolve and redeposit on larger ones to reduce total surface energy; precipitation reactions, where dissolved ions recombine to form a new solid phase; and leaching in hydrometallurgy, where valuable metals are selectively dissolved from ores. The simulation is also conceptually linked to corrosion and electrochemical dissolution of metals.

Drug formulators engineer crystal habit, particle size distribution, and coating materials to achieve targeted dissolution profiles — for example, enteric coatings that resist dissolution in stomach acid but dissolve rapidly in the intestine. In semiconductor fabrication, selective wet etching dissolves specific crystal planes of silicon or gallium arsenide to pattern nanoscale features. Hydrometallurgical leaching tanks use agitation and temperature control to maximize the dissolution rate of copper or gold minerals from crushed ore, applying exactly the same Fick-diffusion principles shown in this simulation.

Researchers are using in situ atomic force microscopy (AFM) and synchrotron X-ray techniques to watch individual crystal faces dissolve at the nanometer scale in real time, revealing step-retreat and etch-pit dynamics that bulk models like the shrinking core cannot capture. Molecular dynamics simulations now resolve the detachment energy of individual ions from kink sites on the crystal surface. In pharmaceutical science, machine learning models are being trained to predict intrinsic dissolution rates from crystal structure alone, aiming to speed up early-stage drug candidate screening without physical experiments.