Исследование комбинированной Френкеля-Тонкса и конвективной неустойчивости Текст научной статьи по специальности «Физика»

Алексей Юрьевич Чирков родился в 1976 г, окончил МГТУ им. Н.Э. Баумана в 2000 г. Канд. техн. наук, доцент кафедры "Теплофизика" МГТУ им. Н.Э. Баумана. Автор около 50 научных работ в области физики плазмы.

A.Yu. Chirkov (b. 1976) graduated from the Bauman Moscow State Technical University in 2000. Ph. D. (Eng.), assoc. professor of "Thermal Physics" department of the Bauman Moscow State Technical University. Author of about 50 publications in the field of plasma physics.

УДК 532.59:537.29

И. Н. Алиев, С. О. Ю р ч е н к о, Е. В. Назарова

ИССЛЕДОВАНИЕ КОМБИНИРОВАННОЙ ФРЕНКЕЛЯ-ТОНКСА И КОНВЕКТИВНОЙ НЕУСТОЙЧИВОСТИ

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

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

©

Первая задача — это случай однородного вертикального электрического поля в сочетании с "подогреваемым дном". В такой задаче представляет интерес взаимное влияние неустойчивости заряженной поверхности (неустойчивость Френкеля-Тонкса) [1] и конвективной неустойчивости [2].

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

Вторая задача — это сочетание подогрева плоского слоя жидкости сверху с вертикальным электрическим полем. Здесь может иметь место существенное изменение картины волновых движений в отличие от аналогичного случая в отсутствие поля, так как при подогреве сверху все нормальные возмущения затухают со временем, но система имеет колебательные моды собственных движений [2], как в случае волн на поверхности раздела двух жидкостей различной плотности в поле сил тяжести, когда верхняя жидкость движется с постоянной фоновой скоростью (неустойчивость Кельвина-Гельмгольца) [3]. Представляет интерес возможность управления началом режима неустойчивости Френкеля-Тонкса при помощи вертикального нагрева сверху, хотя вопрос о взаимодействии факторов, приводящих к монотонной и колебательной неустойчивости, в данном случае остается открытым.

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

При указанных допущениях справедлива линеаризованная система уравнений, состоящая из уравнений Навье-Стокса, непрерывности, теплопереноса в приближении Буссинеска-Обербека [3], уравнения Пуассона [1] для электрического потенциала с соответствующими

граничными условиями:

г > к : Дф = 0;

ди V® к > г > 0: — = - двТ + "Ди; д£ р

дТ

— + иУТ = хДТ, и = 0; г = 0 : иж = = 0, Т = 0; (1)

г = к : ф - 4па£ = 0, и - =0, Т - = 0; е , \ с \ е о ди* , дих п

-рд£ + р + рс + рЕ = 2П"др "дХ + ~д7 = 0;

г ^ : Уф = -4паег,

где ф — потенциал электрического поля, и — скорость точек жидкости; р, Т — отклонения значений давления и температуры от равновесных, р — плотность жидкости; в — температурный коэффициент линейного расширения жидкости; V, ц — кинематическая и динамическая вязкость; £ — отклонение точек поверхности от равновесного положения;

X — температуропроводность; а — поверхностная плотность электри-

с е

ческого заряда; р , рЕ — капиллярное и электростатическое давление; ег — орт в вертикальном положении.

Капиллярное давление над искривленной поверхностью жидкости [3]

1 1 \ л , Л д2 д2

д- + Д) - 7Д-£, Д1 = + д?,

где Я1, Д2 — радиусы кривизны поверхности; 7 - коэффициент по-

1 2 д2

верхностного натяжения; в случае плоской волны лапласиан Д1 = ——.

дх2

Для синусоидальной волны потенциал электрического поля и электростатическое давление с учетом граничных условий имеют вид [1]

ф = 4па £ exp [—k (z — h)] — 4na (z — h);

= —2na2 — 4na2k£,

- — 8П (уФ)

г=Н

где к — волновое число.

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

/ д \ ди

