Document 854847

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
М атематика. Механика. Информатика
2011
Вып.1(5)
УДК 539.3
Решение трехмерных задач теплопроводности
методом фиктивных канонических областей
С. Л. Гладкий 1 , Л. Н. Ясницкий 2
Пермский государственный технический университет, Россия, 614990, Пермь, Комсомольский пр., 29
Пермский государственный университет, Россия, 614990, Пермь, ул. Букирева, 15
1
spikey@mail.ru; (342) 271-61-68; 2yasn@psu.ru,
1
2
Показана возможность применения метода фиктивных канонических областей (ФКО) для
решения задач стационарной теплопроводности в сложных трехмерных областях. Рассмотрены преимущества метода ФКО. Разработана программа, реализующая метод ФКО. Применение метода продемонстрировано на примере решения объемной задачи с высокой точностью.
Ключевые слова: краевая задача; каноническая область; численно-аналитический метод.
1. Введение
ния только линейных краевых задач математической физики. Решение нелинейных задач
– предмет дальнейших исследований.
Решение краевых задач является одним
из ведущих направлений математического
моделирования. Первые аналитические методы решения краевых задач были предложены
в работах Ж.Л. Д’Аламбера и Ж.Б.Ж. Фурье и
с успехом применялись для решения задач в
простых областях. Однако современные инженерные задачи требуют получения решений краевых задач в трехмерных конструкциях сложной конфигурации. Считается, что
применение аналитических методов для таких
задач невозможно или крайне затруднительно. Поэтому большинство расчетов сложных
элементов инженерных конструкций выполняется численными методами. Тем не менее
применение аналитических методов расчетов
на современных высокопроизводительных
компьютерах может быть эффективным для
задач, где первостепенным фактором является
точность и надежность получаемых решений.
Один из таких методов – метод фиктивных
канонических областей (ФКО). Данный метод
является, по сути, методикой выбора базисных функций к методу Треффтца и позволяет
получать высокоточные решения для областей сложной конфигурации. Следует отметить, что метод ФКО применяется для реше-
2. Метод фиктивных канонических
областей
Пусть требуется решить линейную краевую задачу: найти функцию U x  , удовлетворяющую в пределах некоторого тела
D  R3
дифференциальному уравнению в
частных производных
L U x   Rx , x  D ,
(1)
и на поверхности S тела D – граничным
условиям
(2)
B U x  B* x, x  S ,
где L и B – заданные линейные дифференциальные операторы с постоянными коэффициентами; R x  и B* x  – заданные функции
координат x .
Согласно методу ФКО [1], решение краевой задачи ищется в виде конечной суммы
N
U x   U R x    c n u n x  ,
(3)
n 1
где U R x  – любое частное решение уравнения (1), c n – постоянные коэффициенты,
© С. Л. Гладкий, Л. Н. Ясницкий, 2011
41
С. Л. Гладкий, Л. Н. Ясницкий
u n x  – базисные функции, каждая из которых тождественно удовлетворяет однородному уравнению
L U x   0 .
(4)
Базисные функции выбираются следующим образом. Существуют такие области,
называемые каноническими, для которых известны общие решения однородного уравнения (4) (например, полученные методом разделения переменных Фурье) в виде бесконечного ряда
 

 
U V ς V   c nV u nV ς V , ς V V .
преимущество позволяет существенно сократить время на создание объемной модели и ее
подготовку к решению.
3. Решение трехмерных задач
теплопроводности
В стационарной задаче теплопроводности искомой функцией является температура,
которая при отсутствии тепловых источников
удовлетворяет уравнению Лапласа [4]
2 T  0 ,
(6)
где T – температура; 2 – оператор Лапласа.
На поверхности тела могут быть заданы
условия первого, второго или третьего рода.
Для применения метода ФКО к решению трехмерных стационарных задач теплопроводности необходимы общие решения
уравнения Лапласа для различных объемных
канонических областей – сферы, цилиндра и
т.д. Данные решения могут быть получены
методом разделения переменных Фурье в различных системах координат (СК). В работе
[5] приведены общие решения уравнения
Лапласа в сферической и цилиндрической СК.
В сферических координатах r ,  ,  
оператор Лапласа имеет следующий вид:
(5)
n 0
Такие решения являются общими в том
смысле, что подбором коэффициентов cVn из
них можно выделить частные решения, удовлетворяющие достаточно произвольным краевым условиям на границе канонической области V  R 3 с любой точностью. Под достаточно произвольными понимаются граничные
условия, не имеющие разрывов, изломов и
т.д. Согласно методу ФКО исходное тело
D погружается в пересечение нескольких
канонических областей V1 V2  ... . Решение краевой задачи для исходного тела D
определяется в виде суммы решений, относящихся к каноническим областям. В каждом
таком решении ограничивают число слагаемых. Неизвестные коэффициенты c nV находятся из условия приближенного удовлетворения краевым условиям (например, по методу наименьших квадратов).
Таким образом, метод ФКО является
приближенным аналитическим методом. Решение краевой задачи, полученное методом
ФКО, тождественно удовлетворяет дифференциальному уравнению и приближенно –
краевым условиям.
Главным преимуществом метода ФКО
является возможность простой и надежной
оценки точности полученных решений [1–3].
Это обусловлено тем, что дифференциальное
уравнение задачи удовлетворяется тождественно и для оценки точности достаточно выполнить анализ погрешности удовлетворения
граничным условиям. Также к преимуществам
метода ФКО можно отнести снижение размерности исходной задачи на единицу, поскольку
для решения объемной задачи требуется информация только о поверхности тела. Данное
2
1   2  
1
  2

