Представлена численная модель, описывающая динамику плотности частиц плазмы при филаментации фемтосекундного излучения в воздухе. Результаты моделирования хорошо согласуются с экспериментальными данными. Исследованы процессы, приводящие к накачке излучательных уровней N2 и N2 +. Модель предсказывает резкое падение электронной температуры и плотности в течение 1 нс. Для первой положительной системы азота наблюдается превышение населенности верхнего излучательного уровня над населенностью нижнего в течение 550 пс.