A two-dimensional particle-in-cell/Monte Carlo simulation code is developed and used to examine the electrodynamics of a finite-size dust (particulate) plasma. The model that can be applied to a dusty plasma bounded by surfaces or to a dust cloud embedded in an unbounded plasma accounts for the dynamics of plasma and dust particles, includes the collision absorption of plasma species by the dust with the use of a multistep Monte Carlo algorithm, and solves Poisson's equation for the electric fields. Charging of isolated particulates is compared with estimates based on orbital motion-limited current collection. The charge and potential structure of finite-length dust clouds embedded in a plasma are examined, and the results are compared to previous one-dimensional solutions. The average dust charge is found to decrease with increasing dust density. Low-density dust clouds do not perturb the potential and their average charge approaches that of isolated particulates. In the case of dense dust clouds the potential forms a sheathlike structure while the charge forms a double layer with the negative charge residing on the outer boundary of the cloud. Dust clouds initialized in the sheath of a surface biased to a potential are also examined. It is found that the dusty plasma electrodynamics is affected by the applied surface potential, plasma drift energy, and dust density. Collective effects in high-density clouds alter the sheath, which can expand or collapse. It is found that particulates can be charged negatively or positively depending on the drift energy and the applied surface potential. ¿ American Geophysical Union 1994