Alternance ortimization method in the coefficient inverse heat conduction problem

Abstract


This paper deals with the extremely formulated inverse heat conduction problem (IHCP); the goal is to define the variable coefficient of thermal conductivity. On the basis of error of uniform approximation of resulting temperature distribution to required distribution, the problem is formulated as parametrical optimization problem which corresponds to nominally correct formulation of problem. The value to be defined is considered as sought optimal control action. The solution of optimization problem is based on the alternance properties of sought optimal solutions and allows to obtain optimal values of sought parameters vector with stated structure of control action. Test example shows the possibility of practical use of the proposed method to solve the considered type of problems.

Full Text

Решение коэффициентных ОЗТ, позволяющих идентифицировать коэффициенты уравнения теплопроводности, по сегодняшний день остается одной из актуальных задач инженерной теплофизики. Особое практическое значение имеют ОЗТ в нелинейной постановке, когда теплофизические характеристики нагреваемых изделий зависят от температуры. Наиболее общий универсальный подход к решению нелинейных обратных задач [1-3] состоит в применении методов регуляризации решений [4, 5] на основе численного моделирования [6]. Решение ОЗТ в вариационной постановке с использованием условно-регулярных методов [7-9] в ряде случаев позволяет получить явную аналитическую форму решения или свести задачу к системе алгебраических уравнений. В настоящей статье для решения коэффициентной ОЗТ в экстремальной постановке [1, 10], предусматривающей минимизацию функционала, являющегося мерой соответствия рассчитанного и экспериментального температурного состояния, предлагается использовать известные специальные альтернансные свойства искомых оптимальных решений [11-13]. Температурное поле описывается в зависимости от пространственной переменной и времени одномерным нелинейным уравнением нестационарной теплопроводности с зависящими от температуры теплофизическими характеристиками (удельной теплоемкостью , плотностью материала и теплопроводностью при стандартном допущении ) (1) и соответствующими краевыми условиями. Коэффициент температуропроводности в (1) является функцией температуры. ОЗТ заключается в восстановлении зависимости на основе экспериментальных данных , полученных в результате измерения температуры в некоторой фиксированной точке . Сформулируем экстремальную постановку рассматриваемой нелинейной коэффициентной ОЗТ. По заданному температурному распределению требуется определить зависимость , минимизирующую невязку между и точным аналитическим или приближенным численным решением краевой задачи (1), соответствующим искомой функции . Для оценки этой невязки используем ошибку равномерного приближения результирующего температурного поля к требуемому [11-13] на заданном временном интервале . Будем считать, что искомая зависимость заведомо определена в заданном классе функций переменной с точностью до выбора некоторого вектора постоянных параметров, например в виде, диктуемом известными аналитическими решениями нелинейной задачи теплопроводности, и, следовательно, задана известной функцией переменной . В этом случае решение ОЗТ может быть сведено к задаче оптимального управления, в которой вектор параметров принимается за искомое управляющее воздействие. Задача оптимального управления теперь может быть сформулирована следующим образом. Для объекта (1) с соответствующими краевыми условиями необходимо найти вектор параметров зависимости , обеспечивающих на заданном интервале выполнение минимаксного соотношения . (2) Решение полученной задачи математического программирования (2) может быть осуществлено альтернансным методом, основанным на использовании специальных альтернансных свойств разности . В соответствии с данными свойствами [11-13] для разности на интервале достигаются знакочередующиеся максимальные по абсолютной величине значения в точках , число которых на единицу превышает число искомых параметров. На основании указанных свойств составляется замкнутая система соотношений для предельных разностей температур в этих точках относительно всех неизвестных, дополненная условиями существования экстремума функции в точках , не совпадающих с границами интервала: (3) При известной форме кривой невязки на интервале наблюдения система соотношений (3) трансформируется к системе уравнений, решение которой позволяет найти оптимальные значения параметров зависимости , минимизирующие функционал (2). Для подтверждения работоспособности предложенного метода рассматривалась нелинейная задача теплопроводности (1) в полуограниченной области изменения пространственной переменой и во времени , дополненная соответствующими краевыми условиями вида (4) Аналитическое решение нелинейной задачи (1), (4) для коэффициента температуропроводности вида (5) и вида (6) получены в [14] и приведены в [15]. При решении ОЗТ температура поверхности нагреваемого тела полагается известной, а определению подлежит переменный коэффициент , структура которого имеет вид (5). Экспериментальное температурное распределение было получено как решение прямой задачи (1), (4) при зависимости коэффициента температуропроводности в виде (6). Задача восстановления коэффициента сводится к задаче параметрической оптимизации вида (2) относительно двух искомых параметров , для определения значений которых составляется при система вида (3) из трех соотношений для предельных значений в точках относительно неизвестных , с дополнительными условиями существования экстремума функции в точках и . Применение предложенного метода для решения рассматриваемой ОЗТ при численных значениях позволило получить значения искомых параметров при структуре вида (5). Приведенные результаты (рис. 1, 2) показывают удовлетворительную точность идентификации нелинейного коэффициента температуропроводности при максимальной погрешности, достигаемой в начальный момент идентификации, не превышающей 4 %. Рис. 1. Восстановление коэффициента температуропроводности: 1 - истинное значение ; 2 - его аппроксимация t а б Рис. 2. Результаты решения ОЗТ: а - погрешность аппроксимации коэффициента температуропроводности; б - ошибка приближения температурного поля Решение коэффициентных ОЗТ, сформулированных в экстремальной условно-корректной постановке, при оценке невязки между результирующим температурным полем и требуемым на основе ошибки равномерного приближения может быть сведено к решению задачи математического программирования. Для ее решения может быть использован специальный альтернансный метод оптимизации, сводящийся к решению систем уравнений, основанных на альтернансных свойствах оптимальных решений для погрешности приближения температур относительно параметров идентифицируемого воздействия.

