Комментарии 9
Это всё хорошо, но где сравнение с экспериментом?
Интересно, а как покажет температурное поле тепловизор
Вам следовало бы указать, с одной стороны, на ограниченность аналитической модели, а, с другой, на пользу аналитических моделей для общего понимания протекающих процессов.
У вас plot_surface рисует температуру как высоту над плоскостью (r, z) – но это график функции двух переменных, а не поле в пространстве. А поле-то осесимметричное, значит T(r, z) уже задаёт его целиком: если раскрасить сам разрез через contourf/pcolormesh, а ось Z вернуть под настоящую координату, получится реальная геометрия – где именно в пространстве горячо. Заодно dblquad в двойном цикле по 2500 точкам пересчитывает один и тот же интеграл заново; подынтегральное выражение чётно по phi, так что интеграл от 0 до pi с умножением на 2 плюс векторизация ускорят всё на порядок.
Джентльмены, мой коммент про сравнение с экспериментом был шуточным. Ясное дело, что это учебно-иллюстративная задача по урматфизу (первая часть статьи), а во второй части ценность именно в демонстрации возможностей построения таких картинок в Питоне.
---
Про моделирование конфорки может и можно написать действительно интересную статью, но надо чтоб там был какой-то изюм и неочевидности.
Вот на поле скоростей конвекции было бы интересно посмотреть. Но это гидродинамика, а не скучное уравнение теплопроводности. Кстати, а какая тут модель отвода тепла?

Температурное поле круглой конфорки. Визуализация