Multicomponent reactive flows arise in many practical applications such as combustion, atmospheric modelling, astrophysics, chemical reactions, mathematical biology etc. The objective of this work is to develop a rigorous mathematical theory based on the principles of continuum mechanics. The global existence of weak solutions will be discussed. Further discussion on the combined quasineutral and inviscid limit of the corresponding Navier-Stokes-Poisson system will be presented.