УДК 517.929+517.95
DOI: 10.18698/23 09-3684-2024-3 -140167
Аналитические решения нелинейных уравнений с запаздыванием, используемых при математическом моделировании процессов переноса
© В.Г. Сорокин
Институт проблем механики им. А.Ю. Ишлинского РАН, Москва, 119526, Россия
Рассматриваются нелинейные уравнения переноса с постоянным запаздыванием. Во введении дается краткий обзор публикаций, в которых исследуются математические модели переноса с запаздыванием и разрабатываются численные методы решения соответствующих задач. В основных разделах статьи описывается более сорока уравнений переноса с постоянным запаздыванием и различными коэффициентами переноса, которые допускают точные аналитические решения. Кинетические функции всех рассматриваемых уравнений содержат свободные параметры или произвольные функции. Получены решения с аддитивным, мультипликативным, обобщенным и функциональным разделением переменных, а также решения типа бегущей волны и автомодельные решения. Многие решения выражаются через элементарные функции. Для некоторых типов уравнений сформулированы теоремы о «размножении» решений. Описанные уравнения и их решения могут быть использованы для оценки точности численных методов интегрирования соответствующих нелинейных задач переноса с запаздыванием.
Ключевые слова: нелинейные уравнения переноса с запаздыванием, уравнения в частных производных с запаздыванием, аналитические решения, решения с обобщенным и функциональным разделением переменных, автомодельные решения
Введение. Уравнения переноса с запаздыванием применяются для описания конвективного переноса некоторой субстанции, когда рассматриваемый процесс, явление или система обладают свойством последействия, то есть зависимости текущего состояния от значений искомых величин в некоторый прошлый момент времени. Простейшее уравнение переноса с постоянным запаздыванием имеет вид
и( + аих = / (и, н), н = и(х, ?-г), (1)
где и = и( х, t) — искомая функция, н = н(х, t) — функция с запаздыванием, /(и, н) — кинетическая функция, а >0 — коэффициент переноса, имеющий физический смысл скорости, г >0 — постоянное запаздывание.
Замечание 1. В литературе часто можно встретить уравнения переноса с запаздыванием по х (см., например, [1-3]):
и( + аих = / (и, н), н = и( х -г, {). (2)
Так как уравнения (1) и (2) симметричны относительно аргумен-
тов х и t (с точностью до переобозначения константы а и кинетической функции / ), в следующих разделах для простоты и определенности рассматриваются только уравнения вида (1) с запаздыванием по t.
Помимо уравнений переноса с одним постоянным запаздыванием (1), в литературе также встречаются уравнения с несколькими постоянными запаздываниями, с пропорциональным запаздыванием, одновременно с пропорциональным запаздыванием по одной переменной и пропорциональным запаздыванием по другой переменной, и др. Рассмотрим далее некоторые примеры уравнений переноса с такими запаздываниями, которые применяются в биологии и медицине.
В [4] с помощью линейного уравнения переноса с двумя запаздываниями по х (формально, с запаздыванием и с опережением) описывается процесс возбуждения нейронов. Рассматриваемое уравнение имеет вид:
щ + а(х, t)их = Ь(х, t)(щ ~ и) + с(х, t)(щ - и), (3)
где и = и(х, t) — вероятность того, что в момент времени t деполяризация мембраны нейрона меньше х; щ = и( х — ^, t), щ = и( х + г2, t); а( х, t), Ь( х, ?), с( х, t) — достаточно гладкие функции своих аргументов. Положительные постоянные тх и г2 позволяют учесть увеличение
или уменьшение вероятности деполяризации нейрона в случае, соответственно, возбуждающего или тормозного импульса. Статья [4] посвящена разработке разностной схемы для решения соответствующей начально-краевой задачи, анализу ее сходимости и устойчивости.
Линейное уравнение с пропорциональным запаздыванием по х рассматривается в [5] для описания роста и деления клеток, структурированных по размеру:
щ + аих = р2Ьщ — Ьи — /ли, щ = и(рх, t),
где и = и( х, t) — распределение плотности клеток размером х в момент времени t. Данное уравнение моделирует динамику популяции клеток, растущих со скоростью а > 0 и делящихся со скоростью Ь >0 на р >1 дочерних клеток размером х / р (обычно р = 2 ). Предполагается также, что клетки гибнут со скоростью л > 0 . Путем сведения рассматриваемой начально-краевой задачи к последовательности более простых задач Коши в [5] получено аналитическое решение начально-краевой задачи для произвольных начальных данных.
В [6] рассматривается модель популяции клеток с одновременным учетом пролиферации и созревания, которая основана на нелинейном уравнении переноса с постоянным запаздыванием по времени и пропорциональным по х:
ut + = cw(1 — w) — bu, w = u(px, t—r), (4)
где u = u(x, t) — плотность пролиферирующих клеток зрелости x в момент времени t, a(x) >0 — скорость созревания (обычно a(x) = rx, r >0 — некоторая константа), b >0 и c >0 — некоторые константы, r >0 — время запаздывания, p >0 — коэффициент пропорциональности. Переменная x называется «переменной созревания» (от англ. "maturation variable") и отвечает за уровень зрелости клетки. В данной модели учитывается, что динамика текущей популяции клеток зависит как от популяции в некоторый момент времени в прошлом, так и от популяции клеток другого уровня зрелости. В [6] изучается существование и локальная устойчивость решений различного рода для различных начальных данных и значений параметров уравнения. Похожие модели рассматриваются также в [7, 8]. Более сложное нелинейное уравнение с произвольным аргументом по x и постоянным запаздыванием по времени:
u + a(x)ux = f (t, u, w), w = u(^(x), t - r),
где 0 < x ) < x, изучается в [9].
Замечание 2. Модели, основанные на одиночных уравнениях переноса без запаздывания, рассматриваются в [10], а на системах таких уравнений — в [11].
Замечание 3. Некоторые модели, основанные на реакционно-диффузионных уравнениях с запаздыванием (и других уравнениях математической физики второго порядка с запаздыванием), можно найти в [12].
Большой объем литературы посвящен разработке и исследованиям численных методов решения уравнений переноса с запаздыванием. Численные методы для уравнений переноса с постоянным запаздыванием рассматриваются, например, в [1-3]. Одношаговые разностные методы для уравнений переноса с запаздыванием по обеим переменным разрабатываются в [13]. В [14] рассматриваются разностные схемы для линейного уравнения переноса с постоянными и переменными запаздыванием и опережением по x вида (3). Предложены две разностные схемы, одна из которых условно устойчива, а другая — безусловно устойчива. В [15] разрабатывается численный метод для уравнения переноса с пропорциональным запаздыванием по x и постоянным запаздыванием по времени вида (4). Монография [16] посвящена разностным численным методам решения уравнений в частных производных с запаздыванием, в том числе уравнениям переноса с запаздыванием. Предложен подход, заключающийся в использовании аналогов методов, известных для уравнений без запаздывания, и интерполяции с заданными свойствами для учета запаздывания.
Точные аналитические решения уравнений в частных производных с запаздыванием используются в формулировках тестовых задач, применяемых для проверки адекватности и оценки точности численных и приближенных аналитических методов интегрирования таких уравнений. Особенно полезны уравнения, которые содержат произвольные функции и допускают точные решения со свободными параметрами. Такие уравнения и точные решения позволяют легко формулировать большой набор тестовых задач, которые в том числе дают возможность наблюдать и сравнивать работу методов в обычных и особых случаях, например, в устойчивых и неустойчивых областях. Изучение точных решений также способствует лучшему пониманию рассматриваемых процессов и явлений.
Полная групповая классификация уравнений математической физики с запаздыванием реакционно-диффузионного типа и типа Клейна — Гордона проводилась в [17-19]. В результате было получено некоторое количество уравнений, допускающих точные аналитические решения, однако общее число таких уравнений было сравнительно небольшим. Более продуктивным оказался метод функциональных связей, разработанный в [20, 21]. Он позволяет строить точные решения для нелинейных уравнений в частных производных с запаздыванием, содержащих произвольные функции. С его помощью были получены точные решения для нескольких сотен уравнений в частных производных с запаздыванием (см., например, [22-27]). В [28, 29] предложен метод структурной аналогии решений, который позволяет строить точные решения уравнений в частных производных с запаздыванием на основе точных решений более простых уравнений без запаздывания. Точные решения уравнений математической физики с запаздыванием различного вида (в том числе с пространственно анизотропным временным запаздыванием) получены этим методом в [28-33]. Результаты последних лет по точным решениям уравнений математической физики второго порядка с постоянным, пропорциональным и переменным запаздыванием обобщены и структурированно изложены в монографии [12].
Все указанные выше публикации и монографии содержат точные решения уравнений в частных производных второго и старших порядков с запаздыванием. Практически отсутствуют публикации по точным аналитическим решениям нелинейных уравнений первого порядка с запаздыванием. На данный момент существует только одна работа [34]. В ней проводится полная групповая классификация уравнения переноса с запаздыванием вида
и( + иих = ¥ (и, н), н = и(х, t -г). Получено четыре уравнения, которые допускают точные решения:
одно уравнение с произвольной функцией, допускающее решение типа бегущей волны, одно линейное уравнение и два уравнения со степенными нелинейностями. Полученные решения содержат произвольные функции и свободные параметры.
В статье далее термин «точное аналитическое решение нелинейного уравнения переноса с запаздыванием», или просто «точное решение», применяется в случаях, когда решение выражается:
1. в элементарных функциях,
2. через неопределенные интегралы,
3. через решения уравнений переноса без запаздывания,
4. через решения ОДУ (или систем ОДУ) с запаздыванием или без запаздывания.
Допустимы также комбинации решений из пп. (1) - (4).
Замечание 4. Точным решениям уравнений математической физики первого, второго и более старших порядков без запаздывания посвящена книга [35], которая содержит более четырех с половиной тысяч уравнений, допускающих точные решения.
Далее рассматриваются нелинейные уравнения переноса с запаздыванием:
и( + аих = ¥ (и, н), н = и(х, t -г),
где и = и(х, t), а >0 — постоянный (а в некоторых случаях переменный а = а(х) > 0 , но не зависящий от и ) коэффициент переноса, ¥ (и, н) — кинетическая функция, г >0 — постоянное запаздывание. Рассмотрено восемнадцать таких уравнений. Для трех уравнений сформулированы теоремы о размножении решений.
Рассматриваются также нелинейные уравнения переноса с запаздыванием и зависящим от искомой функции коэффициентом переноса:
и + Н (и)их = ¥ (и, н), н = и( х, t - г),
где Н (и) — коэффициент переноса степенного, экспоненциального или логарифмического вида. Рассмотрено двадцать семь таких уравнений
Точные аналитические решения всех рассмотренных уравнений содержат свободные параметры; сами уравнения — свободные параметры и/или произвольные функции. Данное обстоятельство делает рассматриваемые уравнения и их точные решения подходящими для составления обширного каталога тестовых задач, который можно затем использовать для всестороннего тестирования соответствующих численных методов. Описанные точные решения могут быть также использованы при изучении качественных особенностей рассматриваемых и родственных уравнений с запаздыванием.
Точные аналитические решения уравнений переноса с запаздыванием и постоянным коэффициентом переноса.
Уравнение 1. Нелинейное уравнение переноса с запаздыванием
и + аих = и(Ъ 1п и + с 1п w + ё), w = и(х, t-г), (5)
допускает точное решение с функциональным разделением переменных
и (х, t) = exp
%п¥п (t)
п=1
где N — любое натуральное число, а функции уп = у ^) описываются нелинейной системой ОДУ с постоянным запаздыванием:
у'0 = Ъу - ау + су + ё, ¥[ = Ъуг ~^ + \)ау/м +су/г, ¡ = 1,...,М-1, у/'м=Ьу/м+су7м.
Здесь = цг1 (V - т), 1 = 0,..., N.
Уравнение 2. Рассмотрим нелинейное уравнение переноса с запаздыванием
и +аих = и/ ^ / и), w = и( х, t -г). (6)
1. Уравнение (6) допускает точное решение с мультипликативным разделением переменных
и = еЛху(0, (7)
где Л — произвольная постоянная, а функция у = у/(() удовлетворяет ОДУ с запаздыванием
у'+ аЛу-у / (у /у) = 0, у = у(t -г).
2. Уравнение (6) имеет также решение, более общее, чем (7):
и = exp(аx + р )в( г), г = Лх + р,
где а, Р, Л, у — произвольные постоянные , а функция в = в(г) удовлетворяет ОДУ с запаздыванием
(аЛ + у)в' + (аа + Р)в-в/(е'Ртв /в) = 0, в = в(г -о), о = уг.
3. Уравнение (6) допускает точное решение вида
и = еЛу(х, t), у(х, t) = у(х, t - г),
где Л — произвольная постоянная, а у(х, t) — г - периодическое решение линейной задачи для уравнения переноса:
V + аух = х, ^) = х, £ - г), (8)
где Ъ = / (е~Яг) -Я.
Для удобства обозначим общее решение задачи (8) как
V = V(х,Ъ) . Оно имеет вид:
ад
V (х, Г; Ъ) = еЪх/а А 0С8(Д/ - ух) + Вп вт(Д/ - ух)], (9)
п=0
2«п 2«п
Дп =-, 7 =-, (10)
г аг
где Аи, В — произвольные постоянные, для которых ряд (9)-(10) и его производные (V) и (V ) сходятся (сходимость можно обеспечить, если положить Аи = В =0 при п > Ы, где N — произвольное
натуральное число).
4. Уравнение (6) допускает также точное решение вида
и = eЯív(х, £), v(х, £) = х, £ - г),
где Я — произвольная постоянная, а v(х, £) — г - апериодическое решение линейной задачи для уравнения переноса:
V + = Ы, v(х, £) = -v(х, £ - г), (11)
где Ъ = / (-е-Яг) -Я.
Общее решение задачи (11), которое для удобства обозначим
V = V(х,Ъ) , имеет вид:
ад
VI (х, Ъ) = еЪх/а £[А совД - ух) + Вп ът(Дп1 - упх)1 (12)
п=0
Дп=, у=«шп, (,э)
г аг
где А, Вп — произвольные постоянные, для которых ряд (12)-(13) и его производные (У2 ) и (У2) сходятся (сходимость можно обеспечить, если положить Аи = В =0 при п > Ы, где N — произвольное натуральное число).
Уравнение 3. Нелинейное уравнение переноса с запаздыванием
и + аих = Ъи 1п и + и/ ^ / и), w = и(х, £-г), (14)
допускает решение с мультипликативным разделением переменных
и = ехр(С/х/а +Я/ Ъ)И>),
где С — произвольная постоянная, а функция у = у(£) удовлетворяет ОДУ с постоянным запаздыванием
у' = у[Я + Ъ 1пу + / {у/у)], у = у(£-г).
Здесь Я — произвольная постоянная.
Уравнение 4. Рассмотрим нелинейное уравнение переноса с запаздыванием
щ + аих = -Ъи 1пи + и/(w/ик), w = и(х,£-г), к > 0. (15)
1. Уравнение (15) при Ъ = (1п к)/г допускает решение с мультипликативным разделением переменных
и = ехр( Ае~ы )р( х),
где А — произвольная постоянная, а функция р = р(х) описывается нелинейным ОДУ
ар' + Ър1пр = р/ (рр- ).
2. Уравнение (15) при Ъ = (1пк)/г допускает другое решение с мультипликативным разделением переменных
и = ехр( Ахе~ы )у(£),
где А — произвольная постоянная, а функция у = у(£) описывается нелинейным ОДУ с запаздыванием
у' + у[Ъ 1пу-/ (у/ук ) + Аае~ы ] = 0, у = у(£-г).
Уравнение 5. Рассмотрим нелинейное уравнение переноса с запаздыванием
щ + аих = / (и - w), w = и(х, £-г), (16)
где / (... ) — произвольная функция.
1. Уравнение (16) имеет точное решение с аддитивным разделением переменных, линейное по х :
и = Ах + у(г), (17)
где А — произвольная постоянная, а функция у = у(£) описывается ОДУ первого порядка с постоянным запаздыванием
у' = /(у-у) - Аа, у = у(£-г). (18)
ОДУ с запаздыванием (18) имеет линейное по £ частное решение у(£) = Я£ + В, где В — произвольная постоянная, а Я — корень алгебраического (трансцендентного) уравнения Аа + Я- / (гЯ) = 0.
2. Уравнение (16) имеет также решение, более общее, чем (17):
и = Сх + С^ + в(г), г = Рх + уt, (19)
где С, С, Р, У — произвольные постоянные, а функция в(г) удовлетворяет ОДУ первого порядка с постоянным запаздыванием
(ар + у)в = /(в-в+ С2г) -Са - С, в = в(г -о), о = уг.
Значениям С — С — 0 в (19) соответствует решение типа бегущей волны.
Уравнение 6. Рассмотрим нелинейное уравнение переноса с запаздыванием
и = С/х/а + ), (20)
которое при Ъ = 0 переходит в уравнение (16).
1. Уравнение (20) имеет точное решение с аддитивным разделением переменных
и = С/х/а + ), (21)
где С — произвольная постоянная, а функция у = у/^) удовлетворяет ОДУ с запаздыванием
у' = Ъу + /(у-у), у = у(-г).
2. Уравнение (20) допускает также решение, более общее, чем (21):
и = СеЫ а + в( г), г = Рх + уt,
где С — произвольная постоянная, а функция р = (р(х) удовлетворяет ОДУ с запаздыванием
ар = Ър + / (СгеЬх1 а ) - СеЬх1" •
Теорема 1 (о нелинейной суперпозиции решений). Пусть щ (х, t) — некоторое решение нелинейного уравнения с запаздыванием (20), а функция у = V(х,V,Ъ) является любым г- периодическим решением линейного уравнения переноса (8). Тогда сумма
и = и0 (х, ^ + V (х, t; Ъ)
также является решением уравнения (20). Общий вид функции V (х, t; Ъ) задается формулами (9)—(10).
Уравнение 7. Нелинейное уравнение переноса с запаздыванием
и + аих = Ъи + /(и - kw), w = и(х, t - г), к >0, (22)
допускает точное решение с обобщенным разделением переменных
и = С ехр[с£ + (Ь - с)х / а] + у(х), с =11п к,
т
где С — произвольная постоянная, а функция у = у(х) удовлетворяет нелинейному ОДУ
ау' = Ьу + f ()), ) = (1 - k )у.
Теорема 2 (обобщает теорему 1). Пусть и0(x, Г) — некоторое решение нелинейного уравнения с запаздыванием (22), а функция
V = V (х, Р, Ь) является любым т- периодическим решением линейного уравнения переноса (8). Тогда сумма
и = щ (х, {) + ес¥1 (х, V; Ь - с), с = 11п к,
т
также является решением уравнения (22). Общий вид функции
V (х, Р, Ь) задается формулами (9)—(10).
Уравнение 8. Рассмотрим нелинейное уравнение переноса с запаздыванием
и + аих = Ьи + f (и + км>), w = и(х, t -т), к >0, (23)
которое отличается от уравнения (22) знаком перед параметром к в аргументе кинетической функции. Справедлива следующая теорема.
Теорема 3. Пусть и0 (х, {) — некоторое решение нелинейного уравнения с запаздыванием (23), а функция V = V(х,Р,Ь) — любое т- апериодическое решение линейного уравнения переноса (11). Тогда сумма
и = щ (х, Р) + ес^2 (х, Р, Ь - с), с = 11п к,
т
также является решением уравнения (23). Общий вид функции V(х, t;Ь) задается формулами (12)—(13).
Уравнение 9. Рассмотрим более сложное нелинейное уравнение переноса с запаздыванием
и + аих = ^ (и - w) + wg (и - w) + к(и - w), w = и( х, t - т), (24)
где / (...), g (...), к(...) — произвольные функции (одну из двух функций f или g можно без ограничения общности положить равной нулю).
1. Уравнение (24) допускает решение с обобщенным разделением переменных
и = ф(х) + у (х), (25)
где функции р = р(х) и у = у(х) удовлетворяют системе ОДУ:
аф = р[/(() + ё (гр^
ау' = У[/ (^Р) + ё (()]" Р[1 + т§ (()] + Н(гр)-
2. Уравнение (24) допускает более сложное, чем (25), решение с обобщенным разделением переменных
N 2 тгп
и = £[р» (х) С08(Д/) + У (х) )] + в( х) + £( х), = —,
И=1 ^
где N — любое натуральное число, а функции рп = рп (х), уп = уп (х), в = в(х) и £ = £( х) определяются из системы ОДУ:
С = Рп [/(в) + ё (в)]-Ру, аУ = Уп[/(в) + ё (гв)] + ^пРп, ав' = в[/(гв) + ё (гв)], а£ ' = £[/(гв) + ё (гв)] - в[1 + гё (гв)] + И(гв).
Отметим, что третье нелинейное ОДУ допускает тривиальное решение в = 0: в этом случае остальные уравнения становятся линейными ОДУ с постоянными коэффициентами.
Уравнение 10. Рассмотрим нелинейное уравнение переноса с запаздыванием, которое обобщает уравнение (22):
и + аих = и/(и - км) + мё(и - км) + И(и - км),
(26)
м = и(х,^-г), к > 0.
1. Уравнение (26) допускает точное решение с обобщенным разделением переменных
и = ер( х) + у (х), с = — 1п к, (27)
г
где функции р = р(х) и у = у( х) удовлетворяют системе ОДУ:
аР' = Р[/(Л) +1 ё (Л) - c], Л = (1 - к )У, к
аУ' = у[/(л)+ё (л)]+КЛУ
2. Уравнение (26) допускает более сложное, чем (27), решение с обобщенным разделением переменных
и = ес |в(х) + £[рп (х) соь(Рп1) + ¥п (х) вт(Р;)]| + £(х),
1 2лп
c = ln k, ßn =-,
г г
где N — любое натуральное число, а функции (рп = (рп (х), = (х), в = в(х) и % = %(x) определяются из системы ОДУ:
«рП = % U(л)+1 g(л) - c] -ßnVn, л = (1 - k)%, k
<w'„ =Wn[f (Л) +1 g(Л) - c] + ßnPn, k
ав' = e[f (л) +1 g (л) - c], k
а% ' = %[f (Л) + g (Л)] + Кл). Уравнение 11. Нелинейное уравнение переноса с запаздыванием
ut + аих = uf (и + kw) + wg (и + kw) + h(u + kw),
(28)
w = u(x,t-г), k > 0,
которое обобщает уравнение (23), допускает решение с обобщенным разделением переменных
N
и = ect (х) cos(ßnt) + ¥п (х) sin(ßnt)] + % (х),
n х
n=1
. _ ж(2п-1)
с = 1п к, рп =-,
Т т
где N — любое натуральное число, а функции ри = ри (х), = уп (х) и % = х) определяются из системы ОДУ:
арп = (Рп [/(Л) -1Я (Л) - с] -Р„¥„, Л = (1 + к^ к
а¥'п = ¥п[/ л)-1Я (Л) - с] + л(п, к
' = %[/(л) + Я (л)] + л(л).
Уравнение 12. Нелинейное уравнение переноса с запаздыванием и + аих = и/ (и2 + м?) + мя (и2 + м?), м = и( х, £ -т), (29) допускает точное решение с обобщенным разделением переменных:
и = Рп (х)со*Ш) + Уп(фшСдд » = 0+1+2.....
2т
где функции рп = ри (х), уп = уп (х) определяются из системы ОДУ:
арп = Рп/(Л) - (-1)>й£ (Л)-Рп¥п, V = + ¥1n, а¥'п =¥п/(Л) + (-1)" Рп% (П) + РпРп •
Уравнение 13. Рассмотрим нелинейное уравнение переноса с запаздыванием
и + аих = и/ (и - км, м / и), м = и(х, t-г), к > 0, (30) где / (у, у ) — произвольная функция двух аргументов.
и = и (г), z = ¡х + Л, Уравнение (30) допускает точное решение типа бегущей волны
(аР + Я)и' = и/ (и - кЖ,Ж / и), Ж = и (г -ст), а = Лг. (31)
Уравнение (31) при ¡3 = —
а
' 1Л
—/10,
У V
-Л
имеет частное решение
и (г) = С^2, У = -1п к.
ст
Точные аналитические решения уравнений переноса с запаздыванием и переменным по х коэффициентом переноса.
Уравнение 14. Нелинейное уравнение переноса с запаздыванием
и + а(х)их = Ь(х)/(и - м), м = и(х, t - г), (32)
коэффициенты которого зависят от х, допускает точное решение с аддитивным разделением переменных в замкнутом виде:
\
йх.
и=t+![ Ой[ / (г)Ь( х) -11
Уравнение 15. Нелинейное уравнение переноса с запаздыванием
и + а(х)их = Ь(х)и/(м / и), м = и(х, t - г), (33)
имеет точное решение с мультипликативным разделением переменных
и = вМф( х),
где Л — произвольная постоянная, а функция р = р(х) описывается линейным ОДУ
а( х)р + [Л - Ь( х) / (е~Лг )]р = 0.
Уравнение 16. Нелинейное уравнение переноса с запаздыванием
и + а(х)их = Ь(х)и + и/ (м / и), м = и(х, г-т), (34)
имеет точное решение с мультипликативным разделением переменных
и = р(х)у(г),
где функции р(х) и у = у (г) описываются соответственно линейным ОДУ и нелинейным ОДУ с запаздыванием:
а( х)р' (х) + [1- Ь( х)]р( х) = 0, ^ ^
у' = у[1 + /(у /у)], у = у(г -т),
1 — произвольная постоянная.
Нелинейное ОДУ с запаздыванием (35) допускает частные решения вида у = Лв(, где Л — произвольная постоянная, ( определяется из алгебраического (трансцендентного) уравнения
( = 1 + /(в-(Т).
Уравнение 17. Нелинейное уравнение переноса с запаздыванием, которое имеет достаточно общий вид:
и + а(х)их = Ь(х)/(и, м), м = и(х, г - т), (36)
допускает точное решение с функциональным разделением перемен-
Ь( х) -1
ных
, Л
и = и (г), г = г +П
йх,
а( х)
где функция и = и (г) определяется из нелинейного ОДУ первого порядка с запаздыванием
и' = / (и ,Ж), Ж = и( г -т).
Уравнение 18. Нелинейное уравнение переноса с запаздыванием и + а(х)их = и/(х, и - м) + я(х, и - м), м = и(х, г -т), (37) допускает точное решение с обобщенным разделением переменных
и = р(х)г + у(х),
где функции р = р(х) и у = у(х) описываются ОДУ первого порядка
а( х)р' (х) = р( х) / ( х, тр( х) ), а( х)у ' (х) = у( х)/ ( х, тр( х) ) + я (х, тр( х) ) - р(х).
Точные аналитические решения уравнений переноса с запаздыванием и зависящим от искомой функции и коэффициентом переноса.
Уравнение 19. Нелинейное уравнение переноса с запаздыванием и + (ои + а0 )их = Ьи + Ьм, м = и(х, t - г), (38)
допускает точное решение с обобщенным разделением переменных
и = щ (t) + щ (t) х,
где функции щ = щ (t) и щ = щ (t) описываются системой ОДУ с запаздыванием
щ = -ащщ2 -а0щ2 + Ьщ— + Ь2щ—, = -г),
щ2 = -а—щ2 + Ь\¥2 + Ь2щ2 , ¥2 = ¥2 ( - г). Уравнение 20. Нелинейное уравнение переноса с запаздыванием и + (аи + а0 )и = ки2 + Ьи + Ь2м, м = и(х, £-г), (39)
допускает точное решение с обобщенным разделением переменных
и = щ (^ + щ ^) ехр(кх / о ),
где функции щ = щ (t) и щ = щ (t) описываются системой ОДУ с запаздыванием
щ = + кщ—2 + Ь2щ1, щ1 = - г),
^2
Ь1
V а1 У
щ + кщщ + Ь2щ, щ = щ ^ - г).
Уравнение 21. Нелинейное уравнение переноса со степенными нелинейностями и запаздыванием
и + = ЬмИ+1 + си + км1-" + ши1~пм?п, м = м(х, г - г), (40)
которое зависит от шести параметров, допускает точное решение с функциональным разделением переменных
и = [р(£ ) ехр(Ьпх / а) + ",
где функции р = р(£ ) и щ = ) описываются системой ОДУ с запаздыванием
у' = кп + спщ + Ьпщ2 + тпщ, щ = - г), р' = спр + Ьпрщ + тпр, р = р(£ -г).
Уравнение 22. Нелинейное уравнение переноса с экспоненциальными нелинейностями и запаздыванием
щ + ав1иих = Ьв1и + с + кв~1и + тв1(м,-и), м = и(х, г-т), (41)
которое зависит от шести параметров, допускает точное решение с функциональным разделением переменных
и = —1п[Лехр((г + Ь1х/а)-у], у = —(с + т(с + т)2 -4Ьк), 1 2Ь
где Л — произвольная постоянная, ( — корень трансцендентного уравнения ( — т1в~(т = (с - Ьу)1.
Уравнение 23. Нелинейное уравнение переноса с запаздыванием
щ + ауГ112их = Ьи112 + /(и112 - м1/2), м = щ(х,г - т), (42)
которое зависит от произвольной функции / (...), допускает точное решение с функциональным разделением переменных:
и = [р( х)г + у( х)]2,
где р = р( х) и у = у( х) описываются системой ОДУ
2ар' = Ьр- 2р2,
2ау ' = Ьу- 2ру + / (тр). Уравнение (43) является уравнением Бернулли и имеет решение
Ь
--х Л
р = Ь(Лв 2а + 2) ,
где Л — произвольная постоянная.
Уравнение 24. Нелинейное уравнение переноса с запаздыванием
щ + аиких = и/ (м / и), м = и(х, г -т), (44)
имеет точное решение с мультипликативным разделением переменных
и = к (Ьх / а + Л)1/ку (г ),
где Л — произвольная постоянная, а функция у = у (г) удовлетворяет нелинейному ОДУ с запаздыванием
щ + аиких = Ьик+1 + и/(м / и), м = щ(х, г - т).
Уравнение 25. Рассмотрим нелинейное уравнение переноса с запаздыванием
щ + аиких = Ьик+1 + и/(м / и), м = и(х, г -т). (45)
1. Уравнение (45) допускает точное решение с мультипликативным разделением переменных
и = р( х)у(г),
где функция у = у (г) описывается нелинейным ОДУ с запаздыванием
у' = 1ук+1 + у/(у/у), у = у(г-т),
а функция р = р(х) удовлетворяет ОДУ Бернулли
ар ' = Ър-1рк+1. (46)
Здесь 1 — произвольная постоянная.
ОДУ Бернулли (46) имеет общее решение
р = Ь к [1 + Л ехр(Ькх / а)]к,
где Л — произвольная постоянная.
2. Уравнение (45) допускает другое точное решение с мультипликативным разделением переменных
и = ехр(Ьх / а)у (г), где функция у = у (г) описывается нелинейным ОДУ с запаздыванием
у ' = у/(у УХ у = у(г-т). (47)
3. Уравнение (45) допускает еще одно точное решение с мультипликативным разделением переменных
1
и = ехр(Ь(к +1) х / а) к+1 у (г),
где у = у(г) описывается нелинейным ОДУ с запаздыванием (47). Уравнение 26. Нелинейное уравнение переноса с запаздыванием
щ + аиких = ¿ик-п+: + и1-п/(ип - мп), м = и(х, г - т), (48)
допускает точное решение с функциональным разделением переменных
и =
1
» Ьп . . Л +—х + у(г) а
где Л — произвольная постоянная, а функция у = у (г) удовлетворяет нелинейному ОДУ с запаздыванием
у' = п/(у-уХ у = у(г -т).
Уравнение 27. Рассмотрим нелинейное уравнение переноса с запаздыванием
и + аеЛиих = /(и - м), м = и(х, t - г). (49)
1. Уравнение (49) допускает точное решение с аддитивным разделением переменных
и =11п( Ах + Б) + щ(1),
Л
где А и Б — произвольные постоянные, а функция щ = ) удовлетворяет нелинейному ОДУ с запаздыванием
_ аА _
щ = /(щ-щ —геЛщ, щ = щ(1 -г).
Л
2. Уравнение (49) допускает также точное решение вида
и = —1п( Ах + Б) + 0(г), г = С +1 + — 1п( Ах + Б), Л А
где А , Б, С — произвольные постоянные, а функция 0 = 0(г) удовлетворяет нелинейному ОДУ с запаздыванием
(1 + авЛ0)0' = /(0-0) - °АеЛ0, 0 = 0(г-г).
Уравнение 28. Нелинейное уравнение переноса с запаздыванием и + аемих = ЬеЛи + /(и - м) м = и(х, Г - г), (50)
допускает точное решение с аддитивным разделением переменных:
Ь
и = — х + ), а
где функция щ = ) описывается нелинейным ОДУ с запаздыванием
щ = /Д-w), щ = Щ(г-г)• Уравнение 29. Нелинейное уравнение переноса с запаздыванием
и + аеЛиих = Ье(Л-У)и + е~уи /(еуи -еум) м = и(х,t-г), (51)
допускает точное решение с функциональным разделением переменных
и = — 1п
У
Ьу
— х + щ(£ ) а
где функция у = у (г) описывается нелинейным ОДУ с запаздыванием
¥' = 7Ду-у), ¥ = ¥(г~т). Уравнение 30. Нелинейное уравнение переноса с запаздыванием щ + (а 1п и + Ь)их = си 1п и + иД (м / и), м = и(х, г-т), (52)
имеет точное решение с мультипликативным разделением переменных
и = еы ау(г),
где у = у (г) удовлетворяет нелинейному ОДУ с запаздыванием у' = у[Д (у /у) - Ьс / а], у = у(г -т). Уравнение 31. Нелинейное уравнение переноса с запаздыванием
иг + Д'(и)их = а1Д(и) + а2Д(м) + а3 + Ь ^Д(и^- Д,
Д (и) (53)
м = и (х, г -т),
где Д (... ) — произвольная функция, допускает точное решение в неявном виде
Д (и) = Лерх+1г--, р = а + а2е1т,
а + а
где Л — произвольная постоянная, а константа 1 определяется из трансцендентного уравнения 1 - Ь + Ье~1т = 0.
Уравнение 32. Нелинейное уравнение переноса с запаздыванием
\ЬД (и ) + ЬД (м) + Ь1
и, + Д (и)их = а[Д (и) - Д (м)] +\^() ^ ( ) 31, ....
Д (и) (54)
м = и (х, г - т),
допускает точное решение в неявном виде
Д(и) = Лерх+1г --Ь-, р = а(1 - е-1т),
Ь + Ъ2
где Л — произвольная постоянная, а константа 1 определяется из трансцендентного уравнения 1- Ь - Ь2е~1т = 0.
Уравнение 33. Нелинейное уравнение переноса с запаздыванием
щ + аи-1,2их = Д(и112 -м1/2) + uV2g{uV2 -м1/2), м = и(х,г-т), (55)
где / (...) и g (...) — произвольные функции, допускает точное решение с функциональным разделением переменных
и = [р( х)г + щ( х)]2,
где функции р = р(х) и щ = щ( х) определяются из системы ОДУ
2ар' = рg[гр] - 2р2, 2ащ' = / [гр] + уg[гр] - 2рщ
Уравнение 34. Нелинейное уравнение переноса с запаздыванием
щ + аиких = и/ (м / и) + uk+1g (м / и), м = и(х, г-г), (56)
допускает точное решение вида
и = АеРх, Р = g(e~-г) / а,
где константа Л определяется из алгебраического (трансцендентного) уравнения Л = / (е~Лг).
Уравнение 35. Нелинейное уравнение переноса с запаздыванием
и + аиких = /(ик+1 - м^1) + u-kg(uk+1 - мк+г), м = и(х, Г-г), (57)
допускает точное решение вида
и = (А + Рх + Л)кТТ, р = — / (Лг),
а
где А — произвольная постоянная, а константа Л определяется из алгебраического (трансцендентного) уравнения Л = (1 + к) g (Лг).
Уравнение 36. Нелинейное уравнение переноса с запаздыванием
и + аеЛ"их = /(и -м) + eЛug(и -м), м = и(х,г -г), (58) допускает точное решение типа бегущей волны
и = 0( г), г = Рх + Лг, где функция 0( г) определяется из ОДУ с запаздыванием
(Л + ареЛ0)0' = / (0-0) + eЛ0g(0-0), 0 = 0(г -с), ст = Л г. Уравнение 37. Нелинейное уравнение переноса с запаздыванием и + аеЛиих = /(еЛи - еЛ) + e-Лug(еЛи - еЛ), м = и(х, г - г), (59) допускает точное решение вида
1 Л
и = - 1п( А + ух + рг), у = -/ (Рг), Л а
где Л — произвольная постоянная, а Р определяется из алгебраического (трансцендентного) уравнения р = 1g (Рт).
Уравнение 38. Нелинейное уравнение переноса с запаздыванием
и, + [а^и)+Ь]и, = с.g (и) + с g (м) + Д [g (и) - g (м)],
g (и) (60)
м = и (х, г -т),
где Д (...) и g (...) — произвольные функции, допускает точное решение в неявном виде
g (и) = ^^ х + у(г), а
где функция у = у (г) удовлетворяет ОДУ с запаздыванием
у' = с2 (у - у) - Ь(С1 + С2) + /(у - у), а
у = у(г -т).
Уравнение 39. Нелинейное уравнение переноса с запаздыванием щ + а(х)иких = Ь(х)ик+1 + иД (м / и), м = и(х, г-т), (61)
имеет точное решение с мультипликативным разделением переменных
и = ехр ( СЬ(х) йх I^ а(х)
у (г),
где функция у = у (г) удовлетворяет ОДУ с запаздыванием
у' = уДД УХ у = у(г-т). Уравнение 40. Нелинейное уравнение переноса с запаздыванием щ + а(х)е1иих = Ь(х)е1и + Д(и - м), м = и(х, г - т), (62)
допускает точное решение с аддитивным разделением переменных
и = ГЬ( х) йх + у (г), а (х)
где функция у = у (г) удовлетворяет ОДУ с запаздыванием
у' = ДУ-у), у = у(г~т).
Уравнение 41. Нелинейное уравнение переноса с запаздыванием
щ + а/'(щ)их = 8ти) -/И] + м/ (и}~/(")], _
/(и) (63)
ы = и( х, t -г),
где / (...), g (...) и й(...) — произвольные функции, допускает точное решение в неявном виде
/ (и) = Л + ^^ х + а
где Л — произвольная постоянная, а константа 1 удовлетворяет алгебраическому (трансцендентному) уравнению 1 = Н(1г).
Уравнение 42. Нелинейное уравнение переноса с запаздыванием
Щ + а/' (и)их = / (и) g
/ (и)
, /(и) !г ( /(^)Л
/'(и) I/(и)
(64)
ы = и (х, t -г),
допускает точное решение в неявном виде / (и) = Л ехр
х + 1
\ а ,
где Л — произвольная постоянная, а 1 определяется из алгебраического (трансцендентного) уравнения 1 = к(вЯт).
Уравнение 43. Нелинейное уравнение переноса с запаздыванием
щ + g (и )их =—^ ( bg(щ) + к[/(и ) - / (^)]), ы = и t-T), (65) / (и)
допускает точное решение в неявном виде
/ (и) = Ьх + у/^), где функция у = у/^) удовлетворяет ОДУ с запаздыванием
у' = к (у-у), у = у(t-г). Уравнение 44. Нелинейное уравнение переноса с запаздыванием
щ + а(х)/'(и)их = Ь(хМ/(щ) - /(*)] + М/ (щ)- "С(Ы)],
/(Щ) (66)
ы = и(х, t - г),
допускает точное решение в неявном виде
/ (Щ) = Л + g (1г) |Ь(х) йх + 1, -1 а (х)
где Л — произвольная постоянная, а константа X определяется из алгебраического (трансцендентного) уравнения X = Н(Лт).
Уравнение 45. Нелинейное уравнение переноса с запаздыванием
где A — произвольная постоянная, а константа X — корень алгебраического (трансцендентного) уравнения X = h(e~Xr).
Выводы. Рассмотрено сорок пять нелинейных уравнений переноса с запаздыванием, в том числе восемнадцать уравнений с постоянным коэффициентом переноса и двадцать семь уравнений с коэффициентом переноса, зависящим от искомой функции. Все уравнения содержат свободные параметры и/или произвольные функции. Получены новые точные решения с аддитивным, мультипликативным и обобщенным разделением переменных, решения типа бегущей волны и автомодельные решения. Для трех уравнений с постоянным коэффициентом переноса сформулированы теоремы о размножении решений. Все решения выражаются в элементарных функциях, либо через неопределенные интегралы, либо через решения обыкновенных уравнений с запаздыванием или без него, либо через решения уравнений в частных производных без запаздывания. Результаты работы могут быть полезны для анализа рассматриваемых и родственных нелинейных уравнений в частных производных с запаздыванием, а также для тестирования численных и приближенных аналитических методов интегрирования таких уравнений.
Автор благодарит А.Д. Полянина за внимание к статье и полезные замечания.
Работа выполнена по теме государственного задания (№ госрегистрации 124012500440-9).
ЛИТЕРАТУРА
[1] Singh P., Sharma K.K. Numerical solution of first-order hyperbolic partial differential-difference equation with shift. Numerical Methods for Partial Differential Equations, 2010, vol. 26, pp. 107-116.
[2] Singh P., Sharma K.K. Finite difference approximations for the first-order hyperbolic partial differential equation with point-wise delay. arXiv.org, 2022. DOI: https://doi.org/10.48550/arXiv.1012.0974
ut + a(x)f'(u)ux = b(x)f (u)g TTT
If (u)
w = u( x, t -r),
допускает точное решение в неявном виде
' fw V fwh (fM"
f (u) J + f' (u ) h [ f (u) J, (67)
h
[3] Karthick S., Mahendran R., Subburayan V. Method of lines and Runge-Kutta method for solving delayed one dimensional transport equation. Journal of Mathematics and Computer Science, 2023, vol. 28, pp. 270-280.
[4] Sharma K.K., Singh P. Hyperbolic partial differential-difference equation in the mathematical modeling of neuronal firing and its numerical solution. Applied Mathematics and Computation, 2008, vol. 201, pp. 229-238.
[5] Zaidi A.A., Van Brunt B., Wake G.C. Solutions to an advanced functional partial differential equation of the pantograph type. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2015, vol. 471, art. no. 20140947. DOI: http://doi.org/10.1098/rspa.2014.0947
[6] Rey A.D., Mackey M.C. Bifurcations and traveling waves in a delayed partial differential equation. Chaos, 1992, vol. 2, p. 231-244.
[7] Mackey M.C., Rudnicki R. Global stability in a delayed partial differential equation describing cellular replication. Journal of Mathematical Biology, 1994, vol. 33, pp. 89-109.
[8] Dyson J., Villella-Bressan R., Webb G.F. A semilinear transport equation with delays. International Journal of Mathematics and Mathematical Sciences, 2003, vol. 32, pp. 2011-2026.
[9] Mackey M.C., Rudnicki R. A new criterion for the global stability of simultaneous cell replication and maturation processes. Journal of Mathematical Biology, 1999, vol. 38, pp. 195-219.
[10] Rhee H.-K., Aris R., Amundson N.R. First-order partial differential equations, Vol. 1. Theory and application of single equations. Mineola, New York, Dover Publications, 1986, 543 p.
[11] Rhee H.-K., Aris R., Amundson N. R. First-order partial differential equations, Vol. 2. Theory and application of hyperbolic systems of quasilinear equations. Mineola, New York, Dover Publications, 1989, 544 p.
[12] Полянин А.Д., Сорокин В.Г., Журов А.И. Дифференциальные уравнения с запаздыванием: Свойства, методы, решения и модели. Москва, ИПМех РАН, 2022, 463 с.
[13] Kamont Z., Przadka K. Difference methods for first order partial differen-tial-functional equations with initial-boundary conditions. USSR Computational Mathematics and Mathematical Physics, 1991, vol. 31, pp. 1476-1488.
[14] Sampath K., Veerasamy S., Agarwal R.P. Stable difference schemes with interpolation for delayed one-dimensional transport equation. Symmetry, 2022, 14 (5), art. no. 1046. DOI: https://doi.org/10.3390/sym14051046
[15] Solodushkin S.I., Yumanova I.F., De Staelen R.H. First order partial differential equations with time delay and retardation of a state variable. Journal of Computational and Applied Mathematics, 2015, vol. 289, pp. 322-330.
[16] Пименов В.Г. Разностные методы решения уравнений в частных производных с наследственностью. Екатеринбург, Изд-во Уральского университета, 2014, 134 с.
[17] Meleshko S.V., Moyo S. On the complete group classification of the reaction-diffusion equation with a delay. Journal of Mathematical Analysis and Applications, 2008, vol. 338, pp. 448-466. DOI: 10.1016/j.jmaa.2007.04.016
[18] Long F.-S., Meleshko S.V. On the complete group classification of the one-dimensional nonlinear Klein—Gordon equation with a delay. Mathematical Methods in the Applied Sciences, 2016, vol. 39, no. 12, pp. 3255-3270. DOI: 10.1002/mma.3769
[19] Lobo J.Z., Valaulikar Y.S. Group analysis of the one dimensional wave equation with delay. Applied Mathematics and Computation, 2020, vol. 378, art. no. 125193. DOI: 10.1016/j.amc.2020.125193
[20] Polyanin A.D., Zhurov A.I. Functional constraints method for constructing exact solutions to delay reaction-diffusion equations and more complex nonlinear equations. Communications in Nonlinear Science and Numerical Simulation, 2014, vol. 19, pp. 417-430. DOI: https://doi.org/10.1016/j.cnsns.2013.07.017
[21] Polyanin A.D., Zhurov A.I. The functional constraints method: Application to non-linear delay reaction-diffusion equations with varying transfer coefficients. International Journal of Non-Linear Mechanics, 2014, vol. 6, pp. 267-277. DOI: 10.1016/j.ijnonlinmec.2014.09.008
[22] Полянин А.Д., Сорокин В.Г., Вязьмин А.В. Нелинейные реакционно-диффузионные уравнения гиперболического типа с запаздыванием: точные решения, глобальная неустойчивость. Математическое моделирование и численные методы, 2014, № 4, с. 53-73.
[23] Полянин А.Д., Сорокин В.Г. Точные решения нелинейных уравнений в частных производных с переменным запаздыванием типа пантографа. Вестник НИЯУ «МИФИ», 2020, т. 9, № 4, с. 315-328.
[24] Polyanin A.D., Zhurov A.I. Exact separable solutions of delay reaction-diffusion equations and other nonlinear partial functional-differential equations. Communications in Nonlinear Science and Numerical Simulation, 2014, vol. 19, no. 3, pp. 409-416. DOI: 10.1016/j.cnsns.2013.07.019
[25] Polyanin A.D., Zhurov A.I. Generalized and functional separable solutions to nonlinear delay Klein—Gordon equations. Communications in Nonlinear Science and Numerical Simulation, 2014, vol. 19, no. 8, pp. 2676-2689. DOI: 10.1016/j.cnsns.2013.12.021
[26] Polyanin A.D., Zhurov A.I. Nonlinear delay reaction-diffusion equations with varying transfer coefficients: Exact methods and new solutions. Applied Mathematics Letters, 2014, vol. 37, pp. 43-48. DOI: https://doi.org/10.1016/j.aml.2014.05.010
[27] Polyanin A.D., Zhurov A.I. New generalized and functional separable solutions to nonlinear delay reaction-diffusion equations. International Journal of Non-Linear Mechanics, 2014, vol. 59, pp. 16-22. DOI: https://doi.org/10.1016/j.ijnonlinmec.2013.10.008
[28] Полянин А.Д., Сорокин В.Г. Построение точных решений нелинейных уравнений математической физики с запаздыванием с помощью решений более простых уравнений без запаздывания. Вестник НИЯУ «МИФИ», 2020, т. 9, № 2, с. 115-128.
[29] Polyanin A.D., Sorokin V.G. Nonlinear pantograph-type diffusion PDEs: Exact solutions and the principle of analogy. Mathematics, 2021, vol. 9, art. no. 511. DOI: https://doi.org/10.3390/math9050511
[30] Polyanin A.D., Sorokin V.G. A method for constructing exact solutions of nonlinear delay PDEs. Journal of Mathematical Analysis and Applications, 2021, vol. 494, art. no. 124619. DOI: 10.1016/j.jmaa.2020.124619
[31] Polyanin A.D., Sorokin V.G. Construction of exact solutions to nonlinear PDEs with delay using solutions of simpler PDEs without delay. Communications in Nonlinear Science and Numerical Simulation, 2021, vol. 95, art. no. 105634. DOI: 10.1016/j.cnsns.2020.105634
[32] Polyanin A.D., Sorokin V.G. Reductions and exact solutions of nonlinear wave-type PDEs with proportional and more complex delays. Mathematics, 2023, vol. 11, art. no. 516. DOI: https://doi.org/10.3390/math11030516
[33] Polyanin A.D., Sorokin V.G. Exact solutions of reaction-diffusion PDEs with anisotropic time delay. Mathematics, 2023, vol. 11, art. no. 3111. DOI: https://doi.org/10.3390/math11143111
[34] Tanthanuch J. Symmetry analysis of the nonhomogeneous inviscid Burgers equation with delay. Communications in Nonlinear Science and Numerical Simulation, 2012, vol. 17, pp. 4978-4987. DOI: 10.1016/j.cnsns.2012.05.031
[35] Polyanin A.D., Zaitsev V.F. Handbook of Nonlinear Partial Differential Equations, 2nd ed. Boca Raton, CRC Press, 2012, 1912 p.
Статья поступила в редакцию 15.04.2024
Ссылку на эту статью просим оформлять следующим образом: Сорокин В.Г. Аналитические решения нелинейных уравнений с запаздыванием, используемых при математическом моделировании процессов переноса.
Математическое моделирование и численные методы, 2024, № 3, с. 140-167.
Сорокин Всеволод Григорьевич — канд. физ.-мат. наук, научный сотрудник Института проблем механики им. А.Ю. Ишлинского РАН. e-mail: [email protected]
Analytical solutions of nonlinear transport equations with delay
© V.G. Sorokin
Ishlinsky Institute for Problems in Mechanics RAS, Moscow, 119526, Russia
Nonlinear transport equations with constant delay are considered. The introduction provides a brief overview ofpublications that study transport mathematical models with delay and develop numerical methods for solving the corresponding problems. The main sections of the article describe more than forty transport equations with constant delay and various transfer coefficients, which allow exact analytical solutions. The kinetic functions of all considered equations contain free parameters or arbitrary functions. Additive, multiplicative, generalized, and functional separable solutions, as well as traveling-wave and self-similar solutions are obtained. Many solutions are expressed in terms of elementary functions. For some types of equations, theorems on the "multiplication " of solutions are formulated. The described equations and their solutions can be used to evaluate the accuracy of numerical methods for integrating the corresponding nonlinear transport problems with delay.
Keywords: nonlinear transport equations with delay, delay partial differential equations, analytical solutions, generalized and functional separable solutions, self-similar solutions
REFERENCES
[1] Singh P., Sharma K.K. Numerical solution of first-order hyperbolic partial differential-difference equation with shift. Numerical Methods for Partial Differential Equations, 2010, vol. 26, pp. 107-116.
[2] Singh P., Sharma K.K. Finite difference approximations for the first-order hyperbolic partial differential equation with point-wise delay. arXiv.org, 2022. DOI: https://doi.org/10.48550/arXiv.1012.0974
[3] Karthick S., Mahendran R., Subburayan V. Method of lines and Runge-Kutta method for solving delayed one dimensional transport equation. Journal of Mathematics and Computer Science, 2023, vol. 28, pp. 270-280.
[4] Sharma K.K., Singh P. Hyperbolic partial differential-difference equation in the mathematical modeling of neuronal firing and its numerical solution. Applied Mathematics and Computation, 2008, vol. 201, pp. 229-238.
[5] Zaidi A.A., Van Brunt B., Wake G.C. Solutions to an advanced functional partial differential equation of the pantograph type. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2015, vol. 471, art. no. 20140947. DOI: http://doi.org/10.1098/rspa.2014.0947
[6] Rey A.D., Mackey M.C. Bifurcations and traveling waves in a delayed partial differential equation. Chaos, 1992, vol. 2, p. 231-244.
[7] Mackey M.C., Rudnicki R. Global stability in a delayed partial differential equation describing cellular replication. Journal of Mathematical Biology, 1994, vol. 33, pp. 89-109.
[8] Dyson J., Villella-Bressan R., Webb G.F. A semilinear transport equation with delays. International Journal of Mathematics and Mathematical Sciences, 2003, vol. 32, pp. 2011-2026.
[9] Mackey M.C., Rudnicki R. A new criterion for the global stability of simultaneous cell replication and maturation processes. Journal of Mathematical Biology, 1999, vol. 38, pp. 195-219.
[10] Rhee H.-K., Aris R., Amundson N.R. First-order partial differential equations, Vol. 1. Theory and application of single equations. Mineola, New York, Dover Publications, 1986, 543 p.
[11] Rhee H.-K., Aris R., Amundson N. R. First-order partial differential equations, Vol. 2. Theory and application of hyperbolic systems of quasilinear equations. Mineola, New York, Dover Publications, 1989, 544 p.
[12] Polyanin A.D., Sorokin V.G., Zhurov A.I. Differencial'nye uravnenija s za-pazdyvaniem: Svojstva, metody, reshenija i modeli [Delay differential equations: properties, methods, solutions and models]. Moscow, IPMech RAS Publ., 2022, 463 p.
[13] Kamont Z., Przadka K. Difference methods for first order partial differen-tial-functional equations with initial-boundary conditions. USSR Computational Mathematics and Mathematical Physics, 1991, vol. 31, pp. 1476-1488.
[14] Sampath K., Veerasamy S., Agarwal R.P. Stable difference schemes with interpolation for delayed one-dimensional transport equation. Symmetry, 2022, 14 (5), art. no. 1046. DOI: https://doi.org/10.3390/sym14051046
[15] Solodushkin S.I., Yumanova I.F., De Staelen R.H. First order partial differential equations with time delay and retardation of a state variable. Journal of Computational and Applied Mathematics, 2015, vol. 289, pp. 322-330.
[16] Pimenov V.G. Raznostnye metody reshenija uravnenij v chastnyh proizvodnyh s nasledstvennostju [Difference methods for solving partial differential equations with heredity]. Ekaterinburg, Ural University Publ., 2014, 134 p.
[17] Meleshko S.V., Moyo S. On the complete group classification of the reaction-diffusion equation with a delay. Journal of Mathematical Analysis and Applications, 2008, vol. 338, pp. 448-466. DOI: 10.1016/j.jmaa.2007.04.016
[18] Long F.-S., Meleshko S.V. On the complete group classification of the one-dimensional nonlinear Klein—Gordon equation with a delay. Mathematical Methods in the Applied Sciences, 2016, vol. 39, no. 12, pp. 3255-3270. DOI: 10.1002/mma.3769
[19] Lobo J.Z., Valaulikar Y.S. Group analysis of the one dimensional wave equation with delay. Applied Mathematics and Computation, 2020, vol. 378, art. no. 125193. DOI: 10.1016/j.amc.2020.125193
[20] Polyanin A.D., Zhurov A.I. Functional constraints method for constructing exact solutions to delay reaction-diffusion equations and more complex nonlinear equations. Communications in Nonlinear Science and Numerical Simulation, 2014, vol. 19, pp. 417-430. DOI: https://doi.org/10.1016Zj.cnsns.2013.07.017
[21] Polyanin A.D., Zhurov A.I. The functional constraints method: Application to non-linear delay reaction-diffusion equations with varying transfer coefficients. International Journal of Non-Linear Mechanics, 2014, vol. 6, pp. 267-277. DOI: 10.1016/j.ijnonlinmec.2014.09.008
[22] Polyanin A.D., Sorokin V.G., Vyazmin A.V. Nonlinear delay reaction-diffusion equations of hyperbolic type: exact solutions and global instability. Mathematical Modeling and Computational Methods, 2014, no. 4, pp. 53-73.
[23] Polyanin A.D., Sorokin V.G. Exact solutions of nonlinear partial differential equations with pantograph type variable delay. Vestnik Natsional'nogo IssledovateTskogo Yadernogo Universiteta "MIFI", 2020, vol. 9, no. 4, pp. 315-328.
[24] Polyanin A.D., Zhurov A.I. Exact separable solutions of delay reaction-diffusion equations and other nonlinear partial functional-differential equations. Communications in Nonlinear Science and Numerical Simulation, 2014, vol. 19, no. 3, pp. 409-416. DOI: 10.1016/j.cnsns.2013.07.019
[25] Polyanin A.D., Zhurov A.I. Generalized and functional separable solutions to nonlinear delay Klein—Gordon equations. Communications in Nonlinear Science and Numerical Simulation, 2014, vol. 19, no. 8, pp. 2676-2689. DOI: 10.1016/j.cnsns.2013.12.021
[26] Polyanin A.D., Zhurov A.I. Nonlinear delay reaction-diffusion equations with varying transfer coefficients: Exact methods and new solutions. Applied Mathematics Letters, 2014, vol. 37, pp. 43-48. DOI: https://doi.org/10.1016Zj.aml.2014.05.010
[27] Polyanin A.D., Zhurov A.I. New generalized and functional separable solutions to nonlinear delay reaction-diffusion equations. International Journal of Non-Linear Mechanics, 2014, vol. 59, pp. 16-22. DOI: https://doi.org/10.1016/j.ijnonlinmec.2013.10.008
[28] Polyanin A.D., Sorokin V.G. Construction of exact solutions for nonlinear equations of mathematical physics with delay using solutions of simpler equations without delay. Vestnik Natsional'nogo Issledovatel'skogo Yadernogo Universiteta "MIFI", 2020, vol. 9, no. 2, pp. 115-128.
[29] Polyanin A.D., Sorokin V.G. Nonlinear pantograph-type diffusion PDEs: Exact solutions and the principle of analogy. Mathematics, 2021, vol. 9, art. no. 511. DOI: https://doi.org/10.3390/math9050511
[30] Polyanin A.D., Sorokin V.G. A method for constructing exact solutions of nonlinear delay PDEs. Journal of Mathematical Analysis and Applications, 2021, vol. 494, art. no. 124619. DOI: 10.1016/j.jmaa.2020.124619
[31] Polyanin A.D., Sorokin V.G. Construction of exact solutions to nonlinear PDEs with delay using solutions of simpler PDEs without delay. Communications in Nonlinear Science and Numerical Simulation, 2021, vol. 95, art. no. 105634. DOI: 10.1016/j.cnsns.2020.105634
[32] Polyanin A.D., Sorokin V.G. Reductions and exact solutions of nonlinear wave-type PDEs with proportional and more complex delays. Mathematics, 2023, vol. 11, art. no. 516. DOI: https://doi.org/10.3390/math11030516
[33] Polyanin A.D., Sorokin V.G. Exact solutions of reaction-diffusion PDEs with anisotropic time delay. Mathematics, 2023, vol. 11, art. no. 3111. DOI: https://doi.org/10.3390/math11143111
[34] Tanthanuch J. Symmetry analysis of the nonhomogeneous inviscid Burgers equation with delay. Communications in Nonlinear Science and Numerical Simulation, 2012, vol. 17, pp. 4978-4987. DOI: https://doi.org/10.1016/jxnsns.2012.05.031
[35] Polyanin A.D., Zaitsev V.F. Handbook of Nonlinear Partial Differential Equations, 2nd ed. Boca Raton, CRC Press, 2012, 1912 p.
Sorokin V.G., Cand. Sc. (Phys. — Math.), researcher at Ishlinsky Institute for Problems in Mechanics RAS. e-mail: [email protected]