представлена нелинейная математическая модель состояния свободной поверхности гравитационно стекающей вязкой жидкой пленки в процессах тепло-массообмена. Это нелинейное дифференциальное уравнение в частных производных четвертого порядка содержит как пространственные производные, так и производные по времени. Коэффициенты модели включают параметры поверхностного натяжения, термокапиллярных сил, испарения. Осуществлен переход к разностному уравнению — аналогу исходной математической модели состояния свободной поверхности жидкой пленки. Разработаны вычислительные алгоритмы исследования неустойчивости и состояния свободной поверхности волнового течения жидкой пленки для умеренных чисел Рейнольдса. Проведены вычислительные эксперименты по нелинейному развитию возмущений и выявлению неустойчивых режимов течения жидкой пленки воды, в частности, при неустойчивости Марангони. Найдены области неустойчивости жидкой пленки, в которых выделены режимы с максимальным значением инкремента как при свободном стекании, так и при испарении. Рассчитаны критические значения чисел Марангони, при которых происходит разрушение жидкой пленки. Представлены результаты численного моделирования нелинейного развития возмущений и формирования состояния свободной поверхности испаряющейся жидкой пленки, которые могут быть использованы при конструировании или модернизации существующего пленочного оборудования, а также при разработке технологических процессов в жидких пленках.