r

r r  r  r 2 sin 2 ( )  2
2
1

 2
r sin(  ) 
 

 sin(  )
.
 

(7)
Общее решение уравнения Лапласа в сферической СК может быть представлено в следующем виде:
T r ,  ,    c 0  c1 t1 (r ) 

n
4
  c n m i t n m i r ,  ,   
n 1 m 1 i 1

6
  c n 0 i t n 0 i r ,  ,
n 1 i  5
1
t1 ( r )  ,
r
t n m1 r, ,   r n sin( m ) Pnm cos( ) ,
t n m 2 r,  ,    r n cos(m  ) Pnm cos( ) ,
t n m 3 r ,  ,   
42
sin( m  ) P cos( ) 
r
m
n
n 1
,
(8)
Решение трехмерных задач теплопроводности методом фиктивных канонических областей
t n m 4 r ,  ,   
cos(m  ) Pnm cos( ) 
,
r n 1
t n m 6 r , , z   K m  r  sin( m ) cos( z ) ,
tn m 7 r , , z   K m  r  cos(m ) sin(  z ) ,
t n 0 5 r,    r Pn cos( ) ,
n
t n 0 6 r ,   
t n m 8 r , , z   K m  r  cos(m ) cos( z ) ,
Pn cos( ) 
,
r n 1
t n m 9 r , , z   J m  r  sin( m ) sinh(  z ) ,
t n m 10 r , , z   J m  r  sin( m ) cosh( z ) ,
где Pn – полиномы Лежандра, Pnm – присоединенные функции Лежандра [6].
Данное решение является общим для полой сферы, т. е. подбором коэффициентов данного решения можно удовлетворить достаточно произвольным граничным условиям на поверхности полой сферы произвольных внешнего и внутреннего радиуса. Если в данном решении оставить только несингулярные базисные функции, то получим общее решение для
шара, если оставить только сингулярные
функции, получим решение для сферической
полости в бесконечном пространстве.
В
цилиндрических
координатах
r,  , z  оператор Лапласа имеет вид
t n m12 r , , z   J m  r  cos(m ) cosh( z ) ,
t n m 13 r , , z   Ym  r  sin( m ) sinh(  z ) ,
t n m 14 r , , z   Ym  r  sin( m ) cosh( z ) ,
t n m 15 r , , z   Ym  r  cos(m ) sinh(  z ) ,
t n m 16 r , , z   Ym  r  cos(m ) cosh( z ) ,
t n 017 r , z   I 0  r  sin(  z ) ,
t n 019 r , z   K 0  r  sin(  z ) ,
2
t n 0 20 r , z   K 0  r  cos( z ) ,
t n 0 21r , z   J 0  r  sinh(  z ) ,
Общее решение уравнения Лапласа в цилиндрической СК может быть представлено в
следующем виде:
t n 0 22 r , z   J 0  r  cosh( z ) ,
t n 0 23 r , z   Y0  r  sinh(  z ) ,
T r ,  , z   c0  c1 t1 r   c 2 t 2  z 


t n 0 24 r , z   Y0  r  cosh( z ) ,
16
  c n m i t n m i r ,  , z  
t0 m 25 r ,   r m sin( m ) ,
n 1 m 1 i 1

t0 m 26 r ,   r m cos(m ) ,
24
  c n 0 i t n 0 i r , z  
n 1 i 17

t 0 m 27 r ,   r  m sin( m ) ,
28
   c0 m i t 0 m i r ,   
t0 m 28 r ,   r  m cos(m ) ,
m 1 i  25

(9)
t n 018 r , z   I 0  r  cos( z ) ,

1     1 

. (9)
r   2
r r  r  r  2 z 2
2
2 
t n m 11r , , z   J m  r  cos(m ) sinh(  z ) ,
t0 m 29 r , , z   r m sin( m ) z ,
32
   c0 m i t 0 m i r ,  , z ,
