Периодические граничные условия
Материал из Department of Theoretical and Applied Mechanics
Виртуальная лаборатория>Периодические граничные условия
Содержание
Курсовой проект по механике дискретных сред
- разработчик Теницкая Татьяна
- руководитель Кузькин Виталий
Краткое описание
- Метод периодических граничных условий был разработан для решения задач теории жидкостей и плотных газов. Он состоит в том,что вокруг расчетной области строятся ее «образы» с актуальным положением частиц. И частицы «реальной» области взаимодействуют с частицами в «образе», а если частица пересекает границу расчетной области, она появляется с другой стороны.
- В теореме Нетер утверждается, что каждой непрерывной симметрии физической системы соответствует некоторый закон сохранения:
- однородности времени соответствует закон сохранения энергии,
- однородности пространства соответствует закон сохранения импульса,
- изотропии пространства соответствует закон сохранения момента импульса,
- калибровочной симметрии соответствует закон сохранения электрического заряда и т. д.
- Но для классической системы частиц с периодическими условиями сохранение момента импульса нарушается. Этот эффект наглядно проиллюстрирован в данной курсовой работе.
Цель проекта
- Визуализация системы частиц с периодическими граничными условиями.
- Построение графиков зависимости кинетического момента от времени для одной частицы, двух частиц, многих частиц.
Математическая модель
Граничные условия:
если , то
если , то
если , то
если , то
Где x и у - это координаты частицы, а w и h - ширина и длина окна соответственно.
Кинетический момент вычисляется по формуле:
Результаты
Текст программы на языке JavaScript (разработчик Теницкая Татьяна):
Файл "OneBall.js"
1 function main_particle_1() {
2
3 var ctx = canvas_particle_1.getContext("2d");
4 var w = canvas_particle_1.width; //длина
5 var h = canvas_particle_1.height; //высота
6
7
8 var w_1 = w/3;
9 var h_1 = h/3;
10
11 rx = new Array(0, w_1, 2*w_1, 0, w_1, 2*w_1, 0, w_1, 2*w_1);
12 ry = new Array(60, 60, 60, 60 + h_1, 60 + h_1, 60 + h_1, 60 + 2*h_1, 60 + 2*h_1, 60 + 2*h_1);
13
14 // координаты шара
15 var v0 = 1; // скорость шара
16 var r = 7; // радиус шара
17 var dt = 5;
18 var alfa = 45 / 180 * Math.PI;
19 var steps = 0;
20
21 function step()
22 {
23 tick();
24 draw();
25 }
26
27 var vGraph = new TM_graph( // определить график
28 "#vGraph", // на html-элементе #vGraph
29 250, // сколько шагов по оси "x" отображается
30 -1, 1,0.2); // мин. значение оси Y, макс. значение оси Y, шаг по оси Y
31
32 function tick()
33 {
34 for (i = 0; i < rx.length; i++)
35 {
36 vx = v0 * Math.cos(alfa);
37 vy = v0 * Math.sin(alfa);
38 steps+=1;
39
40 rx[i] += vx*dt;
41 ry[i] += vy*dt;
42
43
44 if (rx[i] >= w)
45 {
46 rx[i] = rx[i] - w;
47 }
48
49 if (rx[i] <= 0)
50 {
51 rx[i] = rx[i] + w;
52 }
53
54 if (ry[i] >= h)
55 {
56 ry[i] = ry[i] - h;
57 }
58
59 if (ry[i] <= 0)
60 {
61 ry[i] = ry[i] + h;
62 }
63 }
64
65 var L = (rx[8] * vy - ry[8] * vx)/240.4164;
66
67 console.log(L);
68
69 if (dt % 0.5 == 0) vGraph.graphIter(steps, L); // подать данные на график
70 }
71
72
73
74
75 function draw()
76 {
77 ctx.clearRect(0, 0, w , h ); // очистить экран
78 for (var i = 0; i < rx.length; i++)
79 {
80 var xS = rx[i];
81 var yS = ry[i];
82
83
84 ctx.beginPath();
85 ctx.fillStyle = "#00008B";
86 ctx.arc(xS, yS, r , 0, 2 * Math.PI, false);
87 ctx.closePath();
88 ctx.fill();
89
90 ctx.beginPath(); // начать рисование
91 ctx.fillStyle="#000000";
92 ctx.moveTo(w_1, 0); // переместить "карандаш" в точку
93 ctx.lineTo(w_1, h); // нарисовать "карандашом" линию до точки
94 ctx.stroke();
95
96 ctx.beginPath();
97 ctx.fillStyle="#000000";
98 ctx.moveTo(2*w_1, 0);
99 ctx.lineTo(2*w_1, h);
100 ctx.stroke();
101
102 ctx.beginPath();
103 ctx.fillStyle="#000000";
104 ctx.moveTo(0, h_1);
105 ctx.lineTo(w, h_1);
106 ctx.stroke();
107
108 ctx.beginPath();
109 ctx.fillStyle="#000000";
110 ctx.moveTo(0, 2*h_1);
111 ctx.lineTo(w, 2*h_1);
112 ctx.stroke();
113
114 }
115 }
116
117
118 setInterval(step, 1000/100); // функция step будет запускаться 100 раз в секунду (100 раз / 1000 мс)
119
120 }
Текст программы на языке JavaScript (разработчик Теницкая Татьяна):
Файл "TwoBalls.js"
1 function main_particle_1() {
2
3 var ctx_11 = canvas_particle_11.getContext("2d");
4 var ctx_12 = canvas_particle_12.getContext("2d");
5 var ctx_13 = canvas_particle_13.getContext("2d");
6 var ctx_21 = canvas_particle_21.getContext("2d");
7 var ctx_22 = canvas_particle_22.getContext("2d");
8 var ctx_23 = canvas_particle_23.getContext("2d");
9 var ctx_31 = canvas_particle_31.getContext("2d");
10 var ctx_32 = canvas_particle_32.getContext("2d");
11 var ctx_33 = canvas_particle_33.getContext("2d");
12
13
14 var width = canvas_particle_11.width;
15 var height = canvas_particle_11.height;
16
17 var rx_1, ry_1, rx_2, ry_2; // координаты шара
18 var v0 = 4, vx_1, vy_1, vx_2, vy_2;
19 var r = 7; // радиус шара
20 var dt_1 = 0; dt_2 = 0;
21 var x0_1 = width/2, x0_2 = width/2, y0_1 = 0, y0_2 = height;
22 var steps = 0;
23 var val = 1;
24 var alfa = 45/180*Math.PI;
25
26 function step()
27 {
28 tick();
29 draw();
30 }
31
32 var vGraph = new TM_graph( // определить график
33 "#vGraph", // на html-элементе #vGraph
34 250, // сколько шагов по оси "x" отображается
35 -1,1,0.2);
36
37 function tick()
38 { // вычисление новой позиции шар
39 vx_1 = v0 * Math.cos(alfa);
40 vy_1 = v0 * Math.sin(alfa);
41 vx_2 = -v0 * Math.cos(alfa);
42 vy_2 = -v0 * Math.sin(alfa);
43
44 steps+=1;
45 if (val==1)
46 {
47 rx_1 = vx_1*dt_1 + x0_1 + r - 2.7; //
48 ry_1 = vy_1*dt_1 + y0_1 + r - 2.7;
49 dt_1+= 1;
50 rx_2 = vx_2*dt_2 + x0_2 - r + 2.7; //
51 ry_2 = vy_2*dt_2 + y0_2 - r + 2.7;
52 dt_2+= 1;
53
54
55 }
56 else
57 {
58 vy_1 = -vy_1;
59 vx_1 = -vx_1;
60 rx_1 = vx_1*dt_1 + width - r + 2.7;
61 ry_1 = vy_1*dt_1 + x0_1 - r + 2.7;
62 dt_1+= 1;
63 vy_2 = -vy_2;
64 vx_2 = -vx_2;
65 rx_2 = vx_2*dt_2 - 2.7 + r; //
66 ry_2 = vy_2*dt_2 + x0_2 - 2.7 + r;
67 dt_2+= 1;
68 }
69
70 if (rx_2 < 0 || rx_1 - r> width)
71 {
72 //alert(rx_2);
73 val = 0;
74 dt_1 = 0;
75 dt_2 = 0;
76 }
77
78 if (ry_1 < 0 || ry_2 - r > height)
79 {
80 val = 1;
81 dt_1=0;
82 dt_2 = 0;
83 }
84
85
86 var L = ((rx_1 * vy_1 - ry_1 * vx_1)+(rx_2 * vy_2 - ry_2 * vx_2))/376.2;
87 console.log(L);
88
89 if (dt_1 % 0.5 == 0) vGraph.graphIter(steps, L); // подать данные на график
90
91
92
93 }
94
95 function draw() // рисование шара
96 {
97
98 ctx_11.fillStyle= "#DC143C" ; // цвет закраски
99 ctx_11.clearRect(0, 0, width, height); // очистить экран
100 ctx_11.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
101 ctx_11.fill();
102 ctx_11.beginPath();
103
104 ctx_11.fillStyle = "#00008B";
105 ctx_11.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
106 ctx_11.fill();
107 ctx_11.beginPath();
108
109 ctx_12.fillStyle = "#DC143C"; // цвет закраски
110 ctx_12.clearRect(0, 0, width, height); // очистить экран
111 ctx_12.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
112 ctx_12.fill();
113 ctx_12.beginPath();
114
115 ctx_12.fillStyle = "#00008B";
116 ctx_12.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
117 ctx_12.fill();
118 ctx_12.beginPath();
119
120 ctx_13.fillStyle = "#DC143C"; // цвет закраски
121 ctx_13.clearRect(0, 0, width, height); // очистить экран
122 ctx_13.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
123 ctx_13.fill();
124 ctx_13.beginPath();
125
126 ctx_13.fillStyle = "#00008B";
127 ctx_13.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
128 ctx_13.fill();
129 ctx_13.beginPath();
130
131 ctx_21.fillStyle = "#DC143C"; // цвет закраски
132 ctx_21.clearRect(0, 0, width, height); // очистить экран
133 ctx_21.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
134 ctx_21.fill();
135 ctx_21.beginPath();
136
137 ctx_21.fillStyle = "#00008B";
138 ctx_21.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
139 ctx_21.fill();
140 ctx_21.beginPath();
141
142 ctx_22.fillStyle = "#DC143C"; // цвет закраски
143 ctx_22.clearRect(0, 0, width, height); // очистить экран
144 ctx_22.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
145 ctx_22.fill();
146 ctx_22.beginPath();
147
148 ctx_22.fillStyle = "#00008B";
149 ctx_22.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
150 ctx_22.fill();
151 ctx_22.beginPath();
152
153 ctx_23.fillStyle = "#DC143C"; // цвет закраски
154 ctx_23.clearRect(0, 0, width, height); // очистить экран
155 ctx_23.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
156 ctx_23.fill();
157 ctx_23.beginPath();
158
159 ctx_23.fillStyle = "#00008B";
160 ctx_23.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
161 ctx_23.fill();
162 ctx_23.beginPath();
163
164 ctx_31.fillStyle = "#DC143C"; // цвет закраски
165 ctx_31.clearRect(0, 0, width, height); // очистить экран
166 ctx_31.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
167 ctx_31.fill();
168 ctx_31.beginPath();
169
170 ctx_31.fillStyle = "#00008B";
171 ctx_31.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
172 ctx_31.fill();
173 ctx_31.beginPath();
174
175 ctx_32.fillStyle = "#DC143C"; // цвет закраски
176 ctx_32.clearRect(0, 0, width, height); // очистить экран
177 ctx_32.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
178 ctx_32.fill();
179 ctx_32.beginPath();
180
181 ctx_32.fillStyle = "#00008B";
182 ctx_32.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
183 ctx_32.fill();
184 ctx_32.beginPath();
185
186 ctx_33.fillStyle = "#DC143C"; // цвет закраски
187 ctx_33.clearRect(0, 0, width, height); // очистить экран
188 ctx_33.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
189 ctx_33.fill();
190 ctx_33.beginPath();
191
192 ctx_33.fillStyle = "#00008B";
193 ctx_33.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
194 ctx_33.fill();
195 ctx_33.beginPath();
196 }
197
198 setInterval(step, 1000/30); // функция step будет запускаться 30 раз в секунду (30 раз / 1000 мс)
199
200 }
