Уважаемые форумчане,

У кого есть опыт работы с функцией mesh?

Написан код, посвященный решению задачи теплопроводности: рассматривается нагрев воды лазерным импульсом, решается задача теплопроводности.

Задается начальное распределение температуры среды, граничные условия и изучается распространение тепла с течением времени.

Изначально код работал с функцией generateMesh, задающей сетку с фиксированным размером ячейки, и при этом позволял выполнять расчеты для нескольких лазерных импульсов.

При этом конечное распределение температуры после 1-го лазерного импульса использовалось как начальное распределение для 2-го и т.д. Временной профиль теплового источника в коде задается переменной “A”.

С реализацией в 2021 версии Matlab сетки с переменным размером ячейки появилась возможность сократить время расчетов, увеличив шаг сетки ближе к стенкам области расчета (подальше от лазерного источника). Как только я заменил в коде задание сетки с равномерным шагом на неравномерный, код, к сожалению, перестал работать, в случае, если число импульсов более 1.

Расчет для нескольких импульсов реализован в цикле, в конце которого я запоминаю распределение температуры в узлах сетки в переменную Last_U. При попытке подставить это распределение в качестве начального для второго импульса выскакивает ошибка.

Я, очевидно, не понимаю какой-то тонкости работы с узлами сетки, не понимаю, на что ругается код и как эту ошибку исправить.

"Unable to perform assignment because the size of the left side is
14424-by-1 and the size of the right side is 1-by-15111".

Помогите, пожалуйста, исправить ошибку.

В приложении архив с файлами:

For_forum_uniform.m - рабочий файл с равномерным шагом сетки.

For_forum_nonuniform.m - неработающий файл с неравномерным шагом сетки.

A_profile.mat - файл с переменной "А" временного профиля теплового источника.

Фильтрация

Тип

Сортировка

как задать базисную функцию от acos((x).^(2))-(x).^(3)

Теги

Ответы Подписаться на ответы

Популярные посты

Темы

Популярные теги

Ближайшие события

Ответы