Abstract. In this article, we propose a model describing the transport of trapped particles in a surface potential. The potential confines particles close to the surface, increasing the number of surface collisions. First, we consider the case of noncharged particles. From a kinetic description, we rigorously derive a two dimensional Boltzmann equation. In the case of charged particles we introduce the coupling with the Poisson equation. We perform a formal asymptotic analysis which leads to a two dimensional Boltzmann equation coupled with a three dimensional Poisson equation. We illustrate the charged particle model with some numerical simulations of a gas discharge on a satellite solar array. We use a particle in cell (P.I.C.) scheme that is a particle discretization for the Boltzmann equation and a Fourier approximation for the Poisson equation.Key words. Boltzmann equation, Poisson equation, surface collisions, asymptotic model, simulation of gas discharge, particle in cell scheme, particle method AMS subject classifications. 82C70, 82C40, 82C21, 82B40, 82B21, 35J05, 76M28 DOI. 10.1137/050642897 1. Introduction. In this article, we are interested in particles subject to a given external potential in a half space. We propose a mathematical model to describe their transport, and we perform numerical simulations. The starting point of our model is a kinetic model constituted of the Vlasov equation set in the half space. It is completed by some boundary conditions for the description of the surface collisions.In order to reduce the cost of the numerical simulations, it is classical to derive some asymptotic models with a smaller number of variables than the kinetic description. It is possible when the physical conditions allow one to do it, for example when particles are subject to many collisions. The resulting model depends on the considered physical process (see [6] and the references given there).Here we want to establish an asymptotic model when the applied potential confines the particles close to the surface, increasing the number of collisions with it. Furthermore, we assume that the dominant surface collision process is specular; the other collisions are supposed to be only a perturbation. Although it is not true in practical situations we assume that the error is of the same order as the error we make with the asymptotic limit.This work is a continuation of [15], in which the considered problem is the same as here but the dominant surface collisions are supposed diffusive. In [15], the resulting model is a diffusive model in two space dimensions; its variables are the time, the position on the surface, and the total energy of the particles. These kinds of models are called in the literature Fokker-Planck models or "SHE model" (for spherical harmonics expansion) because of its earlier derivation in [27]