WWW.KONF.X-PDF.RU
БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА - Авторефераты, диссертации, конференции
 

Pages:     | 1 || 3 | 4 |

«Кириловский Станислав Викторович УПРАВЛЕНИЕ ВОЗМУЩЕНИЯМИ ГИПЕРЗВУКОВОГО ВЯЗКОГО УДАРНОГО СЛОЯ С УЧЕТОМ РЕАЛЬНЫХ СВОЙСТВ ГАЗА 01.02.05 – Механика жидкости, газа и плазмы ДИССЕРТАЦИЯ на ...»

-- [ Страница 2 ] --

–  –  –

Рис.1.11. Изолинии мгновенных пульсаций плотности (а,б,в,г) и графики среднеквадратичных пульсаций плотности (д,е,ж,з) в сечении x=0.9: (а,д) при взаимодействии ВУС с внешними акустическими волнами на спектре частот f1=5кГц, f2=10кГц, f3=15кГц, …, f10=50кГц; (б,е) с возмущениями типа вдув-отсос; (в,ж) с акустическими волнами и возмущениями типа вдув-отсос, запущенными в противофазе;

(г,з) запущенными в синфазе; М=21, =0, x1 x2 = 0.155 0.17 На рис.1.11 приведены изолинии мгновенных пульсаций плотности и графики среднеквадратичных распределений пульсаций плотности в сечении x = 0.9 при взаимодействии ВУС с многочастотными внешними акустическими волнами (рис.1.11а), с многочастотными возмущениями типа вдув-отсос (рис.1.11б), с акустическими волнами и возмущениями типа вдув-отсос, запущенными в противофазе (рис.1.11в), и запущенными в синфазе (рис.1.11г).



Видно, что при противофазном запуске наблюдается существенное подавление возмущений на верхней границе ПС. Соответственно при синфазном запуске получено двукратное увеличение пульсаций на границе ПС.

Полученные таким образом данные по фазам i начальных возмущений вдув-отсос, необходимых для подавления возмущений, генерируемых в ВУС внешними акустическими волнами в спектре частот, могут быть представлены в виде линейной зависимости от периода колебаний возмущений с соответствующей частотой:

i Wi 0.0027 1.23, (7) где Wi безразмерная частота возмущений. Фазы начальных возмущений, необходимые для усиления генерируемых в ВУС возмущений i, могут быть получены из формулы (7) путм добавления половины периода, т.е. i i + На рис.1.12а точками представлена зависимость от частоты начальных фаз возмущений вдув-отсос, необходимых для подавления, подобранных индивидуально для каждой частоты, и рассчитанных по формуле (7). Отличие начальных фаз подобранных индивидуально от фаз, рассчитанных по формуле (7) не превышает 6%.

Полученные амплитуды Bi начальных возмущений типа вдув-отсос, необходимые для подавления возмущений, генерируемых в ВУС внешними акустическими волнами в спектре частот, также могут быть представлены в виде зависимости от начальной амплитуды акустических волн Аi с соответствующей частотой:

Bi 922.7 ( Ai ) 2 35.14 Ai 0.0353, (8) где Аi начальной амплитуды акустических волн. На рис.1.12б приведены зависимости начальных амплитуд возмущений вдув-отсос от частоты, полученные для одночастотных возмущений и с использованием формулы (8).

Погрешность полученной формулы (8) не превышает 13%.

–  –  –

Рис.1.12. Зависимость от частоты фаз (а) и амплитуд (б) возмущений типа вдув-отсос, необходимых для подавления возмущений в ВУС: символы – расчтные данные, для одночастотных возмущений, линия – данные, рассчитанные по формулам (6) и (7) Численные эксперименты по подавлению возмущений в ВУС с использованием данных, полученных по формулам (7) и (8), показали тот же результат, что и эксперименты на основе индивидуального подбора фазы и амплитуды для каждой монохроматической частоты, т.е. существенное уменьшение пульсаций плотности и границе ПС.

–  –  –

10-2 10-2 10-3 10-3 10-4 10-4 10-5 10-5

–  –  –

0.01 0.01 0.01 0.001 0.001 0.001 0.0001 0.0001 0.0001

–  –  –

Рис.1.14. Спектры пульсаций плотности в сечении x=0.9 на границе ПС (y=0.11) при взаимодействии ВУС с акустическими возмущениями в спектре частот: М=21, =0, (а) с

–  –  –

На рис.1.14 изображены спектры пульсаций плотности на границе ПС при взаимодействии ВУС с акустическими возмущениями в спектре частот с суммарной амплитудой А=0.004 (а), А=0.04 (б) и А=0.12 (в). Видно, что с увеличением амплитуды наблюдается генерация возмущений в высокочастотном диапазоне спектра (f 50 кГц), что приводит к перекачке энергии из низкочастотной области в высокочастотную.

По спектральным амплитудам возмущений плотности на границе ПС (y=0.11, в сечении x=0.9) подсчитаны коэффициенты преобразования, равные отношению амплитуды возмущений плотности, генерируемых на верхней границе ПС, к амплитуде начальных возмущений hi Ai / Ai0. На рис.1.15 приведены коэффициенты преобразования для трх случаев: 1 при взаимодействии ВУС с акустическими возмущениями на спектре частот с

–  –  –

частоты f в кГц, (1) воздействие на ВУС акустических возмущений на спектре частот c A = Ai = 0.





04, (линия) - на спектре частот с Ai /10, (2) воздействие монохроматических i 1 акустических волной указанной частоты c Ai, (заштрихованная область) – воздействие многочастотных акустических волн со случайной фазой при A=0.04: М=21, =0 Для проверки влияния фазы начальных многочастотных акустических возмущений производились расчты с набором из десяти случайных фаз, для каждой из десяти частот из многочастотного пакета, c суммарной амплитудой A=0.04. На рис.1.15б заштрихованной областью показана область значений коэффициентов преобразования, полученных спектральным анализом генерируемых в ПС возмущений. Видно, что в этой области лежат коэффициенты преобразования для монохроматического воздействия. Это говорит о том, что нарушение фазовой связи между частотами в пакете, взаимодействующем с ВУС, ослабляет перекачку энергии в высокочастотную область и не уменьшает амплитуды генерируемых в ВУС возмущений на частотах внешнего возмущения.

Также производилось моделирование управления возмущениями, генерируемыми в ВУС внешними акустическими волнами на спектре, состоящем из десяти частот, равных простым числам fi = (5, 11, 17, 23, 29, 31, 37, 41, 43, 47) кГц. Такие частоты выбраны с целью нарушить фазовую связь между отдельными частотами в пакете. Амплитуды для начальных акустических возмущений выбирались согласно экспериментальным данным по амплитудно-частотной зависимости для возмущений в свободном потоке аэродинамической трубы Т-325 (рис.1.9). Амплитуды и фазы для возмущений типа вдув-отсос для подавления и усиления возмущений в ВУС выбирались согласно эмпирическим формулам (7) и (8).

