ОГЛАВЛЕНИЕ
стр.
УСЛОВНЫЕ ОБОЗНАЧЕНИЯ ......................................... 4
ВВЕДЕНИЕ ....................................................6
ГЛАВА I. Применение коротковолнового приближения для расчета распространения трехмерных акустических возмущений в осесимметричных каналах переменного сечения со звукопоглощающими стенками при наличии стационарного потока газа......................................19
1.1. Постановка задачи о распространении волн в каналах переменного сечения с облицованными стенками
и вывод основных соотношений ....................... 19
1.2. Примеры распространения акустических волн в каналах, моделирующих тракты ВРД......................30
ГЛАВА II.Разработка численного метода расчета трехмерных акустических полей в областях сложной конфигурации .....................................................'. . . 39
2.1. Постановка задачи расчета акустических полей в области при наличии стационарного потока газа с учетом излучения волн в открытое пространство . . 39
2.2. Построение алгоритма расчета на основе метода конечных элементов (МКЭ) с использованием вариационного принципа..........................................44
2.3. Исследование эффективности схем МКЭ разного порядка точности при решении задач излучения звука в случае покоящейся среды .............................. 52
2.4. Исследование эффективности схемы ИСЗ второго по-
-з-
рядка точности для расчета акустических полей при наличии однородного потока газа ................. 61
ГЛАВА III.Результаты численных расчетов трехмерных акустических полей в областях сложной конфигурации при наличии неоднородного потока газа ............ 66
3.1. Метод расчета стационарного безвихревого сжимаемого потока газа...........................................66
3.2. Исследование эффективности схемы МКЭ второго порядка для расчета акустических полей при наличии неоднородного стационарного потока газа . . 75
3.3. Результаты иллюстративных расчетов излучения звука из канала, моделирующего воздухозаборник ВРД.....................................................79
ВЫВОДЫ...........................................................83
ЛИТЕРАТУРА......................................................167
- УСЛОВНЫЕ ОБОЗНАЧЕНИЯ
г, $) t В
и ’ ^ го-
Р’?>Г
и^,*,А
М
£ V £
к,Х
V
(т,п)
кр
ч
э
В, Б Т У
4>
- цилиндрическая система координат
- время
- безразмерная осевая координата
- составляющие вектора скорости газа вдоль осей X7Г ? 'д' соответственно
- давление, плотность и показатель адиабаты газа
- осевая, радиальная составляющие вектора скорости газа, плотность и скорость звука в стационарном потоке
- число Маха
- малые параметры, характеризующие градиенты стаци-онарных параметров и амплитуду возмущений
физическая частота возмущений осевое и радиальное волновые числа функция Бесселя и ее аргументы фаза возмущений
номера моды (тангенциальной, радиальной) приведенное осевое волновое число и частота акустическая проводимость амплитуда возмущений коэффициенты отражения коэффициент прохождения потенциал акустической скорости ■ полный потенциал скорости (для расчета стационарного течения)
• якобиан преобразования
-5-
Нижняе индексы
* - характерные значения параметров
0 - значения параметров в "О" приближении
1 - значения параметров в "I" приближении
М - значение на стенке канала
е - значение на конечном элементе
Верхние индексы
- комплексное сопряжение т - транспонирование -1 - обратный оператор
-6-
X-€ «
ВВЕДЕНИЕ
Быстрое развитие в последние десятилетия авиационного транспорта привело к сильному увеличению уровня зашумленности на местности, особенно в окрестности аэропортов. В связи с этим встала задача снижения уровня шума, излучаемого авиационными силовыми установками с воздушно-реактивными двигателями (ЕРД). Разрабатываемые на последующие годы стандарты ICA0 на уровень шума, производимого самолетами на местности, особенно на режимах взлёта и посадки, предусматривает постоянное снижение указанного уровня. Разработка мероприятий, направленных на уменьшение генерируемого и выносимого из силовой установки с ВРД шума, требует развития теории возбуждения и распространения акустических волн в каналах и методов расчета указанных явлений. При этом важно уметь достаточно точно учитывать влияние различных факторов, таких, как наличие и свойства звукопоглощающих конструкций (ЗПК), геометрии каналов и обечайки, режима обтекания и т.д.,На характеристики излучае- 47 мого шума..развитие аналитических подходов, предназначенных vy для решения этой задачи, в некотором смысле ограничено, поскольку с их помощью удается получить решение лишь в ограниченном числе наиболее простых случаев. Вместе с тем, развитие в последние годы численных методов, а также рост быстродействия и памяти ЭВл позволили перейти к прямому численному моделированию процессов распространения акустических волн в каналах ВРД и их излучению в открытое пространство.
Настоящая работа посвящена разработке методов расчета распространения волн в каналах переменной площади поперечного сечения при налички стационарного потока газа и звукопог-
-7-
лощающих конструкций. Б первой главе диссертации для оценки влияния геометрии канала, стационарного потока и свойств звукопоглощающей облицовки на изменение амплитуды распространяющегося по каналу акустического возмущения развивается аналитический метод, являющийся обобщением метода работы [35] на случай, когда стенки канала облицованы ЗПК. Как и в случае каналов с твердыми стенками, в облицованном канале могут существовать участки, где становится существенным отражение акустических волн, связанное с наличием точек поворота ВКБ-реше-ния. Вне указанных участков в рассматриваемом приближении отражение несущественно. Ещё одна особенность полученного решения связана с неоднозначностью волновых чисел как функций комплексной акустической проводимости. Вблизи точек ветвления имеет место взаимодействие различных радиальных мод и одномодовое рассмотрение оказывается несправедливы!.:. В работе показано, что учет особенностей, связанных с точками ветвления, является существенным с точки зрения повышения эффективности шумоглушения. Вместе с тем, в случае, когда площадь канала меняется достаточно резко, разработанный метод не позволяет получать результаты с необходимой точностью, к тог,у же с его помощью невозможно учесть эффекты, связанные с отражением волн от открытого конца канала и их излучением в открытое пространство. В связи с этим,во второй главе для расчета акустических полей, возникающих при распространении и излучении волн в областях при наличии стационарного потока газа, разрабатывается более общий численный метод. Указанный метод основан на вариационном принципе, предложенном в работе [&4 ] и описывающем распространение линейных гармонических по времени акустических возмущений в кусочно-изоэнтропических кусочно-изоэнергетических потоках газа,содержащих поверхности танген-
-8-
циального разрыва типа границы струи. Для поиска стационарной точки функционала (и определения таким образом неизвестных узловых значений комплексного потенциала возмущения) использован метод конечных элементов (ШЗ) с квадратичной аппроксимацией искомой функции на 8-узловоы изопараметрическом элементе. Результирующая система линейных уравнений решалась методом исключения Гаусса.
Для расчета акустических полей, возбуждаемых волной заданной моды, необходимо задавать параметры стационарного потока, такие как число Маха и плотность, в расчетных точках. Вообще говоря, указанные параметры могут быть получены любым способом, однако ясно, что точность их вычисления будет влиять на характеристики излучения. В третьей главе настоящей работы для расчета чисто дозвукового обтекания концевой части канала (например, воздухозаборника) безвихревым изоэнтро-пическим потоком сжимаемого газа был разработан метод расчета, основанный на вариационном цринципе Бейтмена Е'З'і] .
Для аппроксимации неизвестной функции (с помощью потенциала скорости) использовались те же 8-узловые изопараметрические элементы. Итерации по нелинейности осуществлялись с помощью матричного аналога метода Ньютона-Рафсона, что обеспечило их очень быструю сходимость. Сам расчет стационарного течения проводился на сетке, предназначенной для расчета акустики, при этом значения необходимых параметров потока вычислялись непосредственно в узлах этой сетки. Все это позволило отказаться от использования для расчета стационарных параметров программ, основанных на других методах (например, конечно-разностных), избежать достаточно сложных пересчетов параглет-ров с сетки на сетку и связанных с этим дополнительных погреш-
-9-
ностей. Результаты выполненных расчетов акустических полей позволили выявить эффективность разработанного метода применительно к задача!/, излучения волн из каналов переменной площади при наличии стационарного потока газа. Метод позволяет такие рассчитывать и чисто канальные задачи о распространении волн в неоднородных газовых потоках. В качестве иллюстрации возможностей метода и созданного комплекса программ приведены результаты расчетов излучения различных мод из каналов, моделирующих воздухозаборник ВРД, работающий на различных режимах. Проведено сравнение результатов созданного метода с результатами, полученными методом Винера-Хопфа и коротковолновым приближением.
К настоящему времени проблемам разработки методов расчета распространения и излучения волн применительно к задачам аэроакустики посвящено уже значительное число работ (см.,например, обзорную статью [ 13 )• Остановимся сначала на чисто каналь-
ных задачах, в которых исследуется распространение звука по каналам ( например, внутренним трактам ВРД) при наличии ЗПК (или без него) в присутствии стационарного потока газа, без учета излучения волн в открытое пространство. Такие задачи представляют интерес в связи с исследованием акустических характеристик внутренних трактов ВРД, самолетных воздуховодов в системах кондиционирования и т.д., то есть в случаях, когда учет эффектов излучения не обязателен. Так в работах [й-41 рассматривается задача о затухании волн в каналах постоянного поперечного сечения (плоских и осесимметричных) при наличии звукопоглощающих покрытий. Для описания свойств БПК использована импеданеная модель, причем свойства покрытия считались постоянными. Стационарный поток газа отсутствовал. Расчет акустических полей проводился с помощью метода Винера-Хопфа [й-]
-10- V \
Учет влияния однородного потопа газа на затухание мод в канале с ЗПК исследован в работах [5- <31 .Вопросы оптимизации свойств ЗПК, то есть подбора его характеристик таким образом, чтобы обеспечить максимальное затухание выбранной распространяющейся моды, изучались в работах Ц 5 ,7 ] , причем поле в канале возбуждалось либо точечным источником [53 , либо заданной при-
ходящей модой [?1 . Рассмотрение велось в рамках аналитических методов. Несмотря на то, что рассмотренные в этих работах простые модели явления дают информацию о качественном поведении решения и помогают выявить его особенности, для целей практики представляет интерес учет влияния реальной конфигурации газовых трактов Г'ТД (сопла и воздухозаборника), по которым происходит вынос шума из двигателя, на эффективность работы звукопоглощающего покрытия. Это влияние является следствием как изменения по длине канала площади его поперечного сечения, так и существенной неодномерности стационарного потока газа. Кроме того, акустические свойства облицовки могут быть разными в разных сечениях канала (в случае, например, секционных ЗПК, состоящих из секций с разной величиной акустической проводимости), что также необходимо учитывать при расчете затухания звука. Учет упомянутых эффектов приводит к сильному усложнению задачи, поэтому применение аналитических методов для ее решения наталкивается на значительные трудности. Вместе с тем имеются работы [10)343 , 'У посвященные указанной задаче и выполненные в рамках различных аналитических подходов. В основе многих решений задачи о распространении звука в каналах с потоком лежит приближенный подход, предложенный Цянь Сэю-Сэнем [21 ] и развитый в дальней-
9 Г
шем Крокко и Чжен Сёнь-и [221 ♦ При таком подходе течение в \,
канале рассматривается как одномерное, а соответствующая сис-
-II- л
I
и
тема уравнений линеаризуется относительно малых нестационарных возмущений. С помощью полученных таким образом уравнений исследуется поведение решений, гармонически зависящих от времени. Указанный подход использовался, например, в работах [2£, 2&-ъо], где с его помощью были получены и исследованы некоторые частные решения. На основании такого же подхода в‘[*4*} было получено решение в ВКБ -приближении для случая коротких волн. Однако, несмотря на наличие значительного количества результатов, полученных в перечисленных работах [26-341 , все они касаются вопросов распространения плоских волн в одномерных потоках газа.
Из работ, посвященных задаче распространения трехмерных акустических волн в каналах прежде всего следует выделить работы Б.Цинна и Л.Крокко [49,20], в которых рассмотрено распространение возмущений типа радиально-тангенциальных мод в осе-сшметричных каналах с медленно меняющейся площадью поперечного сечения, /шторами были получены соотношения, учитывающие нелинейные члены относительно малых возмущений и определяющие отражение таких возмущений от дозвуковой части сопла Лаваля. Средний стационарный поток в этих работах считался одномерным. Хотя соотношения, определяющие распространение возмущений в работах [19,20] были получены, какие-либо результаты расчетов такого распространения там не приведены. В работе [л?] рассмотрено отражение вращающихся акустических мод при их распространении против потока в канале, мало отличающемся от цилиндрического. При этом двумерность среднего потока учтена через влияние малой нормальной составляющей скорости. В работах [24,27] с помощью метода многих масштабов исследовано распространение нормальных акустических мод в плоских и кольцевых каналах с медленно меняющейся площадью поперечного сечения, причем давление
и плотность в среднем считались постоянными во всем канале ' при наличии градиентов продольной и поперечной скоростей.
В работах теоретически и экспериментально иссле-
дуется отражение трехмерных акустических возмущений от конической дозвуковой части сопла Лаваля. При этом в работах [36 , а?] в теоретическом анализе стационарный поток рассмотрен в одномерном приближении течения в стоке, что является достаточно грубым приближением реального течения, приводит к бесконечно большому градиенту скорости на звуковой линии и, следовательно, к ошибке в вычислении коэффициента отражения. В работе 1>3] проведено прямое сравнение теоретических и экспериментальных результатов расчета распространения трехмерных волн в каналах. Теоретическая модель основывалась на методе многих масштабов, разработанном ранее в работах [24,2?] для
«''V '
плоских и осе сиг,не тричных каналов с потоком. Исследуемый канал состоял из прямолинейного акустически облицованного осесимметричного канала с внутренним телом переменного сечения С с твердыми стенками). Совпадение результатов расчетов и эксперимента в случае отсутствия стационарного потока было хорошим, а в случае, когда присутствовал поток с '»/'= 60 м/сек, расчетный метод давал неправильное распределение фазы возмущения давления по угловой координате, что авторы объясняют неоднородностью свойств ЗПК в тангенциальном направлении.
Г
Некоторое усовершенствование метода [24,2?] сделано в работе [*1^1 » Сводится оно к тому, что вместо переменных,
которыми являются амплитуды скорости, давления и плотности, используют переменные, в которые входит функция осевого волнового числа, описывающая главное изменение возмущения при его распространении вдоль канала. Такая замена переменных позволяет
-13-
увеличить точность расчетов в случаях, когда нет сильного отражения распространяющейся моды и мало её взаимодействие с соседними радиальными модами. Таким образом, применение аналитических методов ограничено случаями распространения плоских волн в каналах переменного сечения, трехмерных мод в слабонеоднородных каналах при наличии потока, но без ЗПК, либо при наличии ЗПК, но в отсутствие потока. Для практики, однако, представляет интерес учет совместного влияния геометрии канала, свойств ЗПК и неоднородного потока газа. При этом приближение слабого изменения свойств ЗПК и геометрии канала по его длине является вполне приемлемым во многих случаях применительно к каналам авиационных ГТД, поэтому представляется целесообразным разработать метод расчета распространения трехмерных акустических волн в каналах с потоком при наличии ЗПК в рамках коротковолнового приближения [38] . Указанному вопросу и посвящена первая глава настоящей диссертации.
Для полноты укажем, что существует уже значительное число работ,[40-61] , посвященных решению указанной задачи численно с помощью метода конечных элементов [40-4?], конечно-разностного метода [46,49] , метода интегральных уравнений [_?6] • В указанных работах также решено значительное число задач о распространении волн в каналах переменной плопрди с учетом как стационарного потока (правда, в одномерном приближении), так и звукопоглощающих облицовок. Однако, помимо того, что численные методы не позволяют получить качественных особенностей решения, они обладают ограничениями другого рода, связанными с быстродействием и объемом памяти современных ЭВМ. Дело в том, что для задач акустики все численные методы в конце концов приводят к необходимости решать системы линейных алгебраических уравнений
-14- V
большого порядка относительно узловых значений неизвестной функции (как правило, комплексной амплитуды потенциала, либо амплитуд давления и скорости). Вместе с тем, рост быстродействия и памяти ЭВМ, совершенствование численных методов, а также необходимость оптимизации акустических свойств реальных устройств, обладающих сложной геометрией, приводят ко все более широкому внедрению численных методов. Наиболее перспективным в настоящее время является метод конечных элементов (МКЭ) в сочетании с вариационным принципом. Помимо легкости описания сложной геометрии объектов, присущей МКЭ вообще, в рамках такой постановки оказывается достаточно легко учитывать различные типы граничных условий.До недавнего времени, однако, упомянутый вариационный принцип существовал только для самых простых случаев течения. В СВЯЗИ С ЭТИМ, В большинстве работ [40-47, 50,
5'4 ,59,60], посвященных расчету распространения волн с помощью МКЭ, использовался либо метод Галеркина, либо метод наименьших квадратов, а в качестве искомых переменных (в случае неоднородного потока газа) брались амплитуды акустической скорости, давления и плотности. Все это существенно увеличивает порядок результирующей системы уравнений и сильно снижает эффективность метода вцелом. Вместе с тем, к настоящему времени с помощью V указанного метода решено уже значительное число в основном чисто канальных задач по расчету распространения плоских и трехмерных волн в каналах переменного сечения с потоком и без него, при наличии ЗПК и в каналах с твердыми стенка!,и.Для моделирования свойств ЗПК использовалась импедансная модель, поскольку она достаточно правильно отражает свойства ЗПК, к тому же граничное условие в виде алгебраической связи искомых величин наиболее просто реализуется в програмках расчета. Стационарное те-
- 15 -
чение в канале считалось одномерным [ 47 ] , квазиодномер-
ным - с учетом вертикальной составляющей скорости £5<э> <£.03 либо рассчитанным по несжимаемой теории с приближенным учетом сжимаемости [693 . Для аппроксимации неизвестных функций используются как линейные треугольные элементы [ ь«э 3 , так и изопараметрические восьмиузловые [44,47] и эрмитовы [45,51] элементы. Таким образом, из проведенного анализа работ следует, что е рамках используемых постановок (с методами Галеркина и наименьших квадратов) численные методы расчета распространения трехмерных волн в каналах в случае больших частот и мод колебаний оказываются недостаточно эффективными из-за больших потребных объемов памяти ЭЖ и времени расчета. Вместе с тем, во многих практических ситуациях допустимо использование приближенных аналитических подходов.
Несколько иначе дело обстоит для задачи расчета излучения трехмерных акустических волн из открытого конца канала. Для решения этой задачи первоначально использовался метод Винера-Хоп-фа[б2] , позволивший решить задачу об излучении звука из круглой полубеоконечной трубы постоянного сечения в отсутствии потока газа. В дальнейшем в рамках того же метода была решена аналогичная задача в случае отличного от нуля числа Маха однородного потока [6Т>] , а также в случае, когда газ в трубе и спутном потоке течет с разными скоростями [643. В работах [е5,ба] была решена задача об излучении заданной тангенциально-радиальной моды из отрезка цилиндрической трубы,(моделирующей мотогондолу) конечной длины при наличии однородного потока газа. Дальнейшее развитие подхода, использованного авторами этих работ, сделано в работах [67, а?] , где решена задача об излучении мод из трубы в спутном потоке, причем учитывалось
- Київ+380960830922