The thermodynamics of a disordered planar vortex array is studied numerically using a new polynomial algorithm which circumvents slow glassy dynamics. Close to the glass transition, the anomalous vortex displacement is found to agree well with the prediction of the renormalization-group theory. Interesting behaviors such as the universal statistics of magnetic susceptibility variations are observed in both the dense and dilute regimes of this mesoscopic vortex system.The behavior of vortices in dirty type-II superconductors has been a subject of intense studies in the last decade [1]. Aside from the obvious technological significance of vortex pinning, understanding the physics of such interacting many-body systems in the presence of quenched disorder is a central theme of modern condensed matter physics. Similarities between the randomly-pinned vortex system and the more familiar mesoscopic electronic systems [2] are highlighted by a recent experimental study of a planar vortex array threaded through a thin crystal of 2H-NbSe 2 by Bolle et al. [3]. Interesting behaviors, including the sampledependent magnetic responses known as "finger prints", have been observed for such a mesoscopic vortex system. The disordered planar vortex array is well studied theoretically [4][5][6][7][8][9]. It is one of the few disorder-dominated systems for which quantitative predictions can be made, including a finite-temperature "vortex glass" phase [5] characterized by anomalous vortex displacements [4], and universal variation of magnetic susceptibility [9]. However, until the work of Bolle et al., there were hardly any experimental studies of this system, with difficulties stemming partly from the weak magnetic signals in such 2d systems. Also, numerical simulations have been limited by the slow glassy dynamics [10], although the availability of special optimization algorithms did lead to the elucidation of the zero-temperature problem in recent years [11]. In this letter, we describe numerical studies of the thermodynamics of the vortex glass via a mapping to a discrete dimer model with quenched disorder. A new polynomial algorithm for the dimer problem circumvents the glassy dynamics and enables us to study large systems at finite temperatures. Our results obtained in the dilute (single-flux-pinning) regime compare well with the experiment by Bolle et al. [3], while those obtained in the collective-pinning regime strongly support the renormalization-group theory of the vortex glass, including its prediction of universal susceptibility variation [9].The Model: The dimer model consists of all complete dimer coverings {D} on a square lattice L as illustrated in Fig. 1(a). The partition function iswhere the sum in the exponential is over all dimers of a given covering, and T d is the dimer temperature. Quenched disorder is introduced via random bond energies ǫ ij , chosen independently and uniformly in the interval (− 1 2 , 1 2 ). The dimer model is related to the planar vortex-line array via the well-known mapping to the solid-on-...