The Flory-Rehner theoretical description of the free energy in a hydrogel swelling model can be broken into two swelling components: the mixing energy and the ionic energy. Conventionally for ionized gels, the ionic energy is characterized as the main contributor to swelling and, therefore, the mixing energy is assumed negligible. However, this assumption is made at the equilibrium state and ignores the dynamics of gel swelling. Here, the influence of the mixing energy on swelling ionized gels is quantified through numerical simulations on sodium polyacrylate using a Mixed Hybrid Finite Element Method. For univalent and divalent solutions, at initial porosities greater than 0.90, the contribution of the mixing energy is negligible. However, at initial porosities less than 0.90, the total swelling pressure is significantly influenced by the mixing energy. Therefore, both ionic and mixing energies are required for the modeling of sodium polyacrylate ionized gel swelling. The numerical model results are in good agreement with the analytical solution as well as experimental swelling tests.
Keywords: Flory–Huggins equation; Flory–Rehner theory; finite deformation; ionized hydrogels; mixing energy; polymer mechanics; sodium polyacrylate; superabsorbent polymer; swelling.