t0 m 30 r , , z   r m cos(m ) z ,
m 1 i  29
t1 r   ln( r ) ,
t 2 z   z ,
t0 m 31r , , z   r  m sin( m ) z ,
t0 m 32 r , , z   r  m cos(m ) z ,
t n m1 r,  , z   I m  r  sin( m  ) sin(  z) ,
t n m 2 r,  , z   I m  r  sin( m  ) cos(  z) ,
где  
t n m 3 r , , z   I m  r  cos(m ) sin(  z ) ,
n
L
, L – период изменения темпера-
туры по оси z , I m , K m , J m , Ym – функции
Бесселя [6].
Данное решение является общим для
полого цилиндра. Если в этом решении оставить только несингулярные базисные функ-
t n m 4 r , , z   I m  r  cos(m ) cos( z ) ,
t n m 5 r , , z   K m  r  sin( m ) sin(  z ) ,
43
С. Л. Гладкий, Л. Н. Ясницкий
ции, то получим общее решение для сплошного цилиндра, если оставить только сингулярные – получим решение для цилиндрической полости в бесконечном пространстве.
Решения, содержащие тригонометрические
функции координаты z , соответствуют решению для длинного цилиндра, гиперболические функции – решению для короткого цилиндра.
Решение задачи было выполнено с помощью разработанной универсальной программы REGIONS MF. Согласно методу ФКО
исходное тело было погружено в пересечение
фиктивных канонических областей – полой
сферы и четырех сферических полостей. Схема погружения представлена на рис.2.
Результаты решения задачи представлены на рис. 3 и рис. 4. Максимальная погрешность удовлетворения граничным условиям по тепловому потоку составила менее
1%. Отметим, что в инженерных расчетах такая погрешность может рассматриваться как
следствие идеализации реальных граничных
условий, поскольку все граничные условия
задаются с некоторой допустимой точностью.
То есть данное решение является точным для
некоторых скорректированных граничных
условий [2], незначительно (менее 1%) отличающихся от исходных.
4. Пример решения задачи
Рассмотрим следующую стационарную
задачу теплопроводности. На рис. 1 изображено исходное трехмерное цилиндрическое
тело с отверстием в боковой стенке. На поверхности тела заданы условия третьего рода:
на внешней и торцевой поверхностях температура окружающей среды равна 0, на внутренней поверхности и поверхности отверстия
– 1. Коэффициент теплоотдачи на всей поверхности задан равным 1 (все значения приведены в безразмерном виде).
Рис. 3. Распределение температуры
Рис. 1. Исходное цилиндрическое тело
с отверстием
Рис. 4. Распределение теплового потока
5. Заключение
Рис. 2. Схема погружения исходного тела в пересечение ФКО – полой сферы V1
В работе рассмотрено применение метода фиктивных канонических областей для
решения трехмерных задач теплопроводности. Приведен пример решения задачи для
и сферических полостей V2 ,V3 ,V4 ,V5
(сечение в середине отверстия)
44
Решение трехмерных задач теплопроводности методом фиктивных канонических областей
объемного тела сложной конфигурации, для
которой получено решение с высокой точностью – менее 1% по тепловому потоку. Таким
образом, показано, что метод ФКО может
применяться для решения современных
сложных инженерных задач. Поскольку метод
ФКО позволяет получать надежные аналитические решения с высокой точностью, его областью применения могут быть задачи ответственного назначения, в которых точность и
надежность являются превалирующими факторами.
2. Гладкий С.Л., Ясницкий Л.Н. Об оценке погрешности метода фиктивных канонических
областей // Известия Академии наук. Механика твердого тела. М., 2002. №6. C.69–75.
3. Гладкий С.Л., Степанов Н.А., Ясницкий
Л.Н. Интеллектуальное моделирование
физических проблем / НИЦ "Регулярная и
хаотическая динамика". Москва–Ижевск,
2006. 200 с.
4. Карслоу Г., Егер Д. Теплопроводность
твердых тел. М.: Наука, 1964. 488 с.
5. Полянин А.Д. Справочник по линейным
уравнениям математической физики. М.:
Физматлит, 2001. 576 с.
6. Справочник по специальным функциям с
формулами, графиками и математическими таблицами / под ред. М.Абрамовиц. М.:
Наука, 1979. 832 с.
Список литературы
1. Ясницкий Л.Н. Метод фиктивных канонических областей в механике сплошных
сред. М.: Наука, 1992. 128 с.
Solution of 3D heat conduction problems with
fictitious canonic region method
S. L.Gladkiy1, L. N. Yasnitskiy2
1
Perm state technical university, Russia, 614990, Perm, Komsomolskiy pr., 29
Perm state university, Russia, 614990, Perm, Bukireva st., 15
1
spikey@mail.ru; (342) 271-61-68
2
yasn@psu.ru
2
The fictitious canonic region (FCR) method application for solving steady thermal conduction problems in complicated 3D structures is shown. The method’s advantages are considered. A computer
program for the FCR method is realized. A high-precision solution of 3D problem is presented as an
example of the FCR method’s application.
Key words: boundary problem; canonic region; numerical-analytical method.
45