а б в г Рис.1.16. Графики среднеквадратичных пульсаций плотности в сечении x=0.9: М=21, =0, x1 x2 = 0.155 0.17, (а) при взаимодействии ВУС с внешними акустическими волнами на спектре частот f1=5кГц, f2=11кГц, f3=17кГц, …, f10=47кГц; (б) с возмущениями типа вдувотсос; (в) с акустическими волнами и возмущениями типа вдув-отсос, запущенными в противофазе; (г) запущенными в синфазе; пунктирная линия граница пограничного слоя На рис.1.16 изображены графики среднеквадратичных пульсаций плотности в сечении при взаимодействии ВУС с внешними x=0.9 акустическими волнами на спектре частот f1=5кГц, f2=11кГц, f3=17кГц, …, f10=47кГц (рис.1.16а), с возмущениями типа вдув-отсос (рис.1.16б), с акустическими волнами и возмущениями типа вдув-отсос, запущенными в противофазе (рис.1.16в) и запущенными в синфазе (рис.1.16г). Также на рисунке пунктирной линией изображена граница пограничного слоя. Видно что, как и в случае с набором частот f1=5кГц, f2=10кГц, f3=15кГц, …, f10=50кГц наблюдается подавление пульсаций плотности в случае противофазного запуска возмущений и усиление пульсаций в случае синфазного запуска. Это говорит о слабом влиянии фазовых связей частот в пакете на эффективность интерференционного метода управления возмущениями. Данное утверждение подкрепляется тем фактом, что нелинейные процессы проходят схожим образом в ВУС (см. главу 1.2) при воздействии как акустических волн, так и возмущений типа вдув-отсос.

1.4. Гиперзвуковое обтекание системы пластина-кососрезный свисток В трубных экспериментах, одним из способов генерации локализованных возмущений в ВУС является использование кососрезного газодинамического свистка. В численном моделировании подобный локализованный источник возмущений моделируется заданием граничного условия для поперечного массового расхода на некотором участке поверхности вблизи передней кромки пластины, как это было показано в предыдущих главах. В работах [60, 80] было показано, что использование таких граничных условий хорошо воспроизводит структуру пульсаций плотности в вязком ударном слое и фазовые скорости волн для случая слабонелинейных возмущений газодинамического свистка.

В отличие от случая слабонелинейных возмущений, моделирование нелинейных пульсаций сталкивается с рядом существенных трудностей. В частности, вблизи передней кромки пластины искажения среднего течения уже нельзя назвать малыми, и они не могут быть корректно воспроизведены введением пульсаций массового расхода на поверхности. Чтобы решить эти проблемы, выполнено параметрическое численное исследование обтекания системы «пластина – кососрезный свисток» путем решения полной нестационарной задачи обтекания этой системы при наличии автоколебаний в полости свистка с помощью пакета ANSYS Fluent 12.1 [81]. Результаты отдельных расчетов сопоставлены с данными экспериментов.

1.4.1. Тестирование пакета ANSYS Fluent Тестирование пакета производилось сравнением ANSYS Fluent расчетных данных, полученных при решении задачи гиперзвукового обтекания пластины, расположенной под нулевым углом атаки, при воздействии внешних акустических волн набегающего потока двумя пакетами (пакет ИТПМ и ANSYS Fluent). Размер расчтной области и параметры сетки брались теми же, что и в предыдущих главах (см. параграф 1.1).

При решении стационарной задачи с помощью пакета ANSYS Fluent на левой и верхней границах расчетной области с помощью граничного условия pressure-far-field (здесь и далее использованы термины пакета Fluent) задавался равномерный гиперзвуковой поток с параметрами M=21, p=0.7682 Па, T=12.68К. На нижней границе, не занятой поверхностью пластины, а также на правой границе ставилось граничное условие pressure-outlet, т.е. решение экстраполировалось изнутри расчетной области. На поверхности пластины задавались граничные условия wall с условием скольжения и скачка температуры. Так как условия скольжения и скачка температуры не предусмотрены в пакете ANSYS Fluent, была написана подпрограмма UserDefined-Function (UDF – модуль, встраиваемый пользователем) [82] для задания ненулевой тангенциальной компоненты скорости на поверхности пластины u и температуры T согласно формулам (3). Величина градиента рассчитывалась как разница значения на поверхности и в центре ячейки над поверхностью, делнная на расстояние от центра ячейки до поверхности.

Также решалась нестационарная задача о взаимодействии ВУС с акустическими возмущениями набегающего потока. Для этого также c помощью специально написанного модуля UDF на верхней и левой границах расчетной области задавалась суперпозиция стационарного течения и плоских монохроматических акустических волн. Здесь приведены формулы, используемые в UDF для задания акустической волны:

–  –  –

волнового вектора, связанного с частотой f дисперсионным соотношением 2f (верхний (нижний) знак соответствует быстрой k RT M cos( ) 1 (медленной) акустической волне). Индекс «» соответствует параметрам набегающего потока, штрих при параметрах обозначает возмущения. После введения акустических возмущений уравнения НавьеСтокса интегрировались до момента выхода нестационарного решения на установившийся периодический режим.

Полученные с помощью пакета ANSYS Fluent данные сопоставлялись с результатами численного моделирования, выполненного с помощью пакета, созданного в ИТПМ [76], на основе решения двумерных уравнений Навье – Стокса с использованием схем сквозного счета высокого порядка точности (см.

главу 1.1).

На рис.1.17 приведено распределение плотности в стационарном течении (а) и среднеквадратичных пульсаций плотности для случая взаимодействия ВУС с медленной акустической волной с частотой f=38.4 кГц (б). Видно, что данные, полученные с использованием пакета ANSYS Fluent, удовлетворительно согласуются с результатами, полученными с помощью пакета ИТПМ и экспериментальными данными, полученными в аэродинамической трубе Т-327А ИТПМ СО РАН [60].

а б Рис.1.17. Распределение плотности в стационарном течении (а) и среднеквадратичных пульсаций плотности при взаимодействии с медленной акустической волной (б) по нормали к пластине в сечении x = 155 мм (а), x = 151 мм (б): М=21, = =0, f=38.4 кГц, А=0.024: 1 – расчеты по программе ИТПМ СО РАН; 2 - расчеты с использованием ANSYS Fluent; 3 экспериментальные данные [60] 1.4.2. Постановка задачи и численное моделирование Далее в рамках пакета ANSYS Fluent 12.1. выполнено численное моделирование двумерной нестационарной задачи с автоколебаниями при гиперзвуковом обтекании системы пластина – кососрезный свисток с учетом геометрических параметров свистка. Расчеты проводились на вычислительных кластерах ИТПМ СО РАН и информационно-вычислительного центра (ИВЦ) НГУ.

Расчетная область представляла собой прямоугольник. На левой и верхней границах расчетной области с помощью граничного условия pressurefar-field задавался равномерный гиперзвуковой поток (М = 21, p = 0.7682 Па, Т=12.68К) направленный вдоль оси x под углом. На нижней и правой границе ставилось граничное условие с параметрами pressure-outlet p = 0.7682 Па, Т=12.68К. В центральной части расчетной области размещалось тело трапециевидной формы, представляющее собой модель пластины (см. рис.1.18). Под пластиной помещалось тело в форме кососрезной каверны, имитирующее газодинамический свисток. На всех поверхностях пластины и свистка задавались граничные условия типа wall с условием прилипания и температурой поверхности Tw = 300K. Длина пластины L в расчетах составляла 100 мм, что соответствовало числу Рейнольдса ReL= 6104;

размеры задней границы полости свистка составляли 12 мм.

