Работа посвящена вопросам реализации разработанной авторами математической модели дыхательной системы человека, предназначенной для прогнозирования возникновения патологий органов дыхания, обусловленных негативным воздействием загрязняющих компонентов атмосферного воздуха. Предложенная модель описывает легкие как упруго-деформируемую насыщенную двухфазную пористую среду, испытывающую большие градиенты перемещений. Поскольку аналитическое решение поставленной существенно нелинейной задачи представляется нереализуемым, для решения предлагается прибегать к численным методам с пошаговыми процедурами. Предложен алгоритм решения связанной задачи фильтрации воздуха в упруго-деформируемой пористой среде. Численное решение нелинейной подзадачи деформирования двухфазной среды легких осуществляется методом конечных элементов, подзадачи фильтрации — методом конечных объемов. Для реализации алгоритма разработан комплекс программ (на языке C++) с применением технологий параллельных вычислений. На основе томографических снимков легких, получаемых с помощью интерактивного программного продукта ITK-SNAP, выполняется восстановление трехмерной формы легких. С использованием пакета ANSYS ICEM CFD строится объемная конечно-элементная сетка. Численное моделирование течения воздуха в легких человека производится для персонализированной трехмерной геометрии. Представлены поля давления газовой фазы в легких человека в различные моменты дыхательного цикла. Разработанную модель в дальнейшем планируется рассматривать как инструмент для определения зон риска развития патологий органов дыхания, обусловленных негативным воздействием аэрогенных факторов среды обитания.