Интерференционные методы регистрации полей тепловых потоков и давления на обтекаемых поверхностях тема автореферата и диссертации по механике, 01.02.05 ВАК РФ
Голубев, Максим Павлович
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Новосибирск
МЕСТО ЗАЩИТЫ
|
||||
2009
ГОД ЗАЩИТЫ
|
|
01.02.05
КОД ВАК РФ
|
||
|
ГОЛУБЕВ МАКСИМ ПАВЛОВИЧ
ИНТЕРФЕРЕНЦИОННЫЕ МЕТОДЫ РЕГИСТРАЦИИ ПОЛЕЙ ТЕПЛОВЫХ ПОТОКОВ И ДАВЛЕНИЯ НА ОБТЕКАЕМЫХ ПОВЕРХНОСТЯХ
01.02.05. Механика жидкости газа и плазмы
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук
1 7 Г"'.': ?.Г:СЗ
Новосибирск, 2009
003476650
Работа выполнена в Институте теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения Российской академии наук
Официальные оппоненты:
доктор физико-математических наук, Миронов Сергей Григорьевич
доктор физико-математических наук,
профессор Маркович Дмитрий Маркович
Ведущая организация: Центральный аэрогидродинамический
институт им. проф. Н. Е. Жуковского (ЦАГИ)
Научный руководитель:
кандидат физико-математических наук, Павлов Александр Алексеевич
Защита состоится « 16 » октября 2009 г. в « 14 » часов на заседании диссертационного совета Д 003.035.02 в Институте теоретической и прикладной механики Сибирского отделения Российской академии наук по адресу: 630090, Новосибирск-90, ул. Институтская, 4/1.
Отзыв на автореферат в двух экземплярах, заверенный печатью учреждения, просьба направлять на имя ученого секретаря диссертационного совета Д 003.035.02.
С диссертацией можно ознакомиться в библиотеке Института теоретической и прикладной механики СО РАН.
Автореферат разослан «_ / ьСЯН^Л^^Л 2009г.
Ученый секретарь диссертационного совета
д. т. н.
Засыпкин И. М.
ОБЩИЕ ПОЛОЖЕНИЯ ДИССЕРТАЦИОННОЙ РАБОТЫ
Актуальность работы. Развитие авиационной и космической техники предъявляет все более высокие требования к количеству и качеству информации, получаемой при изучении аэрофизических процессов. Это ведет к необходимости разработки и внедрения новых, более производительных методов регистрации, что обусловлено следующими факторами:
- любые методы имеют ограничения по чувствительности, точности и диапазону измерений, а также по пространственному и временному разрешению;
-расширение линейки доступных методов позволяет одновременно регистрировать несколько параметров и проводить регистрацию одного параметра альтернативными методиками. Это увеличивает достоверность измерений и служит подтверждением адекватности используемых методов.
-важным является снижение стоимости эксперимента, что требует усовершенствования методов измерения в плане повышения качества и количества получаемой за единицу времени информации.
Оптические методы обладают присущими только им достоинствами, что служит дополнительным стимулом их развития. Это возможность бесконтактной панорамной регистрации интересуемого параметра; возможность получения "мгновенной", по сравнению с характерным временем процесса, информации; высокое пространственное разрешение; возможность в некоторых случаях получения информации, недоступной при использовании других методов.
В настоящее время известен ряд панорамных оптических методов регистрации тепловых потоков и давления на обтекаемых поверхностях и различных способов их реализации. Это регистрация температуры с использованием покрытий на основе термочувствительных красителей, индикаторов плавления, жидкокристаллических (ЖК) покрытий, тепловизионная регистрация поверхностной температуры, регистрация давления с использованием барочувстви-тельных люминесцентных покрытий, регистрации поверхностного трения и давления, основанного на применении эластичных покрытий и др.
Наряду с достоинствами, эти методы обладают рядом недостатков. Например, покрытия на основе термокрасителей и индикаторов плавления, меняющих свою окраску при переходе определенного уровня температуры, позволяют визуализировать только одну или несколько изотерм. Ограничением в использовании ЖК-покрытий является узкий температурный рабочий диапазон. При использовании тепловизоров необходимо наличие информации о коэффициенте серости поверхности. Возможно влияние отраженного поверхностью излучения от более нагретых участков самой модели и конструктивных узлов установки, что особенно проявляется при проведении экспериментов на высокоэнтальпий-ных установках. Регистрация давления с помощью люминесцентных покрытий требует наличия кислорода в газе. Метод фактически регистрирует не давление, а концентрацию кислорода. Чувствительность ограничена 0.5-1 %/кПа. Малая скорость выравнивания концентрации кислорода в матрице ограничивает применение методики при исследовании быстропротекающих процессов. При ис-
пользовании эластичных покрытий на результаты могут оказывать влияние изгибы, перекосы и сдвиги модели. Время обработки занимает 20-30 с на кадр, что не позволяет работать в режиме реального времени.
Эти и иные факторы, ограничивающие точность, динамический диапазон, быстродействие и другие важные метрологические характеристики известных оптических методов, подтверждают актуальность работ по их развитию.
Целью работы являлось теоретическое обоснование, разработка и экспериментальное подтверждение работоспособности новых методов регистрации параметров газовых потоков на обтекаемых поверхностях. Для этого решались следующие задачи.
- Разработка общего подхода регистрации параметров потока на обтекаемых поверхностях с использованием методов оптической интерферометрии.
- Разработка и реализация на аэродинамических установках оптического метода регистрации тепловых потоков на обтекаемых поверхностях.
- Разработка и реализация на аэродинамических установках оптического метода регистрации полей давления на обтекаемых поверхностях.
- Развитие разрабатываемых методик для регистрации параметров на искривленных произвольно ориентируемых поверхностях.
Научная новизна работы заключается в том, что
1. Предложен общий подход для панорамных интерференционных методов регистрации тепловых потоков и полей давления на обтекаемых поверхностях с использованием в качестве датчика слоя прозрачного вещества, изменяющего свои оптические свойства под воздействием регистрируемого параметра.
2. Отработана методика регистрации тепловых потоков. В приближении одномерности и линейности тепловой задачи получены общие соотношения, позволяющие по интерферометрическим данным восстанавливать поверхностное распределение теплового потока.
3. Отработана методика регистрации полей давления. В линейном приближении получены соотношения для отклика барочувствительного покрытия, позволяющие по интерферометрическим данным восстановить распределение давления на исследуемой поверхности.
4. Показана возможность использования разрабатываемых методик на искривленных, произвольно ориентированных поверхностях с использованием интерференции в диффузном свете (фазовая спекл-интерферометрия).
5. Получены экспериментальные результаты, доказывающие перспективность разрабатываемых методик. В частности, впервые визуализировано поле давлений от дозвуковой струи (Яе < 1000), набегающей на стенку, с изменением давления на уровне нескольких миллиметров водяного столба и поле давлений за ударной волной, распространяющейся в канале ударной трубы.
Практическая ценность результатов
Панорамный метод регистрации тепловых потоков позволяет проводить измерения непосредственно тепловых потоков без решения обратной тепловой задачи и может применяться при различных режимах обтекания. Особый инте-
pec представляет при измерениях на высокоэнтальпийных установках импульсного или кратковременного действия. Метод используется для панорамной регистрации полей тепловых потоков на обтекаемых поверхностях на установках ИТПМ СО РАН при до-, сверх- и гиперзвуковых режимах обтекания. Может применяться не только в аэрофизическом эксперименте. С успехом используется для визуализации излучения субмиллиметрового диапазона лазера на свободных электронах.
Панорамный метод регистрации полей давления расширяет возможности измерения и визуализации данного параметра в широком диапазоне режимов обтекания. Для его реализации не требуется специализированных источников лазерного излучения и наличия кислорода в газе. Регистрируется непосредственно давление. Полученные результаты показывают, что минимальный уровень визуализируемых давлений может составлять менее 0.1 мм вод. ст., частотный диапазон регистрируемых пульсаций - до 200 КГц, пространственное разрешение - не хуже 1 мм. В настоящее время иных датчиков с подобными метрологическими характеристиками не существует даже для локальных измерений. Метод позволил впервые визуализировать поле давлений от дозвуковой струи (Re < 1000), набегающей на стенку под углом, с изменением давления на уровне нескольких миллиметров водяного столба и поле давлений за ударной волной, распространяющейся в канале ударной трубы.
Достоверность результатов подтверждается данными тестовых экспериментов, повторяемостью и хорошим согласием результатов оптической диагностики с данными полученными альтернативными методами.
Основные защищаемые положения
- Общий подход реализации панорамных интерференционных методов регистрации параметров на обтекаемых поверхностях с использованием в качестве датчика слоя прозрачного вещества, изменяющего свои оптические свойства под воздействием регистрируемого параметра.
-Аналитические соотношения, позволяющие по интерферометрическим данным восстанавливать поверхностное распределение теплового потока.
- Аналитические соотношения для отклика барочувствительного покрытия, позволяющие по интерферометрическим данным восстановить распределение давления на исследуемой поверхности.
-Результаты, полученные в тестовых экспериментах и при исследовании потоков в аэродинамических трубах, при дозвуковых, сверхзвуковых и гиперзвуковых режимах обтекания, подтверждающие работоспособность и эффективность рассматриваемых методов.
Апробация основных результатов. Основные результаты диссертации представлены в 14 печатных работах, в том числе в учебном пособии, публикациях в журналах и в материалах и трудах научных конференций. Результаты докладывались на семинарах ИТПМ СО РАН, СибНИА и XI - XIII Международных конференциях по методам аэрофизических исследований (Новосибирск, 2004, 2007, 2008); 12th International Symposium on the Flow Visualization
(Goettingen, Germany, 2006); IX. Международной научно-технической конференции "Оптические методы исследования потоков" (Москва, 2007) и др.
Личный вклад автора в работу заключается в активном участии в разработке оборудования, необходимого для проведения тестовых экспериментов и реализации представленных в диссертации методов на аэродинамических установках, в тестовых экспериментах по отработке рассматриваемых методик, в аэро физических экспериментах с их применением, а также в интерпретации и обобщении полученных экспериментальных данных. Непосредственно автором получены основные аналитические соотношения для отклика тепло- и барочув-ствительных покрытий.
Структура работы. Работа состоит из введения, 3-х глав, выводов и списка литературы. Объем диссертации составляет 127 страниц, в том числе 72 иллюстрации, б таблиц, и список литературы из 73 наименований.
СОДЕРЖАНИЕ РАБОТЫ
Во введении, наряду с общими положениями, дан обзор методов регистрации температуры, тепловых потоков и давления на обтекаемой поверхности. Рассматриваются особенности их использования в аэрофизическом эксперименте. Приводятся характеристики используемого в работе оборудования.
В главе 1 рассматривается оптический метод регистрации тепловых потоков. Оптическая схема приведена на рис. 1. На плоскопараллельный слой прозрачного вещества толщиной L подается коллимированный пучок когерентного излучения А о под углом у. Часть излучения отражается от грани Pi-волна Аи а часть, преломившись и пройдя через пластину, отражается от задней грани Р2 - волна А2. Волны A¡ и А2 интерферируют между собой.
В результате действия теплового потока Q происходит изменение температуры по толщине слоя, что приводит к изменению разницы оптических путей AS между волнами А\ и А2 и вида интерференционной картины. По сдвигу полос на интерферограммах, зарегистрированных через известный промежуток времени, можно определить значение Q. В общем виде распределение интенсивности на интерферограмме описывается соотношением.
1(х, у) = /, (х о/) + 12(х, у) + 2^0 ,y)I2(x,y) cos[ ДФ(х, у)]- (1)
Зависимость Acp[j2(x, у)] получена при следующих допущениях.
1. Изменения показателя преломления п не приводят к изменению направления распространения излучения - параксиальное приближение.
2. Направление излучения ортогонально поверхности слоя, т. е. у « 0.
Рис.
1. Схема, поясняющая работу метода и фрагменты типичных интерферограмм.
3. Учитывались только температурные коэффициенты изменения показателя преломления р = дп /дТ и линейного расширения а = дЫ ЬдТ. Влияние упругости и пьезооптических эффектов не учитывалось.
4. Распространение тепла вдоль поверхности не учитывалось.
5. Все тепло, вошедшее в датчик, не выходит из его объема. Показано, что с учетом принятых допущений
где N - номер интерференционной полосы, X - длина волны зондирующего излучения, Ср - теплоемкость, к = ХрСр/2(Р + па) - постоянный для вещества коэффициент. В соотношении (2) температура не присутствует, что исключает необходимость решения тепловой задачи. Измерение сводится к регистрации минимум двух интерферограмм через заданный промежуток времени.
Непосредственно регистрируемой величиной является сдвиг интерференционных полос ДЛГ за время Ал Измерение временных интервалов практически не вносит ошибки в результаты. Более существенной является ошибка в определении Ш. Для качественных интерферограмм точность определения Ш обычно составляет 8А/У я 0,05. Ошибка, связанная с этой величиной,
В эксперименте 8AN может быть значительно больше при шаге полос, сравнимом с пространственным разрешением телекамеры. В работе приводятся соотношения, учитывающие этот фактор. Более сложными для анализа являются связанные с принятыми допущениями погрешности.
1. Перетекание тепла вдоль поверхности. Эффект сказывается в основном в погрешности определения локальных значений теплового потока и ухудшает пространственное разрешение метода.
2. Эффекты упругости и связанные с ними эффекты фотоупругости. В рамках данной работы эти эффекты не рассматриваются.
3. Конечность толщины датчика. При выводе (1) предполагалось, что тепло, вошедшее через рабочую поверхность, остается в объеме датчика, однако это может не выполняться. Для датчика и подложки, изготовленных из одного материала с коэффициентом теплопроводности К при включении теплового потока Q\ в момент ? = 0 тепловой поток Qi{t) через нерабочую поверхность датчика
Для проверки метода проводились измерения теплового потока от струи, распространяющейся вдоль плоской поверхности (рис. 2, а). Струя газа 1, истекающая из сопла 2, распространяется вдоль поверхности датчика 3 из стекла ЛК5 толщиной 10 мм, рамерами 80x80 мм. Тепловой поток регулировался мощностью подогревателя 4, работающего от источника тока 5. Управление расходом производилось редуктором б. Излучение лазера 7 (Я=0.65 мкм) лин-
Q(x,y) = hdN/dt,
(2)
bQ(x,y) = /c-5AAVAi = 0.05-k/At.
Рис. 2. Схема регистрации тепловых потоков оптическим методом (а) и с помощью ДТП (б).
зой 8 с фокусным расстоянием 500 мм формировалось в плоскопараллельный пучок. Интерферограммы в отраженном свете регистрировались ПЗС камерой 9 и сохранялись на компьютере 10. Угол падения у не превышал 10~2.
При тех же параметрах проводились измерения (см. рис. 2,6) с помощью паспортизованного термопарного датчика (ДТП) I, установленного на массивной плоской пластине 2, из алюминиевого сплава, выполняющей роль термостата. Покрытие из оргстекла 3 выполняло роль поверхности. Показания регистрировались вольтметром 5. Размеры ДТП 10x10 мм, толщина 2 мм, чувствительностью 288 Дж/(мВ-м ). В табл. 1 приведены оптические свойства и коэффициент /с, рассчитанный по справочным данным для ЛК5, и измеренный экспериментально для органического стекла (ОС) и стекла фирмы Вакеге, также использовавшихся в экспериментах.
Таблица 1. Оптические характеристики некоторых стекол и кристаллов
Вещество Пс Рс'Ю6, К 1 7 -I а 10 , К р, г/см3 Ср, Дж/(гК) 2 к, Дж/см
ЛК5 1,476 6,2 33,9 2,270 0,674 4,440
ОС 1,492 - 350 1,18 1,48 0,442
ВаЬегэ - - - - - 5,155
Измерения проводились следующим образом. Включались струя воздуха и подогреватель. После выхода подогревателя на стационарный режим осуществлялся быстрый ввод датчика в поток, и производилась регистрация необходимых данных. Без изменения параметров истечения струи проводились последовательно измерения термопарным, оптическим и снова термопарным датчиками. Так как ДТП регистрирует средний по площади тепловой поток, результаты оптических измерений усреднялись по соответствующей площади (рис. 3). Различие результатов измерений оптическим методом и ДТП, не превышало 5%, что является хорошим показателем для подобных измерений.
Рис. 3. Положение ДТП и область интегрирования при оптических измерениях. 1 - сопло, 2 - область усреднения.
г =0.1 с
/ = 3.12 с
Рис. 4. Обтекание цилиндра, а - датчик ЛК5, V/ = 44 Вт; 6 - датчик ОС, V/ = 24,9 Вт.
На рис. 4 показаны интерферо-граммы, полученные при исследовании обтекания струей цилиндра диаметром 5 мм и высотой 15 мм, установленного вертикально на плоской поверхности. Видно, что датчик из ОС обладает более высокой чувствительностью по сравнению с датчиком из ЛК5. Количественные измерения показали, что чувствительность ОС в 10 раз выше, чем ЛК5. ОС обладает также и более высоким пространственным разрешением. Сдвиг полос в областях с градиентами теплового потока для ОС более резок и четче выражен, что объясняется более низким коэффициентом температуропроводности.
Для оценки пространственного разрешения и теплового потока через обратную поверхность датчика, были проведены эксперименты по регистрации растекания тепловой метки от точечного источника тепла - сфокусированного лазерного излучения мощностью 25 мВт для датчика из ОС. Показано, что пространственное разрешение при временах регистрации Д/ < 1 с составляет порядка 1 мм, а глубина распространения тепловой метки при Д/ ~ 2 с не превышает 1 мм.
На рис. 5 представлены результаты применения метода на гиперзвуковой азотной трубе Т-327 для исследования обтекания десятипроцентного дельта-крыла с углом стреловидности %=70° и затупленными передними кромками при М = 21, То = 1700 К, Ро = 84 бар. Длина модели составляла 100 мм. На верхнюю плоскую поверхность наклеивался датчик из стекла фирмы ВАГ^ЕИ^ толщи-
Z мм
Рис. 5. Регистрация тепловых потоков на наветренной стороне делта-крыла.
ной 1,5 мм. Измерения проводились для углов атаки а = 0°, 5° и 10°. На рисунке приведены интерферограммы и профили тепловых потоков по размаху крыла для сечений 10, 20, 30 и 40 мм от носика модели, усредненные по трем сериям экспериментов. Наблюдается хорошее совпадение результатов для разных серий.
Метод применялся также при исследованиях взаимодействия косого скачка уплотнения с турбулентным слоем при М = 3 в сверхзвуковой аэродинамической трубе Т-325М. Так как регистрируется непосредственно количество тепла вошедшего в датчик, метод может использоваться не только при исследовании газовых потоков. Например, в работе [1] с его помощью визуализировалась распределение интенсивности излучения в пучке лазера на свободных электронах (ЛСЭ), работающего в субмиллиметровом диапазоне длин волн, для которого в настоящее время не существует полевых фотоприемников.
Во 2-й главе рассмотрен метод, основанный на применении упругих веществ в качестве материала датчика для регистрации полей давления. При этом АЛ1 = А5(Ар). Оптическая схема аналогична используемой для регистрации тепловых потоков (см. рис. 1). Изменение А£(Д/з) в зависимости от давления, сводится к нахождению компонент вектора деформаций. Для регистрации давления по сдвигу интерференционных полос требуется решение обратной задачи. Рассматривалась задача по нахождению деформаций упругого слоя толщиной Я. Одна из поверхностей слоя жестко закреплена на подложке (исследуемая модель). Ось Ох направлена ортогонально поверхности датчика в сторону подложки. Для верхней поверхности 2=0, для нижней - г = Н.
Задача решалась при следующих допущениях и граничных условиях:
1. Поверхностные нагрузки задаются в виде
Р(х,у) = Р0сов(кх), F(;t,^>) = F0sin(fc*:), где Р- изменение нормального давления, Г - поверхностное трение. При этом задача сводится к двумерной, с зависимостью всех величин только от лиг.
2. Деформации считаются малыми.
3. Объемные силы (гравитация) не учитываются.
4. Деформации на нижней границе равны 0 (условие прилипания).
5. Граничные условия на верхней поверхности деформированного слоя сносятся на границу недеформированного слоя, т. е. на г = 0.
Решалось уравнения равновесия [2]
(1 - 2ст)Ди + gIe^d (сНУ «) = 0, (3)
где и = («х, ит) - вектор деформаций, а силы входят через граничные условия. Найдено аналитическое решение уравнения (3), для изменения оптического пути в слое толщиной Н для нагрузки задаваемой условием 1.
ДВД = №*со5{кх)\ (4)
АБ*(к,Н) = 2«о(с!+с3) + 2а[(еШ- 1)(С1 с2/к) - (е"№- 1)(</3 -с3+ с4/к)], (5) 2 2
где а = (по +2)(ло - \)16щ, щ - коэффициент преломления невозмущенного слоя. Коэффициенты с,- и зависят только от Н, к, модуля Юнга Е и коэффици-
ента Пуассона ст используемого вещества. Зависимость п от плотности р находилась по формуле Лоренц Лорентца. Выражение (5) для AS задает спектральную чувствительность слоя к нагрузке с пространственной частотой к. При этом проведена оценка силы поверхностного трения, влияющего на изменение
_з
AS. Показано, что для пограничного слоя ASp/ASp < 10 , и трением можно пренебречь. Для изменения разности хода оптического пути от давления получено:
AS(x) = Р(х)К(к, Н), где К(к, Н) = AS*(к, Н)/Р0. Для анализа решения удобнее использовать относительное изменение оптического пути, зависящее от безразмерной пространственной частоты кН\
AsikH) = AS*(к, H)/S0 = AS*(к, Щ!(2щН) = Рй К(к, Я)/(2п0Я) (6)
На рис. 6 приведены графики As(kH)IE для разных значений о.
Из соотношения (6) видно, что As(kH) с точностью до постоянного множителя совпадает с передаточной функцией К{к, Н) и зависит только от пространственной частоты и параметров материала. В результате, для гармонической нагрузки может быть решена обратная задача:
Р(х) = AS(x)/K(k, Н). (7) 0.001 0.01 о.1 1 ю loo Ш Используя преобразование Фу-Рис. 6. Отклик упругого слоя на гармоническую на- рье, для нагрузки Р(х, у) будем грузку для различных значений ст при Е= 1. иметь
Щх, у) = у)](с, ц)}(х, у), (8)
где Т, - операторы прямого и обратного преобразования Фурье, при этом 2 2 -1/2
/с = + п ) . Выражения (7), (8) описывают формирование отклика на заданную нагрузку. Используя выражение (8), можно найти передаточную функцию:
Л) = Т [AS(x, у)](Ъ Tbl? [Р(х, у)№ Л)- (9)
Т. е. слой вещества может быть прокалиброван, и для известной передаточной функции по интерферограммам можно найти распределение давления:
(10)
Р{х, у) = Т'Х {Т [AS(x, уЖ, т!)}(х, .у).
Передаточная функция для известных Е и а может быть найдена по соотношению (6). Возможно ее определение и с помощью калибровок. Как показывает анализ решения, для этого достаточно знать ее значение всего в двух точках, в качестве которых удобно брать Кн{0) и в максимуме при кН~ 2. Отклик слоя на одностороннее сжатие в барокамере позволяет с высокой точностью определить значение Кц(0). Для нахождения максимума необходимо опреде-
Рис. 7. Схема регистрации полей давления. 1 - трубка; 2 - ветвь струи с большим расходом; 3 - ветвь струи с малым расходом; 4 - датчик давления; 5 - дренажное отверстие.
лить отклик слоя на известную нагрузку для пространственных частот вблизи кН~ 2. В качестве известной нагрузки в работе при калибровках использовалось давление от ламинарной цилиндрической струи, ортогонально падающей на плоскую поверхность. Для подобных струй распределение давления может быть найдено по известным эмпирическим формулам [3], а также с помощью дренажных измерений.
В работе приводятся результаты экспериментов по регистрации давления от дозвуковой струи воздуха, набегающей под углом 0 на плоскую поверхность (рис. 7). Отношение длины L трубки, из которой выдувается струя, к ее диаметру Lid > 50, что обеспечивает на выходе параболический профиль скорости. Вводятся безразмерные параметры h* = hld-расстояние от оси выходного сечения трубки до стенки и х* = x/d - координата вдоль поверхности от критической точки (максимума давления). Использовался датчик с # = 4 мм.
На рис. 8, а, б показаны интерферограммы, зарегистрированные в процессе эксперимента (поток слева). На рис. 8, в-д даны изображения с вычитанием базового кадра - визуализация в муаровых полосах. Современное оборудование позволяет проводить процедуру вычитания в реальном режиме времени и получать наглядную информацию о структуре течения в темпе эксперимента. В зоне контакта струи с поверхностью наблюдается максимум давления, далее по потоку - область возмущений. При увеличении давления растут поперечные размеры струи, возрастает давление в точке контакта и амплитуда возмущений. Подобные возмущения иначе называют псевдозвуковыми волнами. Они генерируются пульсирующей точкой контакта, и связанное с ними колебание давлений знакопеременно. При увеличении давления струя полностью переходит в нестационарный режим.
При угле 0 = 90° реализуется осесимметричная картина стационарного режима. Зависимость статического давления на поверхности от расстояния до оси г* = rid можно описать эмпирическим законом [3]:
Др(г*)/Лр(0)'
: ехр(-с2г*2)
(И)
-0.1
-vi -0 75
где с = 1.48/?* ' при h* < 6.2 и с = 5Ah* ' при h* > 6.2. Для найденной передаточной функции слоя, по соотношению (11) были просчитаны визуализации поля давления (рис. 9, а). Соответствующие им экспериментальные данные представлены на рис. 9, б. При повышении скорости струи возникают поперечные колебания, что приводит к появлению неустойчивости (рис. 9, в).
В настоящее время не существует панорамных методов, позволяющих в режиме реального времени проводить визуализацию полей давлений и пульса-
Рис. 8. Регистрации полей давления от струи, а - интерферограмма зарегистрированная без струи и 6 - со струей; в-д - визуализация полей давления в муаровых полосах от струи набегающей на поверхность: 0 = 25°, Ь* = 2 (в, г) и 0 = 45°, Ь* = 5.2 (Э); г - графики и АР вдоль оси пристенного течения для струи (изображение г).
ций, характеризующихся малыми величинами порядка 1 мм. вод. ст. С помощью разработанного метода впервые удалось в процессе эксперимента визуализировать динамику потока на поверхности, создаваемого осесимметричными струями с различными конфигурациями падения на стенку.
Рассматриваемые датчики могут быть чувствительны одновременно и к давлению, и к тепловым потокам. На рис. 10. представлен пример визуализации теплового потока от струи воздуха, зарегистрированный датчиком давления. Спустя некоторое время струя была резко выключена, что позволило зарегистрировать интерферограммы, со сдвигом полос, связанным лишь с нагревом датчика. Четко регистрируется нулевая полоса, соответствующая Q= 0. Это говорит о наличии в данной области отрывной зоны. Таким образом, становится возможной регистрация сразу двух параметров; давления и тепловых потоков
Рис. 9. Визуализация поля давления датчиком толщиной 4 мм от нормально падающей струи при постепенном увеличении скорости на срезе (ф = 90°, h* = 4): а - теоретический моделирование; б - экспериментальные данные; в - экспериментальные данные для нестационарного режима и положение среза трубки, на одном датчике. В диссертации даны рекомендации по использованию одного покрытия для одновременной регистрации полей давления и тепловых потоков.
Глава 3 посвящена реализации рассматриваемых методов с использованием диффузного света. Необходимость этого связана с тем, что в некоторых случаях реализация схем, работающих в зеркально отраженном свете, сопряжена с рядом трудностей. Рассматривается две группы факторов, ограничивающих их использо-Рис. 10. Визуализация теплового потока, вание. К первой относятся ограничения,
связанные с конструктивными особенностями аэрофизической установки и геометрическими параметрами исследуемой модели, а также с большой возможной кривизной исследуемой поверхности, когда невозможно зарегистрировать ее изображение в зеркально отраженном свете. В результате применение зеркальных схем ограничено плоскими или слабо искривленными поверхностями с ориентацией, обеспечивающей возможность регистрации отраженного излучения. Ко второй группе относятся случаи, когда невозможно создать оптически совершенное покрытие с поверхностями интерференционного качества. Для регистрации интерферограмм в регулярных или, по крайней мере, в криволинейных, но четко регистрируемых полосах, необходимо выполнение условия 8ASo = 2щ-ЪН < 1/4. Здесь 5H и 5ДS - локальные вариации толщины датчика H и оптического пути, связанные с шероховатостью поверхностей. Использование диффузного зондирующего излучения в некото-
рых случаях позволяет обойти эти проблемы. При этом изображения модулированы спекл-структурой.
Нахождение среднестатистических геометрических характеристик спеклов является довольно сложной задачей. Однако для понимания основных принципов рассматриваемых в работе методов достаточно приближенных соотношений для их характерных размеров:
5 ~ £> = Ш, € ~ 52 / X = Ц?! Я)2 = Х/02 = 5/9. где 5 и € - среднестатистический поперечный и продольный размеры спеклов соответственно, .Г - фокусное расстояние приемного объектива, £> - диаметр зрачка объектива. Отношение 0 = Б* = Д/Р- угловая апертура или относительное отверстие объектива £)* является в то же время угловым диапазоном регистрируемых пространственных частот.
Оптические методы с использованием регистрации изображений модулированных спекл структурой, претерпевающей изменения под действием исследуемого параметра, в настоящее время находят широкое применение в физическом эксперименте. В работе кратко рассматриваются основные принципы оптических методов с использованием спекл технологий. Все эти методы, с некоторым приближением, можно разделить на три основные группы.
Спекл-фотография. Позволяет регистрировать поперечное относительно направления регистрации смещение объекта. Широко используется при регистрации деформаций, полей скорости (Р1У-методы), визуализации газовых потоков растровыми (фоновыми) теневыми методами и др. Возможна регистрация изображений в когерентном и в белом свете. В качестве аналога спеклов при этом выступают изображения присутствующих в потоке частиц, неоднородности структуры светорассеивающей поверхности, элементы растра используемого для пространственной модуляции зондирующего пучка и т. п.
Корреляционная спекл-фотография. Предназначена для регистрации смещений, ортогональных исследуемой поверхности. В данной ситуации, с увеличением смещения происходят все более сильные изменения микроструктуры спеклограммы, что позволяет, используя алгоритмы на основе корреляционного анализа получать информацию о сдвигах.
Фазовая (корреляционная) спекл-интерферометрия.
В рассмотренных выше подходах, понятие оптической интерференции не используется, поэтому в современной литературе эти методы чаще называют методами спекл-фотографии. Однако известны методы, в которых описание изменения спекл-структуры невозможно без использования понятия интерференции. К ним относятся и методы, рассматриваемые в третьей главе диссертации.
Интерференционное соотношение (1) верно и для интерференции волн с хаотической пространственной модуляцией амплитуды и фазы (диффузные волны). При интерференции диффузных волн результирующая волна модулирована спеклами, размеры которых определяются спектром пространственных частот суммарной волны. Пусть в плоскости (х, у) интерферируют две волны с
комплексными амплитудами Л\2(х, у) = а\2(х, >')ехр[Д(/112(х, _у)]. Интенсивность имеет вид
/=/1+/2 + 2(/1/2)Шсо8(Ду), где /у = Д\|/= \|/2 - Если в ходе эксперимента пространственное распределение амплитуд а] и а2 остается постоянным, а разность фаз волнами А\ и Л2 меняется на Дер = Дф[■/(/)], где / параметр характеризующий исследуемый процесс, то интерферограммы (спеклограммы) зарегистрированные с некоторым интервалом времени Дг = ¿2 - можно представить в виде:
1Л = /1 + /2 + 2(/1/2)1/2со5(Дч/), /,2 = Л + /2 + 2(/,/2)Шсо5(Ду + Дф). Из-за хаотичного распределения начальных фаз и амплитуд, невозможно визуально связать Дф с изменением вида спеклограмм. Однако при соответствующей обработке, можно получить информацию о Дф, а значит и об изменении параметра Возможно применение различных алгоритмов обработки изображений. В основе одного из них лежит следующая процедура:
I* = <(1а ~ I,\)2> ~ 1 - соэДф. (12)
В соотношении (12) угловые скобки означают усреднение по массиву спеклов в физически малой области (споте). При этом предполагается, что интенсивности интерферирующих волн 1\ и 12 не претерпевают существенного изменения за время между регистрацией спеклограмм. Это условие сводится к малости изменения разности хода между волнами А\ и Аг по сравнению с продольными размерами спеклов, т. е. ДфЯ/(2я)«Необходимость пространственного усреднения приводит к уменьшению пространственного разрешения и чувствительности метода. Как было показано в работе [5], для регулярных ин-терферограмм потенциальная точность определения сдвига полос может составлять ДА^ = Дф/2я ~ 1 /лК, где К - количество градаций серости используемой телекамеры. При работе со спеклами точность определения сдвига интерференционных полос будет составлять Ш- Дф/2л = 2/К.
Схема регистрации интерферограмм в регулярных полосах с использованием диффузного освещения показана на рис. 11. Датчик 1, нанесенный на поверхность модели 2, освещается излучением от лазера 3 прошедшим через
диффузор 4. Изображение поверхности, модулированное интерференционными полосами, объективом 5 фокусируется в приемной плоскости 6. В фокальной плоскости объектива расположена диафрагма 7 с радиусом отверстия Я. Диффузную световую волну, отраженную от модели, можно представить набором плоских волн распространяющихся под различными углами. Для каждой отдельной волны будет наблюдаться интерференционная картина, вид которой описывается выражением, аналогичным соотношениям для зеркальной схемы. Результирующее
Рис. 11. Схема регистрации интерферограмм при диффузном освещении.
изображение будет являться суммой интерферограмм от отдельных плоских
волн с волновыми векторами, лежащими внутри некоторого телесного угла 2
О = др , определяемого параметрами оптической схемы. Для тонких датчиков и малых изменений угла р эти картины практически совпадают, что позволяет наблюдать полосы равной толщины. В работе приведены материалы, позволяющие определить возможность регистрации интерферограмм в регулярных полосах с использованием диффузного освещения, в зависимости от геометрических параметров используемой схемы.
Иногда не удается нанести чувствительное покрытие с качеством поверхностей, позволяющим регистрировать интерферограммы в регулярных полосах. В этих случаях целесообразно изготавливать одну или обе поверхности датчика диффузно рассеивающими. Интерференционное соотношение (1) остается справедливым и в этом случае и позволяет проводить измерения. На рис. 12 представлены результаты визуализации давления от дозвуковой струи воздуха, распространяющейся вдоль поверхности датчика. Одна из поверхностей датчика (внутренняя) была рассеивающей, а другая отражающей. Давление в максимуме порядка 5 см вод. ст. Для изготовления датчика использовались вещества с большим (рис. 12, а) и малым (рис. 12, б) временами отклика. Алгоритм обработки спеклограмм основан на вычитании по модулю из рабочих изображений базового кадра, зафиксированного до включения исследуемого процесса.
Рис. 12. Поле давления от струи, распространяющейся вдоль поверхности. Давление в максимуме АР ~5 см. вод. ст.
Для подтверждения возможности использования упругих покрытий при исследовании быстропротекающих процессов проведена визуализация ударной волны в ударной трубе. Перепад давления Д.Р=1.6атм, скорость фронта 450 м/с. На окно ударной трубы наносился датчик с диффузной поверхностью с К(0) = 15 см. вод. ст. Результаты визуализации показаны на рис. 13. Фронт волны четко локализуется. Здесь же представлены результаты обтекания цилиндра, установленного поперек канала. Отчетливо наблюдается изгиб фронта волны. Пространственное разрешение схемы не позволило определить величину скачка давления. По оценкам, сдвиг полос на скачке составляет Д7У ~ 100. При тех же параметрах проводилась визуализация теневым методом. С точностью до повторяемости параметров процесса положения волны, регистрируемые обоими методами, совпадают.
и»???" -ч- •
даухз; тшШШШШшШ
В диссертации приводятся результаты регистрации радиационных потоков с использованием диффузного освещения, на примере визуализации изображения в инфракрасном (ИК) излучении изогнутого металлического стержня с приближенной температурой 250 °С. Изображение фокусировалось ИК-объективом на датчик тепловых потоков, нанесенный на подложку, изготовленную из кремния.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ
ш 111к /
^■^Ш&шЯШШШш
Рис. 13. Поле давлений настенке ударной трубы.
1. Предложен общий подход применения методов оптической интерферометрии для панорамной регистрации полей давления и тепловых потоков на обтекаемых поверхностях. В качестве датчиков используются слои прозрачного вещества, изменяющие свои оптические свойства под воздействием регистрируемого параметра.
2. Отработана методика регистрации тепловых потоков. В приближении одномерности и линейности тепловой задачи получены соотношения, позволяющие по
интерферометрическим данным восстанавливать распределение непосредственно тепловых потоков на обтекаемых поверхностях. Метод успешно используется на установках ИТПМ СО РАН.
3. Отработана методика регистрации полей давления. В приближении линейности задачи упругости получены общие соотношения для отклика баро-чувствительного покрытия, позволяющие по интерферометрическим данным восстановить распределение давление на исследуемой поверхности. Потенциальные возможности метода позволяют регистрировать поля изменения давления на уровне Ар < 0.1 мм вод. ст., с быстродействием порядка 5 мкс, что недоступно с применением известных датчиков и методик.
4. Показана возможность применения разработанных методик на искривленных, произвольно ориентированных поверхностях с использованием интерференции в диффузном свете (фазовая спекл интерферометрия).
5. Получены экспериментальные результаты доказывающие перспективность разрабатываемых методик. В частности, впервые визуализировано поле давлений от дозвуковой струи (Яе < 1000) набегающей на стенку под углом с изменением давления на уровне нескольких миллиметров водяного столба и поле давлений за ударной волной, распространяющейся в канале ударной трубы.
Цитируемая литература
1. BolotinV.P., Cherkasskiy V.S. et all. Status of the Novosibirsk Free Electron Laser and First Experiments with High Power Terahertz Radiation // First International Conference on Submillimeter Science and Technology: Ahmedabad 2004. Physical Research Laboratory.
2. Ландау Л.Д., Лифшиц E.M. Теоретическая физика. В 10-ти т. Т. VII. Теория упругости: Учеб. пособие. -М.: Наука. Гл. ред. физ-мат. лит., 1987. - 248 с.
3. Горшков А.Г., Старовойтов Э.И., Тарлаковский Д.В. Теория упругости и пластичности: Учеб. Для вузов. -М.: Физматлит, 2002. -414 с.
4. Юдаев Б.Н., Михайлов М.С., Савин В.К. Теплообмен при взаимодействии струй с преградами. -М.: Машиностроение, 1977. -248 с.
5. Pavlov A., Pavlov Al. A., Golubev М. P. Development of interferemetry methods // 12th International Symposium on Flow Visualization (ISFV-12): Proc. Goettingen, 2006. CD Rom Proceedings, ISBN 0-9533991-8-4, ISFV-12-3.3. Paper 159:1-10 p.
Основные результаты работы представлены в публикациях:
1. Маслов Н. А., Бойко В. М., Голубев М. П., Павлов А. А., Павлов Ал. А., Поздняков Г. А. Оптические приборы: Учебное пособие - Новосиб. гос. ун-т, Новосибирск, 2008. 166 с.
2. Голубев М. П., Павлов А. А., Павлов Ал. А., Шиплкж А. Н. Оптический метод регистрации тепловых потоков // Прикладная механика и техническая физика. 2003. Т. 44, № 4. стр. 174-184.
3. Golubev М. P., Pavlov A. A., Pavlov Al. A. Use of phase spekle interferometry for pressure field registration// Int. Conference on the Methods of Aerophysical Research: Proc. Part 1. Novosibirsk, 2008. p. 57-58.
4. Golubev M. P., Pavlov A. A., Pavlov Al. A. Panoramic optical method of heat-flow registration // 12th International Symposium on Flow Visualization (ISFV-12): Proc. Goettingen, 2006. CD Rom Proceedings, ISBN 0-9533991-8-4, ISFV-12-3.3. Paper 157:1-10 p.
Ответственный за выпуск М.П. Голубев
Подписано в печать 27.08.2009 Формат бумаги 60x 84/16, Усл. печ. л. 1.0, Уч.-изд. л. 1.0, Тираж 100 экз., Заказ № 13
Отпечатано на ризографе ЗАО "ДОКСЕРВИС 630090, Новосибирск, Институтская, 4/1
ВВЕДЕНИЕ
1. Общие положения диссертационной работы
2. Краткий обзор оптических методов
3. Некоторые особенности использования оптических методов в аэрофизическом эксперименте
4. Используемое в работе оборудование
ГЛАВА 1. Интерференционный метод регистрации тепловых потоков
1.1. Интерференция двух монохроматических волн
1.2. Физические основы метода
1.3. Чувствительность и точность метода
1.4. Оптические схемы измерений
1.5. Результаты тестовых измерений
1.6. Пространственное и временное разрешение метода
1.7. Примеры использования метода
1.8. Основные результаты главы
ГЛАВА 2. Панорамный метод регистрации полей давления
2.1. Физические основы метода
2.2. Постановка и решение задачи теории упругости 2.4. Анализ решения
2.4. Временные характеристики
2.5. Определение передаточной функции
2.6. Тестовые эксперименты
2.7. Основные результаты главы
ГЛАВА 3. Фазовая спекл-интерферометрия
3.1. Ограничения применения схем, работающих в зеркально отраженном свете
3.2. Спеклы - основные понятия
3.3. Интерференция диффузных волн - спекл интерферометрия
3.4. Схемы регистрации интерферограмм в регулярных полосах
3.5. Схемы регистрации и примеры использования фазовой спекл интерферометрии
3.6. Основные результаты главы ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ СПИСОК ЛИТЕРАТУРЫ v
ВВЕДЕНИЕ 1гОбщие положения^диссертационной работы
Развитие авиационной и космической техники предъявляет все более высокие требования к количеству и качеству информации, получаемой при изучении трехмерного обтекания летательных аппаратов в широком диапазоне летных параметров. Это ведет к необходимости, наряду с совершенствованием традиционных для аэрофизического эксперимента методов регистрации (измерения с помощью пневмодатчиков, термоанемометров, термопар и др.), внедрения новых или мало распространенных, но более производительных, как количественных, так и качественных методов. Это в полной мере относится и к методам регистрация параметров на обтекаемых поверхностях и, в частности, тепловых потоков и полей давления. Необходимость их развития обусловлена следующими факторами:
• Любые методы измерений имеют ограничения по чувствительности, динамическому диапазону, точности измерений, а также по пространственному и временному разрешению.
• Расширение линейки доступных методов позволяет одновременно регистрировать несколько параметров потока и / или проводить регистрацию одного и того же параметра альтернативными способами. Это увеличивает достоверность значений измеряемой величины и служит достаточно надежным подтверждением адекватности используемых методов.
• Важным является снижение стоимости эксперимента, что требует усовершенствования экспериментальной базы, в том числе и методов измерения, в плане повышения качества и количества получаемой за единицу времени информации.
Указанные факторы важны практически для всех методов диагностики. Однако оптические методы обладают присущими только им достоинствами, что служит дополнительным стимулом их развития. К ним относятся: возможность бесконтактных измерений; возможность панорамной регистрации инте3 ресуемого параметра; возможность получения "мгновенной", по сравнению с - характерным~временем~~исследуемого" процесса, информации; высокое пространственное разрешение; в некоторых случаях возможность получения информации недоступной при использовании других методов.
В настоящее время известен ряд оптических методов регистрации тепловых потоков и давления на обтекаемых поверхностях. Наибольшее распространение получили панорамные методы. Для регистрации температуры широко применяются тепловизоры, а также различные термоиндикаторные покрытия на основе термочувствительных красителей, индикаторов плавления, жидких кристаллов (ЖК покрытия). Для регистрации поверхностного трения используются жидкие пленки, эластичные, люминесцентные и ЖК покрытия. Оптическая регистрация полей давления осуществляется с помощью барочувствитель-ных люминесцентных и эластичных покрытий.
Наряду с достоинствами, данные методы обладают определенными недостатками. Например, термоиндикаторные покрытия на основе термочувствительных красителей и индикаторов плавления, резко меняющих свою окраску при переходе определенного уровня температуры, позволяют визуализировать только одну или, в некоторых случаях, несколько изотерм. Существенным ограничением в использовании ЖК-покрытий является узкий температурный рабочий диапазон. При использовании тепловизоров необходимо наличие информации о коэффициенте серости поверхности, что требует дополнительных калибровочных измерений. Возможно влияние отраженного исследуемой поверхностью излучения от более нагретых участков, как самой модели, так и конструктивных узлов установки, что особенно сильно может проявляться при проведении экспериментов на высокоэнтальпийных установках.
Регистрации давления с помощью люминесцентных покрытий требует наличия кислорода в газе. Метод фактически регистрирует его концентрацию, а не давление газа в потоке. Относительно малая скорость выравнивания концентрации кислорода в матрице ограничивает применение методики при исследовании быстропротекающих процессов. Чувствительность, выражающаяся в процентном изменении интенсивности люминесценции при изменении давле--ния-на -1-кПа,—зачастую~ограничена1ша^ёШе1^075-1 %/кПа и не может быть улучшена.
Эти и другие факторы, ограничивающие точность, динамический диапа- ' зон, быстродействие и другие важные метрологические характеристики известных оптических методов регистрации газодинамических параметров на обтекаемых поверхностях подтверждают актуальность работ по их развитию.
Цель и основные задачи работы
Целью работы являлось разработка новых экспериментальных методов регистрации параметров газовых потоков на обтекаемых поверхностях. Теоретическое обоснование и экспериментальное подтверждение работоспособности разрабатываемых методик.
Для достижения поставленной цели решались следующие задачи:
• Разработка общего подхода регистрации параметров потока на обтекаемых поверхностях с использованием методов оптической интерферометрии.
• Разработка и реализация на аэродинамических установках оптического метода регистрации тепловых потоков на обтекаемых поверхностях.
• Разработка и реализация на аэродинамических установках оптического метода регистрации полей давления на обтекаемых поверхностях.
• Развитие разрабатываемых методик для регистрации параметров на искривленных произвольно ориентируемых поверхностях.
Научная новизна
1. Предложен общий подход для панорамных интерференционных методов регистрации тепловых потоков и полей давления на обтекаемых поверхностях с использованием в качестве датчика слоя прозрачного вещества, изменяющего свои оптические свойства под воздействием регистрируемого параметра.
2. Отработана методика регистрации тепловых потоков. В приближении одномерности и линейности тепловой задачи получены общие соотношения, позволяющие по интерферометрическим данным восстанавливать поверхностное распределение теплового потока. Метод получил название Global Interference-Heat-Flux meter (GIHF)
3. Отработана методика регистрации полей давления. В линейном приближении получены аналитические соотношения для изменения оптической толщины слоя упругого вещества на гармоническую (косинусоидальную) нагрузку, позволяющие по интерферометрическим данным восстановить распределение давление на исследуемой поверхности. Метод получил название Global Interference Surface Pressure meter (GISP)
4. Показана возможность использования разрабатываемых методик на искривленных, произвольно ориентированных поверхностях с использованием интерференции в диффузном свете (фазовая спекл интерферометрия).
5. Получены экспериментальные результаты, доказывающие работоспособность и перспективность разрабатываемых методик.
Практическая г{енность результатов
• Панорамный метод регистрации тепловых потоков (GIHF) позволяет проводить измерения непосредственно тепловых потоков без решения обратной тепловой задачи. Метод может использоваться при различных режимах обтекания. Особый интерес представляет при измерениях на высокоэнтальпийных установках импульсного или кратковременного действия. Метод используется для панорамной регистрации полей тепловых потоков на обтекаемых поверхностях на установках ИТПМ СО РАН при до-, сверх- и гиперзвуковых режимах обтекания. Может применяться не только в аэрофизическом эксперименте. С успехом используется для визуализации излучения субмиллиметрового диапазона лазера на свободных электронах.
• Панорамный метод регистрации полей давления на поверхностях (GISP) существенно расширяет возможности измерения и визуализации данного параметра в широком диапазоне режимов обтекания. При его реализации не требуется специализированных источников лазерного излучения и наличия кислорода в газе. Регистрируется непосредственно давление на поверхности модели. Теоретические и экспериментальные оценки чувствительности и быстродействия показывают, что минимальный уровень визуализируемых давлений может составлять-менее~0г1~мм~водгст77в~частотном диапазоне регистрируемых пульсаций до 200 КГц, с пространственным разрешением < 1 мм. В настоящее время иных датчиков с подобными метрологическими характеристиками не существует даже для локальных измерений. Метод позволил впервые визуализировать поле давлений от дозвуковой струи (Re<1000) набегающей на стенку под углом с изменением давления на уровне нескольких мм вод. столба и поле давлений за ударной волной распространяющейся в канале ударной трубы.
Основные защищаемые положения
• Общий подход реализации панорамных интерференционных методов регистрации параметров на обтекаемых поверхностях с использованием в качестве датчика слоя прозрачного вещества, изменяющего свои оптические свойства под воздействием регистрируемого параметра.
• Аналитические соотношения, позволяющие по интерферометрическим данным восстанавливать поверхностное распределение теплового потока.
• Аналитические соотношения для отклика барочувствительного упругого покрытия, позволяющие по интерферометрическим данным восстанавливать поверхностное распределение давления.
• Результаты тестовых экспериментов и данные, полученные при исследовании потоков в аэродинамических трубах, при дозвуковых, сверхзвуковых и гиперзвуковых режимах обтекания, подтверждающие работоспособность и эффективность рассматриваемых методов.
Достоверность результатов
Достоверность результатов подтверждается данными тестовых экспериментов, повторяемостью и хорошим согласием результатов оптической диагностики с данными полученными альтернативными методами.
Личный вклад автора в работу заключается в активном участии в конструкторской разработке оборудования необходимого для проведения тестовых экспериментов и для реализации на аэродинамических установках, представленных в диссертации оптических методов. Непосредственно автором получены основные- аналитические~соотношения для отклика тепло- и барочувстви-тельных покрытий. Автор принимал непосредственное участие в тестовых экспериментах по отработке рассматриваемых методик и в аэрофизических экспериментах с их применением, а также в интерпретации и обобщении полученных экспериментальных данных.
Апробация основных результатов
Основные результаты работы представлены в учебном пособии, публикациях в журналах и в сборниках статей, в материалах и трудах научных конференций [1-12, 51, 52]. Результаты докладывались на семинарах ИТПМ СО РАН, Сиб.НИА, ЦАГИ и XI- XIII Международных конференциях по методам аэрофизических исследований (Новосибирск 2004, 2007, 2008); 12th International Symposium on the Flow Visualization (Goettingen, Germany, 2006); IX. Международной научно-технической конференции "Оптические методы исследования потоков" (Москва ОМИП-2СЮЗ, 2007, 2009).
Структура работы
Объем диссертации составляет 127 страниц, в том числе 72 иллюстрации, 6 таблиц, и список цитируемой литературы из 73 наименований.
Ниже во введении дается классификация оптических методов, и рассматриваются возможности их использования в аэрофизическом эксперименте. Приведен краткий обзор оптических методов регистрации (визуализации) полей тепловых потоков и давления на поверхности модели. Рассматриваются некоторые особенности аэрофизического эксперимента, влияющие на результаты оптической диагностики интерференционными методами и приводятся краткие характеристики используемых в работе источников света и фотоприемников.
В первой главе изложены теоретические основы панорамного оптического метода регистрации тепловых потоков (GIHF). Приведены результаты экспериментов по определению чувствительности, пространственного и временного разрешения метода и даны примеры его использования.
Во второй главе изложены теоретические основы панорамного оптического ^-методарегистрациидавледиОйГпо^ Приведены результаты экспериментов по определению чувствительности, пространственного и временного разрешения метода и даны примеры его использования.
В третей главе рассматриваются возможности реализации GIHF и GISP методов с использованием диффузного освещения, что позволяет регистрировать параметры на криволинейных, произвольно ориентированных поверхностях. Приведены примеры использования данного подхода.
В заключении изложены основные результаты работы.
ОСНОВНЫЕРЕЗУЛЬТАТЫ РАБОТЫ------ —
1. Предложен общий подход применения методов оптической интерферометрии для панорамной регистрации полей давления и тепловых потоков на обтекаемых поверхностях. В качестве датчиков используются слои прозрачного вещества, изменяющие свои оптические свойства под воздействием регистрируемого параметра.
2. Отработана методика регистрации тепловых потоков. В приближении одномерности и линейности тепловой задачи получены соотношения, позволяющие по интерферометрическим данным восстанавливать распределение непосредственно тепловых потоков на обтекаемых поверхностях. Метод успешно используется на аэродинамических установках ИТПМ СО РАН.
3. Отработана методика регистрации полей давления. В приближении линейности задачи упругости получены общие соотношения для отклика баро-чувствительного покрытия, позволяющие по интерферометрическим данным восстановить распределение давление на исследуемой поверхности. Потенциальные возможности метода позволяют регистрировать поля изменения давления на уровне Ар <0.1 мм. вод. ст., с быстродействием < 5 мкс, что недоступно с применением известных датчиков и методик.
4. Показана возможность применения разработанных методик на искривленных, произвольно ориентированных поверхностях с использованием интерференции в диффузном свете (фазовая спекл интерферометрия).
5. Получены экспериментальные результаты доказывающие перспективность разрабатываемых методик. В частности, впервые визуализировано поле давлений от дозвуковой струи (Re<1000) набегающей на стенку под углом с изменением давления на уровне нескольких мм водяного столба и поле давлений за ударной волной, распространяющейся в канале ударной трубы.
1. А. А. Павлов, М. П. Голубев, Ал. А. Павлов. Оптический метод регистрации тепловых потоков // Препринт № 3-2002, ИТПМ СО РАН, Новосибирск, 2002
2. М. П. Голубев, А. А. Павлов, Ал. А. Павлов, А. Н. Шиплюк. Оптический метод регистрации тепловых потоков // Прикладная механика и техническая физика, Т. 44, № 4, стр. 174-184, 2003
3. М. П. Голубев, А. А. Павлов, Ал. А. Павлов. Оптический метод регистрации полей тепловых потоков // Труды VII международной научно-технической конференции ОМИП-2003, (Москва, 24-27 июня 2003 г.) изд. МЭИ, М., 2003, стр. 476-479
4. М. P. Golubev, A. A. Pavlov, Al. A. Pavlov. Optical method of heat-flow-field registration // Proc. at XII Int. Conference on the Methods of Aerophysical Research. Novosibirsk, 28 June 3 July, 2004. Novosibirsk, Russia, Vol. 2, pp. 157-161
5. M. P. Golubev, A. A. Pavlov, Al. A. Pavlov. Panoramic optical method of heat-flow registration // 12th International Symposium on Flow Visualization (ISFV-12), Goettingen, Germany, September 10-14, 2006. Art № 157
6. Н. А. Маслов, В. М. Бойко, М. П. Голубев, А. А. Павлов, Ал. А. Павлов, ГгАт ПозднЖовГОптические прйборыГ"Учебное пособшГ// Новосиб. гос. ун-т, Новосибирск, 2008. 166 с.
7. В. M. Бойко, A. M. Оришич, А. А. Павлов, В. В. Пикалов. Теоретические основы и методы оптической диагностики в аэрофизическом эксперименте. Учебн. пособие // Новосиб. гос. ун-т. Новосибирск, 2008. 412 с.
8. Горлин С.М. и Слезингер И.И. Аэромеханические измерения. Методы и приборы // М.: Наука, 1964
9. Полежаев Ю.В., Юревич Ф.Б. Тепловая защита // Под ред. А.В. Лыкова. М.: Энергия, 1976
10. А. А. Павлов. Развитие панорамных методов оптической диагностики пространственных течений // Автореферат диссертации на соискание ученой степени к.ф.-м.н., Новосибирск, 1995
11. В. Н. Бражко, Н. А. Ковалева, Г. И. Майкапар. О методе измерения теплового потока с помощью термоиндикаторных покрытий // Ученые записки ЦА-ГИ. 1989, Т. XX, №1
12. В. П. Кондакова. Расчетные материалы для определения коэффициентов тёШоотдачй с помощью термоиндикаторов"// В.П. Кондакова, М.В. Рыжкова // Тр. ЦАГИ, 1970, Вып. 1175. с. 240-256
13. В. Я. Боровой. Течение газа и теплообмен в зонах взаимодействия ударных волн с пограничным слоем // М.: Машиностроение, 1983
14. Таблицы физических величин: Справочник // Под ред. И. К. Кикоина. -М.: Атомиздат, 1976
15. Г. М. Жаркова, А. С. Сонин. Жидкокристаллические композиты // Новосибирск: Наука, 1994
16. Г. М. Жаркова, В. М. Корнилов, В. А. Лебига, С. Г. Миронов, А. А. Павлов. Методы и средства исследований течений в аэрогазодинамическом эксперименте // Теплофизика и аэромеханика, т.4, №3, 1997, стр. 283-294
17. Г. М. Жаркова, В. Н. Коврижина, В. И. Корнилов, А. А. Павлов. Метод регистрации полей температур с помощью жидкокристаллических покрытий // Теплофизика и Аэромеханика, т. 3, № 4, стр. 369-279, 1996 г.
18. А. М. Харитонов Техника и методы аэрофизического эксперимента. 4.2. Методы и средства аэрофизических измерений: учебник // Новосибирск: Изд-во НГТУ, 2007
19. А. Н. Петунин Методы и техника измерений параметров газового потока // М.: Машиностроение, 1996
20. Е. У. Репик, В. К. Кузенков Погрешность измерения статического давления при использовании дренажных отверстий // ИФЖ, 1988, № 6
21. А. Г. Гришков, В. А. Шилов. Многоточечные средства для измерения давления в аэродинамическом эксперименте // Тр. ЦИАМ, 1982, №972, с. 59-65.
22. Волоконно-оптические датчики. Вводный курс для инженеров и научных работников. Под ред. Э. Удца // М.: Техносфера, 2008. 520 с.
23. J. I. Peterson, R. V. Fitzgerard. New technique of surface flow visualization based on oxygen quenching of fluorescence // Review of Scientific Instruments, vol. 51, May 1980, p. 670-671
24. S. Tavoularis. Measurement in Fluid Mechanics // England, Cambrige University Press, 2005
25. T. Liu, B. Campbell, S. Burns, J. Sullivan. Temperature and Pressure-Sensitive Paints in Aerodynamics // Applied Mechanics Reviews, Vol. 50, No. 4, pp. 227-246
26. T. Liu, M. Guille, J. Sullivan. Accuracy of Pressure Sensitive Paint // AIAA Journal, Vol. 39, No. 1
27. J. Bell, B. McLachlan. Image Registration for Pressure-Sensitive Paint Applications // Experiments in Fluids, Vol. 22, No. 1, pp. 78-86
28. Jim Crafton, Sergey Fonov, G. Edward, et al. Simultanious measurements of pressure and deformation on a UCAV in the SARL // 11th International Symposium on Flow Visualization, August 9-12, 2004, University of Notre Dame, Notre Dame, Indiana, USA
29. Egami, Fey, Klein, Sitzmann, Wild, Transition detection on high-lift devices in DNW-KKK by means of temperature sensitive paint // 11th International Symposium on Flow Visualization, August 9-12, 2004, University of Notre Dame, Notre Dame, Indiana, USA
30. S. Fonov, G. Jones, J. Crafton, V. Fonov, L. Goss. The development of optical technique for the measurement of pressure and skin friction // Measurement Science and Technology, Vol.17, pp. -1261-1268, 2006
31. JC А. Васильев. Теневые методы. // M:: Наука, 1968.
32. Д. Холдер, Р. Норт. Теневые методы в аэродинамике // М: Мир, 1966
33. М. М. Скотников. Теневые количественные методы в газовой динамике // М.: Наука, 1976
34. М. М. Скотников. Применение теневого прибора к исследованию сверхзвукового пограничного слоя// Оборонгиз, 195446; G.S. Settles. Schlieren and shadowgraph: techniques: -Visualizing phenomena in transparent media//Berlin: Springer,. 2006
35. В. Г1. Клочков, JI. Ф. Козлов, И: В. Потыкевич, М. С. Соскин. Лазерная анемометрия, дистанционная спектроскопия и интерферометрия // Киев: Нау-кова думка, 1985 . 48. М. Борн, Э. Вольф. Основы оптики // М.: Наука, 1970
36. Г. С. Ландсберг. Оптика // М.: Наука, 197650; Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 6. Гидродинамика // М.: Наука, 1986
37. A. A. Pavlov, Al. A. Pavlov, М. P. Golubev. Development of interferometry methods // 12th International Symposium on Flow Visualization (ISFV-12), Goettin-gen, Germany, September 10-14, 2006. Art №159
38. М. Я. Кругер, В. А. Панов, В. В. Кулагин, Г. В. Погарев, Я. М. Кругер,
39. A. М. Левинзон. Справочник конструктора оптико-механических приборов // Под ред. М.Я. Кругера и В.А. Панова. М.; Л.: МАШГИЗ, 1963
40. А. А. Блистанов, В. С. Бондаренко, Н. В. Переломова, Ф. Н. Стрижевская,
41. B. В. Чкалова, М. П. Шаскольская. Акустические кристаллы: Справочник // Под ред. М.П. Шаскольской. М.: Наука, Гл. ред. физ.-мат. лит., 1982
42. В. М. Золотарев, В. Н. Морозов, Е. В. Смирнова. Оптические постоянные природных и технических сред: Справочник // Л.: Химия, Ленингр. отд., 1984
43. А. Г. Горшков, Э. И. Старовойтов, Д. В. Тарлаковский. Теория упругости и пластичности // Учеб.: Для вузов. М.-.ФИЗМАТЛИТ, 2002. - 414 с
44. А. А. Аскадский, Ю. И. Матвеев. Химическое строение и физические свойства полимеров // М.: ХИМИЯ, 1983. 248 с.
45. Л. Д. Ландау, Е. М. Лифшиц. Теоретическая физика, Т. VII. Теория упругости. Учеб. Пособие // М.: Наука. Гл. ред. физ-мат. лит., 1987, 248 с.
46. В. Л. Бидерман. Теория механических колебаний // М.: Высшая школа, 1980,405 с.
47. Я. Г. Пановко. Введение в теорию механических колебаний // М. : Наука:, 1991 .
48. Б. Н. Юдаев, М. С. Михайлов, В. К. Савин. Теплообмен при взаимодействии струй с преградами // М.: Машиностроение, 1977, 248 с.
49. С. Г. Попов. Измерение воздушных потоков // М.-Л., Гостехиздат 1947
50. Д. И. Блохинцев. Акустика неоднородной движущейся среды, 2-е изд. // М.: Наука, 1981
51. М. Франсон. Оптика спеклов // М: Мир, 1980V
52. Н. А. Фомин. Спекл-интерферометрия газовых потоков // Минск: Наука и техника, 1989
53. X. К. Абен. Интегральная фотоупругость // Таллин: Валгус, 1975. 218 с.
54. И. А. Разумовский. Интерференционно-оптические методы механики деформируемого твердого тела // Учеб. Пособие. М.: Изд-во МГТУ им. Н. Э. Баумана, 2007, 240 с.
55. N. Akino, Т. Kunugi, К. Iehimiya, et al. Improved liquid-crystal thermometry excluding human color sensation // J. Heat Transfer, 1989, Vol. 111, p. 558-565
56. A. Burch, M. J. Tokarski. Production of Multiple Beam Fringes from Photographic Scateres // Optica Acta, 15,101, 1968
57. Ю. И. Островский, M. M. Бутусов, Г. В. Островская. Голографическая интерферометрия // М.: Наука, 1977
58. Ч. Вест. Голографическая интерферометрия // М.: Мир, 1982, стр. 504
59. Р. Кольер, К. Беркхарт, JI. Лин. Оптическая голография // М.: Мир, 1973, стр. 688t! й