[ р— - пД - 2п ДП = (рд - ^ Д- - 4па2^ Д-£.

В случае малых отклонений от равновесных значений нелинейными по скорости слагаемыми в уравнениях (1) можно пренебречь, после чего по аналогии с [2] выбрать единицы измерения:

[x] = [z] = h; [t] = h2/v; [и] = x/h; [p] = pvx/h2; [T] = Ah.

Для элементарных волн возмущения скорости, температуры, давления зависят от времени и горизонтальной координаты как exp [i (kx — ut)], где ш — комплексная частота, k — волновое число. Тогда в безразмерных переменных математическая формулировка уравнений задачи и граничных условий имеют вид

du _ ^^ дТ

— = —Vp + Au + Ra Тег; Pr — — uez = ДТ; divu = 0;

dt dt

z = 0 : ux = uz = 0, T = 0; dux

z = 1: —- — ikuz = 0, uz + iu Pr T = 0;

dz (2)

(-гш + 3k2)

2 д3 Uz

dz dz3

=-k2' 1 +

„ gßAh4 ^ v vy pgh2 ^ 4na2

Ra = ^-, Pr =-, ^ = —^, W = , F =

их ' х' д^3' т '

где Яа — число Рэлея, Рг — число Прандтля.

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

их = и (г) ехр [г (кх — ш£)]; р = Р (г) ехр [г (кх — ш£)];

= V (г) ехр [г (кх — ш£)]; Т = 0 (г) ехр [г (кх — ш£)], где и, V, Р, 0 — амплитудные функции. Отсюда и" — (к2 — гш) и — гкР = 0; V" — (к2 — гш) V — Р' + Яа 0 = 0; 0" — (к2 — гш Рг) 0 + V = 0; V' + гки = 0.

Будем искать систему линейно независимых решений уравнений системы (3) в виде

и (г), V (г), Р (г) ,0 (г) - ехр[(г].

Общее решение системы (3) для каждой функции есть линейная комбинация линейно независимых решений. Подстановка приведен-

(3)

ных решений в систему (3) и требование нетривиальности общего решения приводят к спектральному уравнению вида

с2 - к2 + гш 0 -гк 0

0 С2 - к2 + гш -С Яа

0 1 0 С2 - к2 + гш Рг

ik

С

о

о

= 0.

Алгебраическое уравнение, определяющее спектр значений £ получается после раскрытия детерминанта и элементарных преобразований:

(С2 - к2) (С2 - к2 + гш) (С2 - к2 + гш Рг) + Яа к2 = 0. (4)

Уравнение (4) — бикубическое относительно показателя экспоненты, т.е. спектр значений £ представляет собой три пары комплексно-сопряженных чисел, а искомые амплитудные функции являются линейными комбинациями гиперболических синусов и косинусов.

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

Если ^ — корни спектрального уравнения (4), то скорости и отклонения температур можно искать в виде

их = Zj А ехр г] ехр [г (кх - ш£)]; = -гkAj ехр г] ехр [г (кх - ш£)];

гкА, ехр [С,-г] г./7 в = .2 Зи2_р р ехр [г (кх - ш*)], с/ - к2 + гш Рг

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

Подстановка в граничные условия приводит к однородной системе линейных уравнений относительно констант:

С Aj = 0; ^ Aj = 0;

Aj

j - k2 + гш Pr

= 0;

(с2 + k2) Aj exp [Zj] = 0;

Z2 - k2

Z2 - k2 + гш Pr

Aj exP [Cj] = 0;

(

k2 1

К, [с,2 + (гш - 3k2)] +

F ,

(5)

k +

VW w )

С2 - k2 + гш Pr

V

Aj exp [Zj] = 0.

/

z = 1 : = „ * = 0.

Требование нетривиальности системы линейных уравнений относительно констант Aj дает трансцендентное дисперсионное уравнение c определителем шестого порядка.

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

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

д 2и* dz2

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

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

— фактор поверхностного заряда (причина неустойчивости Френ-келя-Тонкса) — монотонный, так как выйдя из положения равновесия, поверхность жидкости отклоняется вплоть до проявления нелинейных эффектов движения или даже до отрыва капли (экспериментально наблюдаемое появление "конусов Тейлора");

— фактор динамической раскачки верхней движущейся жидкостью (причина неустойчивости Кельвина-Гельмгольца) — колебательный.

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

Монотонные факторы должны приводить к реализации монотонного возмущения, а следовательно, условием критичности числа Рэлея будет

Ra = Rac ^ гш = 0.

Тогда условие критичности числа Рэлея может быть определено с учетом (4) и (5) как условие нетривиальности решения системы

линейных алгебраических уравнений

ZjA = 0; £ Aj = 0;

Aj

j - k2

Aj exp [Zj] = 0;

/

Kj [Zj-3k2] +

= 0; (Zj2 + k2) Aj exp [Zj] = 0;

w) 1