Угол наклона свистка в расчетах изменялся по отношению к пластине с ростом угла атаки для получения автоколебаний в свистке (при = 5-10 =10, при = 15-20 =20, при = 25-30 =30). После установления автоколебаний в полости свистка уравнения интегрировались до установления периодического режима колебаний. При дальнейшем расчете по целому числу периодов колебаний рассчитывался Фурье спектр сигналов в контрольных точках для получения спектральных амплитуд пульсаций давления B.

Использовалась прямоугольная расчетная сетка с характерными размерами ячейки 0.10.1 мм и 0.30.3 мм. Величина шага по времени составляла 210-7 с, и на каждом шаге по времени производилась запись в файл величины давления в 4-х контрольных точках: в центре задней границы полости свистка и на поверхности пластины в сечениях x=0.05L, 0.5L, 0.95L.

Для расчетов использовался 2D density-based солвер, неявная схема 2-го порядка точности по времени и неявная схема 2-го порядка точности по пространству с Roe-FDS методом расщепления конвективных потоков.

а б

Рис.1.18. Поля мгновенного статического давления в расчетной области: а) =15, =20;б) =25, =30

Примеры полученных полей мгновенного статического давления в расчетной области и Фурье спектров возмущений давления в свистке представлены на рис.1.18 и рис.1.19, соответственно. Из рис.1.19 видно, что в полости свистка наблюдается ярко выраженный пик основной частоты колебаний.

Основная доля возмущений, генерируемых кососрезным свистком, распространяется ниже передней кромки модели, тем не менее, некоторая доля возмущений проникает за переднюю кромку пластины и далее развивается вниз по потоку.

а б

–  –  –

Рис.1.19. Фурье спектры пульсаций давления в кососрезном свистке: а) =15, =20;

б) =25, =30 На рис.1.20 представлены Фурье спектры пульсаций давления на поверхности пластины в сечении x=0.05L (вблизи передней кромки) и в сечении x=0.95L. Сравнение их со спектрами на рис.1.19 показывает, что при переходе через переднюю кромку пластины меняется не только интенсивность возмущений давления, но и соотношение между амплитудами основной частоты и гармоник. Дальнейший рост гармоник (в большей степени гармоник более высоких порядков) также наблюдается в конце пластины.

–  –  –

Рис.1.20. Фурье спектры пульсаций давления на поверхности пластины в сечении x=0.05L (а) и x=0.95L (б): =25, =30 Эксперименты по взаимодействию возмущений, создаваемых цилиндрическим кососрезным газодинамическим свистком, с ударным слоем на пластине под углом атаки производились в гиперзвуковой азотной аэродинамической трубе Т-327А ИТПМ СО РАН Мироновым С. Г. и Цырюльниковым И. С. [83]. Пульсации давления в полости свистка и на поверхности пластины в сечении x=0.95L измерялись пъезокерамическими датчиками пульсаций давления. На рис.1.21а изображены зависимости амплитуд пульсаций давления на основной частоте и частоте гармоник на поверхности пластины в сечении x=0.95L от угла атаки, полученные в трубных экспериментах.

Далее проведено параметрическое исследование амплитуд и частот автоколебаний в свистке в зависимости от угла атаки пластины и геометрии свистка. На рис.1.21б изображены зависимости амплитуд пульсаций давления на основной частоте и частоте гармоник на поверхности пластины в сечении x=0.95L от угла атаки, полученные в вычислительных экспериментах.

Спектральные амплитуды пульсаций давления здесь нормированы на амплитуду основной частоты.

–  –  –

,, Рис.1.21. Экспериментальные (а) и расчетные (б) амплитуды пульсаций давления на поверхности пластины, нормированные на амплитуду основной частоты, в зависимости от угла атаки в сечении x=0,95L для основной частоты f0 (1) и первых четырех гармоник (2,3,4,5): 1 - f0= 4,5 кГц, 2 - f1=2 f0, 3 - f2=3 f0, 4 - f3=4 f0, 5 - f4=5 f0, М = 21, ReL = 6104 Сравнение расчетных данных с экспериментальными (рис.1.21) показывает хорошее качественное, и даже количественное соответствие.

Наблюдается рост амплитуд гармоник от угла атаки. Для более точного соответствия эксперименту, вероятно, нужно производить численное моделирование трхмерной задачи.

1.5. Выводы по главе 1 Показано численно, что для управления интенсивностью пульсаций, генерируемых акустическими возмущениями в гиперзвуковом (М 12) вязком ударном слое на пластине, расположенной под малыми углами атаки ( 15), может использоваться интерференционный метод. Этот метод позволяет подавлять или усиливать пульсации, генерируемые в ударном слое внешними возмущениями, путем воздействия на них контролируемыми возмущениями с поверхности пластины типа вдув-отсос с определенной амплитудой и фазой.

Показана возможность применения интерференционного метода управления возмущениями, генерируемыми в ВУС внешними акустическими волнами в спектре частот. Для параметров возмущений типа вдув-отсос, подобранных для усиления/подавления возмущений ВУС, предложена математическая зависимость, позволяющая определять необходимую амплитуду и фазу возмущений типа вдув-отсос в зависимости от амплитуды и частоты внешних акустических волн.

Выполнено численное моделирование нестационарной задачи гиперзвукового обтекания системы «пластина – кососрезный свисток» с учетом геометрических параметров свистка при наличии в нем автоколебаний.

Исследовано влияние угла атаки на возникновение пульсаций в ударном слое пластины, их интенсивность и спектральный состав. Получены данные об амплитудной и фазовой связи пульсаций в кососрезном свистке и в ударном слое на пластине при вариации угла атаки пластины. Установлено что с ростом угла атаки наблюдается рост амплитуд гармоник. Результаты численного моделирования сопоставлены с экспериментальными данными, полученными в аэродинамической трубе Т-327А ИТПМ СО РАН.

Результаты приведенных исследований представлены в [80,83-102].

ГЛАВА 2 Метод пористых звукопоглощающих покрытий

Для управления ламинарно-турбулентным переходом при гиперзвуковых скоростях потока в [103] было предложено использование звукопоглощающих покрытий. Суть предложенного метода сводится к поглощению энергии акустических возмущений с помощью пористого покрытия, что приводит к замедлению роста возмущений второй моды и увеличению длины ламинарного участка течения.

В работе [104] впервые получено экспериментальное подтверждение работоспособности предложенного метода, в которой исследовался ламинарнотурбулентный переход в пограничном слое на остром конусе, половина которого была покрыта пористым звукопоглощающим материалом, в потоке азота и углекислого газа с М 5 под нулевым углом атаки.

Загрузка...
Задержка ламинарно-турбулентного перехода вследствие поглощения акустических возмущений пористым покрытием показана и при бльших числах Маха как с (М= 6) помощью линейной теории устойчивости [105] так и в экспериментальных исследованиях [32-34,105-107] (6 М 8), в том числе с разным углом затупления носовой части конуса [108]. Численное моделирование восприимчивости и развития возмущений на плоской пластине со звукопоглощающим покрытием при М = 6 [109] показало, что покрытие эффективно уменьшает скорость роста второй моды.

