О тороидальных поверхностях вращения с постоянной средней

advertisement
В.И. Ильгисонис, А.А. Сковорода, Е.А. Сорокина
УДК 517.972.6
О ТОРОИДАЛЬНЫХ ПОВЕРХНОСТЯХ ВРАЩЕНИЯ С ПОСТОЯННОЙ
СРЕДНЕЙ КРИВИЗНОЙ
В.И. Ильгисонис1—3, А.А. Сковорода1, Е.А. Сорокина1—3
1НИЦ
«Курчатовский институт», Москва, Россия
университет дружбы народов, Москва, Россия
3Национальный исследовательский ядерный университет «МИФИ», Москва, Россия
2Российский
Показано, что экстремаль объёма, заключённого внутри тороидальной поверхности заданной площади, доставляется поверхностью постоянной средней кривизны. Получено точное решение соответствующей вариационной задачи, и проведён его параметрический анализ в пределах большой и малой средней кривизны. Доказано отсутствие гладких торов постоянной средней
кривизны и продемонстрировано, что искомая экстремальная поверхность имеет, по крайней мере, одно ребро, расположенное
на наружном обводе тора.
Ключевые слова: поверхности постоянной средней кривизны, торы, токамак.
ON THE TOROIDAL SURFACES OF REVOLUTION WITH CONSTANT MEAN
CURVATURES
V.I. Ilgisonis1—3, A.A. Skovoroda1, E.A. Sorokina1—3
1NRC
«Kurchatov Institute», Moscow, Russia
Friendship University of Russia, Moscow, Russia
3National Research Nuclear University «MEPhI», Moscow, Russia
2Peoples'
It is shown that the surface with constant mean curvature encloses the extremal volume among all toroidal surfaces of given area. The
exact solution for the corresponding variational problem is derived, and its parametric analysis is performed in the limits of high and
small mean curvatures. An absence of smooth torus with constant mean curvature is proved, and the extremal surface is demonstrated to
have at least one edge located on the outer side of the torus.
Key words: surfaces with constant mean curvature, tori, tokamak.
ВВЕДЕНИЕ
Среди систем магнитного удержания высокотемпературной плазмы наибольшее распространение получили токамаки — осесимметричные системы, в которых силовые линии магнитного поля
образуют набор тороидально вложенных поверхz
ностей [1], называемых «магнитными поверхностями» (см. рис. 1). Геометрические свойства этих
поверхностей в значительной степени определяют
способность токамака удерживать плазму, а также
устойчивость плазмы по отношению к магнито1
гидродинамическим и дрейфовым возмущениям. В
0
частности, важное значение имеют характеристики
r
кривизны силовых линий и магнитных поверхностей в целом, определяющие дрейфы заряженных
частиц плазмы, направление её макроскопического
5
4
вращения и возможность раскачки ряда неустойчивостей плазмы (см., например, [2—4]). Наряду с
3
2
этим задача оптимизации токамака требует увеличения отношения объёма удерживаемой тороиРис. 1. Сечение плоскостью r—z магнитной конфигурации
типового токамака: 1 — сечения магнитных поверхностей; дальной плазмы к площади, ограничивающей
2 — дивертор; 3 — диверторные пластины; 4 — сепаратри- плазму тороидальной поверхности, через которую
са; 5 — засепаратрисный слой (scrape-off layer)
происходят потери частиц и энергии.
22
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
О тороидальных поверхностях вращения с постоянной средней кривизной
В данной статье показано, что такая оптимизация естественным образом приводит к задаче нахождения тороидальной поверхности с постоянной средней кривизной. Как известно, теорема Александрова
[5] доказывает отсутствие гладких торов без самопересечений с постоянной средней кривизной. Вместе
с тем существуют примеры незамкнутых поверхностей с постоянной средней кривизной, к числу которых относятся нодоиды и ондулоиды [6]. Проведённое далее рассмотрение доказывает, что тороидальные поверхности вращения максимального объёма обладают постоянной средней кривизной почти всюду, имея как минимум одно ребро. Важность таких конфигураций в теории магнитных ловушек обсуждалась в [7]. В частности, применительно к магнитной конфигурации токамака это означает возможности организации дивертора — устройства, собирающего потоки вылетающих из системы частиц плазмы
и предохраняющего тем самым от повреждений стенку камеры токамака (см. рис. 1).
ВАРИАЦИОННАЯ ЗАДАЧА О ТОРЕ МАКСИМАЛЬНОГО ОБЪЁМА
Рассмотрим следующую вариационную задачу: найти кривую вращения z(r) (вращение вокруг оси
z — прямой r = 0 в цилиндрических координатах {r, , z}), однозначно определённую в интервале r  [r1, r2]
и заметающую максимальный объём при заданной величине поверхности S, такую, что z(r1) = z(r2) = 0.
Объём будем отсчитывать от плоскости z = 0. Впоследствии вместо задания величины поверхности можно будет задать значение zm  max ( z (r )). Другими словами, фиксируем три числа: r1, r2 и S (или zm).
r  [ r1,r2]
Искомый объём равен
r2
V   2πrz (r )dr ,
r1
поверхность
r2
S   2πr 1  zr2 ,
r1
где zr = dz/dr.
Решение должно иметь вид экстремали функционала
r2
(
)
J  V  λRS = 2πR3  r z  λ 1  zr2 dr .
(1)
r1
Здесь λ — множитель Лагранжа, и все величины под интегралом обезразмерены на некоторый характерный радиус R:
r
r/R; z
z/R.
Значение R выберем позже. Уравнение Эйлера—Лагранжа для функционала J = J/(2R3) имеет вид
r
z'r
λ
1
z' 2
'
.
(2)
Здесь и далее будем опускать тильды, где это не приводит к недоразумениям, а также нижний индекс r у производной, обозначаемой штрихом. Смысл параметра λ тесно связан с понятием средней кривизны H s поверхности, образованной кривой вращения z(r), определяемой уравнением (2), а именно:
λ
1
.
2H s R
(3)
Напомним, что по определению H s — величина локальная, определяемая как среднее арифметическое главных кривизн поверхности в рассматриваемой точке.
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
23
В.И. Ильгисонис, А.А. Сковорода, Е.А. Сорокина
Интегрируя (2), получаем
z'
2rλ
r2
1
z' 2
= c1,
(4)
где c1 — постоянная интегрирования. Выберем нормировочную величину R так, чтобы значение
rR = r/R = 1 отвечало точке максимума z(r), т.е. z' (rR ) 0. Тогда из (4) следует c1 = 1.
Введём замену переменных z' = –tg, где φ = 0 соответствует этой же точке максимума z(r):
z  r   zm ; zm  z  r  r  1 .
В отсутствие иных экстремумов z(r) на отрезке r  [r1, r2], появление которых с очевидностью противоречило бы нашей вариационной задаче, область r1  r  1 отвечает значениям φ < 0 (z' > 0), а область
π π
,
и для всех φ cos   0.
2 2
1  r  r2 — значениям φ > 0 (z' < 0). Поэтому 
Переписав уравнение (4) в виде
r 2 + 2rsin = 1,
находим
λsin
r
1 λ 2sin 2 .
Поскольку при φ = 0 величина r = 1, двойственность в решении устраняется, т.е.
λsin
r
1 λ 2sin 2.
(5)
При этом должно быть λ < 0, чтобы меньшие значения φ отвечали меньшим значениям r, как это
предполагалось.
Зависимость z(r) находим следующими манипуляциями:
z′ = –tg;
1
dz dr
d d
dz
d
z
λcos
tg;
sin
λsin 1
1 λ 2sin 2
λ 2sincos
1
1 λ sin 
1 λ
2
2
2
;
(F (1 λ 2 )E)
c2 ,
(6)
где эллиптические интегралы F и E берутся как функции аргументов 1 и 2:
α1
arcsin
1 λ 2 sin
1 λ 2sin 2 
; α2
λ
1 λ2
.
(7)
Так как F(0, ) = E(0, ) = 0, то z( = 0) = c2 – λ, т.е. c2 = zm + λ.
ПРИБЛИЖЁННЫЕ РЕШЕНИЯ И ИХ ПАРАМЕТРИЧЕСКИЙ АНАЛИЗ
Построим приближённое решение задачи для случая λ2 << 1. Из системы (5), (6) имеем
r
z
24
λ2
sin 2;
2
λ2
λ(1 cos)

2
1 λsin
zm
sin2
2
(8)
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
О тороидальных поверхностях вращения с постоянной средней кривизной
— с учётом уже описанного выбора постоянной интегрирования. Пренебрежение членами ~λ2 в (8) даёт
решение в виде дуг окружностей, как в изопериметрической задаче. Действительно, при φ = φ 1,2 z = 0 и
cos1, 2
λ
zm
zm
λ
1
λ
zm
.

1
π π
, :
2 2
Отсюда находим граничные значения угла φ в диапазоне 
 1, 2
причём λ
zm
,
λ
arccos 1
(9)
zm . Тогда
sin1, 2
1 cos 1, 2
2
1
2
zm )
1
zm
λ
zm ( 2λ
λ2
r1
1
zm (2 λ
zm );
r2
1
zm (2 λ
zm )
1
λ
zm (2 λ
zm ).
Условия
(10)
позволяют, вернувшись к размерным переменным, по заданным r1 , r2 и zm определить искомые R и λ.
Если же говорить собственно о безразмерном решении, задаваемом (10), то из-за связи r1 и r2 достаточно задать только два безразмерных числовых параметра, в качестве которых удобно выбрать r1 :
0 r1 1 и zm . Тогда величина третьего параметра и значение λ находятся из условий (10).
Соотношения (10) можно также получить из явной формы искомой кривой в переменных r – z. Действительно, в первом порядке разложения по λ из уравнений (8) имеем
λsin 1 r;
λcos zm λ
(11)
z,
так что
(r 1)2
(z
λ)2
zm
λ2 .
(12)
Это же уравнение нетрудно получить, используя только первое из уравнений (11) и связь
cos
(при этом, разумеется, 1
λ
1
r
tg
z'
1 sin 2 
1
2
1 r
λ
λ ). Действительно, интегрируя
r 1
λ
1
1 r

z
λ
λ2
2
1 r
2
(1 r )2
с условием z(r = 1) = zm, получаем
zm
(1 r )2 ,
(13)
откуда немедленно следует (12).
Используем этот же приём для нахождения зависимости z(r) с учётом членов ~λ2 в (8). Из первого
уравнения системы (8) имеем
λsin
1 r
λ
1
sin
2
1 r
1
1
(1 r )
2
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
2
1 r
.
1 r
(14)
25
В.И. Ильгисонис, А.А. Сковорода, Е.А. Сорокина
Заметим, что соотношение (14) работает при r:
1
λ /2
1
λ /2
r
1
λ /2
1
λ /2
,
что практически соответствует ранее (в первом порядке разложения по λ) найденной области
1 λ r 1 λ . Второе уравнение системы (8) даёт искомую зависимость z(r) в форме
z  zm  λ  λ 
2
4  r  1
 r  1
2
2

λ 2 
 r  1   r  1
 r 1 
arcsin 

1 


2
λ
 λ 
 λ 

2

,


(15)
модифицирующей выражение (13). В отличие от (13) кривая (15) уже не является дугой окружности и не
симметрична относительно точки r = 1.
Положив в (15) z = 0, можно найти r1 и r2 . Как уже отмечалось, это возможно не для любого zm.
Безразмерные величины r1 , r2 и zm не вполне произвольны и в общем случае. Связь между ними удобно
проследить в зависимости от величины λ. Очевидно, что поскольку задаваемая (5) зависимость r () монотонна, то r2
r1
2 λ . Если теперь задать значение r1 , то r2 и zm однозначно определятся значением
λ. Действительно, в общем случае из (4) находим
r
z (r )
r1
1 r2
4λ 2 r 2
(r 2 1) 2
dr ,
(16)
откуда немедленно следуют искомые определения
r2
r1
При λ
1 r2
4λ2 r2
(r2 1)2
1
dr
r1
0 соотношения (17) дают r2
(10). При λ
0; zm
1 r2
4λ2 r2
1 1 r1 , zm
λ
(r2 1)2
λ2
(17)
dr.
(1 r1 )2 , что соответствует
соотношения (17) интегрируются в виде
1
r
ln 2
2λ
r1
1 2
(r1
2
r2 2 )
0; zm
1 1 2
(r1 1) ln r1 .
2λ 2
(18)
1
Как видно из (18), при больших λ величина r2 не меняется, а zm убывает как ~ λ .
В обратную сторону при заданной величине r1  (0, 1) решение безразмерной задачи в виде
поверхности постоянной средней кривизны существует лишь для z m < 1 — только в этом случае в соответствии с (17) существует искомое
значение λ. В размерной задаче это отвечает
требованию z m < R, где величина R определяется
0,2
z(r)
0,15
первым из соотношений (17) для заданных r1, 2 .
0,1
Проиллюстрируем сказанное примерами. На
рис. 2—4 показаны кривые (16), рассчитанные
для r1 0,8 и различных значений λ: –0,225
0,05
0
1
1,1
1,2
r
Рис. 2. Образующая поверхности с постоянной средней кривиз0,8
0,9
ной для r1 = 0,8 и λ = –0,225, при этом zm = 0,186, r2 = 1,24
26
(см. рис. 2), –0,5 (см. рис. 3) и –2,0 (см. рис. 4).
Видно, что с ростом λ величина zm убывает, а значение r2 меняется несущественно. Поверхность,
формируемая изображённой на рис. 2 кривой, пе-
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
О тороидальных поверхностях вращения с постоянной средней кривизной
0,015
0,05
0,04
z(r)
z(r)
0,01
0,03
0,02
0,005
0,01
0
0
1,1
1
1,2
r
Рис. 3. Образующая поверхности с постоянной средней кривизной для r1 0,8 и λ = –0,5, при этом zm = 0,045, r 2 = 1,218
0,8
0,8
1
1,1
1,2
r
Рис. 4. Образующая поверхности с постоянной средней кривизной для r1 0,8 и λ = –2,0, при этом zm = 0,011, r 2 = 1,215
0,9
0,9
ресекающей ось z = 0 под прямым углом (1 = –/2), характерна для петель нодоида. В общем случае
кривая (16) пересекает ось z = 0 под углами, отличными от прямого. Таким образом, искомая поверхность имеет, как правило, два ребра.
О ТЕОРЕМЕ АЛЕКСАНДРОВА ПРИМЕНИТЕЛЬНО К СИММЕТРИЧНЫМ ТОРАМ
Полученное нами общее решение (5)—(7) исходной задачи позволяет доказать утверждение об отсутствии гладких тороидальных поверхностей постоянной средней кривизны, вытекающее из теоремы А.Д. Александрова [5]. Проведём сначала рассуждение в классе симметричных относительно
плоскости z = 0 поверхностей. Если бы такая симметричная поверхность существовала, она бы по построению описывалась уравнением (2) с λ = const и, следовательно, формулами (5)—(7). Искомая гладкая кривая вращения пересекала бы ось z = 0 под прямым углом, т.е. при значениях  = /2 и соответственно 1 = /2 (см. (7)). В этих точках из (6) для z = 0 имеем
1
λ
1
2
1
Здесь K (α 2 )
F
1
λ2

, α 2 и E(α 2 )
2
E
( K(α 2 )
(1
λ 2 )E(α 2 ))
λ
zm
0;
(19)
(K (α 2 )
(1
λ )E(α 2 ))
2
λ
zm
0.
π
, α 2 — полные эллиптические интегралы (использована не2
чётность эллиптических интегралов {F, E}( α1, α2 )
{F, E}(α1, α2 )).
Вычитая уравнения (19), приходим к соотношению
K(2) – (1 + 2)E(2) = 0.
(20)
Нетрудно видеть, что условие (20) не может быть выполнено. Действительно, по определению и с
учётом (7)
 /2
K (α 2 )
(1 λ 2 )E(α 2 )
0
 2 (sin 2 x 1)
1 α 22sin 2 x
dx
0
(21)
при λ ≠ 0, что доказывает утверждение.
Неравенство (21) позволяет также заключить, что если кривая z(r), формирующая симметричную
тороидальную поверхность постоянной средней кривизны, пересекает ось z = 0 перпендикулярно, то это
происходит при r < R, а в точке пересечения при r > R на кривой всегда находится «зубик» (см. рис. 2).
Действительно, если данная кривая z(r) пересекает ось z = 0 при  1 = –/2, то в силу (21) второе пересечение, отвечающее зоне r  R, произойдёт при значении  2 = /2, так как z' < 0 при φ > 0. И наоборот,
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
27
В.И. Ильгисонис, А.А. Сковорода, Е.А. Сорокина
если допустить, что кривая z(r) пересекает ось z = 0 при  1 = /2, т.е. под прямым углом, то в силу (21) она
уже не пересечёт ещё раз эту ось ни при каком  2  [–/2, 0).
Вытекающая из этого же неравенства (21) знакоопределённость разности z( = /2) – z( = –/2) позволяет распространить утверждение и на случай несимметричных поверхностей.
ПРИМЕР ЧИСЛЕННОЙ РЕАЛИЗАЦИИ
z, см
Продемонстрируем возможность построения магнитной конфигурации, ограниченной сепаратрисной поверхностью с постоянной средней кривизной, посредством набора витков с током, создающих
полоидальное поле, и сравним указанную конфигурацию со «стандартной», полученной из исходной
поворотом на 90º вокруг магнитной оси. При такой трансформации Х-точка, расположенная в экваториальной плоскости исходной конфигурации, оказывается в более «привычном» положении в нижней части сечения плазменного шнура. Сравнение проведём при одинаковых полном токе, магнитном поле на
оси токамака и среднем (по объёму) давлении плазмы.
Искомая конфигурация, рассчитанная кодом ДИНА [8], показана на рис. 5, здесь же изображены сечения витков с током. С помощью кодов VMEC [9] и JMC [10] были проведены расчёты равновесия и
МГД-устойчивости плазмы в такой конфигурации вплоть до давления плазмы с β = 1%.
Сравнение данной конфигурации со «стан80
дартной» показало, что амплитуды первой гар60
моники вторичных токов в обеих конфигурациях
практически одинаковы, тогда как амплитуда
40
второй гармоники вторичного тока больше в
конфигурации с экваториальной Х-точкой из-за
20
большей треугольности. Конфигурация с сепа0
ратрисной поверхностью с постоянной средней
кривизной оказалась более устойчивой по от–20
ношению к резистивным модам и модам Мер–40
сье, что связано с различиями в глубине и условиях формирования «магнитной ямы». В конфи–60
гурации с нижним расположением Х-точки маг–80
нитная яма образуется только при достаточно
0
50
100
150
200
большой треугольности, тогда как в конфигуr, cм
рации с экваториальной Х-точкой глубокая
Рис. 5. Пример магнитной конфигурации с сепаратрисной поверхностью постоянной средней кривизны: аспектное отношение 2,7, магнитная яма формируется практически по
средний радиус плазмы 0,38 м, эллиптичность 0,73, треугольность всему сечению шнура при небольшой вели0,08, шафрановское смещение 0,4 см, объём плазмы 3 м3
чине шафрановского смещения. С этой точки
зрения «стандартная» конфигурация какимилибо явными преимуществами по сравнению с рассмотренной нами конфигурацией с экваториальной
Х-точкой не обладает.
ЗАКЛЮЧЕНИЕ
Мы продемонстрировали, что поверхность фиксированной площади осесимметричного тора ограничивает максимальный объём в случае, если её средняя кривизна почти всюду постоянна. Такая поверхность должна иметь, по крайней мере, одно ребро, расположенное на внешнем обводе тора. Установлены допустимые диапазоны изменения параметров, характеризующих данную поверхность. Построены точные и приближённые решения соответствующей вариационной задачи, которые могут быть
полезны для создания оптимизированных оболочек различных ёмкостей, а также для расчётов магнитных конфигураций установок токамак. Предъявлен численный пример реализации ограничивающей область плазмы сепаратрисы с постоянной средней кривизной. Отметим, что возможность формирования в
реальном эксперименте сепаратрисы с Х-точкой, расположенной в экваториальной плоскости токамака,
28
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
О тороидальных поверхностях вращения с постоянной средней кривизной
сомнений не вызывает. Примером может служить японский токамак JT-60 с экваториальной Х-точкой,
который продемонстрировал в 1985 г. инженерные преимущества расположения диверторных пластин
на внешнем обводе токамака [11]. Как известно, вскоре эта установка была переделана в «стандартную»
систему JT-60U с нижним расположением Х-точки с целью экспериментального обнаружения улучшенного удержания (H-моды). Возможность реализации Н-моды в исследуемых конфигурациях с экваториальным расположением Х-точки, равно как и другие более «тонкие» свойства таких систем, требует, разумеется, отдельного рассмотрения.
Работа выполнена при частичной поддержке грантами РФФИ (№ 16-02-00640) и Президента РФ для
ведущих научных школ (№ НШ-8243.2016.2).
СПИСОК ЛИТЕРАТУРЫ
1. Соловьёв Л.С., Шафранов В.Д. Замкнутые магнитные конфигурации для удержания плазмы. — В кн.: Вопросы теории
плазмы. Вып. 5. Под ред. М.А. Леонтовича. — М.: Атомиздат, 1967, с. 3—208.
2. Морозов А.И., Соловьёв Л.С. Движение заряженных частиц в электромагнитных полях. — Там же, вып. 2. Под ред.
М.А. Леонтовича. — М.: Госатомиздат, 1963, с. 177—261.
3. Ильгисонис В.И. Классические результаты Б.Б. Кадомцева и вращение плазмы в современных токамаках. — УФН, 2009,
т. 179, с. 790—800.
4. Ilgisonis V.I., Lakhin V.P., Smolyakov A.I., Sorokina E.A. Geodesic acoustic modes and zonal flows in rotating large-aspect-ratio
tokamak plasmas. — Plasma Phys. Control. Fusion, 2011, vol. 53, р. 065008.
5. Александров А.Д. Теоремы единственности для поверхностей «в целом». — Вестник ЛГУ. Сер. математики, механики и
астрономии, 1959, т. 19, с. 5—8.
6. Kenmotsu K. Surfaces with constant mean curvature. — In: Translations of Mathematical Monographs, vol. 221. — Providence:
American Mathematical Society, 2003.
7. Сковорода А.А., Тайманов И.А. О значении средней кривизны в геометрии магнитного поля ловушек для удержания
плазмы. — Физика плазмы, 2010, т. 36, с. 874—878.
8. Khayrutdinov R.R., Lukash V.E. Studies of plasma equilibrium and transport in a tokamak fusion device with the inverse-variable
technique. — J. Comput. Phys., 1993, vol. 109, p. 193—201.
9. Hirshman S.P., Betancourt O. Preconditioned descent algorithm for rapid calculations of magnetohydrodynamic equilibria. —
J. Comput. Phys., 1991, vol. 96, p. 99—109.
10. Nuehrenberg J., Zille R. Theory of Fusion Plasmas. — In: Proc. of the Workshop on Theory of Fusion Plasmas. Varenna, 1987,
vol. EUR 11336 EN. — Bologna: Editrice Compositori, 1988, p. 3—23.
11. JT-60 Team. Recent experiments in JT-60. Plasma physics and controlled nuclear fusion research. — In: Proc. of the Eleventh Intern. Conf. on Plasma Physics. Kyoto, 1986, vol. 1. — Vienna: IAEA, 1987, p. 11—30.
Виктор
Игоревич
Ильгисонис,
зам.
директора по научной работе – главный учёный секретарь, д.ф.-м.н., профессор; НИЦ «Курчатовский
институт»,
123182
Москва, пл. Академика Курчатова 1,
Россия
Александр АлексеЕкатерина
Алексеевич
Сковорода,
евна Сорокина, с.н.с.,
г.н.с.,
д.ф.-м.н.,
к.ф.-м.н.; НИЦ «Курветеран
атомной
чатовский институт»,
энергетики и про123182 Москва, пл.
мышленности; НИЦ
Академика Курчато«Курчатовский инва 1, Россия
ститут»,
123182
Москва, пл. Академика Курчатова 1,
Россия
Skovoroda_AA@nrc
ki.ru
Статья поступила в редакцию 30 декабря 2015 г.
Вопросы атомной науки и техники.
Сер. Термоядерный синтез, 2016, т. 39, вып. 1, с. 22—29.
ВАНТ. Сер. Термоядерный синтез, 2016, т. 39, вып. 1
29
Download