В работе численно исследуется долговременная асимптотика вероятности выживания частиц, диффундирующих в двумерных областях со случайно расположенными поглощающими ловушками. Исследованы две области: фрактальная, построенная на основе кривой Коха, и обычная прямоугольная с неравномерной плотностью посева ловушек, имитирующей проекцию фрактала. Методом Монте-Карло смоделировано движение ансамбля частиц, проанализирована зависимость времени жизни системы и финального числа выживших частиц от количества ловушек. Установлено, что во фрактальной области захват частиц происходит медленнее, что приводит к значительно более высокому уровню выживаемости на больших временах. Показано, что для обеих областей зависимость времени жизни от количества ловушек наилучшим образом аппроксимируется гиперболической функцией, что согласуется с теоретическими результатами.
Идентификаторы и классификаторы
- SCI
- Математика
Исследование диффузионных процессов в прикладных задачах часто приводит к необходимости работать с достаточно сложными двумерными и трёхмерными областями [7–9]. В качестве модели для описания объектов со сложной структурой можно использовать фракталы и их производные. В случае исследования скорости захвата частиц на ловушки в фрактальных областях возникла идея использовать аналогичные области с более простой структурой за счёт сгущения посева ловушек на некоторой части области.
Если у вас возникли вопросы или появились предложения по содержанию статьи, пожалуйста, направляйте их в рамках данной темы.
Список литературы
1. Montrol E. W., Weiss G. H. Random Walks on Lattices. Journal of Mathematical Physics. 1965; 6: 167-175.
2. Балагуров Б. Я., Вакс В. Г. Теория диффузных фазовых переходов // Журнал экспериментальной и теоретической физики. 1973. Т. 65. C. 1600-1604.
3. Soghra Safaverdi, Gerard T. Barkema, Eddy Kunnen, Adam M. Urbanowicz, Christian Maes. Saturation of front propagation in a reaction-diffusion process describing plasma damage in porous low-k materials. Physical Review B. 2011; 83: 245320. DOI: 10.1103/PhysRevB.83.245320 EDN: OLSPEH
4. Kunnen E., Barkema G. T., Maes C., Shamiryan D., Urbanowicz A., Struyf A., Baklanov M. R.Integrated diffusion-recombination model for describing the logarithmic time dependence of plasma damage in porous low-k materials. Microelectron. Eng. 2011; 88: 631-634. EDN: ONJQOX
5. Архинчеев В. Е. Влияние дрейфа на временную асимптотику вероятности выживания частиц в средах с поглощающими ловушками // Журнал экспериментальной и теоретической физики. 2017. Т. 151, вып. 2. С. 322-325. DOI: 10.7868/S0044451017020109 EDN: YJHBXZ
6. Maex K., Baklanov M.R., Shamiryan D., Iacopi F., Brongersma S.H., Yanovitskaya Z. Sh. Low Dielectric Constant Materials for Micro Electronics. Journal of Applied Physics. 2003; 93: 8793-8841. DOI: 10.1063/1.1567460 EDN: LHXFVH
7. Rasadujjaman M., Wang X., Wang Y., Zhang J., Arkhincheev V. E., Baklanov M. R. Analytical Study of Porous Organosilicate Glass Films Prepared from Mixtures of 1,3,5- and 1,3-Alkoxysilylbenzenes. Materials. 2021; 14: 1881. EDN: HYZBEA
8. Bini F., Pica A., Novelli S., Pecci R., Bedini R., Marinozzi A., Marinozzi F. 3D FEM Model to Simulate Brownian Motion inside Trabecular Tissue from Human Femoral Head.Comput. Methods Biomech. Biomed. Eng. Imaging Vis. 2022; 10: 500507.
9. Arkhincheev V. E., Kunnen E., Baklanov M. R. Active species in porous media: random walk and capture in traps. Microelectronic Engineering. 2011; 88; 5: 686689. DOI: 10.1016/j.mee.2010.08.028 EDN: SKUGZF
10. Рязанов Г. В. Случайные блуждания на плоской решетке с ловушками // Теоретическая и математическая физика. 1972. Т. 10. С. 271-277. DOI: 10.1007/BF01090731
11. Балагуров Б. Я., Вакс В. Г. Случайные блуждания частицы по решеткам с ловушками // Журнал экспериментальной и теоретической физики. 1973. Т. 65. C. 1939-1943.
12. Численное моделирование диффузии в средах с ловушками: статистика в одномерном случае / В. Е. Архинчеев, Б. В. Хабитуев, Д. Ф. Дерюгин [и др.] // Вестник Бурятского государственного университета. Математика, информатика. 2024. № 3. С. 31-44. DOI: 10.18101/2304-5728-2024-3-31-44
13. Архинчеев В. Е., Хабитуев Б. В., Мальцев С. П. Влияние дрейфа на временную асимптотику вероятности выживания частиц в средах с поглощающими ловушками: численное моделирование // Журнал экспериментальной и теоретической физики. 2025. Т. 167, № 5. С. 703-710. DOI: 10.31857/S0044451025050086 EDN: VCTHAL
14. Arkhincheev V. E., Khabituev B. V., Maltsev S. P. Numerical Simulation of Capture of Diffusing Particles in Porous Media.Computation. 2025; 13; 4: 82. DOI: 10.3390/computation13040082 EDN: ZZRXMM
Выпуск
Другие статьи выпуска
Результаты, полученные с помощью модели делокализованных атомов, впервые установили в области стеклования линейную температурную зависимость параметра ангармонизма для ряда стеклообразных материалов. Суммарный эффект теплового движения, обусловленного ангармонизмом, заключается в возникновении растягивающего теплового давления на межатомные связи. Внутреннее давление, являющееся основной компонентой теплового давления в жидкостях и твердых телах, представляет собой упругую реакцию решетки на внешние воздействия и рассчитывается как скорость, с которой внутренняя энергия вещества меняется при изменении его объема при неизменной температуре. Используя область упругой деформации в ее верхнем предельном значении, на основе предложенной модели выявлена максимальная величина теплового давления. Выявленное значение максимального теплового давления приводит к предельному растяжению межатомной (межмолекулярной) связи, которое обратно пропорционально параметру ангармонизма. Из уравнений состояния и модели делокализованных атомов следует линейная корреляция между обратным значением квадрата параметра ангармонизма и температуры стеклования стеклообразных твердых тел.
В работе представлены результаты компьютерного моделирования процессов охлаждения систем частиц ксенона при различных термодинамических параметрах. Исследование фокусируется на выявлении взаимосвязи между условиями охлаждения (скорость изменения температуры и приложенное давление) и структурными характеристиками образующихся твёрдых фаз. В ходе моделирования варьировались два ключевых параметра: скорость охлаждения - рассмотрены три порядка величин: 10^8, 10^9 и 10^10 К/с; давление - исследованы значения 0.5, 1.0, 2.0 и 4.0 МПа. Основной научный результат работы заключается в установлении чёткой зависимости типа образующейся структуры от кинетических и термодинамических условий процесса. Показано, что при высоких скоростях охлаждения (порядка 10^10 К/с) формируется стеклообразное твёрдое тело - неупорядоченная структура, в которой частицы «замораживаются» в положении, характерном для жидкого состояния, не успевая выстроить кристаллическую решётку; при меньших скоростях охлаждения (10^8 - 10^9 К/с) происходит образование систем с явными структурными признаками кристалла - упорядоченной решёткой, характерной для твёрдого состояния. Кроме того, исследование продемонстрировало существенное влияние давления на структурные особенности получаемых твёрдых фаз. Изменение давления от 0.5 до 4.0 МПа модифицирует параметры образующейся кристаллической решётки, влияя на межатомные расстояния и координацию частиц.
В работе рассматривается задача наилучшего распределения пациентов между медицинскими учреждениями с учетом их вместимости. Негладкая целевая функция этой задачи представляет собой максимальную нагрузку среди всех медучреждений, которую требуется минимизировать. Для поиска решения доказываются достаточные условия оптимальности, базирующиеся на свойстве неубывания целевой функции. На этой основе построен новый метод потенциальных значений прямого-обратного хода, позволяющий находить решения, которые удовлетворяют достаточным условиям оптимальности. С целью тестирования разработанного метода проведен вычислительный эксперимент на тестовых задачах, моделирующих эпидемии в крупных городах. Результаты эксперимента показывают, что данный метод успешно находит решения поставленных задач за меньшее время, чем ранее разработанный метод потенциальных значений прямого хода и оказывается эффективнее программных пакетов, использующих традиционные подходы глобальной оптимальности.
В работе рассматривается задача минимизации вырожденного квадратичного функционала на траекториях квазилинейной управляемой системы. Особенность задачи состоит в том, что решение задачи существует в пространстве обобщенных функций. Для исходной задачи строится вспомогательная задача. Её решение существует в пространстве абсолютно непрерывных функций. Для ее решения применяются методы разложения по параметру и принцип максимума Л. С. Понтрягина. В результате конечного числа итераций строится субоптимальное управление для вспомогательной задачи. С помощью этого управления конструируется субоптимальное управление для исходной задачи. Это управление имеет импульсные составляющие, сосредоточенные в начальный и конечный моменты. Между этими моментами управление является непрерывной функцией. Построенное управление является программным. С его помощью можно построить позиционное управление, которое в случае воздействия возмущений на систему будет порождать импульсно-скользящий режим.
В настоящей статье рассмотрены взаимосвязанные линейные системы обыкновенных дифференциальных и алгебраических уравнений, которые записаны в векторно-матричном виде. Такие постановки достаточно часто встречаются в важных прикладных задачах энергетики, кинетической химии, биологии, моделировании многозвеньевых систем и других областей. В отечественной и зарубежной литературе их принято называть дифференциально-алгебраическими уравнениями. Предполагается, что число уравнений в системе меньше, чем размерность искомой вектор-функции. Используя результаты теории недоопределенных систем линейных алгебраических уравнений, дано понятие нормального решения для выделенного класса задач. Для их численного решения предложен вариант коллокационно-вариационных разностных схем, основанных на решении задачи математического (квадратичного) программирования специального вида. Приведены результаты численных расчетов нескольких модельных примеров, которые иллюстрируют работоспособность данного подхода.
В статье предложено обоснование к методу Рауса - Ляпунова для нахождения стационарных движений механических автономных консервативных систем. Изложение вначале проведено для автономных гамильтоновых систем на которые наложены дополнительные условия в виде некоторых алгебраических равенств. Их применение осуществляется известным способом Лагранжа неявного учитывания отмеченных условий, который состоит в формировании связки из гамильтониана исходной системы и алгебраической суммы указанных условий, умноженных на неопределенные вещественные множители. Последние в дальнейшем рассматриваются наравне с переменными системы. При новом взгляде на интегрирование гамильтоновых систем как бесконечно малое контактное преобразование последнее по существу свелось к нахождению положений равновесия преобразованной системы с учетом введенных ограничений. Принимая во внимание то обстоятельство, что в консервативных системах гамильтониан состоит из суммы кинетической и потенциальной энергий, этот подход в точности переносится на механические автономные консервативные системы общего вида. Для этого формируется связка, состоящая из алгебраической суммы первых интегралов исходной системы уравнений движения с неопределенными вещественными множителями, участвующими наравне с фазовыми переменными. Условие постоянства решений для стационарного движения сводит задачу их нахождения к обращению в нуль всех частных производных по фазовым переменным и множителям Лагранжа как в гамильтоновых системах, так и механических автономных консервативных системах общего вида. В общем случае установлено несовпадение найденных методом Рауса - Ляпунова стационарных движений с решениями положений равновесия изучаемой системы.
В конце прошлого столетия Дж. Б. Майлз доказал, что если целая функция бесконечного порядка имеет нули, распределенные на конечной системе лучей, то её нижний порядок также равен бесконечности. К. Г. Малютин, М. В. Кабанко и Т. В. Шевцова (2022) распространили результат Майлза на истинно аналитические функции бесконечного порядка относительно классической функции роста r на верхней полуплоскости. В данной работе мы распространяем результат К. Г. Малютина, М. В. Кабанко и Т. В. Шевцовой на пространство истинно аналитических функций на верхней полуплоскости бесконечного порядка относительно модельной функции роста M. Понятие модельной функции роста M, введенное Б. Н. Хабибуллиным, охватывает широкий класс функций. В частности, функции, определяемые модельной, могут иметь бесконечный порядок, а также нулевой порядок в классическом смысле.
В данной работе рассмотрена задача Коши для системы уравнений с правыми частями специального вида и постоянными коэффициентами. С помощью метода дополнительного аргумента определены достаточные условия нелокальной разрешимости задачи Коши для системы уравнений с правыми частями специального вида и постоянными коэффициентами. Доказательство нелокальной разрешимости задачи Коши для системы уравнений с правыми частями специального вида и постоянными коэффициентами опирается на глобальные оценки. Приведен пример задачи Коши для системы уравнений с правыми частями специального вида и постоянными коэффициентами, которая имеет единственное нелокальное решение. Пример показывает, что существует задача Коши для системы уравнений с правыми частями специального вида и постоянными коэффициентами, которая имеет единственное нелокальное решение.
Издательство
- Издательство
- БГУ
- Регион
- Россия, Улан-Удэ
- Почтовый адрес
- 670000, Дальневосточный федеральный округ, Республика Бурятия, г. Улан-Удэ, ул. Смолина, 24 "а"
- Юр. адрес
- 670000, Дальневосточный федеральный округ, Республика Бурятия, г. Улан-Удэ, ул. Смолина, 24 "а"
- ФИО
- Дамдинов Алдар Валерьевич (Руководитель)
- E-mail адрес
- univer@bsu.ru
- Контактный телефон
- +7 (301) 2297170
- Сайт
- https://bsu.ru