В предыдущей главе для пластины под небольшими углами атаки (15) показано, что при воздействии акустических волн внешнего потока в ударном слое генерируются и развиваются энтропийно-вихревые возмущения, управление интенсивностью которых возможно путем введения в ударный слой локальных контролируемых возмущений типа вдув-отсос соответствующей амплитуды и частоты (интерференционный метод управления). Однако при увеличении угла атаки пластины управление интерференционным методом оказывается менее эффективным.

Исследования развития возмущений в ВУС на пластине под углом атаки [110] показали, что при увеличении угла атаки доля энтропийно-вихревых возмущений за УВ резко уменьшается, а возмущения давления растут, при этом наиболее интенсивный рост возмущений давления наблюдается на поверхности пластины. Максимум интенсивности пульсаций давления на поверхности пластины наблюдается в гиперзвуковом пограничном слое и при меньших числах Маха [111-113], при которых преобладают акустические возмущения.

Эти выводы согласуются с результатами расчета по линейной теории взаимодействия возмущений с УВ [73,60] которые показали, что при М = 21 и 15 в ударном слое доминируют энтропийно-вихревые возмущения, а акустические волны экспоненциально затухают, но при углах атаки 15 в ударном слое наряду с энтропийно-вихревыми возмущениями возникают и незатухающие акустические возмущения.

Т.о. при 15° в ВУС доминирует энтропийно-вихревая мода возмущений. В этом случае для управления энтропийно-вихревыми возмущениями в ВУС на пластине можно использовать интерференционный метод управления возмущениями, как это было показано в главе 1 настоящей работы. При увеличении угла атаки пластины эффективность интерференционного метода снижается, и для подавления пульсаций давления на поверхности модели можно использовать метод звукопоглощающих покрытий.

2.1. Постановка задачи В настоящей главе рассматривается развитие возмущений в ВУС на пластине L=100 мм со звукопоглощающим покрытием и на сплошной пластине, расположенной под углами атаки 15° в гиперзвуковом потоке азота (M = 21, ReL = 6104), при взаимодействии с внешними акустическими волнами набегающего потока. Моделирование пористых звукопоглощающих покрытий осуществлялось двумя разными способами: 1-ый с помощью граничного условия на поверхности пластины, предложенного С.А. Гапоновым [114] и усовершенствованного А.В. Фдоровым [115] и моделирующего 2-ой взаимодействие возмущений со звукопоглощающим материалом, и заданием конкретной геометрии и структуры звукопоглощающих вставок [116,117] (геометрическая модель). Во втором случае рассматривалось воздействие на ВУС как медленных монохроматических акустических волн, так и широкополосного шума аэродинамической трубы.

–  –  –

стенке (различие между мгновенным значением давления и давлением, соответствующим решению стационарной задачи до введения возмущений);

A th(h) / Z0 коэффициент проницаемости пористого звукопоглощающего слоя; величина пористости (отношение объема пор к общему объему пористого слоя); Z 0 и характеристический импеданс и константа распространения для отдельной поры, безразмерная круговая частота возмущений. Импеданс пористой звукопоглощающей поверхности Z 1/ A для материала, расположенного на жесткой стенке, равен Z 1 Z 0cth (h).

В случае произвольного звукопоглощающего покрытия импеданс является комплексной величиной, значение которого для конкретного типа покрытий можно определить с помощью теоретических или полуэмпирических зависимостей, исследованных, например в [118,119] для материалов с прямолинейными регулярными порами и хаотической волокнистой структурой.

Однако, в качестве упрощенной тестовой модели, а также для реализации условия полного поглощения звука порами, в данном численном моделировании были применены граничные условия для звукопоглощающего покрытия только с действительным значением импеданса Z.

В ряде случаев импеданс такого типа может быть реализован в эксперименте. В частности, значения импеданса Z с малой комплексной частью могут быть получены для материалов с высокой пористостью при значительном увеличении толщины слоя h или с высоким показателем затухания звука в порах в исследуемом диапазоне частот акустических возмущений, тогда cth(h) 1. С учетом направления оси у в данном случае

–  –  –

Это упрощенное граничное условие на части поверхности пластины со звукопоглощающим материалом было использовано в прямом численном моделировании развития возмущений в ударном слое с помощью программы решения уравнений Навье-Стокса, разработанной в ИТПМ СО РАН [76].

–  –  –

Рис.2.1. Расчетная область для прямого численного моделирования с использованием пакета ANSYS Fluent При численном моделировании с помощью геометрической модели пористости расчетная область (рис.2.1) состояла из двух прямоугольных подобластей. Первая подобласть (I) располагалась над поверхностью пластины и включала в себя набегающий поток, ударную волну и вязкий ударный слой.

Вторая подобласть (II) с верхней границей, совпадающей с поверхностью пластины, располагалась под поверхностью пластины и включала в себя пористую звукопоглощающую вставку. Высота первой расчетной подобласти выбиралась из условия, чтобы идущая с передней кромки головная ударная волна не взаимодействовала с верхней границей расчетной области. Правая (выходная) граница отодвинута от задней кромки пластины, чтобы течение в выходном сечении было полностью сверхзвуковым. Использовалась равномерная расчетная сетка с квадратными ячейками, длина стороны которых равнялась 0.1 мм.

В работе рассматривались звукопоглощающие вставки двух видов: в виде набора полузамкнутых двумерных каналов (рис.2.1а) и в виде регулярной ячеистой структуры, состоящей из поверхностей квадратного сечения, расположенных в шахматном порядке (рис.2.1б). Первый вид звукопоглощающей вставки рассмотрен в данной главе, второй вид подробно рассмотрен в главе 3 настоящей работы. Геометрические размеры пор выбирались из соображений сходства со звукопоглощающими вставками, использованными в трубных экспериментах. В данной главе рассматриваются звукопоглощающие вставки в виде набора полузамкнутых двумерных каналов глубиной 25мм, шириной равной 4мм и 1мм, располагающиеся друг от друга на расстоянии 1мм и 0.25мм, соответственно.

Параметры решателя использовались такие же, как и в предыдущей главе, т.е. 2D density-based решатель, неявная схема 2-го порядка точности по времени и неявная схема 2-го порядка точности по пространству с Roe-FDS методом расщепления конвективных потоков.

Сначала производился расчет стационарного течения, при этом с помощью граничного условия pressure-far-field на левой и верхней границах (1) расчетной области (см. рис.2.1) задавался равномерный гиперзвуковой поток с параметрами M=21, p=0.7682Па, T=12.68К, направленный вдоль оси x под углом =30. На нижней границе, не занятой поверхностью пластины, а также на правой границе ставилось граничное условие pressure-outlet (2). На поверхности пластины и поверхностях каналов и пор задавались граничные условия wall с условием прилипания и температурой поверхности Tw = 300K.

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

Для этого с помощью процедуры User-Defined-Function (UDF) был написан модуль, позволяющий на верхней и левой границах расчетной области задавать суперпозицию стационарных граничных условий и периодических возмущений давления, температуры, числа Маха, продольных и нормальных компонент скорости (см. параграф 1.4.1)

–  –  –

Здесь c T R - местная скорость звука, c T R - скорость звука в набегающем потоке.

На рис.2.2б представлен полученный спектр пульсаций плотности в свободном потоке аэродинамической трубы Т-327А. Далее широкополосный шум потока в предположении прямолинейности фронта идущих волн вводился в расчтную область, и выполнялось численное моделирование взаимодействия широкополосного шума аэродинамической трубы с ВУС на пластине под углом атаки.

