CC BY
15
МАТЕМАТИКА И ИНФОРМАТИКА
УДК 517.9
О. Н. Шабловский
КОЛЕБАНИЯ, РЕЗОНАНСЫ И ВОЛНЫ В НЕЛОКАЛЬНОЙ СРЕДЕ С ИСТОЧНИКАМИ
Получены новые точные решения пространственно нелокального волнового уравнения с источниками. Результаты изложены в терминах теории теплопереноса. Нелокальность задачи определяется членом уравнения, содержащим четвертую производную по пространственной координате. Рассмотрены два вида знакопеременных по отношению к температуре объемных источников энергии. Для источника технического происхождения производная от функции источника по температуре положительна, поскольку «высокие» температуры наблюдаются там, где происходит выделение энергии. Для источника, воздействующего на биологическую ткань, наклон функции источника отрицателен, так как в области «низких» температур происходит выделение тепла. Внешнее воздействие на нелокальную среду моделирует неоднородный по координате нестационарный источник энергии, подробно рассмотрены пять примеров таких источников. Аналитические решения представлены в конечном виде. Выполнено сравнение результатов воздействия монотонного и немонотонного (импульсного) по времени реономных источников. Указаны условия трансзвукового перехода скорости волны при распространении возмущения по неоднородному температурному фону. Изучены резонансные ситуации в системе «среда — источник энергии». Определены границы устойчивости/неустойчивости колебательных процессов. Получен безразмерный критерий, заключающий в себе наклон функции источника и параметр нелокальности среды и оказывающий существенное влияние на свойства корреляции «частота колебаний — параметр затухания».
New exact solutions are obtained for a spatially non-local wave equation with sources. The results are set out in the terms of the heat transfer theory. The non-locality of the problem is determined by the value of the fourth spatial derivative. We considered two types of volume energy sources which are alternating with respect to the temperature. For a technical source the derivative is positive, since «higher» temperatures arise from the energy release. For a biological source the source function is negative inclined, because a biological tissue gives off heat in the region of «lower» temperatures. The external influence on a non-local medium is simulated by spatially nonuniform energy source, and we have considered five types of such sources. The analytical solutions are presented in the finite form. The effects of monotonous and nonmonotonous (impulsive) reonomic sources are compared. The conditions for a transonic transition are indicated for the wave of perturbation in the temperature set. Resonance occurrences in the system «medium — energy source» are studied.
5
© Шабловский О. Н., 2018
Вестник Балтийского федерального университета им. И. Канта. Сер.: Физико-математические и технические науки. 2018. № 4. С. 5-14.
The limits of oscillation stability/instability are determined. We found a di-mensionless criterion including the inclination of the source function and the parameter of the medium non-locality. The criterion affects the correlation «oscillation frequency - fading parameter».
Ключевые слова: волновое уравнение, нелокальность, источник энергии, резонанс.
Keywords: wave equation, nonlocality, energy source, resonance.
Введение
Волновые уравнения с источниками (уравнения Клейна — Гордона) относятся к числу фундаментальных уравнений математической физики. Такие уравнения позволяют моделировать сложные явления в различных областях естествознания. В данной работе для определенности будем говорить о процессах волнового теплопереноса в системе «среда — источник энергии». Волновые задачи являются важным элементом динамической теории неравновесных состояний вещества. В качестве примера укажем гиперболическое уравнение теплопроводности, получаемое с помощью вариационных принципов [1; 2] и учитывающее конечную скорость распространения тепловых возмущений:
д Т д 2 Т
— + Y—2
dt Y dt2
д2 Т
= +q и, (1)
дх
где t — время; x — декартова координата; т = T - T0 есть отклонение температуры Т от ее отсчетного значения T0 = const; с — объемная теплоемкость; А — коэффициент теплопроводности; y — время релаксации теплового потока; qv — мощность внутренних источников и стоков энергии; скорость распространения тепловых возмущений равна w = (А/с у )1/2.
Физические аспекты обоснования уравнения (1) изложены в [3], а современное состояние термодинамической теории неклассических процессов переноса массы, импульса и энергии — в [4]. Здесь же приведены примеры учета нелокальных эффектов в задачах динамики тепловых волн. Частным случаем модели (1) является волновое уравнение
(2)
где x' = x/w, kv = qv/(су); с, у, w — const.
Это уравнение описывает быстрые процессы, в которых волновой механизм переноса тепла преобладает над диффузионным: yd/ dt >> 1. Учет пространственно нелокальных эффектов переноса приводит к повышению порядка уравнения (2):
2- д 4Т
- ЕХ 24тТг = К (т, x', t ). (3)
д2 т д2 т - 2 д4 т dt2 д(х')2 ЕХ д(х')
6
c
Здесь х2 = Xi/w4 , и величина ех 2 есть параметр слабой нелокальности задачи.
Обоснование вывода уравнения (3) в пределе слабой нелокальности и подробная библиография даны в [5]. Современное состояние математических исследований нелокального волнового уравнения представлено в работах [6—9].
Отметим, что dx'/dt = N/w = M есть число Маха, определяющее дозвуковой (M2 < 1) либо сверхзвуковой (M2 > 1)скоростной режим; N = dx/dt.
Далее применяем функцию источника
kv = k\ т + su (x', t), kU ^ const. (4)
Знакопеременный источник kU т моделирует процессы энергообмена в системах различной физической природы. Источник технического происхождения (kU > 0) положителен в области «высоких» температур (т > 0), где происходит подвод тепла, и отрицателен в области «низких» температур (т < 0), где тепло отводится, например, вследствие теплообмена между элементом технического устройства и окружающей средой. Источник, типичный для биологической ткани, обладает свойствами, отличающими его от объектов неживой природы [10]. Такой источник (kU < 0) выполняет уравновешивающую роль компенсатора: в области «высоких» температур (т > 0) идет теплоотвод; в области «низких» температур (т < 0) происходит выделение энергии в биологической ткани [11 — 14]. В работе [15] обсуждены некоторые закономерности влияния знака производной dk0/д т на устойчивость волновых процессов в нелокальной среде.
Неоднородный по координате x' реономный источник s„ (x', t) моделирует внешнее воздействие на нелокальную среду. В качестве примеров пространственно-временной неоднородности источника sD мы рассматриваем апериодические и периодические зависимости по x' и t.
Цель работы: построить новые аналитические решения нелокального уравнения (3) с физически содержательными источниками (4) и проанализировать условия возникновения резонансов в системе «среда — источник энергии».
Апериодическое решение
От аргументов (x', t) перейдем к новым независимым переменным g, h: g = exp(kx'+rt'), h = t'exp[n(kx'+rt')], t' = t +10, k, n, r, t0 — const, t0 > 0.
Источник энергии берем в виде, содержащем экспоненциальную зависимость от координаты и времени:
su (g) = sUgb*, b* = b1 + n; sU, b1 — const; sU * 0. (5)
7
Решение уравнений (3), (4) представляется формулами
т( h)= 6о (g)+h 0i (g), (6)
01(g)= Digbl, Di = si/(2rb*), 00 (g) = Dogb0 + D*gb*; Do, D. — const.
Здесь для показателей степени b* (то есть для bi) и b0 имеем два соотношения:
r2 - k2 =((/b2)+ ex2 k4b.2, (7)
ex2 k4b02 +(k2 - r 2 ) + kU = 0.
Следовательно решение допускает два вица варианта значения b2, а именно:
1) b02 = ki/(k4b*2), ki e > 0; 2) b02 = b*2 . (8)
Физическое истолкование данного решения: тепловое поле между неподвижной стенкой x' = 0 и волной возмущения x' = x'w(t), распространяющейся по неоднородному фону т0 = T0(x'), x' > 0. Фронт волны hgbi = Hi = const обладает такими свойствами:
Mw = dx'w/dt = (-r/k)- [1/(kb.t')], t'>0, (9)
x'w = ТГ!nH7T-jt', Hi = t0exp(b*rt0), (10)
b.k t' k
f H \b0/b. D H т w - T(x'w , t') = D0 [Hi] + ^f1 + Di Hi, (11)
(ku )w = kU T w+ ( Hi/t'), k < 0, r > 0, n < 0, b0 < 0, bi < 0.
Две функции x' = x'w(t'), т0 = Tw(t') в (10), (11) дают параметрическое представление неоднородного температурного фона т = ^(x '); t' — параметр; здесь dx0/dx' < 0; фронт волны движется с положительным ускорением вправо, в сторону положительных значений x.
Входные параметры задачи: kU, sU, e, k, n, b*, D0, D*. Из формулы (7) следует, что если r2 > 0, то в решении (6) нужно взять
r = [k2 + ((/b*2)+ ex2 kb] >0,
и это решение — апериодическое по времени. Обсудим свойства скорости волны возмущения. Пусть в начальном состоянии = 0) = 0, то есть 1/ = - тЪ* > 0. В установившемся ^ да) по времени режиме имеем
М^=- т/к > 0. Значит, звуковая точка = £1) = 1 достигается, если (- т/к) > 1. В результате простых вычислений находим
ti = t„/ (mw- i)> 0, M;
CO > i
w
8
Вывод: трансзвуковой переход «дозвук ^ сверхзвук» происходит тем быстрее, чем больше тепловое число Маха в установившемся сверхзвуковом режиме.
Два периодических решения
Рассмотрим два варианта трансформирования решения (6): переход к зависимостям, содержащим колебания по времени t'; преобразование к решению, обладающему периодичностью по координате X.
I. Если входные параметры такие, что r2 < 0, то нужно взять
r = ir*, r* = [-k2 -(/b*2)-ex2 k4b2f >0, kU < 0, e<0.
Выделяя действительную часть решения, находим
s0 = s1 exp(kb*x')cos b*r*t, X > 0, t > 0, s1 Ф 0, t = D0 exp(kb0x')cosb0r*t + D* exp(kb*x')cosb*r*t + + (s(¡2b*n)exp(kb*x'))sin b*r*t, k < 0, b0 > 0, b* > 0.
Мнимая часть решения записывается аналогичным образом. Данное решение содержит свободные параметры: k1, e, s1, k, b*, D0, D*; константа n не входит явно в эти формулы. Периодический по времени источник (13) локализован в правой окрестности X = 0 и, согласно (14), возбуждает неустойчивость по резонансному типу: tsinb*r*t. В зависимости от способа подсчета b0 (см. (8)) решение (14) содержит колебания на одной частоте b*r* либо на двух частотах: b*r* и bür*. В случае D0 = 0, D* ф 0 решение (14) определяет резонансную ситуацию, при которой температура и реономный источник sD колеблются на одной частоте. Зависимость (b*r*)2 от квадрата параметра затухания (kb*)2 немонотонная, имеет минимум при (kb*)2 = - 1/(2ex2). В случае D0 ф 0 имеем решение, содержащее колебания на частоте b0r*, отличающейся от частоты колебаний реономного источника. Из (8) ясно, что связь между квадратами параметров затухания (kbe)2 и (kb*)2 определяется величиной k^/(ex2).
Исходное ограничение r2 < 0, при котором существует решение (14), выполнено для всех (kb*)2 > 0, если 1 < 4kl ex2. Если же 1 > 4kl ex2, то должны быть удовлетворены неравенства (kb* )2 > (kb* либо 0 <(kb* )2 <(kb* )2, где
(kb* )2,2 =( 1 ±yl 1 - 4kl ex2 )Д- 2 ex 2 ),
причем знакам «+» и «-» соответствуют индексы 1 и 2.
Вывод: в структуре данного решения безразмерный комплекс k1 ex2 характеризует мультипликативное взаимодействие источника k¿ тс нелокальной средой, и его числовое значение влияет на свойства корреляции «частота колебаний — параметр затухания».
(12)
(13)
(14)
9
II. В решении (6) сделаем замену к = П. Тогда, выделив действительную часть решения, получаем
= б1 ехр(тЪ*£)со81Ъ*х', х' > 0, £ > 0, (15)
т = О0 ехр(тЪ^ )со81Ъ0 х'+О* exp(тЪ.í )со81Ъ* х'+ (16)
+ (/2тЪ*)ехр(тЪ*£)соя 1Ъ*х'.
Аналогично варианту I здесь допускаются два значения константы
Ъ2:
1) Ъ2 = кЦ(ех214Ъ*2), к\ е > 0; 2) Ъ2 = Ъ*2 . (17)
Выбор знаков следующий: Ъ0 > 0, Ъ* > 0, т < 0, I > 0, к^ > 0, е > 0; величина т2 подсчитывается по формуле (7), в которой к2 = - I2; вид начального температурного поля т^ = 0, х') определяется выбором произвольных постоянных О0, О*; £0 = 0.
Внешний для нелокальной среды источник монотонно релакси-рует во времени и является периодическим по координате X (см. (15)). Решение (16) имеет всплеск амплитуды, который определяется функцией В(£) = íexp(тЪ*í). Эта функция имеет максимум при £ = = 1/(- тЪ*), В(£т) = 1/(- тЪ*е). Если тЪ* — 0, то есть 1т —> да, В(£т) —> да, то наблюдается неограниченный рост амплитуды колебаний в решении (16). Обсудим связь между параметром затухания тЪ* и частотой 1Ъ *:
(тЪ* )2 = ех2 (Ъ*2) - /2Ъ*2 + к1.
Если 1 < 4кЦ, ех2, то при всех 12 Ъ*2 > 0 будет (тЪ*)2 > 0, и процесс устойчив. Если 1 > 4кЦ, ех2, то /2Ъ? = (Ъ* ) есть нижняя граница устойчивости для области высоких частот /2Ъ*2 > (Ъ* ); вместе с тем /2Ъ*2 = (ъ* )2 есть верхняя граница устойчивости для области низких частот 0 < /2Ъ2 < (ъ* )2 • Здесь
(1Ъ*)2,2 =( 1 ±д/1 -4ки ех2)/(ех2),
причем знакам «+» и «-» соответствуют нижние индексы 1 и 2.
Вывод: чем дольше релаксирует (затухает) источник тем большую амплитуду колебаний приобретает температура (16).
Корреляция «параметр затухания — частота» и резонанс
Рассмотрим источник энергии
к0 = к! т + кр^, р1 = const.
(18)
Решение нелокального волнового уравнения (3) имеет вид
т = D0 gb0 + D* gn + h 01, (19)
01 = pU/S1, s1 = (nr )2 - (nk)2 - ех 2 (nk)4 - kU, (20)
1 / 2
D* = - 2nrpU/s1,
где b2 следует вычислять из второго соотношения в (7); оценки параметров, при которых существует положительный корень b02 > 0 , здесь не приводятся; Do, n — произвольные постоянные; x' > 0, t > 0.
Пусть x' = 0 — неподвижная стенка, h = hw = const — волна, распространяющаяся вправо по неоднородному фону т = т0^'):
x'w (t') = (1/ nk )ln(h^t')- (rt'l k), hw = t0 exp(nrt0) > 0, (21)
M-(22)
x'w (t = 0) = 0; r>0, k<0, n<0, b0<0, 1+nrt0<0; тw(t') = D0(h^i')b^" + (D* h^/i')+ hw 01. (23)
В этом решении функция т0^') представлена параметрически двумя зависимостями (21), (23). Для скорости волны (22) переход «дозвук — сверхзвук» обладает такими же свойствами, как в уже изученном решении (10), (11), потому что формулы (22) и (9) имеют одинаковую структуру.
Неустойчивость решения (19) возникает, если s1 = 0 (см. (20)):
(nr )2 = ех 2 (nk )4 + (nk )2 + kU > 0. (24)
Если правая часть этого выражения отрицательная, то решение всюду устойчивое, si > 0. Здесь nr и nk являются параметрами неоднородности по времени и координате. При kU > 0, е > 0 для каждого (nk)2 > 0 существует (nr)2 (см. (24)), дающее неустойчивость. При kU < 0, е < 0 решение устойчивое в случае 1 < 4kU ех2 при всех (nk)2>0; другой вариант устойчивости: 1 > 4k1 ех2, 0 < (nk)2 < (nk^ либо (nk)2 > (nk)2 . В интервале 0 < (nk) < (nk)2 < (nk)2 для каждого (nk)2 существует значение (nr)2, дающее неограниченно большую амплитуду (см. (24)). Если kU < 0, е > 0, то решение устойчивое для 0 <(nk)2 < (nk )2, а для (nk)2 > (nk)2 существуют неустойчивые решения. Если kU > 0, е < 0, то решение устойчивое при (nk )2 > (nk)2, а в интервале 0 < (nk)2 < (nk)2 существуют неустойчивые решения. Отметим, что в представленных здесь вариантах важны знаки величин е и kU, числовое значение kU ех2 имеет количественное влияние на эти оценки. Границы упоминавшихся интервалов подсчитываются по формуле
11
(nk )2,2 =(-1 ±д/1 - 4k1 ех 2 )/( ех 2).
12
Вывод: реономные источники (5) и (18) изменяются монотонно и немонотонно по времени, что приводит к существенным различиям в условиях устойчивости / неустойчивости теплового поля.
Трансформируем решение (18) — (20), взяв к = И и выделяя действительную часть:
= р1* exp(nтt )соя п1х', (25)
т = О0 exp(тЪ0t )со8 /Ъ0х'+О* exp(nтt )со8 п1х'+ (26)
+ 01 í exp(nтí)cos п1х',
01 = р£/«2 , О* =- 2птр\1 в2 = (пт)2 + (п1 )2 - ех2 (п1 )4 - к1, / > 0, п > 0, т < 0, Ъ0 > 0, £0 = 0.
Напомним, что при замене к = й меняется вид второго соотношения в (7). Именно это соотношение несет информацию о взаимодействии источника кЦ, т с нелокальной средой.
Пространственно-периодический источник (25) отличается от (15) немонотонным изменением во времени. Анализ поведения решения (26) дает следующие результаты. Неограниченная амплитуда колебаний получается при э2 = 0:
(пт )2 = ех2 (п1 )4 - (п1 )2 + к1 > 0 . (27)
Если правая часть этого выражения отрицательная, то решение устойчивое, э2 > 0. Пусть к1 > 0, е > 0: если 1 < 4к1 ех2, то для каждого (пТ)2 > 0 существует (см. (27)) резонансное значение (пт)2; если 1 > 4 к1 ех2 , то имеем устойчивость при (п/) < (п/)2 < (п/), а значения (пт)2, дающие неустойчивость, существуют на интервалах 0 < (п/)2 < (п/) и (п/)2 > (п/)2. Пусть к1 < 0, е < 0: решение устойчивое при всех (п/)2 > 0. Если к0 < 0, е > 0, то решение устойчивое для 0 < (п/)2 < (п/)2, а для (п/)2 > (п/) существуют резонансные значения (пт)2, которые подсчи-тываются по формуле (27). Если к1 > 0, е < 0, то решение устойчивое при (п/)2 > (п/)2, а в интервале 0 < (п/)2 < (п/)2 существуют неустойчивые решения. Границы названных интервалов определяются выражением
(п/)22 =( 1 ±у1 1 -4к1 ех2 )/(2ех2).
Вывод: пространственно-периодические решения (25), (26) и (15), (16) при О0 ф 0 содержат колебания на собственной частоте /Ъ0, отличающейся от частоты колебаний источника различия в нестационарных свойствах этого источника приводят к различиям в условиях возникновения резонансов.
Заключение
Получены точные аналитические решения нелокального волнового уравнения (3) с источниками (4) вида (5), (13), (15), (18) и (25). Построены волны возмущения, содержащие трансзвуковой переход «дозвук — сверхзвук» по отношению к скорости w. Проанализированы условия возникновения резонансных ситуаций. Установлена важная роль безразмерного комплекса k1 ех 2 при оценке границ устойчивости / неустойчивости колебаний. В частном случае е = 0 (классическое волновое уравнение) поведение представленных здесь решений определяется прежде всего знаком параметра источника k1. Анализ этих простых вариантов здесь не приводится.
Список литературы
1. Глазунов Ю. Т. Вариационный принцип явлений взаимосвязанного тепло-и массопереноса, учитывающий конечную скорость распространения возмущений // Инженерно-физический журнал. 1981. Т. 40, № 1. С. 134—138.
2. Яворский Н. И. Вариационный принцип для вязкой теплопроводной жидкости с релаксацией // Известия АН СССР. Механика жидкости и газа. 1986. № 3. С. 3 — 10.
3. Никитенко Н. И. Проблемы радиационной теории тепло- и массоперено-са в твердых и жидких средах // Инженерно-физический журнал. 2000. Т. 73, № 4. С. 851 — 859.
4. Jou D., Casas-Vazquez J., Lebon J. Extended Irreversible Thermodynamics. Berlin ; Heidelberg, 2001.
5. Алфимов Г. А. Нелокальное уравнение синус-Гордона: решения типа «кинк» в пределе слабой нелокальности // Нелинейная динамика. 2009. Т. 5, № 4. С. 585—602.
6. Свешников А. Г., Альшин А. Б., Корпусов М. О., Плетнер Ю. Д. Линейные и нелинейные уравнения соболевского типа. М., 2007.
7. Мамчуев М. О. Необходимые нелокальные условия для диффузионно-волнового уравнения // Вестник Самарского государственного технического университета. Естественнонаучная серия. 2014. № 7. С. 45 — 59.
8. Керефов М. А., Геккиева С. Х. Первая краевая задача для неоднородного нелокального волнового уравнения // Вестник Бурятского государственного университета. Математика, информатика. 2016. № 4. С. 76—86.
9. Дрегля А. И., Сидоров Н. А. Идентификация динамики внешней силы при моделировании колебаний // Известия Иркутского государственного университета. Сер. Математика. 2017. Т. 19. С. 105 — 112.
10. Pennes H. H. Analysis of Tissue and Arterial Temperature in the Resting Human Forearm // J. Appl. Phisiol. 1948. Vol. 1, iss. 2. P. 93 — 122.
11. Tung M.M., Trujillo М., Lopez Molina J. A. et al. Modelling the Heating of Biological Tissue Based on the Hyperbolic Heat Transfer Equation // Math. Comput. Model. 2009. Vol. 50, iss. 5—6. P. 665—672.
12. Yang C.-Y. Boundary Estimation of Hyperbolic Bio-Heat Conduction // Int. J. Heat Mass Transfer. 2011. Vol. 54, iss. 11 — 12. P. 2506 — 2513.
13. Lin S.-Y., Chou T.-M. Numerical Analysis of the Pennes Bioheat Transfer Equation on Skin Surface // 2015 Third International Conference on Robot, Vision and Signal Processing (RVSP). IEEE, 2015. P. 71 — 74.
13
14. Mochnacki B., Ciesielski M., Piasecka-Belhayat A. Numerical Solution of the Bio-Heat Transfer Equation with Uncertain Parameters Using the Sensitivity Analysis Method // Defect Diffus. Forum. 2017. Vol. 379. P. 39-47.
15. Шабловский О. Н. Нелокальность и возникновение резонансов в динамике волн // Фундаментальные физико-математические проблемы и моделирование технико-технологических систем. 2017. Вып. 18. С. 125 — 138.
14
Об авторе
Олег Никифорович Шабловский — д-р физ.-мат. наук, проф., Гомельский государственный технический университет им. П. О. Сухого, Беларусь. E-mail: [email protected]
The author
Prof. Oleg N. Shablovsky, Sukhoi State Technical University of Gomel, Belarus. E-mail: [email protected]
