Редактирование: Совершенствование алгоритмов численного моделирования в методе динамики частиц
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 1: | Строка 1: | ||
− | == | + | == Задача == |
+ | * ... | ||
+ | |||
+ | '''Цели''': | ||
* Создание программы на базе существующего кафедрального программного пакета по молекулярно-динамическому расчёту для | * Создание программы на базе существующего кафедрального программного пакета по молекулярно-динамическому расчёту для | ||
** Молекулярно-динамическое моделирования кипящего слоя | ** Молекулярно-динамическое моделирования кипящего слоя | ||
Строка 9: | Строка 12: | ||
'''Замечание''': | '''Замечание''': | ||
* Ниже приведены приёмы и методы, уже запрограммированные и активно отлаживающиеся в настоящее время. | * Ниже приведены приёмы и методы, уже запрограммированные и активно отлаживающиеся в настоящее время. | ||
− | |||
− | |||
− | |||
− | |||
− | |||
== Построение модифицированного метода Рунге-Кутты 4 порядка == | == Построение модифицированного метода Рунге-Кутты 4 порядка == | ||
Строка 116: | Строка 114: | ||
Проблема ошибок такого рода решается как правило заменой ЭВМ на более мощную. Однако это не единственный способ решения проблемы. Существует ряд алгоритмических приёмов, позволяющих добиться меньшей машинной погрешности, работая на прежней ЭВМ. Один из наиболее эффективных приёмов - обезразмеривание. | Проблема ошибок такого рода решается как правило заменой ЭВМ на более мощную. Однако это не единственный способ решения проблемы. Существует ряд алгоритмических приёмов, позволяющих добиться меньшей машинной погрешности, работая на прежней ЭВМ. Один из наиболее эффективных приёмов - обезразмеривание. | ||
− | Обезразмеривание как таковое возникло намного раньше ЭВМ и вовсе не из желания сократить машинную ошибку, о которой тогда и знать не знали. Первоначально обезразмеривание было исключительно математическим приёмом, позволяющим перейти от физической задачи (в которой присутствовали физические величины, обладающие смыслом и размерностью) к абстрактной математической задаче (в которой величины - просто некоторые отвлечённые параметры функции), что | + | Обезразмеривание как таковое возникло намного раньше ЭВМ и вовсе не из желания сократить машинную ошибку, о которой тогда и знать не знали. Первоначально обезразмеривание было исключительно математическим приёмом, позволяющим перейти от физической задачи (в которой присутствовали физические величины, обладающие смыслом и размерностью) к абстрактной математической задаче (в которой величины - просто некоторые отвлечённые параметры функции), что давао возможность применять формальные математическое методы решения уравнений, например метод малого параметра. |
Обезразмеривание позволяет перейти от физических величин, которые могут иметь очень большие (модуль Юнга) или очень малые (постоянная Планка) по модулю значения к их безразмерным аналогам, которые порядка единиц и в крайнем случае десятков. Понятно что второй случай более предпочтителен в вычислительном плане, ибо при оперировании данными безразмерными величинами происходит гораздо более медленное накопление ошибки. | Обезразмеривание позволяет перейти от физических величин, которые могут иметь очень большие (модуль Юнга) или очень малые (постоянная Планка) по модулю значения к их безразмерным аналогам, которые порядка единиц и в крайнем случае десятков. Понятно что второй случай более предпочтителен в вычислительном плане, ибо при оперировании данными безразмерными величинами происходит гораздо более медленное накопление ошибки. | ||
Строка 168: | Строка 166: | ||
<math> f(r_{ij} )=- \frac{\partial u}{\partial r_{ij}} </math> | <math> f(r_{ij} )=- \frac{\partial u}{\partial r_{ij}} </math> | ||
− | Если в качестве потенциала взаимодействия U взять потенциал | + | Если в качестве потенциала взаимодействия U взять потенциал Леннарда-Джонса |
<math> u(r_{ij} )=d[(\frac{a}{r_{ij}} )^{12}-2(\frac{a}{r_{ij}} )^6 ]=d[(\frac{1}{R_{ij}} )^{12}-2( \frac{1}{R_{ij}} )^6 ] </math> | <math> u(r_{ij} )=d[(\frac{a}{r_{ij}} )^{12}-2(\frac{a}{r_{ij}} )^6 ]=d[(\frac{1}{R_{ij}} )^{12}-2( \frac{1}{R_{ij}} )^6 ] </math> | ||
Строка 197: | Строка 195: | ||
... | ... | ||
− | |||
− | |||
− |