–  –  –

Рис.2.2. Пульсации давления от времени (а) и спектр пульсаций плотности (б) в свободном потоке аэродинамической трубы Т-327А ИТПМ СО РАН Для задания направления распространения возмущений под углом по отношению к пластине, т.е. в направлении распространения потока ( = ), расчтная область строилась из двух подобластей (рис.2.3). Первая подобласть I представляла собой прямоугольник и была подобна расчтной области, используемой в предыдущих главах. Вторая подобласть II – это развернутая на угол по отношению к поверхности модели часть прямоугольника, в которой задается распространение возмущений под углом =30°. С помощью разработанного модуля UDF на левой верхней границе второй подобласти Рис.2.3. Расчтная область в задаче моделирования взаимодействия ВУС с широкополосным шумом аэродинамической трубы Т-327А ИТПМ СО РАН:

1 входная граница с граничным условием набегающего потока, 2 границы с условием симметрии, 3 поверхность пластины, 4 выходная граница, 5 вязкий ударный слой на пластине (рис.2.3, граница 1) ставилось нестационарное граничное условие набегающего потока на основе данных по пульсациям плотности в аэродинамической трубе Т-327А. На боковых гранях второй подобласти (рис.2.3, граница 2) ставилось условие симметрии. На нижней границе первой подобласти, не занятой поверхностью пластины, а также на правой границе первой подобласти (рис.2.3, граница 4) ставилось граничное условие На pressure-outlet.

поверхности пластины (рис.2.3, граница 3) задавалось граничное условие wall с условием прилипания и температурой Tw = 300K.

При решении задачи величина статического давления записывалась на каждом шаге по времени (t = 1е-7 с) в некоторых точках расчтной области, а именно в набегающем потоке (x = 0.03м; y = 0.015м) и на поверхности пластины в координатах x=0.03 м, 0.040625 м, 0.051875 м, 0.064375 м,

0.076875 м, 0.089375 м, 0.095 м. Далее производился Фурье анализ записанных данных.

2.2. Влияние звукопоглощающих покрытий на развитие возмущений Воздействие акустических возмущений набегающего потока на ударный слой на сплошной пластине, расположенной под углом атаки =30, рассматривалось с помощью пакета 12.1 и пакета, ANSYS Fluent разработанного в ИТПМ СО РАН. На рис.2.4 приведены среднеквадратичные пульсации давления на поверхности сплошной пластины в точках x/L = 0.3, 0.5,

0.77 и 0.95 при различных частотах акустических возмущений. Видно согласие данных, полученных с помощью пакета ИТПМ и пакета ANSYS Fluent. Также видно, что результаты для случая воздействия широкополосного шума согласуются с результатами для монохроматических акустических волн.

На рис.2.5 приведены распределения среднеквадратичных амплитуд пульсаций давления на поверхности пластины для разных частот при воздействии монохроматических акустических волн и широкополосного шума.

Видно, хорошее совпадение расчтных данных.

–  –  –

Рис.2.4. Среднеквадратичные пульсации давления на поверхности пластины в сечениях x/L = (а) – 0.3; (б) – 0.52; (в) – 0.77; (г) – 0.95 для различных частот возмущений; (1,2) – медленные монохроматические акустические волны; (3) – широкополосный шум; (1) – пакет ИТПМ СО РАН; (2,3) – пакет ANSYS Fluent

–  –  –

Рис.2.5. Распределения среднеквадратичных пульсации давления на поверхности пластины для различных частот монохроматических акустических волн (линии) и широкополосного шума (точки): f = (а) – 20 кГц; (б) – 40 кГц; (в) – 60 кГц; (г) – 80 кГц Далее рассматривалось обтекание пластины, расположенной под углом атаки =30, со звукопоглощающим покрытием в виде каналов на интервале от

–  –  –

поля мгновенных и среднеквадратичных пульсаций давления для частоты 20 кГц как для сплошной пластины, так и для пластины со звукопоглощающей вставкой. Видно, что звукопоглощающая поверхность оказывает существенное влияние на пульсации давления в ударном слое. Внутри каналов вставки распространяются затухающие акустические волны (рис.2.6в) и, отражаясь от заглушенного торца канала, образуют стоячую волну (рис.2.6г).

а б

–  –  –

На рис.2.7 представлены также распределения по координате x среднеквадратичных пульсаций давления на поверхности сплошной пластины и пластины со звукопоглощающим покрытием. Видно, что над поверхностью, расположенной сразу за звукопоглощающим покрытием, наблюдается ослабление амплитуды пульсаций давления. Вследствие присутствия волн сжатия на рбрах пор, структура пульсаций давления над звукопоглощающим покрытием имеет пилообразный вид. Для удобства анализа данных здесь и далее величина пульсаций давления над звукопоглощающим покрытием усреднялась по каждой поре (линия 3 на рис.2.7).

–  –  –

Рис.2.7. Среднеквадратичные пульсации давления на поверхности сплошной пластины (1) и пластины со звукопоглощающей вставкой (2); величина пульсаций, усредннная по каждой поре (3): газ азот, M=21, ReL= 6104, L=0.1м, Tw/T0= 0.25, T0=1200К, A=0.02, f=20кГц, =30, = 4мм На рис.2.8 кривой 1 приведены экспериментальные спектры значений амплитуд пульсаций давления, нормированные на значения амплитуды на сплошной поверхности, в сечении х = 0.95 со вставкой из полузамкнутых цилиндров с шириной каналов = 4мм на пластине, расположенной под углом атаки 30°. Трубные эксперименты проведены Мироновым С.Г. и Цырюльниковым И. С. на трубе Т-327А ИТПМ СО РАН [120]. С помощью пьезодатчика измерялись пульсации давления на поверхности сплошной пластины p’w(solid) и пластины со звукопоглощающей вставкой p’w(porous) в точке x = 0.95L.

Символами 2 (рис.2.8) показаны расчетные данные, полученные с помощью пакета ИТПМ с использованием граничного условия полного поглощения акустических возмущений (см. параграф 2.1.1). Данные численного моделирования показывают качественное совпадение с результатами эксперимента на частотах f 20 кГц. На частотах f 20 кГц длина волны внешних акустических возмущений сравнима с длинной звукопоглощающего участка, что не позволяет говорить о правомерности использования граничного условия (9) для моделирования звукопоглощающего покрытия.

Также на рис.2.8 приведены расчетные данные, полученные при воздействии внешних монохроматических акустических волн (символы 3), и широкополосного шума (символы 4) на ударный слой на пластине с заданием

–  –  –

Рис.2.8. Эффективность поглощения пульсаций давления в ударном слое звукопоглощающими покрытиями в сечении x=0.95L на поверхности пластины, ширина каналов пор = (а) – 4мм и (б) – 1мм: (1) – экспериментальные данные, данные численного моделирования (2) – с граничным условием на поверхности пластины, (3) – с заданием геометрии звукопоглощающей вставки и воздействием монохроматической акустической волны, (4) с заданием геометрии звукопоглощающей вставки и воздействием широкополосного шума

–  –  –

