In this paper, we consider a quasi-neutral plasma expanding in the vacuum gap separating two electrodes. During the expansion, some particles are emitted from the plasmavacuum interface and form a beam in the vacuum. Starting from the two-fluid full EulerPoisson model, we derive an asymptotic model. This asymptotic model consists of a quasi-neutral model in the plasma region, a Child-Langmuir law in the beam region and connection relations at the plasma-beam interface. For this model, we propose a numerical scheme which accounts for the motion of the plasma-beam interface and is much more efficient than the resolution of the original two-fluid Euler-Poisson problem. We demonstrate the efficiency of the model by means of numerical simulations for two different one-dimensional test cases.