Периодические граничные условия — различия между версиями

Материал из Department of Theoretical and Applied Mechanics
Перейти к: навигация, поиск
Строка 13: Строка 13:
 
Граничные условия:
 
Граничные условия:
  
если <math> x > w </math>
+
если <math> x > w </math>,
то <math>\{x = x - w;\} </math>
+
то <math> x = x - w </math>
  
<math> если (x < 0) </math>
+
если <math> x < 0 </math>,
<math>\ то {x = x + w;\} </math>
+
то <math>x = x + w </math>
  
 
<math> если (y > h) </math>
 
<math> если (y > h) </math>

Версия 11:48, 22 декабря 2015

Виртуальная лаборатория>Периодические граничные условия

Краткое описание
Метод периодических граничных условий был разработан для решения задач теории жидкостей и плотных газов.Он состоит в том,что вокруг расчетной области строятся ее «образы» с актуальным положением частиц.И частицы «реальной» области взаимодействуют с частицами в «образе». А если частица пересекает границу расчетной области, она появляется с другой стороны.В теореме Нетер утверждается, что каждой непрерывной симметрии физической системы соответствует некоторый закон сохранения:
однородности времени соответствует закон сохранения энергии,
однородности пространства соответствует закон сохранения импульса,
изотропии пространства соответствует закон сохранения момента импульса,
калибровочной симметрии соответствует закон сохранения электрического заряда и т. д.
Но для классической системы частиц с периодическими условиями сохранение момента импульса нарушается.
Цель работы:

Визуализация системы частиц с периодическими граничными условиями.Построение графиков зависимости кинетического момента от времени для одной частицы, двух частиц, многих частиц.
Граничные условия:

если [math] x \gt w [/math], то [math] x = x - w [/math]

если [math] x \lt 0 [/math], то [math]x = x + w [/math]

[math] если (y \gt h) [/math] [math]\{y = y - h;\} [/math]

[math] еслт (y \lt 0) [/math] [math]\{y = y + h;\} [/math]

Где x и у - это координаты частицы, а w и h - размеры окна.

Так же высчитывается кинетический момент по формуле:

[math]L(t) = \sum_{i\in\wedge(t)} r_i\times mV_i [/math]