The simulation of chemically reacting flows is a basic tool in natural scie
nces as well as engineering sciences to understand and predict complex flow
phenomena (e.g., concentrations of salt in oceans or crystal growth in sem
iconductor industries, see e.g. [7, 19]). The objective of this paper is tw
o-fold: First, a first-order time-splitting scheme is presented that allows
for efficient parallelization of the related quantities in each time-step.
This scheme is based on the decoupled computation of the new velocity-fiel
d and pressure iterates by means of Chorin's projection method. Second, a t
horough analysis of this scheme is given that leads to optimal error statem
ents which apply to general flow situations.