Abstract. We analyze a model of polyelectrolyte gels that was proposed by the authors in previous work. We first demonstrate that the model can be derived using Onsager's variational principle, a general procedure for obtaining equations in soft condensed matter physics. The model is shown to have a unique steady state under the assumption that a suitably defined mechanical energy density satisfies a convexity condition. We then perform a detailed study of the stability of the steady state in the spatially one-dimensional case, obtaining bounds on the relaxation rate. Numerical simulations for the spatially one-dimensional problem are presented, confirming the analytical calculations on stability.1. Introduction. Gels are crosslinked, three dimensional polymer networks that absorb solvent and swell without dissolution [21,22,16,39,2]. In this paper, we study polyelectrolyte gels, in which the polymer network carries charge, and delivers counterions to the solvent. An important feature of many polyelectrolyte gels is that they undergo large and often discontinuous volume changes (called the volume phase transition) in response to various environmental parameters, including pH, temperature and ionic composition [14,37,38]. Some physically interesting characteristics of the volume phase transition include robust hysteresis, coexistence of phases and complicated transient dynamics [12,37]. These volume changes are at the basis of numerous applications of gels to artificial devices [3,31,29,7,13,19] and are thought to underlie certain physiological processes [46,41,37,29]. It is thus of both practical and theoretical interest to develop a dynamic model of polyelectrolyte gels.There have been numerous modeling studies of gels. Studies using purely mechanical models include [36,40] in which static problems are addressed, and [5,44,35,9,10,24,25] in which dynamic problems are the focus. Models of polyelectrolyte gels include the static models of [37,20,34] and the dynamic models of [17,15,18,26,28,43,42,45,27,6,1,20]. In this paper, we initiate a detailed study of a dynamic PDE model of polyelectrolyte gels introduced in [30]. Our model is distinguished from previous models in its careful derivation of the boundary conditions at the gel-fluid interface as well as its satisfaction of a free energy identity, including interfacial terms. We show that our model can be derived from Onsager's variational principle which is a systematic way of deriving equations in soft condensed matter systems [10,33]. The one-dimensional stability calculation is a generalization of the classical work by [40] in which the neutral gel case is studied. We prove the uniqueness of the steady state solution and study the stability as the model being restricted to one spatial dimension. Finally, we have successfully simulated our polyelectrolyte gel model in the one spatial dimension. The simulated equilibrium solution and exponential decay rate both match our analytical calculation.The paper is organized as follows. In section 2, we apply Onsager'...