Aj exp [Zj] = 0;

k2 1-

Vwk+wj

Zj2 j

k2

/

Zi,2 = k2 - Ra1/3 k2/3; Z|,4,5,a = k2 + 1 (1 ± ¿V3) Ra1/3 k2/3, j = 1... 6.

Классическая неустойчивость Френкеля-Тонкса наступает при значении параметра ^ = 2.

На рис. 1 и 2 приведены кривые зависимости критического числа Рэлея от волнового числа для случая отсутствия электрического поля, что совпадает с рассмотренным в работе [2], а также в случае докри-тического и критического поля с параметрами ^ =1 и ^ = 2. Из проведенного численного расчета следует, что электрическое поле снижает значения критических чисел Рэлея. При достижении критического значения напряженности электрического поля на волнах с к =

Rac F=0

0,5 1,0 1,5 2,0 2,5 3,0

Rac F= 1

0,5 1,0 1,5 2,0 2,5 3,0

Рис. 1. Графики критических значений чисел Рэлея в отсутствие поля (Р = 0) и в случае докритического поля (Р = 1)

Rac F=2

1200 п--т-,—

I ■ ... I ■ ... I ■ ... I к 0,5 1,0 1,5 2,0 2,5 3,0

Рис. 2. Графики критических чисел Рэлея в случае критического поля (Р = 2)

начинается неустойчивость Френкеля-Тонкса, а волны с близкими значениями волнового числа, устойчивые в отсутствие подогрева, могут возбуждаться даже при малых значениях чисел Рэлея.

Анализ спектра амплитудных функций возмущений. Как уже было отмечено, спектральное уравнение (4) — кубическое относительно £2, а следовательно, всегда имеет точное решение в радикалах:

,0 гш (Рг +1) . _ „

С2 = к2--^3-^ + Жг, г = 1, 2, 3;

X =

(Vp3 + о2 - о) / - p

\J'p3 + q2 - q

1 p - (Vp3 + o2 - o) 7 г V3 p + (Vp3 + q2 - o) 7 (6)

X2,3 =--! =— ±--! =—;

у Vp3 + o2 - o у Vp3 + o2 - o

1Pr3 +1 2

p = 9pr+Tш;

Rak2 гшз _ iWti / 1 o = — - г^( Pr +1)( Pr - 2) Pr - -

Выражениями (6) определяется спектр собственных амплитудных функций линеаризованных уравнений в приближении Буссинеска. Оно может быть переписано в виде

у3 + гш (Pr + 1) у2 - ш2 Pry + Ra k2 = 0,

(7)

У = Z2 - k2.

Пусть у!, y2, y3 — корни уравнения (4), тогда справедливы соотношения

