Процессы управления и устойчивость



жүктеу 30.48 Mb.
Pdf просмотр
бет18/57
Дата27.12.2016
өлшемі30.48 Mb.
1   ...   14   15   16   17   18   19   20   21   ...   57

ζ

) dZ


ζ

dζ+


M

T,u


G(Z

T

, X



T

) (T, Z


T

) dZ


T

,

(5)



141

где g и G — некоторые интегрируемые по ζ, Z и дифференцируемые

по Z и X функции, характеризующие качество пучка.

Вариацию функционала (5) можно записать в виде

δI = −


T

0 M


ζ,u

ψ



x

(ζ, Z


ζ

)∆

u



h

1

(ζ, Z



ζ

, X(ζ, Z


ζ

), u)+


z



(ζ, Z

ζ

)∆



u

f

1



(ζ, Z

ζ

, u)



(ζ, Z

ζ

) dZ



ζ

dζ,


(6)

где Z


ζ

≡ Z(ζ). При этом сопряженные функции ψ

z

, ψ


x

удовлетворя-

ют системе уравнений



x

= −ψ



x

∂h



1

(ζ, Z, X, u)

∂X

+

M



ζ,u

∂h

2



(X, Z, X , Z )

∂X

(ζ, Z ) dZ



M



ζ,u

ψ



x

(ζ, Z )


∂h

2

(X, Z, X , Z )



∂X

(ζ, Z ) dZ +

∂g(ζ, Z, X)

∂X

,



z



= −ψ


z

∂f



1

(ζ, Z, u)

∂Z

+

M



ζ,u

∂f

2



(Z, Z )

∂Z

(ζ, Z ) dZ



M



ζ,u

ψ



z

(ζ, Z )


∂f

2

(Z, Z )



∂Z

(ζ, Z ) dZ −

− ψ



x



∂h

1

(ζ, Z, X, u)



∂Z

+

M



ζ,u

∂h

2



(X, Z, X , Z )

∂Z

(ζ, Z ) dZ



M



ζ,u

ψ



x

(ζ, Z )


∂h

2

(X, Z, X , Z )



∂Z

(ζ, Z ) dZ +

∂g(ζ, Z, X)

∂Z

.



(7)

и конечным условиям

ψ



z



(T, Z

T

) = −



∂G(Z, X)

∂Z

Z=Z



T

, X=X(T,Z

T

)

,



142

ψ

x



(T, Z

T

) = −



∂G(Z, X)

∂X

Z=Z



T

, X=X(T,Z

T

)

.



Для численного решения задачи управления построим в фазовом

пространстве продольного движения сетку так, чтобы она покрыва-

ла множество M

0

. Это может быть сетка вида



z

i

j



= z

i

0



+ jh

i

,



j = 0, m

i

,



i = 1, n.

(8)


Будем обозначать через z

j

1



,...,j

n

(ζ) точки траекторий, начинающихся



при ζ = 0 в узлах сетки:

z

j



1

,...,j


n

(0) = (z


1

j1

, . . . , z



n

jn

)



.

Аналогично, обозначим соответствующие значения X(ζ, Z), (ζ, Z),



ψ

x

(ζ, Z),



ψ

z

(ζ, Z)



через

x

j



1

,...,j


n

(ζ),


j

1

,...,j



n

(ζ),


ψ

z

j1,...,jn



(ζ),

ψ

z



j1,...,jn

(ζ) соответственно.

Для вычисления интегралов, входящих в правую часть уравне-

ния динамики (4), в уравнения (7), функционал (5), выражение для

вариации функционала (6) будем использовать известные квадра-

турные формулы.

Для простоты рассмотрим случай, когда область M

0

представля-



ет собой n-мерный параллелепипед. Другие случаи можно свести к

обсуждаемому, если положить

j

1

,...,j



n

= 0 в узлах, не принадлежа-

щих области M

0

. Выберем сетку в виде (8) с узлами, лежащими на



гранях параллелепипеда. Обозначим весовые коэффициенты квад-

ратурной формулы через γ

j1,...,jn

.

Используя то, что дифференциальная форма плотности являет-



ся интегральным инвариантом, квадратурную формулу для вычис-

ления функционала I(u) можно привести к виду

I(u)

T

0 j



1

,...,j


n

γ

j



1

,...,j


n

g(ζ, z


j

1

,...,j



n

(ζ), x


j

1

,...,j



n

(ζ))


j

1

,...,j



n

(0) dζ+


+

j

1



,...,j

n

γ



j

1

,...,j



n

G(z


j

1

,...,j



n

(T ), x


j

1

,...,j



n

(T ))


j

1

,...,j



n

(0).


143

К аналогичному виду приводятся формулы для вариации функцио-

нала (6), а также уравнения динамики и уравнения для сопряженных

