Modeling Electronegative Plasma Discharges