У1 + У2 + Уз = -s (Pr +

У1У2 + У1У3 + У2У3 = s2 Pr;

2 (8)

У1У2У3 = Ra k2;

s = гш.

Уравнения (8) позволяют сделать ряд важных выводов о спектре значений Z. Например, требование, чтобы значения Z были действительными (или одно решение Z2 действительное, а два других — комплексно-сопряженные) для любых волн, приводит к выводу о том, что величина s = гш является действительной (известный результат для случая подогреваемой снизу жидкости, т. е. положительных значений числа Рэлея). Возможные значения Z-параметра и их связь с видом амплитудных функций приведены в таблице.

Таблица

Z2 s = iw Vi Вид амплитудной функции

Z2 > 0 s е R, s < 0, Ra > 0, У > 0 Монотонная, экспоненциально затухающая или возрастающая по глубине

Z2 < о s е R, s > 0, Ra < 0, У < 0 Синусоидальная, не затухающая по глубине

Z2 е с s е с, Ra < 0 y е с Немонотонная, экспоненциально затухающая или возрастающая по амплитуде синусоида

Примечание. К — множество действительных чисел; С — множество комплексных чисел.

Перепишем выражения (6) с использованием s-параметра

С2 = - s ( Pr +1)

1 Pr3 + 1

3

+ Xi, i = 1, 2, 3;

P =

2

9 Pr + 1

(9)

Ra k2 s3 , . , .

q = ~T + 27 (Pr +1)( Pr - 2)

1

Pr - -

-2

Действительное значение первого корня достигается при р3 + + д2 ^ 0, что равносильно двум неравенствам:

Rak2

Rak2

2 s3

9 Pr + 1

1 Pr3 + 1 ^ (Pr +1)( Pr - 2)(Pr - 2

£ - -

1 Pr3 + 1

9 Pr + 1

3/2

27

(Pr + 1) (Pr - 2)

Pr - -2

(10)

27

Второй и третий корни действительны при

Rak2

1 Pr3 + 1 = 9 Pr +1

3/2

(Pr + 1) (Pr - 2)

Pr - -

-2

27

(11)

Пара кривых (11) приведена на рис.3. В точках, принадлежащих области р3 + д2 < 0, все три решения комплексные; кривые соответствуют случаю трех действительных корней, причем два из них совпадают, т.е. имеет место вырождение шести линейно независимых амплитудных функций в четыре. Остальные точки отвечают случаю одного действительного и двух комплексно-сопряженных корней.

Ra к2

2s3

0,2

0,1

-0,1

-0,2

рЧ, :f >0

1 1 1 1 р: '+<?2< ..... :0

3,5 \ 2 : 1 1 1 1 1 1,5

р'+с f> 0 \

к

Pr

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

Следует подчеркнуть, что явление бифуркации в задаче о распространении волн, которое выражалось в появлении дополнительных решений в линейной

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

Если действительные значения принимает первая пара решений , тогда решения уравнения

Рис.3. Кривые вырождения и области чисто мнимых, действительных и комплексных решений

Z? = k2 -

s (Pr +1) 2p

3

+

1/3

q

= 0

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

Ra

2k4

9 Pr +1 2 Pr3 + 1

'k?

s

Pr + 1'

3

3

(Pr + 1) (Pr

2)( Pr - 2

27

k 2)

3

(12)

/

Поверхность, определяемая уравнением (12), приведена на рис.4. Для точек, расположенных ниже поверхности (12), решение положительно, что соответствует монотонным функциям профиля, точки над поверхностью (12) определяют функции профиля в синусоидальном виде. При этом действительная часть пары комплексно-сопряженных решений будет иметь соответственно положительное значение над поверхностью и отрицательное — под поверхностью.

Рис. 4. Поверхность =0

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

а2 ;2 в (Рг +1) 2р (2 = к2--^-- +

3

q

1/3'

Z2 =0 - k2 = 3 ( Pr +1 ± 2 ^

,2 _,2 « (РГ +1) Р

Z2;3 = k

,1/3 :

(13)

k2 1

Z2;3 = 0 - - = 3 | РГ +1 ±

Pr3 + 1 Pr + 1 1 '

q = ±s3

1 Pr3 + 1\

3/2

Рг + 1 ) '

На рис.5 приведены кривые, определяемые формулами (13) при равенстве значений С2 нулю, обозначены области положительных и отрицательных значений (2. Ответить на вопрос, каким волнам соответствует определенный амплитудный профиль по глубине, можно имея дисперсионное соотношение, а также зная значение числа Прандтля.

3

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

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

при больших отрицательных числах Рэлея такой механизм не обнаружен, что связано с постоянством плотности в выражении для безразмерного Ш-параметра. Это требует отдельного рассмотрения. Известный вопрос о взаимодействии монотонных и колебательных мод при этом остается открытым.

СПИСОК ЛИТЕРАТУРЫ

1. Л а н д а у Л. Д., Л и ф ш и ц Е. М. Электродинамика сплошных сред. - М.: Физматлит, 2003. - 656 с.

2. Гершуни Г. З., Жуховицкий Е. М. Конвективная неустойчивость несжимаемой жидкости. - М.: Наука, 1972. - 392 с.

3.Ландау Л. Д., Л и ф ш и ц Е. М. Гидродинамика. - М.: Физматлит, 2003. -736 с.

4. А н т о н ю к П. Н. Дисперсионные уравнения для плоской капиллярно-гравитационной волны на свободной поверхности вязкой несжимаемой жидкости // ДАН СССР. - 1986. - Т. 286, № 6. - С. 1324-1328.

3

2,5 2

1,5 1

0,5 0

-0,5 3

2,5 2

1,5 1

0,5 0

-0,5

ръ + q2 >0

$?>0

^>0

'5' 1 i ц s" : 1 ---- I ,5 3

к2_ s

p3+q2> 0

Рг

^2,3 >0

^2,3 <0

^2,3 >0

Рг

12 3 4

Рис. 5. Области положительных и отри-

5. А л и е в И. Н., Ф и л л и п о в А. В. О волнах, распространяющихся по плоской поверхности вязкой проводящей жидкости в электрическом поле // Магнитная гидродинамика. - 1989. - № 4. - С. 94-98.

Статья поступила в редакцию 29.10.2007 Исмаил Новрузович Алиев родился в 1945 г., окончил в 1969 г. Московский инженерно-физический институт (МИФИ). Д-р физ.-мат. наук, профессор кафедры "Физика" МГТУ им. Н.Э. Баумана. Автор более 60 научных работ в различных областях физики.

I.N. Aliev (b. 1945) graduated from the Moscow Engineering and Physical Institute in 1969. D. Sc. (Phys.-Math.), professor of "Physics" department of the Bauman Moscow State Technical University. Author of over 60 publications in the various fields of physics.

С.О. Юрченко родился в 1985г. — студент МГТУ им. Н.Э.Баумана. Автор более 10 научных работ.

S.O. Yurchenko (b. 1985) — student of the Bauman Moscow State Technical University. Author of more than 10 publications.

Екатерина Викторовна Назарова родилась в 1982 г., окончила МГТУ им. Н.Э. Баумана в 2005 г. Соискатель кафедры физики МГТУ им. Н.Э. Баумана.

Ye.V. Nazarova (b. 1982) graduated from the Bauman Moscow State Technical University in 2005. Competitor for a scientific degree in "Physics" department of the Bauman Moscow State Technical University.

В издательстве МГТУ им. Н.Э. Баумана в 2007 г. вышла в свет книга

Недашковский В.М. Хартов В.Я.

Микроконтроллеры AVR. Практикум для начинающих. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2007. - 240 с.

Практикум содержит материалы для изучения микроконтроллеров AVR с архитектурой RISC. Рассмотрены необходимые инструментальные средства — интегрированный пакет для разработки и отладки программ AVR Studio 4, стартовый набор разработчика STK500. Предложен комплект программ для изучения функциональных возможностей микроконтроллеров (27 учебных проектов). Тематика охватывает практически все аспекты архитектуры микроконтроллеров: работу портов, таймеров, арифметическую обработку данных, организацию ввода/вывода по параллельным и последовательным (UART, SPI и I2C) каналам связи, устройств обработки аналоговых сигналов, системы прерывания. Базовые программы могут быть использованы в качестве основы для обучения и самостоятельного программирования на языке Ассемблер AVR в курсовом и дипломном проектировании.

Материалы книги автор использует в учебном процессе в МГТУ им. Н.Э. Баумана.

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

По вопросам приобретения обращаться по тел. 499-263-60-45; e-mail: press@bmstu.ru