The linear double-diffusivity (D-D) model of Aifantis, comprising two coupled Fick-type partial differential equations and a mass exchange term connecting the diffusivities, is a paradigm in modeling mass transport in inhomogeneous media, e.g., fissures or fractures. Uncoupling of these equations led to a higher order partial differential equation that reproduced the non-classical transport terms, analyzed independently through Barenblatt’s pseudoparabolic equation and the Cahn–Hilliard spinodal decomposition equation. In the present article, we study transport in a nonlinearly coupled D-D model and determine the regime-switching of the associated diffusive processes using a revised formulation of the celebrated Lux method that combines forward Fourier transform with a Laplace transform followed by an Inverse Fourier transform of the governing reaction–diffusion (R–D) equations. This new formulation has key application possibilities in a wide range of non-equilibrium biological and financial systems by approximating closed-form analytical solutions of nonlinear models.