Описан способ построения математических моделей взаимосвязанных процессов в пористых средах, исходя из сложной многокомпонентной системы. Показаны возможности пакета FreeFem ++ при решении соответствующих краевых задач для систем квазилинейных параболических уравнений в областях с подвижными границами методом конечных элементов с использованием распараллеливания вычислений.