Рис.2.9. Среднеквадратичные пульсации давления на поверхности сплошной пластины (1) и пластины со звукопоглощающим покрытием = 4мм (2) при воздействии медленной акустической волны A=0.02: (а) f = 10кГц; (б) 20кГц; (в) 30кГц, газ азот, M=21, ReL = 6104, L=0.1м, Tw/T0= 0.25, T0=1200К, =30 Анализ расчетных данных по амплитудам пульсаций давления на поверхности пластины в зависимости от частоты (рис.2.9) показал, что в ударном слое над звукопоглощающей поверхностью наблюдаются их вариации в продольном направлении, т.е. формируется периодическая структура пульсаций давления, напоминающая стоячую волну (линия 2 на рис.2.9).

Оказалось, что в случае, когда положение точки измерения пульсаций давления приходится на минимум этой структуры (по аналогии со стоячей волной узел), то амплитуда пульсаций за звукопоглощающей вставкой минимальна (рис.2.9б). В том же случае, когда положение точки измерения пульсаций давления близко к максимуму этой структуры (по аналогии со стоячей волной пучность), то амплитуда пульсаций максимальна (рис.2.9а). С увеличением частоты внешних возмущений вариации амплитуды в этой структуре ослабляются по продольной координате (рис.2.9в).

К образованию вариаций амплитуд пульсаций давления над звукопоглощающим покрытием приводит изменение амплитудно-фазового баланса возмущений ВУС, вызванного звукопоглощающим покрытием. На формирование пульсаций в ВУС оказывает влияние наличие незатухающего, прошедшего за УВ акустического возмущения, распространяющегося к поверхности пластины, отражение волн от поверхности, распространение затухающих акустических мод от носика модели и колебаний самой УВ под влиянием внешних и внутренних возмущений. Наличие звукопоглощающего покрытия приводит к значительному ослаблению отражнных от поверхности волн.

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

–  –  –

как ширина пор, длина звукопоглощающего участка, начальное x1por и por конечное x2 положение звукопоглощающего покрытия. Во всех расчетах коэффициент пористости был одинаковым и равным 80%.

–  –  –

Рис.2.10. Распределение среднеквадратичных пульсаций давления на поверхности пластины без звукопоглощающей вставки (1) и со вставкой с размером пор = 4мм–(2); = 1мм–(3):

M=21, ReL= 6104, L=0.1м, Tw/T0=0.25, T0=1200К, A=0.02, f=20кГц, =30 Исследование влияния ширины пор на структуру пульсаций давления на поверхности пластины показало, что ширина пор в диапазоне от 1мм до 4мм качественно не влияет на структуру пульсаций вблизи поверхности пластины при f 20кГц. На рис.2.10 представлены распределения среднеквадратичных пульсаций давления на поверхности сплошной пластины и пластины со звукопоглощающей вставкой с шириной пор 1мм и 4 мм. Видно хорошее совпадение формы и положения экстремумов вариаций над звукопоглощающей вставкой.

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

На рис.2.11 представлены распределения среднеквадратичных пульсаций давления на поверхности сплошной пластины и пластины со звукопоглощающей вставкой в интервале от x1por = 40 мм и x1por = 10 мм до por x2 = 90 мм с шириной пор 1мм.

–  –  –

Рис.2.11. График среднеквадратичных пульсаций давления на поверхности сплошной пластины (1) и с разным положением начала звукопоглощающего покрытия, (2) – x1por =40мм;

(3) – x1por =10мм; M=21, ReL = 6104, L=0.1м, Tw/T0 = 0.25, T0=1200К, A=0.02, f=20кГц, =30, = 1мм, x2 =90мм por

–  –  –

por x2 = 240 мм с шириной пор 1мм. Видно, что влияние звукопоглощающего покрытия на амплитуду пульсаций имеет локальный характер, т.е. вниз по потоку степень влияния уменьшается, и амплитуда пульсаций нарастает.

Удлинение звукопоглощающего покрытия вниз по потоку не изменяет структуру и амплитуду пульсаций на начальном участке покрытия (до х = 90мм), но уменьшает амплитуду пульсаций вниз по потоку (х 90мм).

2.4. Комбинированное управление возмущениями ВУС с помощью интерференционного метода и метода звукопоглощающих покрытий В предыдущей главе показано, что при углах атаки 15 энтропийновихревыми возмущениями в ВУС можно эффективно управлять вводимыми контролируемыми возмущениями с поверхности (интерференционный метод).

При увеличении угла атаки ( 15) в ударном слое генерируются незатухающие акустические возмущения, которыми, как было показано выше, можно управлять с помощью метода звукопоглощающих покрытий. В данном параграфе рассмотрена возможность комбинированного применения двух этих методов при углах атаки 15 20 для управления как энтропийновихревыми, так и акустическими возмущениями в ВУС одновременно.

Численное моделирование развития возмущений при одновременном применении интерференционного метода и метода звукопоглощающих покрытий для управления возмущений в ВУС выполнено с помощью пакета, разработанного в ИТПМ СО РАН, с использованием граничного условия полного подавления акустических возмущений (см. параграф 2.1).

