Нa основі аналізу та систематизації підходів до розрахунку теплових процесів
при безконтактній лазерній деформації розроблено математичні моделі температурного поля пластини за умови дії на неї двох лазерів у вигляді квазілінійного нестаціонарного тривимірного рівняння теплопровідності, які
орієнтовано на ефективне числове розв'язування. Розроблено алгоритми розрахунку розподілу щільності потужності лазерного променю. Одержав подальшого розвитку метод побудови адаптивних сіток для використання з методами першого порядку точності за часом другого порядку за просторовими координатами, який забезпечує скорочення часу розрахунку. Запропоновано використовувати сплайни Коханека-Бартелса. Здійснено розпаралелювання адаптивого методу розв'язування рівняння теплопровідності з використанням викликів до системного АРІ для скорочення часу розрахунку. Розроблені методи реалізовані практично в програмній системі, яка використовується у виробництві та навчальному процесі.