CC BY
21
УДК 623.464.34
МЕТОДИКА ЧИСЛЕННОГО РАСЧЁТА ГАЗОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СТРУЙНЫХ ЭЛЕМЕНТОВ ПНЕВМОАВТОМАТИКИ СИСТЕМ УПРАВЛЕНИЯ СЛОЖНЫМИ ТЕХНИЧЕСКИМИ СИСТЕМАМИ
© 2017 В.Н. Самсонов, В.Г. Шахов, И.С. Барманов
Самарский национальный исследовательский университет имени академика С.П. Королёва
Статья поступила в редакцию 10.07.2017
В работе представлены результаты моделирования течения газа в каналах струйных элементов пневмоавтоматики на примере двухвходового усилителя. Разработана математическая модель расчета параметров течения газа методом конечных объёмов, базирующимся на численном решении систем дифференциальных уравнений в частных производных Навье-Стокса, осреднён-ных по методу Рейнольдса. Данный метод реализуется с помощью программного комплекса ЛиБуБ СРХ. Для замыкания этих уравнений применяется модель турбулентности, сформулированная в соответствии с предположением Буссинеска, согласно которому рейнольдсовы напряжения могут быть связаны со скоростью средней деформации через турбулентную вязкость. С помощью математической модели выполнена серия численных расчетов, на основании результатов которых построены характеристики струйного двухвходового усилителя.
Ключевые слова: управление сложными техническими системами, струйная техника, пневмоавтоматика, струйный усилитель, математическая модель, численные расчёты.
Пневмоустройства, построенные на использовании струйных элементов, имеют широкое применение в промышленности, например, в расходомерах газов, датчиках скорости, устройствах автоматики и других пневмоустройствах [1-4]. Данные устройства имеют перспективы развития, в связи с чем совершенствование конструкции и улучшение характеристик струйных элементов является весьма актуальным. Добиться этого можно, используя современные вычислительные продукты.
В работе приведены результаты исследований и апробации математической модели для выполнения расчета параметров потока газа в каналах элементов струйной техники миниатюрных размеров. Для выполнения численных расчетов разработана математическая модель струйного устройства. Представлены результаты моделирования течения газа в каналах струйного двухвходового усилителя, который предназначен для усиления входных сигналов по давлению и расходу и выполнения логических операций «отрицание» и «повторение».
Рассматривается задача течения газа в каналах элементов струйного двухвходового усилителя. Моделирование процессов течения газа
Самсонов Владимир Николаевич, доктор технических наук, профессор кафедры основ конструирования машин. E-mail: [email protected] Шахов Валентин Гаврилович, кандидат технических наук, профессор кафедры конструкции и проектирования летательных аппаратов.
Барманов Ильдар Сергеевич, кандидат технических наук, доцент кафедры основ конструирования машин. E-mail: [email protected]
внутри каналов струйных устройств проведено методом конечных объёмов, базирующимся на численном решении систем дифференциальных уравнений в частных производных Навье-Стокса, осреднённых по методу Рейнольдса. Уравнения Навье-Стокса включает уравнение сохранения массы (уравнение неразрывности), уравнение сохранения количества движения и уравнение сохранения энергии. Данный метод реализуется с помощью программного комплекса ЛЫБУБ СРХ. Для замыкания этих уравнений применяется модель турбулентности, сформулированная в соответствии с предположением Буссинеска, согласно которому рейнольдсовы напряжения могут быть связаны со скоростью средней деформации через турбулентную вязкость.
Уравнение неразрывности записывается в виде:
др+^(ри] ) = 0.
дг дхгу 3
Уравнение сохранения количества движения [5, 6]:
(pUUj)
др
--\
dxi
2 дЦ_
3 Meff ~дХ~
для 1 = 1,2,3.
Проекции и. и и. соответствуют компо-
1 3
нентам их, иу, и2 вектора скорости потока; р - плотность; 8. - символ Кронекера. Параметр Р* представляет собой сумму давления
Р и величины (2/3)р41к, где к - турбулентная кинетическая энергия. Оба этих слагаемых соответствуют нормальным напряжениям. Эффективная вязкость Ме^г определяется суммой ламинарной и турбулентной вязкостей + ^. Уравнение сохранения энергии представляется следующим образом:
| Р+ ¿(1 ) =
д1 дх,
дх,
(
дТ
дк
дх Рг дх
Л
+ -
дк тт
ди
1 У
ди
А
дх,
+
дх. дх,.
- и 41
3 дх 1
м
где Н = к + 0,5иЦ] + к - полная энтальпия; к - энтальпия; Т - температура; Рг1 - турбулентное число Прандтля. Также в систему уравнений включено уравнение состояния. Во всех уравнениях по повторяющимся индексам производится суммирование [7, 8]. Для замыкания уравнений неразрывности, сохранения количества движения и сохранения энергии в исследованиях настоящей работы использовалась гибридная двухзональная модель Ментера (модель переноса сдвиговых напряжений) [9].
Уравнения переноса для характеристик турбулентности в соответствии с гибридной двухзональной моделью турбулентности Ментера (ББТ) имеют вид [10]:
д д
-(рк)+—(ои) =
д дх,
=Рк - 0*рак +
дх.
I \ дк
д д д —(рсо)+—(ра>и\) = ар£>2 - 0ра2 + — х д1 дх,. дх,.
/ ч да
дх.
+
2(1-Рх)раШ2 дк да
с
дх1 дх1
В уравнении переноса для к функция Рк , используемая в модели ББТ для предотвращения нарастания турбулентности в застойных областях, записывается следующим образом [11]:
Рк = ш1п(Рк ,10 • £)к), где Ьк =0рка , а -скорость диссипации кинетической энергии, 0* = 0,09. В двузональной модели Ментера применяются две системы констант, составляющие которых отмечены индексами 1 и 2:
3 5
ак1 = 0,85; аШ1 = 0 = —; Ч = 9 ^к2 =1; ^«2 = 0,856; 02 = 0,0828; «2 = 0,44.
Первая система калибрована по пристеночным течениям, а вторая обладает высокой степенью применения для свободных сдвиговых слоёв. Первая система констант используется во внутренней области пограничного слоя и соответствует к -а модели Уилкокса, а вторая - в области пограничного слоя и соответствует стандартной к - е модели. В модели Ментера используется линейная комбинация этих двух систем констант с гладким переключателем вида ср=+ (1 - F)^2, где ^ - связующая функция. Вблизи стенки в значительной части пограничного слоя эта функция должна быть близка к единице. Однако, по мере перехода от стенки и приближения к границе пограничного слоя она стремится к нулю, чтобы обеспечить независимость от внешних условий, характерную для к - е модели. Связующая функция определяется по формуле:
Fl = Л(а^);
arg1 = ш1Й
max
4к 500у
Л
у «
СВка = шах
0,09ау
1 дк да
4Р°«2к
СВкшу2 Л
2рст
«2
а дх1 дх1
10
-10
где arg1 - переменная, СОка - положительная часть перекрёстных диффузионных членов в уравнении переноса.
По мере удаления от стенки а^ ^ 0. Внутри пограничного слоя первый член представляет отношение масштаба турбулентности к расстоянию от стенки и равен 2,5 в логарифмическом подслое. Он исчезает при приближении к границе подслоя. Второй член записан так, чтобы F1=\ в пределах подслоя (когда исключается использование двухпараметрической дис-сипативной модели турбулентности). При этом параметр 1/у2а является константой вблизи стенки, но в логарифмической зоне он стремится к нулю. Третий аргумент в формуле переменной а^ предотвращает вырожденность решения исходной модели к - а при малых значениях а набегающего потока.
В целях наилучшего моделирования ограниченных стенками потоков, в отличие от модели к - а, в модель ББТ введена расширенная формулировка вихревой вязкости ак
У( =-, где вместо абсолютной
шах(а1а, SF2)
величины завихренности потока О, применяемой в ранних версиях модели Ментера, используется абсолютная величина скорости сдвиговых деформаций £ = (2£г7£г7 )1/2.
Функция ¥ имеет вид:
^ = Ш(а^2);
а^2 = тах
4к 500у
0,09^у у2ю
где arg2 - переменная. Для внешнего слоя константы не изменяются, а для внутреннего слоя добавляется константа а1=0,31.
На основе описанных математических моделей проведён расчёт газодинамических характеристик струйного двухвходового усилителя, трёхмерная геометрическая модель которого представлена на рис. 1. На рис. 2 изображён общий вид расчётной сетки. Для моделирования вязкого взаимодействия стенки каналов и газа предусмотрен градиентный рост толщины ячеек с поверхности канала.
На входах в каналы питания задана величина скорости потока. Скорость на входе в канал питания нижней платы струйного двухвходового усилителя равна 16,5 м/с, а на входе в канал питания верхней платы 14,5 м/с. Величины скоростей подобраны таким образом, чтобы дав-
ление газа на входе в каналы питания соответствовало заданной величине 20 кПа.
Управляющие сигналы, поступающие на нижнюю плату устройства, моделируются граничным условием скорость на входе. Величина скорости потока на входе в управляющие каналы переменная и подбирается из условия соответствия давления газа в канале заданному диапазону. Необходимость перехода от давления к скорости потока при назначении граничных условий на входах в каналы питания и управления диктуется требованиями устойчивого сходящегося итерационного процесса интегрирования задачи.
На входах в дренажные отверстия присвоено условие равенства статического давления, на выходе - атмосферному давлению. На выходе из каналов, где формируются выходные сигналы, граничное условие обеспечивает равенство статического давления потока на выходе из канала атмосферному давлению, подбирая потребный для этого расход. Стенкам каналов логического устройства присвоено условие непротекания. Данное граничное условие позволяет также мо-
а)
б)
Рис. 1. Основные составные части (а) и геометрическая модель внутренних каналов (б)
струйного двухвходового усилителя
Рис. 2. Расчетная сетка струйного двухвходового усилителя
Таблица 1. Результаты расчета параметров потока в каналах усилителя при увеличении давления в канале управления Х1
Давление в каналах управления, Па
Поле скоростей
Поле давлений
Рх1=0, Рх2=0.
Рх1=590, Рх2=0.
Рх1=3416, Рх2=0.
Рх1=4814, Рх2=0.
Таблица 2. Результаты расчета параметров потока в каналах усилителя при уменьшении давления в канале управления Х1
Давление в каналах управления, Па
Поле скоростей
Поле давлений
Рх1=7792, Рх2=0.
Рх1=2502, Рх2=0.
Рх1=1169, Рх2=0.
Рх1=590, Рх2=0.
делировать вязкое взаимодействие газа со стенками каналов.
При расчётах были приняты следующие параметры среды: газ - воздух, сжимаемый, вязкий, температура воздуха 25 °С, динамическая вязкость 1,831 •Ю-5 Пахс, начальная плотность 1,18 кг/м3. режим теплообмена - изотермический, атмосферное давление 101325 Па. Под термином давление далее подразумевается статическое давление, являющееся величиной превышения давления над атмосферным.
Некоторые результаты численного расчета параметров потока в каналах струйного двух-входового усилителя при увеличении и снижении давления в канале управления Х1 приведены в таблицах 1 и 2 соответственно.
Аналогичным образом были проведены расчеты по управляющему каналу Х2. На основании полученных результатов построены характеристики переключения двухвходового усилителя для каналов Х1 и Х2 (рис. 3).
в канале Хх до 590 Па и в канале Х2 до 1152 Па. Внешний вид полученных характеристик струйного двухвходового усилителя качественно совпадает с теоретическими характеристиками. Это позволяет утверждать об адекватности разработанной математической модели и возможностях её применения для расчета параметров потока газа в каналах струйных элементов.
ВЫВОДЫ
Разработана математическая модель и методика численного расчёта параметров потока газа в малых каналах струйных элементов с использованием программных продуктов течения газа в каналах струйного двухвходового усилителя.
На основе математического моделирования методами вычислительной газодинамики получены картины течения газа внутри каналов устройства на различных режимах, включая поля скоростей и давлений.
e-t
г 1
Х2 падение давления
1
1
1 1 / /
л /
б)
Рис. 3. График зависимости относительного расхода воздуха в канале У2 от относительного давления в канале питания Х1 (а) и Х2 (б)
Анализ результатов моделирования позволяет утверждать, что рассматриваемый струйный элемент способен функционировать согласно заданной схеме логических состояний при включении/отключении соответствующих управляющих сигналов. При отсутствии давления на управляющих сигналах струя питания на нижней плате усилителя отклонена влево, что позволяет формировать соответствующий управляющий сигнал на верхней плате усилителя для формирования выходного сигнала в канале Уг При повышении давления в управляющем канале Хх до 4815 Па и в канале Х2 до 4384 Па происходит переключение питающей струи в нижней плате усилителя, что влечет за собой переключение выходного сигнала на верхней плате усилителя, создавая расход в канале У2. Обратное переключение двухвходового усилителя происходит при снижении давления
Определены рабочие характеристики струйного двухвходового усилителя, качественно подтверждающие логику работы устройства, что свидетельствует о высокой достоверности применяемой модели течения газа в каналах струйного устройства.
СПИСОК ЛИТЕРАТУРЫ
1. Касимов А.М. Попов А.И. Расходомеры с нулевым перепадом давления и компенсационные измерители скорости // Датчики и системы. М.: Сенси-дат плюс. 2004. № 4. С. 27-31.
2. Ванин В.А. Родина А.А., Колодин А.Н. Управление автоматической линии гальванопокрытий средствами струйной техники // Материалы III международной научно-практической конференции «Виртуальное моделирование, прототипирование и промышленный дизайн». Тамбов: Тамбовский государственный технический университет. 2016, С. 283-287.
3. Многопоточный эжектор и новое направление для развития струйной техники / Ю.А. Сазонов, А.В. Деговцов, Е.С. Казакова, К.И. Клименко // Территория Нефтегаз. М.: ЗАО «Камелот Паблишинг». 2012. № 4. С. 75-77.
4. Лаврентьев В.В., Метляев Д.Д. Исследование возможности использования элементов струйной техники во взрывателях // Молодежный научно-технический вестник. М.: Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет). 2015. № 12. С. 14.
5. Флетчер К. Вычислительные методы в динамике жидкостей: в 2-х. т. Т. 1. М.: Мир, 1991. 504 с.
6. Флетчер К. Вычислительные методы в динамике жидкостей: в 2-х. т. Т. 2. М.: Мир, 1991. 554 с.
7. Андерсон Д., Таннехилл Дж, Плетчер Р. Вычислительная гидромеханика и теплообмен: в 2-х. т. Т. 1. Пер. с англ. М.: Мир, 1990. 384 с.
8. Андерсон Д., Таннехилл Дж, Плетчер Р. Вычислительная гидромеханика и теплообмен: в 2-х. т. Т. 2. Пер. с англ. М.: Мир, 1990. 392 с.
9. Практика применения и особенности современных методов расчёта аэродинамических характеристик летательных аппаратов на основе решений уравнений Навье-Стокса / В.В. Вождаев, Л.Л. Теперин, С.Л. Чернышев // Труды центрального аэрогидродинамического института им. Проф. Н.Е. Жуковского, вып. 2740. М.: издательский отдел ЦАГИ, 2014. 63 с.
10. MenterF.R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications // AIAA Journal, 1994, V. 32, N 8.
11. Ten year of Industrial Experience with the SST-Turbulence Model / Turbulence, Heat and mass transfer 4 / F. Menter, M. Kuntz, R. Lantry. Edited by Hanjalic K., Nagano Y., Tummers M. USA, CT, Redding: Begell House Inc. 2003.
TECHNIQUE OF NUMERICAL CALCULATION OF GASDYNAMIC CHARACTERISTICS OF JET ELEMENTS OF PNEUMOAUTOMATIC EQUIPMENT OF CONTROL SYSTEMS OF DIFFICULT TECHNICAL SYSTEMS
© 2017 V.N. Samsonov, V.G. Shakhov, I.S. Barmanov
Samara National Research University named after Academician S.P. Korolyov
The paper presents the results of modeling of gas flow in the channels of the fluidic elements of pneumoautomatics - a dual input amplifier. The mathematical model of calculation of parameters of gas flow with finite elements method based on the numerical solution of systems of differential equations Navier-Stokes equations, averaged by Reynolds's method. This method is implemented using the software Ansys CFX. To close these equations a turbulence model, are formulated in accordance with the Boussinesq assumption, according to which Reynolds's tensions can be related to the speed of the average strain through the turbulent viscosity. With the help of mathematical model a series of numerical calculations, on the basis of which the characteristics of the jet of a dual input amplifier. Keywords: Management of difficult technical systems, jet equipment, pneumoautomatic equipment, jet amplifier, mathematical model, numerical calculations
Vladimir Samsonov, Doctor of Technical Science, Professor at the Machine Building Design Department. E-mail: [email protected]
Valentin Shakhov, Candidate of Technical Science, Professor at the Aircraft Construction and Design Department. Ildar Barmanov, Candidate of Technics, Associate Professor at the Machine Building Design Department. E-mail: [email protected]
