Equation (1) is also applied to solvent diffusion in polymers but more recent studies have demonstrated different diffusive behaviors for imbibition. Due to effects of polymer physical properties and polymer/solvent interactions, Alfrey and coworkers have proposed a semi-empirical classification based on polymer relaxation rates, [6] identifying the Fickian and the non-Fickian regimes. According to this model, the amount of solvent absorbed by a polymer membrane follows a power law and the power exponent is related to the transport mechanism.The transport mechanism is related to the bulk properties of the polymer, in particular the physical state. In the rubbery state, when the glass transition temperature T g is above the working temperature T w , molecules show high mobility and it facilitates solvent penetration across the polymer matrix. On the other hand, glassy polymers, with T g < T w , are more rigid and solvent absorption is less easy. [7] In some polymers, the glassy-to-rubbery transition is physically similar to that occurring in hydrated matrices. Water acts like a plasticizer and affects diffusion properties within the polymer. This plasticizer effect leads to nonlinear diffusion, and the diffusion coefficient D can show a dependency on concentration. [8] The identification of the interface between swollen and unswollen regions, and its evolution with time, can help in the evaluation of nonlinear diffusion effects. This internal boundary was experimentally tracked in gelatin beads, and its evolution with time significantly differed from that of the external bead boundary. [9] Mathematical models involving the movement of the internal boundary were developed in food science to study grain hydration [10] while in polymer membranes this study was simplified with a front-fixing transformation [11] and a comprehensive analysis of this aspect is still lacking.When soft materials imbibe a large amount of water their volume can increase. The swelling phenomenon is a kinetic process and it involves mass transport and mechanical properties. [12] In systems with variable volume, concentration profiles are described with the more advanced advection-diffusion model [13] that takes the moving front velocity v into accountNumerical models of swelling, in which the boundary of the domain needs to be calculated as part of the solution, are commonly called Stefan problems. The Stefan problem can describe Absorption Water diffusion and swelling of polymer films can be modeled by nonlinear diffusion equations with fixed or moving boundaries, investigating the Fick and the Stefan problems. To form a better view of the dynamics of the problems, this paper reports the mathematical modeling of water diffusion with linear and nonlinear diffusion coefficients by using the finite element method, and the description of boundary movements calculated with a fully implicit upwind scheme. Experimental results are used to support theoretical data. It is found that porosity strongly affects diffusion properties and proposed models de...