Рассматривается функция Фокса с четырьмя параметрами, которая возникает в теории вырождающихся дифференциальных уравнений с частными производными дробного порядка. В терминах указанной функции были ранее записаны явные решения первой и второй краевых задач в полуполосе для уравнения с оператором Бесселя, действующим по пространственной переменной, и дробной производной по времени.
Для рассматриваемой функции в случае зависимости двух параметров из четырех в работе получена формула преобразования Лапласа, которая выражается через специальную функцию Макдональда. Также получены формулы интегральных преобразований, выражающиеся через обобщенную функцию Райта и более общую H-функцию Фокса.
Вспомогательным средством для доказательства полученных формул является интеграл Меллина–Барнса, с помощью которого записывается рассматриваемая специальная функция. Сходимость несобственных интегралов при этом следует из асимптотических оценок, также приведенных в работе.
Показано, что при частных значениях из формулы преобразования Лапласа следуют известные формулы преобразований экспоненциальной функции и функции Райта со степенными множителями.
Идентификаторы и классификаторы
Некоторые свойства функции (1), такие как представление через контурный интеграл, асимптотические свойства, формулы дифференцирования и интегрирования, рекуррентные соотношения, рассмотрены в работах [9–13]. Отметим при этом, что основные свойства функции (1), такие как, например, представление через контурный интеграл Меллина–Барнса, асимптотические свойства, разложение в степенные ряды, следуют из свойств более общей
-функции. Некоторые интегральные преобразования
-функции Фокса исследованы в работах [1–3].
Среди более поздних работ, посвященных интегральным преобразованиям с различными специальными функциями гипергеометрического типа в ядрах, отметим, например, работы [14–23].
В работах [24–26] развиты методы операторов преобразования для эллиптических и параболических уравнений с операторами Бесселя.
Список литературы
- Прудников А. П., Брычков Ю. А., Маричев О. И. Интегралы и ряды. Т. 3: Дополнительные главы. М.: Физматлит, 2003. 708 с.
- Kilbas A. A., Saigo M. H-Transforms. Theory and Applications / Analytical Methods and Special Functions. vol. 9. Boca Raton, FL: Chapman & Hall/CRC, 2004. xii+389 pp.
- Mathai A. M., Saxena R. K., Haubold H. J. The H-Function. Theory and Applications. Dordrecht: Springer, 2010. xiv+268 pp. DOI: https://doi.org/10.1007/978-1-4419-0916-9.
- Нахушев А. М. Дробное исчисление и его применение. М.: Физматлит, 2003. 271 с.
- Хуштова Ф. Г. Первая краевая задача в полуполосе для уравнения параболического типа с оператором Бесселя и производной Римана–Лиувилля // Матем. заметки, 2016. Т. 99, №6. С. 921–928. EDN: WGAQIH. DOI: https://doi.org/10.4213/mzm10759.
- Кузнецов Д. С. Специальные функции. М.: Высш. шк., 1962. 248 с.
- Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. Т. 1: Гипергеометрическая функция. Функции Лежандра. М.: Наука, 1965. 294 с.
- Хуштова Ф. Г. Вторая краевая задача в полуполосе для уравнения параболического типа с оператором Бесселя и частной производной Римана–Лиувилля // Матем. заметки, 2018. Т. 103, №3. С. 460–470. EDN: YSWJYG. DOI: https://doi.org/10.4213/mzm10986.
- Хуштова Ф. Г. Формулы дифференцирования и формула автотрансформации для одного частного случая функции Фокса // Доклады Адыгской (Черкесской) Международной академии наук, 2020. Т. 20, №4. С. 15–18. EDN: DKAMMT. DOI: https://doi.org/10.47928/1726-9946-2020-20-4-15-18.
- Хуштова Ф. Г. О некоторых свойствах одной специальной функции // Доклады Адыгской (Черкесской) Международной академии наук, 2022. Т. 22, №2. С. 34–40. EDN: LITQCZ. DOI: https://doi.org/10.47928/1726-9946-2022-22-2-34-40.
- Хуштова Ф. Г. Об интегральном представлении Меллина–Барнса одной специальной функции // Известия Кабардино-Балкарского научного центра РАН, 2022. №6. С. 19–27. EDN: TXVTRD. DOI: https://doi.org/10.35330/1991-6639-2022-6-110-19-27.
- Хуштова Ф. Г. Некоторые формулы дробного интегрирования от одной функции Фокса с четырьмя параметрами // Доклады Адыгской (Черкесской) Международной академии наук, 2022. Т. 22, №4. С. 29–38. EDN: NUYVKX. DOI: https://doi.org/10.47928/1726-9946-2022-22-4-29-38.
- Хуштова Ф. Г. К свойствам одной функции Фокса // Вестн. КРАУНЦ. Физ.-мат. науки, 2023. Т. 42, №1. С. 140–149. EDN: FXXPSA. DOI: https://doi.org/10.26117/2079-6641-2023-42-1-140-149.
- Ворошилов А. А. Дробное дифференцирование типа Эрдейи–Кобера H-функции Фокса // Вестн. Гродненск. гос. ун-та им. Янки Купалы. Сер. 2. Мат. Физ. Информ., вычисл. техн. и управл., 2012. Т. 2, №129. С. 11–20. EDN: TSVCDL.
- Авсиевич А. В., Авсиевич В. В. Преобразование Лапласа в системах автоматического управления дробного порядка // Наука и образование транспорту, 2013. №1. С. 195–199. EDN: SJGJKR.
- Авсиевич А. В. Преобразование Лапласа специальных функций Райта // Вестник транспорта Поволжья, 2013. №6. С. 50–52. EDN: RVKGWX.
- Заикина С. М. Обобщённое интегральное преобразование Лапласа и его применение к решению некоторых интегральных уравнений // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2014. №1. С. 19–24. EDN: TFGEOL. DOI: https://doi.org/10.14498/vsgtu1265.
- Qureshi M. I., Kabra D. K., Baboo M. S. Laplace transforms of multiple hypergeometric functions using Mellin–Barnes type contour integration // Asia Pac. J. Math., 2015. vol. 2, no. 2. pp. 94–107.
- Скоромник О. В. Интегральные преобразования с вырожденной гипергеометрической функцией Куммера и нормированной функцией Бесселя в ядрах и интегральные уравнения первого рода в пространстве суммируемых функций // Вестн. Полоцкого гос. ун-та. Сер. С. Фундамент. науки, 2016. №12. С. 104–110. EDN: XRFOMX.
- Karp D., Prilepkina E. G. Applications of the Stieltjes and Laplace transform representations of the hypergeometric functions // Integral Transforms Spec. Funct., 2017. vol. 28, no. 10. pp. 710–731. DOI: https://doi.org/10.1080/10652469.2017.1351964.
- Скоромник О. В. Двумерное интегральное преобразование с модифицированной H-функцией в пространстве суммируемых функций // Вестн. Полоцкого гос. ун-та. Сер. С. Фундамент. науки, 2018. №4. С. 187–193. EDN: UXBAMJ.
- Папкович М. В., Скоромник О. В. Двумерное интегральное преобразование с G-функцией Мейера в ядре в пространстве суммируемых функций // Вестн. Полоцкого гос. ун-та. Сер. С. Фундамент. науки, 2019. №4. С. 131–136. EDN: HFPVNO.
- Mohammed A. O., Rakha M. A., Awad M. M., Rathie A. K. On several new Laplace transforms of generalized hypergeometric functions
- // Bol. Soc. Parana. Mat. (3), 2021. vol. 39, no. 4. pp. 97–109. DOI: https://doi.org/10.5269/bspm.42207.
- Катрахов В. В., Ситник С. М. Метод операторов преобразования и краевые задачи для сингулярных эллиптических уравнений // Современная математика. Фундаментальные направления, 2018. Т. 64, №2. С. 211–426. EDN: AXVBAI. DOI: https://doi.org/10.22363/2413-3639-2018-64-2-211-426.
- Ситник С. М., Шишкина Э. Л. Метод операторов преобразования для дифференциальных уравнений с операторами Бесселя. М.: Физматлит, 2019. 221 с. EDN: YGUEZW.
- Transmutation Operators and Applications / Trends in Mathematics / eds. V. V. Kravchenko, S. M. Sitnik. Cham: Birkhäuser, 2020. xvii+686 pp. DOI: https://doi.org/10.1007/978-3-030-35914-0.
- Лебедев Н. Н. Специальные функции и их приложения. М.: Физматлит, 1963. 358 с.
- Маричев О. И. Метод вычисления интегралов от специальных функций (теория и таблицы формул). Мн.: Наука и техника, 1978. 312 с.
- Wright E. M. On the coefficients of power series having exponential singularities // J. Lond. Math. Soc., 1933. vol. s1-8, no. 1. pp. 71–79. DOI: https://doi.org/10.1112/jlms/s1-8.1.71.
- 30. Wright E. M. The generalized Bessel function of order greater than one // Q. J. Math., 1940. vol. os-11, no. 1. pp. 36–48. DOI: https://doi.org/10.1093/qmath/os-11.1.36.
- Прудников А. П., Брычков Ю. А., Маричев О. И. Интегралы и ряды. Т. 1: Элементарные функции. М.: Физматлит, 2002. 632 с.
- Псху А. В. Уравнения в частных производных дробного порядка. М.: Наука, 2005. 199 с. EDN: QJPLZX.
Выпуск
Другие статьи выпуска
Рассматривается напряженно-деформированное состояние и определяется время до разрушения составного растягиваемого стержня при ползучести в условии воздействия на него активной окружающей среды. Стержень состоит из трех частей, расположенных симметрично по толщине. Принято дополнительное условие: все части составного стержня жестко, без проскальзывания соединены между собой. Ползучесть каждой части стержня описывается степенной моделью с различными параметрами. Для определения времени до разрушения используется кинетическое уравнение, которое описывает накопление повреждений в процессе ползучести. Для каждой части стержня принят одинаковый вид кинетического уравнения, но накопление повреждений происходит под действием напряжений, различных для каждой части стержня. Влияние активной среды определяется диффузионным проникновением ее элементов в материал стержня. Используется приближенный метод решения уравнения диффузии, основанный на введении диффузионного фронта. Анализируется распределение напряжений во времени при условии проникновения активной среды в разные части стержня с различными коэффициентами диффузии. В результате исследований показано, что отношение констант в определяющих соотношениях ползучести частей стержня влияет на характер накопления поврежденности и распределения напряжений, а следовательно, имеется влияние на очередность разрушения частей составного стержня. С ростом показателей степеней в определяющих и кинетических соотношениях время до разрушения составного стержня увеличивается. Определена зависимость времени до разрушения от соотношения коэффициентов диффузии активной среды в части стержня.
Рассмотрена коэффициентная обратная задача определения геометрических параметров продольного прямоугольного паза по собственным частотам изгибных колебаний прямоугольного стержня. Предполагается, что паз проходит не по всей длине, а от определенной точки до правого конца. Для решения задачи стержень с продольным пазом моделируется в виде двух стержней, причем первый не имеет паз, а второй имеет.
В месте соединения используются условия сопряжения, в которых приравниваются величины прогибов, углов поворота, изгибающие моменты и перерезывающие силы. Исследованы закономерности поведения собственных частот изгибных колебаний при изменении длины паза. Предложен метод решения, позволяющий определять искомые параметры по конечному числу собственных значений изгибных колебаний. Показано, что решение однозначно в случае использования частотных спектров относительно взаимно перпендикулярных осей.
Рассматриваются математические модели конвекции–диффузии–реакции, которые относятся к моделям тепломассопереноса и применяются при исследовании природных и техногенных процессов. Для данного класса моделей актуальной является задача идентификации как параметров самой модели, так и входящих в нее граничных условий по результатам измерений значений искомой функции в отдельных точках рассматриваемой области. Задачу усложняет наличие неполных измерений, искаженных случайными помехами.
Решение заключается в разработке комбинированного двухэтапного метода идентификации, основанного на последовательном применении метода минимизации критерия идентификации безградиентного типа и рекуррентного метода оценивания неизвестных входных сигналов. Для применения указанных методов выполняется переход от исходной модели, описываемой уравнениями в частных производных, к дискретной линейной стохастической модели в пространстве состояний, в которой неизвестные граничные условия рассматриваются как неизвестные входные сигналы.
В результате построены новые дискретные линейные стохастические модели конвекции–диффузии–реакции для трех разных типов граничных условий. Предложена общая схема процесса параметрической идентификации, включающая двухэтапную идентификацию неизвестных параметров математической модели и идентификацию неизвестных граничных условий.
Для проверки работоспособности предложенного метода построены компьютерные модели конвекции–диффузии–реакции и выполнена реализация всех алгоритмов на языке MATLAB. Проведена серия вычислительных экспериментов, результаты которых показали, что разработанная двухэтапная комбинированная схема позволяет идентифицировать параметры исходной модели, значения функций, входящих в граничные условия, а также вычислить по неполным зашумленным измерениям оценки функции, описывающей процесс конвекции–диффузии–реакции.
Полученные результаты могут быть использованы не только при исследовании процессов тепломассопереноса, но также при решении задач идентификации параметров моделей дискретных стохастических систем с неизвестными входными сигналами и при наличии случайных помех.
Массоперенос в электромембранных системах в режимах интенсивного тока сопровождается возникновением дополнительных механизмов переноса, которые существенно влияют на эффективность их функционирования. Согласно современным представлениям для разбавленных растворов электролитов среди таких механизмов особенно важными являются электроконвекция и реакции диссоциации/рекомбинации молекул воды. Эти процессы оказывают противоположное действие на эффективность электромембранных технологий.
В исследованиях мембранных систем активно применяются математические модели, учитывающие влияние указанных механизмов, однако они обычно описывают только потенциодинамический режим, при котором устанавливается скачок потенциала в системе. Для интерпретации обширной базы экспериментальных данных по гальванодинамическому режиму (при фиксированной плотности тока) также необходимы инструменты теоретического анализа.
Цель данной работы заключается в разработке математической модели массопереноса в слое раствора электролита у ионообменной мембраны с учетом электроконвекции и диссоциации воды в гальванодинамическом режиме. Модель основана на системе связанных уравнений Нернста—Планка—Пуассона—Навье—Стокса, дополненной новым гальванодинамическим граничным условием для потенциала.
С использованием разработанной модели впервые рассчитаны хронопотенциограммы мембранной системы с учетом влияния как электроконвекции, так и реакции диссоциации/рекомбинации молекул воды. Результаты показали, что отношение концентрации продуктов диссоциации воды к концентрации ионов соли определяет баланс эффектов электроконвекции и диссоциации.
Рассмотрены следующие варианты соотношения эффектов электроконвекции и диссоциации молекул воды:
значимое влияние на массоперенос оказывает электроконвекция, в то время как влияние диссоциации воды минимально;
электроконвекция и диссоциация существенно влияют на процессы переноса:
образование дополнительных носителей заряда в результате диссоциации молекул воды снижает скачок потенциала в слое раствора электролита, что уменьшает интенсивность электроконвекции, в то время как развитие электроконвекции, в свою очередь, замедляет процесс диссоциации;
продукты интенсивной диссоциации молекул воды тормозят развитие электроконвекции.
Вопрос исследования гравитационного поля тел сложной формы (не относящийся к шарообразной) представляет большой интерес для геофизики, астрофизики, математической физики и других областей. Статья состоит из двух частей. В первой части представлен краткий литературный обзор различных методов расчета потенциала гравитационного поля однородного куба в рамках классической механики: получение аналитического решения; как частный случай задачи нахождения гравитационного поля полиэдра; методом конечных элементов; методом мультипольного разложения. Более подробно проанализирован метод расчета потенциала гравитационного поля однородного куба с помощью аналитического решения и мультипольного разложения. Во второй части статьи описан релятивистский случай гравитационного поля однородного куба в рамках постньютоновского формализма в первом и втором приближении. Данный метод расчета выбран по причине чрезвычайной сложности получение решения с помощью уравнений Эйнштейна. Ранее подобные задачи для тел с формой куба не рассматривались. Для решения задачи выбрана физическая модель — координатный равновесный куб, заполненный несжимаемой жидкостью с нулевой скоростью и постоянной плотностью. Получены релятивистские поправки для временной и пространственной координаты. Получен точный аналитический вид этих поправок для области вне куба, а также компоненты метрического тензора. Дано краткое сравнение полученных результатов для релятивистского случая с результатами классического ньютоновского случая. Для области внутри куба решение получено с помощью численных методов. Полученные результаты с достаточной точностью определяют параметры гравитационного поля для однородного куба, рассмотренного в рамках релятивистского подхода. Основное приложение этой задачи в рамках релятивистской физики относится к области математической физики (или, шире, математики).
Работа посвящена дальнейшему исследованию и построению конструктивных методов последовательной параметрической оптимизации неизвестных характеристик нестационарных процессов технологической теплофизики на компактном множестве непрерывных и непрерывно-дифференцируемых функций. Предложенная методика распространяет разработанный алгоритмически точный метод решения на многомерную постановку обратных задач технологической теплофизики и позволяет отыскать физически обоснованную идентифицируемую характеристику на последовательно сходящихся компактных множествах.
В качестве объекта исследования рассматривается двумерное осесимметричное тело канонической формы, где искомой функцией является сосредоточенная величина мощности внутренних теплоисточников. Задача сформулирована в равномерной метрике оценивания температурного отклонения расчетного состояния от экспериментального. В качестве математической модели рассматриваемого объекта используется его модальное описание, на основе которого проведена редукция исходной обратной задачи теплопроводности, сформулированной в экстремальной постановке, к задаче оптимального управления.
Использование предварительной параметризации искомой характеристики процесса приводит к ее представлению в форме кусочно-параболических функций, конкретизируемых в рамках выбранной структуры с помощью вектора параметров. Количество учитываемых параметров определяет точное представление идентифицируемой величины, а их значения отыскиваются в результате решения полученной задачи параметрической оптимизации. Для решения полученной задачи математического программирования относительно оптимальных значений вектора параметров используются, аналогично одномерному случаю, альтернансные свойства искомых экстремалей, в результате чего задача сводится к замкнутой системе соотношений.
Полученные результаты демонстрируют эффективность распространения конструктивного метода последовательной параметрической оптимизации, опробованного на одномерных обратных задачах теплопроводности, на решение двумерных задач с использованием их модального представления. Увеличение числа параметров решений, образующих кусочно-параболическую форму искомой зависимости, приводит к уменьшению погрешности восстановления как искомой сосредоточенной функции, так и пространственно-временного температурного поля во всей двумерной области определения пространственных переменных.
Представлен численный метод расчета полей остаточных напряжений в поверхностно упрочненном призматическом образце с несквозной V-образной трещиной, базирующийся на упругопластическом решении задачи. По полученным результатам проведен подробный анализ распределений остаточных напряжений вблизи дефекта по нескольким контурам. Определено, что при глубине трещины 0.3 мм практически все изучаемые компоненты остаточных напряжений сжатия имеют бóльшие (по модулю) значения, чем при глубине 0.1 мм, либо равные значения.
Предлагается новый алгоритм нахождения приближенных симметрий для дробно дифференциальных уравнений с производными типа Римана–Лиувилля и Герасимова–Капуто, порядок которых близок к целому. Алгоритм основан на разложении дробной производной в ряд по малому параметру, выделяемому из порядка дробного дифференцирования. В линейном приближении такое разложение содержит нелокальный интегро-дифференциальный оператор с логарифмическим ядром.
В результате исходное дробно-дифференциальное уравнение приближается интегро-дифференциальным уравнением с малым параметром, для которого могут быть найдены приближенные симметрии. Доказывается теорема о виде продолжения однопараметрической группы точечных преобразований на новую переменную, порождаемую нелокальным оператором, входящим в разложение дробной производной. Знание такого продолжения позволяет применить к рассматриваемому уравнению приближенный критерий инвариантности.
Предлагаемый алгоритм иллюстрируется на задаче нахождения приближенных симметрий для нелинейного дробно-дифференциального уравнения фильтрации субдиффузионного типа. Показано, что размерность алгебры приближенных симметрий такого уравнения оказывается существенно больше размерности алгебры точных симметрий, что открывает возможность построения большого числа приближенно инвариантных решений. Также на примере линейного дробно-дифференциального уравнения субдиффузии показывается, что алгоритм дает принципиальную возможность находить нелокальные приближенные симметрии определенного вида.
Финансовая система является важной составляющей в регулировании глобальных экономических процессов, поскольку обеспечение безопасности или контроль финансовой системы или рынка является ключом к стабилизации экономики.
Целью данного исследования является выяснение, насколько приближенные аналитические решения, полученные с помощью метода остаточного степенного ряда и метода разложения Эльзаки для дробной нелинейной финансовой модели, соответствуют экономической теории. Здесь понятие дробной производной используется в смысле производной Капуто.
Полученные численные результаты показывают, как приближенные решения реагируют на изменения процентной ставки, инвестиционного спроса и индекса цен. Оба метода показали результаты, согласующиеся с экономической теорией. Это означает, что исследователи могут использовать эти два метода для решения различных задач, связанных с дробными нелинейными моделями в финансовых системах.
Предлагается новый гибридный численный метод с использованием производной Капуто для решения нелинейного дробного уравнения Льенара — метод дифференциального преобразования Халуты. Доказана теорема сходимости данного метода при определенных условиях.
Метод дифференциального преобразования Халуты представляет собой полуаналитическую технику, объединяющую два мощных подхода: метод преобразования Халуты и метод дифференциального преобразования. Основное преимущество этого метода заключается в том, что он позволяет очень быстро находить решения и не требует линеаризации, возмущения или каких-либо других предположений. Предложенный метод подробно описан, а его эффективность продемонстрирована на двух числовых примерах. Результаты вычислений хорошо согласуются с точными решениями, что подтверждает надежность и эффективность предложенного подхода.
Издательство
- Издательство
- СамГТУ
- Регион
- Россия, Самара
- Почтовый адрес
- 443100, Самарская обл, г Самара, Октябрьский р-н, ул Молодогвардейская, д 244
- Юр. адрес
- 443100, Самарская обл, г Самара, Октябрьский р-н, ул Молодогвардейская, д 244
- ФИО
- Быков Дмитрий Евгеньевич (РЕКТОР)
- E-mail адрес
- rector@samgtu.ru
- Контактный телефон
- +7 (846) 2784300
- Сайт
- https://samgtu.ru/