Скачать 1.34 Mb.
|
ГЛАВА 4 РАЗРАБОТКА АНАЛОГО-ЦИФРОВОГО КОМПЛЕКСА ДЛЯ ПАРАЛЛЕЛЬНОГО МОДЕЛИРОВАНИЯ ФИЗИЧЕСКИХ И МАТЕМАТИЧЕСКИХ ПРОЦЕССОВ 4.1 Назначение аналого-цифрового комплекса Одним из самых перспективных методов изучения различных процессов, особенно сложной природы, например, когда гидродинамическая картина исследуемого процесса ощутимо влияет на процессы переноса тепла и/или массы, является системный анализ. При реализации данного метода исследуемый реальный процесс представляется в виде гак называемой физико-химической системы, в состав которой входят элементарные процессы, протекающие одновременно и оказывающие взаимное влияние друг на друга. При системном анализе технологических процессов обычно прибегают к их моделированию. К настоящему времени сложилось два основных вида моделирования, которые находятся в непрерывном взаимодействии и взаимно дополняют друг друга: физическое и математическое. В основе физического моделирования лежит теория подобия и размерности. Оно заключается в экспериментальном исследовании процессов на моделях. При этом исследуется явление не в натуре, а на модели большего или, чаще всего, меньшего масштаба, обычно в специальных лабораторных условиях. Физическое моделирование - единственный надежный способ получения систематических экспериментальных данных о характере и особенностях действительных процессов, которые могут быть положены в основу проектирования новых установок, расчета их характеристик, оптимизации конструктивных и режимных параметров. Особенно это относится к сложнейшим процессам термо-и газодинамики и тепломассообмена, протекающим в энергетических установках. Результаты физического моделирования обычно служат той экспериментальной основой, на которой строится и без которой невозможно математическое моделирование. Двумя недостатками, ограничивающими область применения этого вида моделирования, являются, во-первых, невозможность получения прогнозирующих зависимостей исследуемых процессов, и, во-вторых, его высокая стоимость. Математическое моделирование базируется на знании термодинамических и теплофизических свойств рабочих веществ; уравнений материального и теплового балансов и тепломассопереноса; особенностей действия процессов тепломассообмена; конструктивных особенностей элементов машин и других факторов. Математической моделью реального, физического объекта является математический; объект в виде системы уравнений, дополненной начальными, граничными условиями и описаниями, которые определяют области и диапазоны изменения физических характеристик и технологических параметров, поставленный в соответствие изучаемому физическому объекту, способный замещать оригинал на определенных этапах познавательного процесса и давать при его изучении новую информацию о реальном объекте. Основными преимуществами математического моделирования перед физическим являются следующие: - применение математической модели на этапах замысла и проектирования машин и установок дает возможность заранее определить их параметры и характеристики, что позволяет практически исключить вынужденное применение дорогостоящего метода проб и ошибок, и связанные с этим затраты на реализацию схем и конструкций, которые могут оказаться нерациональными; - математическая модель - чрезвычайно гибкий объект, позволяющий воспроизводить любые как реальные, так и гипотетические ситуации. Поэтому при моделировании появляется возможность имитировать и исследовать особенности работы машин и установок в самых различных условиях, в том числе таких, которые практически невозможно реализовать в опытах, например, сверхвысоких или сверхнизких температур или давлений. Благодаря этому уменьшается потребность в сложном лабораторном оборудовании и в эксплуатационных испытаниях машин; - работу некоторых очень сложных или опасных, с точки зрения их взрывов и пожароопасности или токсичности, систем невозможно проанализировать в необходимом объеме ни с помощью лабораторных, ни с помощью натурных экспериментов. В этом случае применение математической модели представляет собой единственный способ решения этих задач; - математическое моделирование позволяет использовать новейшие методы извлечения информации из экспериментальных данных. Особое значение приобрел в последнее время метод интенсификации, позволяющий оценить ненаблюдаемые в процессе проведения опытов параметры или выбрать наилучшую структуру математического описания зависимости между отдельными факторами. Математические исследования вскрывают природу явлений и тем самым позволяют установить количественные и качественные зависимости изучаемых процессов, но математическое моделирование невозможно без знания уточненных параметров переноса и их зависимости от всех входящих факторов. В практике даже ведущих мировых производителей теплообменной техники, как правило, используется комбинированный способ моделирования, при котором отдельные данные, полученные путем математического моделирования, проверяются с помощью эксперимента и наоборот. Для реализации приведённых рекомендаций предлагается использование аналого-цифрового комплекса [34] Задачей, на решение которой направлен АЦК, является создании эффективного способа моделирования процесса газификации остатков жидкого ракетного топлива в баках отделяющейся части ступени РКН и устройства на основе физической и математической моделей, путем повышения точности, достоверности, верификации и валидации получаемых результатов [25,35] Основной проблемой при разработке АЦК являются: - определение структуры АЦК; - выбор сложности математической и физической моделей; - организация функционирования (внутренних связей) внутри АЦК, в том числе и синхронизация процессов математического и физического моделирования; - разработка критериев для оценки эффективности моделируемого процесса (влияние параметров процесса на конечные характеристики, например влияния состава газовой смеси на удельный импульс и т.д.). Возможность использования моделей различной сложности: - различные физические модели, например, элементы емкости бака (ЭММУ) и топливный отсек реального бака, модельный компонент (вода) и реальный компонент (керосин); - различные математические модели, например, уравнение Навье - Стокса и метод теплового баланса; - оценка на ранних этапах исследования возможных направлений конструирования бортовых систем, например, системы газификации остатков топлива в баках [36] Указанная задача решается тем, что процесс физического и математического моделирования проводят одновременно, при этом, часть параметров исследуемого процесса, например, давление, температуру берут из физической модели процесса и используют в математической модели, синхронность процессов физического и математического моделирования обеспечивают за счёт введения «математического» времени, т.е. введением задержки в процесс математического моделирования, обеспечивающей синхронность, используя математическую модель процесса, осуществляют пересчёт параметров и рассчитывают различные критерии и показатели, например, текущие тепловые потери, прогнозируемое время завершения процесса газификации и на основании этих критериев и показателей и изменяют параметры физического процесса, например, температуру теплоносителя, массовый секундный расход. Конкретные виды физической и математической моделей определяют из условия совместимости по исследуемым параметрам, скоростям протекания процессов, например, использование уравнения Навье - Стокса, решаемого в пакете ANSYS требует на 1 секунду реального процесса в среднем 300-600 секунд вычислительного времени. Таким образом, использование таких сложных математических моделей в составе АЦК нереально. Рисунок 4.1- принципиальная схема АЦК: 1- компрессор, 2- теплоэлектрического нагревателя, 3- дроссель, 4,5,6,8,9,10 – датчики, 7- ЭММУ, 11- контроллер, 12,13 – АЦП с интегрированной системой верификации данных, 14 – управляющий комплекс, 15 – исходные данные, 16 – математическая модель, 17 – блок обработки данных, 18 сервер совместного расчета параметров и критериев. Последовательность работы АЦК (рисунок 4.1): -запуск стенда (запуск компрессора 1 или забор теплоносителя из ресивера, теплоэлектрического нагревателя 2, открытие дросселя 3 максимальное, для наиболее быстрого прогрева системы); - до выхода на рабочий режим (достижения заданной температуры на выходе из теплоэлектрического нагревателя 2 и выходе из экспериментальной емкости) осуществляется измерение только температуры внутри экспериментальной емкости 7; - после выхода стенда на рабочий режим,датчики давления 4,10, температуры 5, 8 (вход и выход экспериментальной емкости 7), расхода 6 на входе в экспериментальную емкость 7 и процентного состава газа 9, подают данные на контроллер11. Контроллер 11 обрабатывает данные и посылает их на аналого-цифровой преобразователь с интегрированной системой верификации данных 12 для перевода показаний датчиков 4-6, 8-10в цифровой. Управляющий комплекс13, совмещенный с хранилищем базы данных (показание всех датчиков), накапливает получаемые данные и задает различные варианты рабочего давления подачи теплоносителя через дроссель3 в магистрали давления, температуру теплоэлектрического нагревателя 2, согласно заранее заданной программе проведения эксперимента 14 (циклограммы секундного расхода теплоносителя и его температура). Управляющий комплекс13 контролирует поступающие измерения с датчиков 4-6, 8-10. После проверки достоверности поступающих данных из датчиков управляющий комплекс 13 привязывает данных ко времени и структурирование в виде таблиц. Таблицы полученных параметров поступают в блок совместного расчета параметров и критериев для обработки сервером совместного расчета параметров и критериев17 [37] Для синхронизации времени физического и математического экспериментов вводиться дополнительные циклы в программу расчета математической модели (для выбранных математических моделей процесс моделирования идет в несколько раз быстрее, чем физический процесс). В процессе работы сервер 17, с определенным временным ожиданием, определяемым в процессе эксперимента, и зависимым от сложности математической модели, передает данные в блок обработки данных 16, который используя математическую модель процесса 15, рассчитывает параметры и критерии, заданные в процессе написания программы обработки данных и математической модели. Данные вносятся в базу данных для сравнительного анализа и моделирование изменений текущих условий протекания физического процесса (изменение параметров во время эксперимента, например температура теплоносителя и его массовый расход) [38] 4.2 Математические модели АЦК Математическая модель процесса испарения в общем виде одинакова как для ЭМУ №1, ЭМУ №2, так и для топливного бака, однако наличие различных устройств, их количества приводит к необходимости конкретного описания каждого объекта. На рисунке 4.2 приведена общая схема взаимодействия тепловых потоков физических объектов (ТН, стенки бака, газ наддува, испаряемая жидкость), участвующих в процессе теплообмена внутри ёмкости (бака), которая иллюстрирует достаточно сложный механизм теплообмена. Рисунок 4.2 – Тепловые потоки и их составляющие (конвективные, кондуктивные, излучение) между участниками взаимодействия На рисунке 4.3 приведено фото ЭМУ№1, а на рисунке 4.4 приведено фото ЭМУ №2 Рисунок 4.3 - Экспериментальная модельной ёмкости №1 для исследования процессов испарения модельных жидкостей (вода, спиртовая смесь) и керосина при параметрах ТН: избыточное давление в ёмкости 0 – 1,5 атм., массовый секундный расход 100 – 500 л/мин.; углы ввода теплоносителя 00 , 300 , 450; температура теплоносителя 1000 С – 1500 С; скорость входа ТН ….м/с; состав теплоносителя – воздух; состав газовой подушки – воздух. Рисунок 4.4 - Экспериментальная модельной ёмкости №2 в вакуумной камере для исследования процессов испарения модельных жидкостей (вода, спиртовая смесь) и керосина в вакуумную камеру. Диапазоны давления в вакуумной камере 0,1 – 1,0 атм. Производительность форвакуумного насоса для откачки газифицированной смеси до 500 л/мин. Абсолютное давление в ёмкости 0,1 – 3,5 атм., массовый секундный расход 100 – 500 л/мин.; углы ввода теплоносителя 00 , 300 ; температура теплоносителя 1000 С – 1500 С; скорость входа ТН ….м/с. Состав теплоносителя – воздух, азот; состав газовой подушки – воздух, гелий. Исходная система уравнений для исследования процесса испарения жидкости взята в соответствии с работой [14], при этом она преобразована с учётом использования не химического механизма получения тепла, а за счёт подачи горячих газов, т.е. термодинамического источника тепла. Изменение массы продуктов испарения в объеме бака можно записать в следующем виде: , (4.1) где: - массовый расход ТН; - массовая скорость испарения жидкости; где - коэффициент массоотдачи определяется экспериментально [24]; - парциальное давление испарившейся жидкости у поверхности жидкости, принимаемое равным давлению насыщенных паров при температуре поверхности (имеются табличные значения[24]); - текущее парциальное давление испарившейся жидкости в свободном объёме бака; - площадь испарения жидкости, является граничным условием, например, «зеркало»; - массовый расход жидкости уходящей через дренажное отверстие, где - коэффициент расхода; - площадь выходного сечения дренажного отверстия; - давление в смеси газов в баке (газ наддува испарившаяся жидкость и ТН); - температура газовой смеси в баке; - газодинамическая функция. Начальные условия в τ=0 , где - масса воздуха в баке. Уравнение изменения плотности газа в емкости имеет вид: , (4.2) где - объем бака. Начальное условие - плотность воздуха. Уравнение изменения давления в газовой подушке в соответствии с [14]: , (4.3) где - суммарное поступление тепловой энергии в емкость; - энтальпия ТН; - энтальпия смеси газов (ТН, воздух, испарившаяся жидкость) на выходе из бака; - адиабатическая постоянная. Температура газовой подушки в емкости определяется по формуле Менделеева –Клапейрона [39]: (4.4) где R – газовая постоянная смеси газов. Распределении тепловой энергии (уравнение теплового баланса) имеет вид: , (4.5) где - энергия расходующаяся на нагрев газа, - конвективная составляющая энергии, где - коэффициент теплоотдачи ТН, является функцией от температуры и давления (); - площадь стенок емкости, в которую заключен газ; - температура ТН - температура газа; - лучистая составляющая энергии, где - постоянная Стефана-Больцмана - степень черноты горячего воздуха, в перспективе предполагается использовать зависимости от различных газов; Температура газа находится из уравнения[39]: , (4.6) где - теплоёмкость воздуха; - масса воздуха в емкости. где - конвективная составляющая теплообмена между жидкостью и теплоносителем; где - площадь поверхности жидкости; - температура теплоносителя; - температура жидкости; - лучистая составляющая теплообмена между жидкостью и теплоносителем; где - степень черноты жидкости; - конвективная составляющая теплообмена между жидкостью и газом; где - коэффициент теплоотдачи газа, находится экспериментально; - лучистая составляющая теплообмена между жидкостью и газом; - конвективная составляющая теплообмена между жидкостью и пластиной; где - температура пластины; - коэффициент теплоотдачи жидкости, находится экспериментально; - лучистая составляющая теплообмена между жидкостью и пластиной; Температуру жидкости находим из уравнения: , (4.7) где - теплоёмкость жидкости; - масса жидкости. где - конвективная составляющая теплообмена между пластиной и газом; где - площадь поверхности пластины; - лучистая составляющая теплообмена между пластиной и газом; где - степень черноты пластины; - конвективная составляющая теплообмена между жидкостью и пластиной; - лучистая составляющая теплообмена между жидкостью и пластиной; - конвективная составляющая теплообмена между пластиной и теплоносителем; - лучистая составляющая теплообмена между пластиной и теплоносителем; Температуру пластины находим из уравнения: , (4.8) где - теплоёмкость пластины; - масса пластины. где - конвективная составляющая теплообмена между стеклянной поверхностью стенки и газом; где - площадь стеклянной поверхности стенки; - температура стеклянной поверхности стенки; - лучистая составляющая теплообмена между стеклянной поверхностью стенки и газом; где - степень черноты стеклянной поверхности стенки; - конвективная составляющая теплообмена между стеклянной поверхностью стенки и теплоносителя; - лучистая составляющая теплообмена между стальной поверхностью стенки и теплоносителем; - конвективная составляющая теплообмена между стальной поверхностью стенки и газом; где - площадь стальной поверхности стенки; - температура стальной поверхности стенки; - лучистая составляющая теплообмена между стальной поверхностью стенки и газом; где - степень черноты стальной поверхности стенки; - конвективная составляющая теплообмена между стальной поверхностью стенки и теплоносителя; - лучистая составляющая теплообмена между стальной поверхностью стенки и теплоносителем; - лучистая составляющая теплообмена с стальной поверхности стенки с внешней средой. - лучистая составляющая теплообмена с стеклянной стенки с внешней средой. Температура стеклянной стенки находим из уравнения: , (4.9) где - теплоёмкость стеклянной стенки; - масса стеклянной стенки. Температура стальной стенки находим из уравнения: , (4.10) где - теплоёмкость стельной стенки; - масса стальной пластины. Начальные условия для жидкости принимаются в виде «зеркала». Значения теплофизических констант для ЭМУ принимаются в следующими: - температура теплоносителя 373-423 К; - массовый секундный расход 100-300 л/мин; - температура в ЭММУ 293°К; - масса газифицируемой жидкости 50 гр. На рисунке 4.5 приведена блок-схема алгоритма расчёта параметров процесса испарения жидкости. Рисунок 4.5. Блок-схема алгоритма расчета параметров процесса испарения жидкости. В приведённой системе уравнений проблемными вопросами является определение коэффициентов теплоотдачи, которые являются функциями конкретных условий взаимодействия и аналитическое определение их невозможно. 4.3 Физические модели АЦК Физическая составляющая АЦК представляет собой набор измерительно-контрольной аппаратуры. Таблица 4.1 Датчики контроля и регистрации
Рисунок 4.6 - Схема подключения датчиков и управления пневмосистемой Программное обеспечение для АЦК используется как стороннее так уникальные программные продукты разработанные специально для данной системы. В настоящее время осуществляется тонкая настройка и отработка алгоритмов работы АЦК на виртуальной модели, на которой сигналы измерительной аппаратуры генерируются математической моделью. Целью данного этапа является отработка синхронности процесса расчета математической модели и процесса протекания физического моделирования. 4.4 Выводы по главе 4 На основании проведённых в главе исследований получены следующие основные результаты 1. Рассмотрен состав и назначение аналого-цифрового комплекса для параллельного математического моделирования процессов, исследуемых на стенде. Показаны возможности АЦК для решения задач валидации и верификации исследуемых на стенде процессов. 2. Приведена математическая модель для исследования процесса тепло-и массообмена, алгоритм её интегрирования, определены возможные физические измеряемые параметры (температура, давление, массовый секундный расход теплоносителя) и места их использования в алгоритме интегрирования для использования в математическом моделировании. Введены критерии, расчёт которых подтверждает (или опровергает) достоверность получаемых результатов. 3. Представлена физическая модель АЦК в составе экспериментального стенда. Сформированы основные требования к функционированию АЦК, в частности, разработки алгоритма синхронизация математического и физического времени. 4. Показана возможность повышения достоверности и точности получаемых результатов математического моделирования за счёт использования в математическом моделировании физических измерений, учитывающих нелинейные эффекты, которые затруднено описать в упрощённой математической модели, например, углы ввода теплоносителя, введение ультразвукового акустического взаимодействия с различной ориентацией зон воздействия и т.д. |
Федеральное государственное бюджетное образовательное учреждение... Министерство образования и науки российской федерации федеральное государственное бюджетное образовательное учреждение высшего профессионального... |
Федеральное государственное бюджетное образовательное учреждение... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
||
Федеральное государственное бюджетное образовательное учреждение... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
Федеральное государственное бюджетное образовательное учреждение... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
||
Федеральное государственное бюджетное образовательное учреждение... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
Отчет о самообследовании деятельности федерального государственного... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
||
Отчет о самообследовании деятельности федерального государственного... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
Федеральное агентство воздушного транспорта Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Московский государственный... |
||
«донской государственный технический университет» (дгту) Кафедра «Иностранные языки» Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
Программ а стратегического развития Федерального государственного... Полное наименование вуза – федеральное государственное бюджетное образовательное учреждение высшего профессионального образования... |
||
Федеральное государственное образовательное учреждение высшего профессионального... Федеральное государственное образовательное учреждение высшего профессионального образования |
Рекомендации по приведению в соответствие мероприятий и планов реализации Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Нижегородский государственный... |
||
Федеральное государственное бюджетное образовательное учреждение... Сущность и этапы организационно управленческого анализа (на примере ОАО "вымпелком") |
Самарский государственный технический университет Федеральное государственное бюджетное образовательное учреждение высшего образования |
||
Федеральное государственное бюджетное образовательное учреждение... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
Федеральное государственное бюджетное образовательное учреждение... Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования |
Поиск |