функций вдоль характеристических линий, начинающихся в узлах

сетки:


δI(u)

T

t



0

j

1



,...,j

n

γ



j

1

,...,j



n

×

× ψ



xj

1



,...,j

n

(ζ)∆



u

h

1



(ζ, z

j

1



,...,j

n

(ζ), x



j

1

,...,j



n

(ζ), u)+


zj



1

,...,j


n

(ζ)∆


u

f

1



(ζ, z

j

1



,...,j

n

(ζ), u)



j

1

,...,j



n

(0) dt;


xj

1



,...,j

n



= −ψ

xj



1

,...,j


n

{

∂h



1

ζ, Z, X, u

∂X

Z=z


j

1

,...,j



n (ζ)

+

+



j1,...,jn

γ

j1,...,jn



∂h

2

X, Z, X , Z



∂X

Z=z


j1,...,jn

(ζ)


Z =z

j1,...,jn

j(ζ)

j

1



,...,j

n

j(0)}−



j

1



,...,j

n



j

1

,...,j



n

xj



1

,...,j


n

j

∂h



2

X, Z, X , Z

∂X

+

+



∂g ζ, Z, X

∂X



Z=z

j1,...,jn

(ζ)

;



zj

1

,...,j



n

= −ψ



zj

1



,...,j

n

{



∂f

1

ζ, Z, u



∂Z

Z=z


j1,...,jn

(ζ)


+

+

j



1

,...,j


n

j



1

,...,j


n

j

∂f



2

Z, Z


∂Z

Z=z


j1,...,jn

(ζ)


Z =z

j1,...,jn

j(ζ)

j

1



,...,j

n

j(0)}−



j

1



,...,j

n



j

1

,...,j



n

zj



1

,...,j


n

j

∂f



2

Z, Z


∂Z

Z=z


j1,...,jn

(ζ)


Z =z

j1,...,jn

j(ζ)

j

1



,...,j

n

j(0)−



144

−ψ

xj

1



,...,j

n

{



∂h

1

ζ, Z, X, u



∂Z

Z=z


j1,...,jn

(ζ)


+

+

j



1

,...,j


n

j



1

,...,j


n

j

∂h



2

X, Z, X , Z

∂Z

Z=z


j1,...,jn

(ζ)


Z =z

j1,...,jn

j(ζ)

j

1



,...,j

n

j(0)}−



j

1



,...,j

n



j

1

,...,j



n

xj



1

,...,j


n

j

∂h



2

X, Z, X , Z

∂Z

Z=z


j1,...,jn

(ζ)


Z =z

j1,...,jn

j(ζ)

j

1



,...,j

n

j(0)+



+



∂g ζ, Z, X

∂Z



Z=z



j

1

,...,j



n (t)

.

Конечные условия для ψ



xj

1

,...,j



n

и ψ


zj

1

,...,j



n

приобретают вид

ψ



zj



1

,...,j


n

(T ) = −


∂G Z, X

∂Z

Z=z



j1,...,jn

(T )


,

ψ



xj

1

,...,j



n

(T ) = −


∂G Z, X

∂X

Z=z



j1,...,jn

(T )


.

Наиболее простой вид полученные выражения примут, если при-

менить квадратурную формулу с одинаковыми коэффициентами

γ

j



1

,...,j


n

для всех узлов. В этом случае интегралы, входящие в эти

формулы, аппроксимируются суммами по узлам от подынтеграль-

ных выражений.

Литература

1. Овсянников Д.А., Дривотин О.И. Моделирование интенсивных

пучков заряженных частиц. СПб.: Изд-во СПбГУ, 2003. 172 c.

145


Дривотин О.И., Семёнов В.В.

Санкт-Петербургский государственный университет

Исследование плотности заряженных частиц

при моделировании пучка

в продольном магнитном поле

Одной из математических моделей, описывающих ансамбль вза-

имодействующих частиц, является уравнение Власова. Его решения

называются самосогласованными распределениями. Известны раз-

личные аналитические решения уравнения Власова для пучка в про-

дольном магнитном поле. При этом рассматривается фазовое про-

странство поперечного движения, и его размерность равна 4.

Целью работы является сравнение плотностей распределения ча-

стиц, получаемых в результате моделирования, с плотностями для

этих аналитических решений.

Уравнения характеристических линий для уравнения Власова,

которые совпадают с уравнениями динамики частиц в самосогласо-

ванном поле в безразмерных переменных имеют вид [1]

q = −(1 − α)q +

M

2

q



3

,

(1)



где α — безразмерный параметр, характеризующий интенсивность

пучка, α ∈ [0, 1], q — безразмерная радиальная координата, q ∈ [0, 1],

M — безразмерная угловая компонента канонически сопряжённого

