La validité de théorie quasi-linéaire (QL) décrivant la saturation de l?instabilité faisceau-plasma a été un sujet très controversé depuis plusieurs décennies dans la communauté des physiciens des plasmas. La question de la validité de la théorie QL est importante car les estimations quasi-linéaire se retrouvent dans toutes les branches de la physique des plasmas (fusion magnétique et inertielle, gyrocinétique, interaction laser-plasma, ...) tant du point de vue des développements analytiques que des simulations numériques. Cette question, est reconsidérée d?une manière nouvelle tant du point de vue analytique que par des simulations numériques qui bénéficient à la fois de la puissance des ordinateurs massivement parallèles et du développement des codes Vlasov (méthode semi- Lagrangienne) réputés pour leur précision et leur faible diffusion. Les simulations numériques du modèle auto-cohérent Vlasov-ondes montrent que la théorie QL reste valide dans le régime fortement chaotique à condition de la considérer d?un point de vue statistique. Cependant on montre aussi qu?il existe un régime non-QL avant la saturation. Par ailleurs, on montre analytiquement l?absence de couplage de modes dans le régime de saturation de l?instabilité où un plateau est présent dans la queue de la fonction de distribution des particules. Cela infirme plusieurs travaux analytiques essayant de contredire la validité de théorie QL dans le régime non linéaire fortement chaotique de l?instabilité faisceau-plasma.