Abstract. Many practically relevant polymers undergoing desorption change from the rubbery (saturated) to the glassy (nearly dry) state. The dynamics of such systems cannot be described by the simple Fickian diffusion equation due to viscoelastic effects. The mathematical model solved numerically is a set of two coupled PDEs for concentration and stress. Asymptotic solutions are presented for a moving boundary-value problem for the two states in the short-time limit. The solutions exhibit desorption overshoot, where the penetrant concentration in the interior is less than that on the surface. In addition, it is shown that if the underlying time scale of the equations is ignored when postulating boundary conditions, nonphysical solutions can result. 1. Introduction. Over the past few decades, much experimental and theoretical work has been devoted to the study of polymer-penetrant systems. In particular, the desorption of penetrants from saturated polymer matrices has been examined due to its wide industrial applicability. One unusual feature of such systems is the change in the polymer from a rubbery state when it is nearly saturated to a glassy state when it is nearly dry. As part of the drying process, a glassy skin often develops at the exposed surface of a polymer whose properties are significantly different from the rest of the polymer-penetrant solution There are many different theories for why the skinning process occurs, including phase separation [17], crystallization [18], and diffusion-induced convection [19]. Nevertheless, for the systems we wish to study, most scientists agree that one important factor is a viscoelastic stress in the polymer entanglement network, which can be as important to the transport process as the well-understood Fickian dynamics [20], [21], [22]. The size of this stress is related to the relaxation time of the viscoelastic polymer matrix. In the glassy skin, the relaxation time is finite, so the stress is an important effect, but in the rubbery region the relaxation time is nearly zero [15], [20], [23]. Nevertheless, we will show that in order for the mathematical model to yield physically meaningful results, at some level the short relaxation time in the rubber must also be taken into account.Numerical and analytical solutions are derived here for model equations for the