импульса, сохраняющегося в силу теоремы Буша.

Уравнение (1) имеет интеграл

H = q


2

+ (1 − α)q

2

+

M



2

q

2



.

Для рассматриваемых самосогласованных распределений интегралы

M и H лежат внутри множества

2ω|M | < H ≤ M

2

+ 1 − α,


M ∈ (−

1 − α,



1 − α).


(2)

146


В работе [1] рассматривалось распределение Капчинского – Вла-

димирского. В настоящей работе рассмотрим распределение с посто-

янной фазовой плотностью внутри множества допустимых значений

(2). Для такого распределения плотность распределения в простран-

стве (M , H) равна

ρ

M H



=

λ

1 − α + M



2

− H


,

(3)


где λ — положительный безразмерный параметр.

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

рассматривался как состоящий из макрочастиц, каждая из которых

представляет бесконечно тонкое кольцо радиуса q с центром на оси

пучка и в плоскости, ортогональной к оси пучка. Для описания ди-

намики таких частиц использовалось уравнение

q = −q +

M

2



q

3

+



α

q

N



q

,

где N



q

— доля макрочастиц, находящихся внутри кольца радиуса q.

При моделировании оценка плотности распределения производи-

лась на основе различных статистик [1]. В качестве одной из таких

статистик предлагается брать расстояние от заданной до ближайшей

точки в пространстве интегралов M , H. Обозначим эту величину че-

рез S

1

. В качестве второй статистики S



2

предполагается брать сред-

нее арифметическое для расстояний от заданной точки до первой и

до второй ближайших к ней. Аналогично, в качестве S

n

берётся сред-



нее арифметическое расстояний от данной точки до n ближайших к

ней. Связь между эмпирическим значением плотности и статистика-

ми выражается формулой

ρ

M H



= k

i

S



i

−2

,



i = 1, 20,

(4)


где черта обозначает усреднение по всем частицам, а k

i

— коэффици-



енты, соответствующие статистике S

i

. Значения k



i

были определены

ранее после ряда численных экспериментов.

147


В данной работе моделирование проводилось при количестве ча-

стиц, равном 10000. Отклонение плотности ρ

M H

от теоретического



значения ρ

M H


на разных этапах моделирования (при z = 0, z = 5,

z = 10, где z — продольная координата) оценивалось с помощью

отношения

∆ =


ρ

M H


ρ

M H


.

(5)


Результаты моделирования приведены ниже в таблице.

Таблица. Значения ∆

i

, при z = 0, z = 5, z = 10



S

i

k



i

z=0


z=5

z=10


1

0,254


4,620

4,560


4,929

2

0,396



1,838

1,836


1,806

3

0,620



1,326

1,298


1,272

4

0,895



1,202

1,190


1,153

5

1,198



1,149

1,145


1,112

6

1,517



1,129

1,125


1,095

7

1,842



1,117

1,112


1,086

8

2,172



1,109

1,103


1,081

9

2,506



1,104

1,097


1,076

10

2,838



1,098

1,091


1,072

11

3,174



1,095

1,087


1,068

12

3,511



1,092

1,083


1,066

13

3,851



1,090

1,081


1,064

14

4,185



1,087

1,077


1,061

15

4,521



1,084

1,074


1,058

16

4,857



1,081

1,071


1,055

17

5,194



1,079

1,068


1,053

18

5,528



1,077

1,066


1,051

19

5,861



1,074

1,063


1,048

20

6,198



1,072

1,061


1,046

21

6,539



1,071

1,060


1,045

22

6,878



1,070

1,058


1,044

Проводилась оценка того, насколько хорошо эмпирическая фор-

мула (4) описывает теоретическую формулу (3). В ходе моделирова-

ния было установлено, что значения плотности, полученные по эм-

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

ским значениям. Отношение (5) между ρ

M H

и ρ


M H

оказалось близко

к единице, и с увеличением номера i статистики S

i

оно стремится



к единице. При этом с ростом количества частиц значение, давае-

мое формулой (5), также стремится к единице. Другими словами,

148


эмпирическое значение (4) стремится к теоретическому (3) с ростом

числа частиц. Кроме того, представленные результаты показывают,

что в процессе интегрирования отношение (5) примерно сохраняет

свое значение, что позволяет говорить о технической устойчивости

рассматриваемых распределений.

В заключении хотелось бы отметить, что, вычисляя плотность

по эмпирической формуле (4), можно оценивать качество алгоритма

задания начального распределения частиц, а также изучать поведе-

ние частиц на разных этапах моделирования, сравнивая полученные

результаты с аналитическим решением уравнения Власова.

Таким образом методика, разработанная первоначально для рас-

пределения Капчинского – Владимирского, применима и для других

распределений.

Литература