На рис.2.13 показаны среднеквадратичные амплитуды пульсаций плотности (рис.2.13а) и пульсаций давления (рис.2.13б) на поверхности пластины. В случае применения только интерференционного метода запуска начальных возмущений) видно существенное (противофазного уменьшение пульсаций плотности и давления на поверхности пластины (линия 3). При помещении звукопоглощающего покрытия в область роста

–  –  –

управления дает двукратное снижение пульсаций по сравнению с применением только интерференционного метода и четырхкратное снижение пульсаций по сравнению с возмущениями, генерируемыми внешними акустическими волнами (линия 1).

–  –  –

Рис.2.13. Амплитуды среднеквадратичных пульсаций плотности (а) и давления (б) на поверхности пластины: М=21, ReL = 1,44105, L=0.24м, =15°, f=38 кГц, А=0.003, B=0.73:

1 – при воздействии внешней монохроматической акустической волны, 2 – при воздействии возмущения вдув-отсос, 3 – интерференционный метод воздействия; 4 – комбинация интерференционного метода и метода звукопоглощающих покрытий Т.о. для углов атаки 15 20, где в ВУС (M = 21, ReL = 1,44105) развиваются одновременно энтропийно-вихревые и акустические возмущения, следует применять комбинированный метод управления для достижения наибольшей эффективности подавления возмущений в ВУС.

–  –  –

Исследовано взаимодействие вязкого ударного слоя на пластине с внешними монохроматическими акустическими возмущениями медленной моды и с широкополосным акустическим шумом аэродинамической трубы Т-327А.

Продемонстрирована высокая эффективность поглощения пульсаций давления (до 50%) звукопоглощающим покрытием на частотах более 15 кГц, однако в области частот менее 15 кГц получено некоторое усиление пульсаций давления по сравнению с пластиной без звукопоглощающих материалов, связанное с возникновением на этих частотах в ударном слое структуры типа продольной стоячей волны. Расчтные данные по влиянию пористых звукопоглощающих материалов на пульсации давления на поверхности пластины согласуются с экспериментальными данными, полученными в аэродинамической трубе Т-327А.

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

Положение начала звукопоглощающего участка значительно влияет на вариации пульсаций давления на поверхности пластины в отличие от положения конца звукопоглощающей вставки. Обнаружено, что варьирование ширины пор в звукопоглощающей вставке длиной 50 мм в пределах 1 мм 4 мм не влияет на е эффективность при частотах возмущений больше 20 кГц.

Впервые численно реализован комбинированный метод управления интенсивностью возмущений в гиперзвуковом ударном слое на пластине под углами атаки свыше 15: одновременное подавление энтропийно-вихревых пульсаций на верхней границе пограничного слоя интерференционным методом и пульсаций давления на поверхности модели с помощью метода звукопоглощающих покрытий.

Результаты приведенных исследований представлены в [120-126].

ГЛАВА 3 Влияние реальных свойств газа на управление развитием возмущений в гиперзвуковых течениях Термин «реальный газ» в аэродинамике обычно используется для того, чтобы подчеркнуть отличие термодинамических и кинетических свойств реальной газовой среды от термодинамических свойств ее простейшей модели

– совершенного газа с постоянным показателем адиабаты. Причины такого отличия могут быть разные:

высокая плотность газа, благодаря чему уравнение состояния реального газа отличается от уравнения состояния совершенного газа из-за действия межмолекулярных сил;

влияние физико-химических процессов, протекающих в газе при высоких температурах: возбуждение колебательных степеней свободы, диссоциация, химические реакции и т.д.;

конденсация, ионизация и т.д.

Различные физико-химические процессы в газах протекают с конечными скоростями и характеризуются своими временами релаксации (или длинами релаксации). В зависимости от соотношения времен релаксации и характерного газодинамического времени течения газа процессы могут быть равновесными (или близкими к равновесным), неравновесными и замороженными. Для различных гиперзвуковых аппаратов и различных траекторий их движения могут существовать различные типы течений.

При входе летательных аппаратов в атмосферу Земли с орбитальными скоростями равновесная температура воздуха за прямой ударной волной превышает 6000К, а характерные длины релаксации изменяются от многих километров на больших высотах до долей миллиметра вблизи Земли, так что в процессе спуска могут встретиться все типы неравновесных течений.

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

Результаты расчетов, представленные в предыдущих параграфах, получены для «холодного» гиперзвукового ударного слоя, а именно при сравнительно небольшой температуре торможения Т0 ~ 1000К, когда температура азота в набегающем потоке не превышает 13К, а максимальная температура в ударном слое меньше 400К. В этих условиях газ не проявляет выше перечисленных реальных свойств. В условиях реального гиперзвукового полета в высокотемпературном газе вблизи поверхности летательного аппарата идут физические процессы, включающие возбуждение различных степеней свободы молекул, диссоциацию и ионизацию молекул, которые оказывают существенное воздействие, как на среднее течение, так и на устойчивость и ламинарно-турбулентный переход в гиперзвуковых сдвиговых слоях [52,53,127].

В данной главе представлены результаты исследования влияния колебательных степеней свободы молекул СО2 на среднее течение и развитие возмущений в вязком ударном слое на пластине, обтекаемой потоком воздуха, углекислого газа и их смесями, с учетом колебательной неравновесности.

Выбор СО2 в качестве рабочего газа определялся низкой характерной температурой возбуждения колебательных состояний молекул. Сравнение расчетных данных с результатами трубных экспериментов, проведнных в аэродинамической трубе ИТ-302М ИТПМ СО РАН, показало, что влияние реальных свойств газа невозможно без учета процессов, идущих в тракте аэродинамической трубы, где формируется поток, влияющий как на характеристики среднего течения, так и на развитие возмущений. Поэтому в схему моделирования включены процессы, идущие в набегающем потоке (в сопловом тракте аэродинамической трубы).

Процессы возбуждения и неравновесности колебательных степеней свободы, оказывая влияние на развитие возмущений, также могут влиять на эффективность методов управления возмущениями ВУС. В данной главе исследована эффективность метода пористых звукопоглощающих покрытий в условиях высокотемпературных термически неравновесных течений.

3.1. Постановка задачи Задача обтекания пластины гиперзвуковым потоком газа с учетом реальных свойств рассматривается в рамках модели совершенного газа.



Pages:     | 1 || 3 | 4 |
 
Похожие работы:

«Игнатенко Евгений Александрович МЕТОДИКА РАССЛЕДОВАНИЯ НЕЗАКОННОЙ ПЕРЕСЫЛКИ НАРКОТИЧЕСКИХ СРЕДСТВ Специальность: 12.00.12 – «Криминалистика; судебно-экспертная деятельность; оперативно-розыскная деятельность» Диссертация на соискание ученой степени кандидата юридических наук Научный руководитель: доктор юридических наук, доцент П.В....»

«Горбунов Юрий Вадимович ИСПОЛЬЗОВАНИЕ ВУЗОВСКИХ НАУЧНЫХ РАЗРАБОТОК ПРИ ФОРМИРОВАНИИ МЕХАНИЗМА УСТОЙЧИВОГО РАЗВИТИЯ ПРОМЫШЛЕННЫХ ПРЕДПРИЯТИЙ Диссертация на соискание ученой степени кандидата экономических наук Специальность 08.00.05 – Экономика и управление народным хозяйством (экономика, организация и управление предприятиями,...»

«ГОРПИНЧЕНКО Ксения Николаевна ОРГАНИЗАЦИОННО-ЭКОНОМИЧЕСКИЙ МЕХАНИЗМ УПРАВЛЕНИЯ ИННОВАЦИОННЫМ ПРОЦЕССОМ: ТЕОРИЯ, МЕТОДОЛОГИЯ И ПРАКТИКА (на примере зернового производства) Специальность 08.00.05 – экономика и управление народным хозяйством: управление инновациями ДИССЕРТАЦИЯ на соискание учной степени доктора экономических наук...»

«Ботнарюк Марина Владимировна Организационно-экономический механизм повышения конкурентоспособности морских транспортных узлов на принципах маркетинга взаимодействия Специальность 08.00.05 «Экономика и управление народным хозяйством (маркетинг)» Диссертация на соискание ученой степени доктора экономических наук Научный...»

«Смирнова Елена Юрьевна Свойства корковых нейронов и механизм обработки информации о цвете в первичной зрительной коре 03.01.02 Биофизика Диссертация на соискание ученой степени кандидата физико-математических наук Научный руководитель: доктор физико-математических наук, Чижов Антон Вадимович Санкт-Петербург – 2015 ОГЛАВЛЕНИЕ 1. ВВЕДЕНИЕ 1.1 Актуальность...»

«ЧЖАН ГОФАН ВЛИЯНИЕ РАЗГРУЗОЧНЫХ ПРОБ НА БИОМЕХАНИЧЕСКИЕ ПАРАМЕТРЫ ГЛАЗА ПРИ ПЕРВИЧНОЙ ОТКРЫТОУГОЛЬНОЙ ГЛАУКОМЕ 14.01.07 глазные болезни Диссертация на соискание ученой степени кандидата медицинских наук Научный руководитель: Д.м.н., Макашова Надежда Васильевна М о с к в а – 2016 ОГЛАВЛЕНИЕ Список сокращений ВВЕДЕНИЕ ГЛАВА I. Обзор литературы Биомеханика склеры. 1. Терминология: понятия биомеханики, ригидности и...»

«Летнер Оксана Никитична ИССЛЕДОВАНИЕ ОСОБЕННОСТЕЙ ДИНАМИКИ АСТЕРОИДОВ, СБЛИЖАЮЩИХСЯ С ЗЕМЛЕЙ Специальность 01.03.01 – астрометрия и небесная механика Диссертация на соискание ученой степени кандидата физико–математических наук Научный руководитель доцент, к.ф.-м.н. Л.Е. Быкова Томск – 2015 СОДЕРЖАНИЕ...»

«ИЗ ФОНДОВ РОССИЙСКОЙ ГОСУДАРСТВЕННОЙ БИБЛИОТЕКИ Карницкая, Элла Николаевна Формирование экономического механизма развития здравоохранения региона в условиях социально­ориентированной рыночной среды Москва Российская государственная библиотека diss.rsl.ru Карницкая, Элла Николаевна Формирование экономического механизма развития здравоохранения региона в условиях социально­ориентированной рыночной среды : [Электронный ресурс] : Дис. . канд. экон. наук : 08.00.05. ­...»

«АРТЕМЬЕВ АНДРЕЙ БОРИСОВИЧ Коррупция в механизме функционирования государства (теоретико-правовое исследование в рамках эволюционного подхода) Специальность 12.00.01 – теория и история права и государства; история учений о праве и государстве ДИССЕРТАЦИЯ на соискание ученой степени доктора юридических наук Научный консультант: доктор юридических наук профессор С.А.КОМАРОВ...»

«Бекежанова Виктория Бахытовна УСТОЙЧИВОСТЬ НЕИЗОТЕРМИЧЕСКИХ ЖИДКОСТЕЙ В РАЗЛИЧНЫХ МОДЕЛЯХ КОНВЕКЦИИ 01.02.05 механика жидкости, газа и плазмы Диссертация на соискание ученой степени доктора физико-математических наук Научный консультант доктор физико-математических наук, профессор В. К. Андреев Красноярск 2015 ОГЛАВЛЕНИЕ ВВЕДЕНИЕ Глава 1....»

«АБРАМОВ АЛЕКСЕЙ ЮРЬЕВИЧ СОВЕРШЕНСТВОВАНИЕ МЕХАНИЗМА РЕГУЛИРОВАНИЯ ОБОРОТА НАРКОТИЧЕСКИХ СРЕДСТВ И ПСИХОТРОПНЫХ ВЕЩЕСТВ В РОССИЙСКОЙ ФЕДЕРАЦИИ 14.02.03 – Общественное здоровье и здравоохранение Диссертация на соискание ученой степени доктора медицинских наук Научный консультант Доктор медицинских наук, профессор Заслуженный деятель науки РФ Михайлова Ю.В. Москва 2015 год СОДЕРЖАНИЕ Введение Глава I. Незаконный оборот наркотиков и наркомания глобальные проблемы современности...»

«Ращектаев Александр Сергеевич Фармако-клиническое обоснование применения «Геприма для кошек» при жировом гепатозе 06.02.03 – Ветеринарная фармакология с токсикологией диссертация на соискание учной степени кандидата ветеринарных наук Научный руководитель: доктор ветеринарных наук, доцент Щербаков П.Н. Троицк – 2015 Оглавление Перечень сокращений в диссертации ВВЕДЕНИЕ Обзор литературы 1. 1.1 Гепатопротекторы....»

«Павлов Александр Борисович ИССЛЕДОВАНИЕ И ОПТИМИЗАЦИЯ СИСТЕМ ЭЛЕКТРОПОДОГРЕВА ПОЛИМЕРНЫХ ТРУБОПРОВОДОВ 05.09.01 — Электромеханика и электрические аппараты Диссертация на соискание ученой степени кандидата технических наук Научный руководитель доктор технических наук профессор Плохов И.В. Псков 2015 ОГЛАВЛЕНИЕ ВВЕДЕНИЕ ГЛАВА 1:...»

«СТЕПАНЕНКО Сергей Владимирович ПРОГНОЗ ДЕФОРМАЦИЙ ГРУНТОВОГО МАССИВА ПРИ СТРОИТЕЛЬСТВЕ ПОЛУЗАГЛУБЛЕННЫХ ПОДЗЕМНЫХ СООРУЖЕНИЙ СПОСОБОМ «СТЕНА В ГРУНТЕ» Специальность 25.00.20 Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика...»

«ГРАЧЕВ Николай Николаевич РАЗВИТИЕ СИСТЕМЫ УПРАВЛЕНИЯ ОХРАНОЙ ТРУДА В СЕЛЬСКОМ ХОЗЯЙСТВЕ Специальность 08.00.05 – Экономика и управление народным хозяйством (экономика труда) Диссертация на соискание ученой степени доктора экономических наук Рязань – СОДЕРЖАНИЕ Стр. СОДЕРЖАНИЕ...»

«Максимов Роман Александрович МЕХАНИЗМ ДЕЙСТВИЯ ПРАВА В ЧРЕЗВЫЧАЙНЫХ СИТУАЦИЯХ (Общетеоретический аспект) Специальность 12.00.01 – теория и история права и государства; история учений о праве и государстве Диссертация на соискание ученой степени кандидата юридических наук Научный руководитель – доктор юридических наук, доцент Фомин...»

«Дундуков Михаил Юрьевич РАЗВЕДКА В ГОСУДАРСТВЕННОМ МЕХАНИЗМЕ США (ИСТОРИКО-ПРАВОВОЙ АСПЕКТ) Диссертация на соискание ученой степени доктора юридических наук Специальность: 12.00.01 — теория и история права и государства; история учений о праве и государстве Научный консультант: доктор юридических наук, профессор Томсинов Владимир Алексеевич МОСКВА ВВЕДЕНИЕ Глава 1. РАЗВИТИЕ РАЗВЕДЫВАТЕЛЬНОЙ ДЕЯТЕЛЬНОСТИ В США (КОНЕЦ...»

«ПАЛАНКОЕВ ИБРАГИМ МАГОМЕДОВИЧ Обоснование параметров технологии проходки шахтных стволов в искусственно замороженных породах Специальности: 25.00.22«Геотехнология (подземная, открытая и строительная)» 25.00.20 «Геомеханика, разрушение горных пород,...»

«Карыев Леонид Геннадьевич ЗАКОНОМЕРНОСТИ И МЕХАНИЗМЫ ПРОЦЕССОВ, ОБУСЛОВЛЕННЫХ ТЕРМОЭЛЕКТРИЧЕСКИМ ВОЗДЕЙСТВИЕМ, В ИОННЫХ КРИСТАЛЛАХ 01.04.07 – Физика конденсированного состояния Диссертация на соискание ученой степени доктора физико-математических наук Научный консультант: доктор физико-математических наук, профессор В.А. Фдоров Тамбов 2015 Автор выражает...»

«ГОРПИНЧЕНКО Ксения Николаевна ОРГАНИЗАЦИОННО-ЭКОНОМИЧЕСКИЙ МЕХАНИЗМ УПРАВЛЕНИЯ ИННОВАЦИОННЫМ ПРОЦЕССОМ: ТЕОРИЯ, МЕТОДОЛОГИЯ И ПРАКТИКА (на примере зернового производства) Специальность 08.00.05 – экономика и управление народным хозяйством: управление инновациями ДИССЕРТАЦИЯ на соискание учной степени доктора экономических наук...»









 
2016 www.konf.x-pdf.ru - «Бесплатная электронная библиотека - Авторефераты, диссертации, конференции»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.