Numerical simulation of active flow control using plasma actuators is dependent upon the development of models accounting for the effects of the actuators on the flow. These can be obtained either by using experimental results coupled with strong assumptions regarding the plasma force distribution, or through numerical simulation. The objective of this paper is to investigate the characteristics of DBD plasma source terms obtained using direct numerical simulation. Using ONERA's in-house plasma solver COPAIER, we propose an insight of a complete chain of numerical models, ranging from the plasma discharge description to the computation of a DBD-induced ionic wind.