1. Семенов В.В. Моделирование самосогласованных распределений

для пучка заряженных частиц в продольном магнитном поле //

Процессы управления и устойчивость: Труды 35-й научной кон-

ференции аспирантов и студентов / Под ред. Н.В. Смирнова, В.Н.

Старкова. – СПб.: Изд-во СПбГУ, 2004. С. 262–267.

2. Дривотин О.И., Овсянников Д.А. Самосогласованные распределе-

ния для пучков заряженных частиц. СПб.: Изд-во СПбГУ, 2001.

108 с.


149

Зимин А.В.

Санкт-Петербургский государственный университет

О построении вейвлет-разложения

на периодической неравномерной сетке

1

Рекомендовано к публикации профессором Демьяновичем Ю.К.



Имеется большое количество работ, в которых строится вейвлет-

ное разложение на равномерной сетке (см. [1–3]). Цель данной рабо-

ты состоит в том, чтобы построить вейвлетное разложение телеско-

пической системы пространств B-сплайнов на неравномерной сетке

[4,5].

1. Некоторые обозначения. Пусть Z — множество всех целых



чисел, а N — натуральное число; будем обозначать 0 : N — мно-

жество целых чисел от 0 до N . На вещественной оси R рассмотрим

сетку X

def


=

{x

j



}

j∈Z


,

X :


. . . < x

−1

< x

0

< x

1

< . . . ,

для которой lim

j→−∞


x

j

= −∞, lim



j→+∞

x

j



= +∞. Промежутки

I

j



= [x

j

, x



j+1

) называются элементарными промежутками сетки X.

Полиномиальный B-сплайн ω

B

j



второй степени (j ∈ Z) определя-

ется формулами

ω

B

j



(t) = (t − x

j

)



2

(x

j+1



− x

j

)



−1

(x

j+2



− x

j

)



−1

,

t ∈ I



j

,

ω



B

j

(t) = (x



j+2

− x


j

)

−1



(x

j+2


− x

j+1


)

−1

(x



j+3

− x


j+1

)

−1



×

× (x


j

−x

j+2



− x

j+3


+ x

j+1


)t

2

− 2(x



j+1

x

j



− x

j+2


x

j+3


)t+

+x

j



x

j+1


x

j+3


− x

j

x



j+2

x

j+3



+

+ x


j

x

j+1



x

j+2


− x

j+1


x

j+2


x

j+3


,

t ∈ I


j+1

,

ω



B

j

(t) = (t − x



j+3

)

2



(x

j+3


− x

j+2


)

−1

(x



j+3

− x


j+1

)

−1



,

t ∈ I


j+2

,

ω



B

j

(t) = 0 при t /



∈ [x

j

, x



j+3

] (так что supp ω

B

j

= [x



j

, x


j+3

]).


(1)

Пусть C


1

— линейное пространство непрерывно дифференци-

руемых функций u(t) с периодом 1, u(t+1) ≡ u(t), t ∈ R. Рассмотрим

1

Работа выполнена при финансовой поддержке РФФИ (гранты 04-01-00692 и



04-01-00026)

150


1-периодические функции Ω

B

j



(t), j = 0, N − 1, определяемые на про-

межутке [0, 1) формулами

B

j



(t) ≡ ω

B

j



(t), t ∈ [0,1),

j = 0, N − 3,

B

N −2



(t) ≡ ω

B

N −2



(t), t ∈ [x

N −2


, 1),

B



N −2

(t) ≡ ω


B

−2

(t),



t ∈ [0, x

N −2


),

B



N −1

(t) ≡ ω


B

N −1


(t), t ∈ [x

N −1


, 1),

B



N −1

(t) ≡ ω


B

−1

(t),



t ∈ [0, x

N −1


).

(2)


Отсюда следует, что Ω

B

j



∈ C

1

.



Каталог: bitstream
bitstream -> С. торайғыровтың публицистикасы
bitstream -> Қ атыстық сын есімдердің лексикалық тіркесімділігі
bitstream -> «Қазақ» газетіндегі ұлт-азаттық көтеріліс туралы мақалалардың маңызы
bitstream -> Бегімбай К. М. Дизайн және бейнелеу өнері салаларына
bitstream -> Грам м атические форм ы вы раж ения ком паративны Х отнош ений в научно-популярной лингвистической литературе
bitstream -> СҮттің Қоректік сапасы және сүттегі микроорганизмдер
bitstream -> Қазақ халқының ою-өрнектерінің Қолданылуы
bitstream -> Л. Н. Гумилев атындағы ЕҰу хабаршысы


Поделитесь с Вашими друзьями:
1   ...   14   15   16   17   18   19   20   21   ...   57


©emirsaba.org 2019
әкімшілігінің қараңыз

    Басты бет