About the authors

Anna N Diligenskaya

Samara State Technical University

244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
(Ph.D. (Techn.)), Associate Professor

References

  1. Алифанов О.М. Обратные задачи теплообмена. - М.: Машиностроение, 1988. - 280 с.
  2. Алифанов О.М. Идентификация процессов теплообмена летательных аппаратов. - М.: Машиностроение, 1979. - 216 с.
  3. Мацевитый Ю.М. Обратные задачи теплопроводности. В 2-х т. - Киев: Наукова Думка, 2002. - 408 с.
  4. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. - М.: Наука, 1979. - 288 с.
  5. Морозов В.А. Регулярные методы решения некорректно поставленных задач. - М.: Наука, 1987. - 240 с.
  6. Самарский А.А., Вабищевич П.Н. Численные методы решения обратных задач математической физики. - М.: Изд-во ЛКИ, 2009. - 480 с.
  7. Иванов В.К., Васин В.В., Танана В.П. Теория линейных некорректных задач и ее приложения. - М.: Наука, 1978. - 206 с.
  8. Лаврентьев М.М., Романов В.Г., Шишатский С.П. Некорректные задачи математической физики и анализа. - М.: Наука, 1980. - 288 с.
  9. Лионс Ж.-Л. Некоторые методы решения нелинейных краевых задач. - М.: Мир, 1972. - 587 с.
  10. Алифанов О.М., Артюхин Е.А., Румянцев С.В. Экстремальные методы решения некорректных задач. - М.: Наука, 1988. - 288 с.
  11. Плешивцева Ю.Э., Рапопорт Э.Я. Метод последовательной параметризации управляющих воздействий в краевых задачах оптимального управления системами с распределенными параметрами // Известия РАН. ТиСУ. - 2009. - № 3. - С. 22-33.
  12. Рапопорт Э.Я., Плешивцева Ю.Э. Алгоритмически точный метод параметрической оптимизации в краевых задачах оптимального управления системами с распределенными параметрами // Автометрия. - 2009. - Т. 45. - № 5. - С. 103-112.
  13. Рапопорт Э.Я., Плешивцева Ю.Э. Специальные методы оптимизации в обратных задачах теплопроводности // Известия РАН. Энергетика. - 2002. - № 5. - С. 144-155.
  14. Fujita H. Text. Res. Jr., 22. - p. 757, 1952.
  15. Лыков А.В. Теория теплопроводности. - М.: Высшая школа, 1967. - 600 с.

Statistics

Views

Abstract - 42

PDF (Russian) - 13

Cited-By


Article Metrics

Metrics Loading ...

PlumX

Dimensions

Refbacks

  • There are currently no refbacks.

Copyright (c) 2014 Samara State Technical University

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies