CC BY
58
моделирования могут служить основой для проведения оптимизации геометрических и теплофизических параметров капиллярного насоса и системы охлаждения в целом.
Библиографический список
1. Двухфазный контур терморегулирования - перспективное направление совершенствования СТР космических аппаратов / С. М. Беднов, П. Д. Вежневец, Ю. М. Лукоянов и др. //Полет. 2003. № 3. С. 37-40.
2. Development of a Two-Phase Capillary Pumped Heat Transport for Spacecraft Central Thermal Bus / T. Hoang,
M. Brown, R. Baldauff, and S. Cummings //Proceeding of2003 Space Technology and Applications International Forum, February 2003. Albuquerque, New Mexico, 2003. P. 49-54.
3. Чи, С. Тепловые трубы. Теория и практика / С. Чи. М. Машиностроение, 1981.
4. Поляев, В. М. Гидродинамика и теплообмен в пористых элементах конструкций летательных аппаратов /
В. М. Поляев, В. А. Майоров, JI. JI. Васильев. М.: Машиностроение. 1988.
5. Патанкар, С. Численные методы решения задач теплообмена и динамики жидкости / С. Патанкар. М.: Энер-шатомиздат, 1984.
E. N. Vasilev, A. A. Dekterev
MATHEMATICAL SIMULATION OFTHE HEAT AND MASS TRANSFER PROCESSES IN THE TWO-PHASE THERMO CONTROL CIRCUITS WITH CAPILLARY PUMP
Heat and mass transfer processes in the two-phase thermo control circuits with capillary pump are considered. Two-dimensional unsteady mathematical model of the hydrodynamic and thermo physics phenomena in the two-phase flow moved in porous media of the capillary pump is based on the numerical solution heat conduction andfiltration equations. Temperature, pressure, velocity, liquid and vapors phase concentrations fields are simulation results.
Keywords: heat transfer, numerical simulation, XML technology, spacecraft.
УДК:534.322
А. С. Орлов
КОНЕЧНО-ЭЛЕМЕНТНОЕ МОДЕЛИРОВАНИЕ ДИФФУЗНОГО ПОЛЯ В АКУСТИЧЕСКИХ РЕВЕРБЕРАЦИОННЫХ КАМЕРАХ БОЛЬШОГО ОБЬЕМА
Рассмотрено построение конечно-элементной модели диффузного поля для акустическихреверберационных камер большого объема. Модель построения с использованием трехмерных акустических элементов и различного вида демпфирования.
Ключевые слова: моделирование, диффузное поле, акустические расчеты, реверберационные камеры.
При проведении наземной экспериментальной отработки космических аппаратов (КА) испытания на акустические воздействия в частотном диапазоне до 4...6 кГц являются общепринятой процедурой. Цель данного вида испытаний состоит в подтверждении прочности конструкции КА и способности аппаратуры выполнять свои функции по время и после акустического воздействия. В то же время это один из самых дорогостоящих видов испытаний. Уменьшить затраты на такие испытания можно за счет предварительного численного моделирования акустического воздействия на КА. В настоящее время сформировалось несколько общепризнанных подходов к такому анализу: замена акустического поля статическим или пульсирующим давлением и решение собственно сопряженной задачи, рассматривающей взаимодействие акустической среды с конструкцией [1]. В работе [2] на примере расчета сотовой панели проведено сравнение этих трех подходов и показано, что результаты, совпадающие с экспериментальными данными, дает только реше-
ние сопряженной задачи «акустическая среда-конструкция».
В акустических расчетах сложных конструкций используются два основных метода: конечно-элементное моделирование (КЭМ) [3] и статистический энергетический анализ (СЭА). Использование КЭМ предъявляет высокие требования к ресурсам (высокопроизводительные компьютеры, подробная модель среды и объекта исследования и т. д.) и наиболее часто применяется для акустического анализа в низкочастотной области. Для решения акустических задач в области средних и высоких частот применяется СЭА в котором прогноз основан на обширной статистической базе, позволяющей оценить реакцию в любом элементе, присутствующем в базе, при прохождении потока энергии, создаваемого акустическим воздействием и, как следствие, в любой системе, состоящей из таких элементов. Этот метод не требователен к вычислительным ресурсам, но чем сложней система, тем больше должна быть база данных, чтобы она смогла обеспе-
чить достоверный результат расчета. Применяется также и комбинированный подход, когда анализ в низкочастотной области проводится с использованием КЭМ, а в области средних и высоких частот - СЭ А.
Постоянно увеличивающиеся производительность и ресурсы используемых компьютерных систем позволяют строить достаточно подробные модели КА и среды. Поэтому основным методом акустического анализа КА становится метод конечно-элементного моделирования. В динамических расчетах с использованием методов КЭМ существует два подхода: расчет с предварительным анализом форм и частот модели и расчет с использованием прямого интегрирования уравнений динамики конструкций и среды по времени [4]. Более приемлемым для рассматриваемых задач является прямое интегрирование уравнений, так как оно не требует линеаризации задачи и учета сотен тысяч форм и частот для сложных систем.
В любом случае первым этапом решения задачи является построение модели диффузного акустического поля, создаваемого в реверберационной камере при испытаниях. Согласно определению, диффузным звуковым полем является поле, все направления потоков энергии волн которого равновероятны, а плотность акустической энергии такого поля по всему объему постоянна [5].
Распространение волн в акустической среде описывается волновыми уравнениями вида
1 5 ы
Дм—-—-L(u) = f(r,t),
С ot
(1)
клика (неявное или явное интегрирование) как результат изменений малого давления, производимого в газе (жидкости). Часто в связанных конструкционно-жидкостных задачах задание давления в жидкости (даже в линейной постановке) достаточно для получения корректного отклика конструкции.
Уравнение равновесия для малых колебаний сжимаемой, адиабатической жидкости, с зависимой от скорости инерцией и демпфированием представлено в пакете АВА()и8 в виде
др , „ , ди/ , „ ч д2и/
= 0,
(2)
где р - избыточное давление в жидкости; г/- смещение частицы жидкости; х - пространственные координаты частицы; с}- плотность жидкости; у - объемное сопротивление (сила на единицу объема и скорости); © и / - переменные, от которых могут зависеть р^и у.
Обычно рассматривается достаточность точности стационарных жидкостных скоростей до 0,1 М.
Основными характеристиками поведения жидкости (газа) считаются вязкость, линейность, сжимаемость и т. д., так что давление в уравнении (2) может быть представлено выражением (3):
р = -КЛх,%.)—uf
f ' дх
(3)
где Дм - оператор Лапласа; с - константа, характеризующая свойство среды; t - время; Ь(и) - дополнительные члены, учитывающие диссипацию и дисперсию сред, в которых распространяются волны; fir, t) -функция, характеризующая внешние воздействия.
Задачи, связанные с распространением волн в линейных диспергирующих и недиспергирующих средах, с отражением и преломлением волн на границах раздела однородных сред, сводятся к решению уравнений типа (1) с различными граничными условиями. Следует заметить, что в нелинейных задачах нарушается принцип суперпозиции и происходит взаимодействие волн различных частот. При этом характер протекания волновых процессов существенно зависит от соотношения дисперсионных и нелинейных свойств этих процессов [13].
Большинство современных пакетов КЭМ и СЭА предлагают решение акустических задач в линейной и нелинейной постановке с широким диапазоном задания граничных условий и воздействий. Одним из таких пакетов является пакет АВ AQU S.
Пакет КЭМ ABAQUS предлагает набор конечных и граничных элементов для моделирования газожидкостной среды, подверженной малым давлениям, и граничные условия для связи этих акустических элементов с конструкцией модели. Например, анализ стационарного (линейного) гармонического отклика может быть выполнен для связанной структурно-акустической системы. Стационарная процедура основана на прямом решении связанных комплексных гармонических уравнений. Конечные элементы акустического поля могут использоваться в анализе процедуры нелинейного от-
где кг модуль всестороннего сжатия жидкости.
Акустическая среда в пакете АВАСЮЗ может иметь зависимый от скорости расход энергии, создаваемый вязкостью жидкости или пористым сопротивлением материала. К тому же основные граничные условия обеспечиваются акустической средой, включающей импеданс либо реакцию границ.
Возможными граничными условиями на поверхности акустической среды являются следующие:
- давление в граничных узлах (давление также может задаваться в любых узлах модели);
- замена сосредоточенной нагрузки в граничных узлах заданием градиентов давления на единицу плотности акустической среды;
- связь меяеду акустической средой и конструкцией, которая осуществляется либо процедурами, основанными на стыковки поверхностей, либо с помощью связывающих элементов А81;
- импеданс среды, характеризуемый поглощающей границей между акустической средой и твердой поверхностью, либо вибрацией конструкции, либо излучением в некую внешнюю бесконечную область;
- давление в единичном объеме акустической среды как результат приходящей волны.
Сопротивление потока и свойства поглощающих поверхностей при анализе стационарного отклика могут быть функциями частоты, но принимаются постоянными в процедуре прямого интегрирования уравнений.
Таким образом, пакет конечно-элементного моделирования АВАСЮЗ позволяет решать широкий класс задач акустики, в том числе и задачи акустического нагружения конструкций.
Для описания распространения волн по конструкции обычно достаточно уравнений плоских волн, для описа-
нияже акустической среды должны использоваться уравнения сферических волн.
Прямое интегрирование волновых уравнений предполагает задание внешних воздействий в виде временной функции. Задание режимов испытаний производится, как правило, в виде зависимости уровня звукового давления, дБ, от частоты. В этом случае для моделирования условий испытаний производится замена частотного представления внешнего воздействия временной зависимостью.
Разобьем исходный спектр нагрузок на части из условия обеспечения четырех точек на октаву по оси частот. Для каждой точки на оси частот определяется уровень давления, после чего формируется функция нагружения, изменяющего во времени, в виде суммы синусов соответствующей амплитуды и частоты.
Пересчет давления выполняется следующим образом [7]:
V
Р = Р0-10®, (4)
где 1\ = 2 ■ 10 * Па: 17— численное значение уровня звукового давления, дБ; Р - искомое значение уровня звукового давления, Па.
Значение фазы выбирается по случайному закону. Затем суммарное значение звукового давления приводится к требуемому (расчетному) уровню давления:
Рх = 1£Р< •8ш(ш,Г + фД (5)
i
где / - суммарное значение звукового давления, Па; Р - уровень давления на соответствующей частоте, Па; ш. - набор частот взят из условия 4 частоты на октаву; ср - фаза, значение которой принимается по случайному закону; ? - время.
Полученное давление воспроизводится локальными источниками звука. Специальные демпфирующие элементы вводятся как на стенки камеры, так и в акустическую среду. Сначала данная технология была апробирована при построении КЭМ акустической камеры небольшого объема (30 м3) [2].
Переход к построению КЭ модели акустической среды в камере большого объема связан с появлением дополнительных трудностей: во-первых, появляются ограничения, связанные с ресурсами компьютерной системы и возникает необходимость менять размеры акустических элементов; во-вторых, необходимо варьирование за счет демпфирования по объему камеры и введения в модель, помимо демпфирующих, еще и отражающих элементов. Но в целом технология построения КЭМ сохраняется.
Рассмотрим далее процедуру построения конечно-эле-ментной модели диффузного поля для акустичсской ревер-берационной камеры РК660, представляющей собой прямоугольный параллелепипед размером 6,94 х 8,7 х 11м. Погрешность задания давления находится в пределах ±3 дБ.
Представим стандартный график уровня звукового давления (суммарныйуровень 142 дБ), задаваемый в акустической камере (рис. 1). Данное воздействие заменялось по формуле (5) временным воздействием.
Конечно-элементарная модель акустической камеры РК660 (рис. 2) состоит из 51 185 узлов, 288 347 элементов (в том числе 1 048 демпфирующих элементов, равномерно распределенных по всему объему камеры, и 1 873 отражающих элементов, установленных рядами на стенках
камеры). Объемный модуль акустической среды принят равным 1.424/'.'5. Согласно определению [8], расчет времени реверберации проводится по анализу затухания амплитуды импульсов различной длительности на 60 дБ.
Октавные уровни
Рис. 1. Уровень звукового давления (суммарный уровень 142 дБ)
Расстояние между узлами конечно-элементной (КЭ) сетки (0,26 м) выбирается исходя из возможности получения корректного результата на существующих вычислительных мощностях. Корректность результата определяется отсутствием трендов в графиках давления, амплитуды которых эквивалентны входному воздействию с погрешностью, не превышающей ± 3 дБ во всех точках акустической камеры. Геометрия узлов и элементов камеры была построена в пакете РАТЯ АМ, а затем экспортирована в файл Сеотейулпр формата АВА(|Щ*. Аналогичную процедуру можно выполнять в пакете N А5ТН.АМ.
Узлы КЭ-сетки, расположенные в углах КЭМ акустической камеры, зафиксированы по шести степеням свободы. Количество и местоположение отражающих и демпфирующих элементов, а также некоторые параметры файлов описания определяются результатами процесса верификации КЭМ акустической камеры.
Рис. 2. КЭМ акустической камеры РК660 (стрелка показывает точку приложения источника акустического воздействия)
Так как верификация - это процесс, исследующий качество численного приближения модели, используемой как основа для предсказания, то при верификации выполняется варьирование различными параметрами модели:
- распределением отражающих элементов;
- распределением демпфирующих элементов;
- определением времени реверберации на различных частотах;
- параметрами акустической среды;
- расположением узлов КЭ-сетки для контроля давления в объеме акустической камеры.
Формализовать процесс верификации практически невозможно. Временные затраты для получения приемлемого результата определяются в первую очередь опытом и квалификацией расчетчика. Следует также добавить, что вопросы валидации в данной статье не рассматриваются, так как уравнения (2), (3) достаточно корректно описывают происходящие в акустической камере процессы, т. е. решение задач верификации фактически включает и процедуру валидации.
Всего было проведено более 30 расчетов и получено свыше 500 графиков давлений по различным контрольным точкам (рис. 3). Помимо данных точек, контроль давления проводится в углах акустической камеры, на поверхностях стенок камеры и выборочно в других (практически произвольных) точках камеры. При этом
Рис. 3. Распределение контрольных точек внутреннего объема камеры в зоне возможной установки объекта испытаний
[хЮ3]
0.20
ей
§
■§ 0.00 ІІ Е 2 LL
-0.20
0,00 0.40 0.80 1.20 1.60 2.00 2.40 2.80
Time
Рис. 5. Давление в центре акустической камеры (завершение верификации)
внутри зоны установки объекта испытаний обеспечивается выполнение требований по уровням звукового давления (252 Па), а в других точках - приблизительное равенство давлений для симметричных точек при выполнении требований по погрешностям, укладывающимся в пределы ±3 дБ. Время реверберации должно соответствовать экспериментальным значениям, полученным в данной камере (например, на частоте 1 000 Гц время реверберации составляет 4,6 с).
На начальном этапе верификации количество демпфирующих и отражающих элементов минимально (один отражающий элемент расположен в центре камеры, а демпфирующие элементы находятся в углах камеры). На графике давления (рис. 4) присутствует достаточно большой положительный тренд, при этом значения давлений существенно отличаются от необходимых значений в 252 Па. На завершающем этапе верификации (рис. 5) тренд отсутствует (график симметричен относительно 0), уровень давления соответствует требуемым значениям. Разброс давлений в контрольных точках укладывается в необходимые пределы (±3 дБ). Корректность разработанной модели также подтверждается графиком (рис. 6). Ве-
Time
Рис. 4. Давление в центре акустической камеры (начальный этап верификации)
1x103]
0.00 0.40 0.80 1.20 1.60 2.00 2.40 2.80
Time
Рис. 6. Давление в центре стенки акустической камеры
личина давления у стенки акустической камеры примерно вдвое выше, чем в центре (эффект наложения падающей и отраженной волн). Расчетное время реверберации на частоте 1 ООО Гц составляет около 4,1с.
Таким образом, разработана модель КЭМ ревербе-рационной акустической камеры большого объема, позволяющая в дальнейшем проводить анализ поведения конструкций при акустическом воздействии.
Библиографический список
1. Felippa, С. A. Partitioned analysis of coupled system/
С. A. Felippa, К. C. Park // Chapter 3 of Computational Methods for Transient Analysis, ed. by T. Belytschko and T. J. R.Hughes. North-Holland, Amsterdam, 1983. P. 157-219.
2. Орлов, А. С. Анализ нагружения сотовой панели акустическими воздействиями / А. С. Орлов // Эффек-
тивность сотовых конструкций в изделиях авиационно-космической техники : материалы II Между на р. науч,-практ. конф. Днепропетровск, 2007. С. 142-150.
3. Бате, К. Численные методы анализа и метод конечных элементов / К. Бате, Е. Вилсон. М.: Стройиздат, 1982.
4. Шимкович, Д. Г. Расчет конструкций в MSC/ NASTRAN for Windows / Д. Г. Шимкович. М.: ДМК-Пресс, 2001.
5. Лепендин, JI. Ф. Акустика / JI. Ф. Лепендин. М. : Высш. ж., 1978.
6. Виноградова, М. Б. Теория волн / М. Б. Виноградова, О. В. Руденко, А. П. Сухоруков. М.: Наука, 1979.
7. Справочник по технической акустике: пер. с нем. / под ред. М. ХеклаиХ. А. Мюллера. JI.: Судостроение, 1980.
8. Лопашев, Д. 3. Методы измерения и нормирование шумовых характеристик / Д. 3. Лопашев, Г. Л. Осипов, Е. Н. Федосеева. М.: Изд-во стандартов, 1983.
A. S. Orlov
FINITE ELEMENT MODELING OF DIFFUSIVE FIELD IN BIG VOLUME ACOUSTIC REVERBERANT CAMERAS
In this paper it is considered creation offinite element model of diffusive field for big volume acoustic reverberant cameras. Model created by usage 3d acoustic elements and different types of damping.
Keywords: modeling, diffusive field, acoustic calculation, reverberation chambers.
УДК539.3
H. В. Молокова
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ И АНАЛИЗ ПРОЦЕССА ЗАГРЯЗНЕНИЯ ПОРИСТОЙ СРЕДЫ
Приведено численное решение математической модели двухфазной фильтрации, учитывающей движение углеводородных загрязнителей и воздуха в пористом грунте. Модель включает систему уравнений в частных производных с дополнительными условиями. В число дифференциальных уравнений входит уравнение баланса массы в элементе пористой среды - уравнение неразрывности, а также дифференциальные уравнения движения. Для замыкания системы вводятся уравнения состояния рассматриваемого загрязнителя и среды. Начальные и граничные условия соответствуют фильтрационному процессу, начиная с поверхности грунта и начальной стадии разлива загрязнителя. Проводится сравнительный анализ результатов математического моделирования с экспериментами.
Ключевые слова: модель двухфазной фильтрации, геофильтрация, математическое моделирование, порис-
тая среда, углеводородный загрязнитель.
Геофильтрационная задача о движении углеводородного загрязнителя и воздуха с учетом гравитационного влияния в трехмерной постановке описывается системой уравнений:
ds
т-------Ь
dt
д(р и х) + д(р иу) + d{puz)
дх
ду
dz
= 0,
5(1-5)
т—-------- +
dt
д(р и х) + д(риу) + d(puz)
дх
ду
dz
(1)
= 0.
^ _____________*■
Щ =~ k— (V^-p^), (2, а)
h
Is ( сЛ _,
и2= ~к——(Ур2~ P2g)- (2,6)
М-2
Начальные условия:
-при/= 0: 5j=52 ,s = s0(x,y),p =р0 наГр
s = 0,p=p наГ,. (3)
’Г Г атм 2 v 7
- начальное давлениерх = рxgH распределено по
