Аннотация:
Получено решение типа контрастной структуры в системе уравнений теории эффекта бареттирования, включающей нелинейное сингулярно возмущенное параболическое уравнение и дополнительное нелокальное интегральное соотношение. Доказано существование решения с внутренним переходным слоем, построено асимптотическое приближение этого решения. Получены оценки основных физических параметров модели, совпадающие с экспериментальными данными, а также с оценками, полученными ранее другими методами.
Ключевые слова:
сингулярно возмущенные параболические уравнения, уравнения реакция-диффузия, контрастные структуры, внутренние слои, асимптотические методы, эффект бареттирования.
Одной из актуальных задач теории сингулярных возмущений в настоящее время является исследование нелинейных сингулярно возмущенных уравнений в частных производных, решения которых имеют пограничные и внутренние слои. Такие уравнения представляют интерес как в качественной теории дифференциальных уравнений, так и во многих прикладных задачах, например в математических моделях химической кинетики, синергетики, нелинейной теории волн, биофизики и других областей физики [1]. В настоящей статье рассматривается математическая модель, которая описывает известный в электротехнике эффект бареттирования [2]–[4], возникающий в системе железо+водород.
Указанные выше модели и процессы описываются нелинейными параболическими уравнениями с малыми параметрами при производных. Решения таких задач могут содержать узкие области быстрого изменения параметров: пограничные или внутренние переходные слои (контрастные структуры) различных типов – стационарные или движущиеся фронты (см. [5]–[7] и ссылки в этих работах). Причиной образования переходных слоев (контрастных структур) в сингулярно возмущенных моделях типа реакция-диффузия-адвекция может служить выполнение условия баланса реакции в некоторой точке или на некоторой кривой, лежащей в области рассмотрения [7], [8], либо баланса адвекции [9], либо баланса адвекции и реакции [10].
В работах [11]–[13] был рассмотрен так называемый критический случай для нелинейных сингулярно возмущенных уравнений реакция-диффузия-адвекция, когда условие баланса реакции выполняется тождественно, т. е. в любой точке области. В этом случае процедура определения точки локализации внутреннего переходного слоя несколько усложняется, и главный член асимптотики положения внутреннего слоя определяется в процессе построения высших членов асимптотики. Более того, в пространственно однородном одномерном случае с помощью методов, применяемых в работах [11]–[13], невозможно определить положение внутреннего переходного слоя.
В рассматриваемой в настоящей статье модели для нелинейного сингулярно возмущенного уравнения реакция-диффузия также реализуется пространственно однородный критический случай. Однако наличие дополнительного нелокального интегрального условия понижает степень автономности задачи по пространственной переменной и позволяет определить положение внутреннего переходного слоя. Доказано существование решения с внутренним переходным слоем, построено асимптотическое приближение этого решения. Получены оценки основных физических параметров модели, совпадающие с экспериментальными данными, а также с оценками, полученными ранее другими методами.
Статья содержит пять разделов, в которых сформулирована физическая модель и математическая постановка задачи, проведено построение асимптотического приближения решения с внутренним переходным слоем и доказано существование решения типа контрастной структуры с асимптотикой указанного вида. В разделе 5 проведено обсуждение полученных результатов и сравнение параметров модели, рассчитанных на основе построенного асимптотического приближения, с экспериментальными данными.
2. Формулировка модели и постановка задачи
Эффект бареттирования заключается в способности железной проволоки, находящейся в атмосфере водорода и являющейся элементом электрической цепи, стабилизировать ток в этой цепи в некотором диапазоне приложенных напряжений [2], [3].
Пусть тонкая железная нить длины $l$ помещается в водородную среду и включена в цепь источника регулируемого напряжения. Математически данная система описывается следующими уравнениями:
$$
\begin{equation}
e =J \int_0^l r(T)\, d y.
\end{equation}
\tag{2}
$$
Уравнение (1) выражает условие баланса тепла на нити, здесь $T$ и $T_0$ – температуры нити и окружающей среды, $c$ – удельная теплоемкость, $\rho$ – плотность, $\lambda$ – коэффициент теплопроводности, $d$ – диаметр нити, $\alpha$ – коэффициент теплоотдачи, $J$ – ток, текущий в цепи, $\tau(T)$ – сопротивление единицы длины нити. Зависимость $\tau(T)$ для железной проволоки, находящейся в среде водорода, нелинейна и имеет характерный $S$-образный вид [3], [4], причем при некоторых значениях параметров $\alpha, d, J$ возможно, что однородное стационарное уравнение (1) (т. е. при $\partial T/\partial t=0$, $\partial T/\partial y=0$) имеет три изолированных корня $T_1, T_2, T_3$.
Интегральное соотношение (2) есть закон Ома для участка цепи, $e$ – напряжение на концах нити, $\int_0^{l} r(T(y))\, d y$ – полное сопротивление участка. Величина напряжения $e$ предполагается постоянной во времени ($e$ может быть задано различными значениями из некоторого диапазона).
Будем считать контакты на концах нити теплоизолированными:
Заметим, что для набора параметров модели, реализованных в эксперименте, описанном в работе [2], $\varepsilon^2 \sim 10^{-4}$, т. е. можно считать $\varepsilon$ малым параметром.
Далее будет рассмотрен только стационарный случай, для которого в работах [2] и [14] приведен ряд результатов, касающихся теоретического исследования задачи (1)–(3), однако некоторые вопросы остались невыясненными. В частности, был получен только главный член формального асимптотического приближения и не было доказано существование решения, близкого к построенной асимптотике. В настоящей работе проведено уточнение асимптотики для этой задачи, а также доказано существование решения с построенной асимптотикой и получены оценки точности асимптотического приближения.
3. Формальная асимптотика решения
Для построения асимптотического приближения используется подход, предложенный в [14], основанный на алгоритме [5], а для обоснования асимптотики применяется метод сшивания (см., например, [6]).
Рассмотрим следующую интегро-дифференциальную систему уравнений:
где $V>0$ – заданная постоянная, а $\varepsilon>0$ – малый параметр.
Будем предполагать выполненными следующие условия.
Условие 1. Пусть $R(u)$ – достаточно гладкая положительная функция в области $D_u$ – рассматриваемой области определения переменной $u$.
Условие 2. Пусть вырожденное уравнение $f(u,j):=u-j^2 R(u)=0$ при $u\in D_u$ и $j\in D_j$ (где $D_j$ – некоторый интервал изменения $j$) имеет три упорядоченных корня $\varphi^{(-)}(j)<\varphi^{(0)}(j)<\varphi^{(+)}(j)$, причем выполнены неравенства
Для определенности будем рассматривать переход от корня $\varphi^{(-)}$ к корню $\varphi^{(+)}$. Положение внутреннего переходного слоя $x^*(\varepsilon)$ определим как точку пересечения решения и корня $\varphi^{(0)}(j)$ вырожденного уравнения, т. е. положим
Вид асимптотики обусловлен автономностью уравнения (5) и граничными условиями Неймана. Заметим, что внутренний переходный слой представлен только членом нулевого порядка $Q_0^{(\pm)}(\xi,j)$. C учетом граничных условий Неймана легко видеть, что в асимптотике полностью отсутствует пограничная часть, описывающая пограничные слои на концах отрезка $[0,1]$.
Согласно алгоритму [5] функции $Q_0^{(\pm)}$ определяются как решения задач
Процедура $C^1$-сшивания асимптотики справа и слева от точки $x^*(\varepsilon)$ приводит к уравнению для определения главного члена асимптотического приближения тока $j$. Действительно, имеем
Как показано ниже, неравенство (14) обеспечивает локализацию переходного слоя внутри интервала $(0;1)$. В силу условия 2 неравенство (14) может быть также записано в виде
где $n$ – произвольное натуральное число. Такая асимптотическая оценка следует из экспоненциального затухания функций $Q_0$ (см., например, [5]).
Приравняв к нулю коэффициент при $\varepsilon^0$ в равенстве (17), получим выражение для главного члена асимптотического приближения положения внутреннего переходного слоя $x_0$:
Продолжая процесс, можно показать, что все коэффициенты более высокого порядка по $\varepsilon$, т. е. $x_2, x_3,\dots$, описывающие асимптотическое приближение положения внутреннего слоя $x^*(\varepsilon)$, равны нулю.
Итак, формальная асимптотика задачи (5)–(7) полностью построена.
4. Обоснование асимптотического приближения
Очевидно, что функция $U(x,\varepsilon;j_0,x_1^*),$ где $x_1^*(\varepsilon ):=x_0+\varepsilon x_1$, является дважды непрерывно дифференцируемой на интервале $(0;1)$ и вместе с $j=j_0$ удовлетворяет точно уравнению (5), т. е.
Эта функция удовлетворяет также граничным условиям (7) и интегральному уравнению (6) по невязке произвольного натурального порядка по $\varepsilon$, т. е. для любого натурального $n$ выполнено
Пусть $j_{\gamma}(\varepsilon):=j_0+\varepsilon^m \gamma$, $x_{\delta}(\varepsilon):=x_0+\varepsilon x_1+\varepsilon^m\delta$, где $m$ – произвольное натуральное число, а $\delta$ и $\gamma$ – некоторые постоянные, не зависящие от $\varepsilon$.
Слева от переходного слоя рассматривается вспомогательная задача
Известно, что при выполнении условий 1–3 для произвольных фиксированных $m$, $\gamma$ и $\delta$ задачи (23) и (24) имеют соответственно решения $u^{(-)}(x,\varepsilon;j_{\gamma}(\varepsilon),x_{\delta}(\varepsilon))$ и $u^{(+)}(x,\varepsilon;j_{\gamma}(\varepsilon),x_{\delta}(\varepsilon))$ с пограничным слоем вблизи точки $x=x_{\delta}(\varepsilon)$, для которых имеют место асимптотические представления
Тогда дважды непрерывно дифференцируемая функция $u(x,\varepsilon;j_{\gamma^*}(\varepsilon),x_{\delta^*(\gamma^*)}(\varepsilon))$ и величина $ j_{\gamma^*}(\varepsilon)$ будут удовлетворять задаче (5)–(7).
Для доказательства существования $\delta^*$ достаточно заметить, что для произвольного $\gamma$ найдется достаточно большое $\delta_0>0$ такое, что при достаточно малых $\varepsilon$ при любых натуральных $n$ выполняются неравенства
Знак в последних двух неравенствах определяется условием 2. В силу непрерывной и монотонной зависимости величины $M[u(x,\varepsilon;j_{\gamma}(\varepsilon),x_{\delta}(\varepsilon)),j_{\gamma}(\varepsilon)]$ от $\delta$ существует такая непрерывная функция $\delta=\delta^*(\gamma)$, что равенство (27) выполнено.
Из представления (26) и условий 1–3 следуют неравенства, справедливые для произвольного положительного $\gamma_0$ и достаточно малых $\varepsilon$:
В силу непрерывной зависимости функции $\partial{u(x,\varepsilon;j_{\gamma}(\varepsilon),x_{\delta^*(\gamma)}(\varepsilon))}/\partial x$ от $\gamma$ из последних неравенств следует, что существует $\gamma^*$, удовлетворяющая равенству (28).
Таким образом, доказана следующая теорема.
Теорема 1. Пусть выполнены условия 1–3. Тогда при достаточно малом $\varepsilon >0$ существуют функция $u(x,\varepsilon)$ и величина $j(\varepsilon)$, удовлетворяющие задаче (5)–(7), для которых при произвольном натуральном $m$ справедливы оценки
Замечание 1. Нетрудно показать, что ток $j(\varepsilon)$, удовлетворяющий задаче (5)–(7), зависит от напряжения $V$ и от малого параметра $\varepsilon$, т. е. является некоторой функцией $j=J(V,\varepsilon)$, не обращающейся тождественно в постоянную.
Предположим противное, тогда $(\partial J/\partial V) \equiv 0$. Это предположение противоречит тождеству
Таким образом, величина $O(\varepsilon^m)$, фигурирующая в приведенном в теореме соотношении (32), является функцией, зависящей от $V$ и, очевидно, от $\varepsilon$, причем эта функция не обращается тождественно в постоянную.
5. Заключение
Асимптотический подход применен для доказательства существования решения интегро-дифференциальной системы (5)–(7). Известные результаты, сформулированные в работах [1] и [2], получили строгое математическое обоснование.
Во-первых, видно, что данная система способна стабилизировать ток $j$ в некотором диапазоне изменения напряжения $V$ (эффект бареттирования).
Действительно, если заменить в (6) $V$ на $V_1$, то это не повлияет на значение тока $j_0$, определяемое в процессе построения асимптотики, а изменит только положение точки перехода в нулевом порядке
и, следовательно, изменит и общее сопротивление проволоки $\int _0^1 R(u)\,dx$.
Из теоремы 1 и замечания 1 следует, что ток $j$ все-таки зависит от напряжения $V$, но чрезвычайно слабо: эта зависимость имеет сверхстепенной характер относительно малого параметра $\varepsilon$.
Во-вторых, соотношение (16) позволяет оценить ширину зоны бареттирования, т. е. того диапазона, в котором ток можно считать (с учетом предыдущего замечания) не зависящим от приложенного напряжения,
Эта формула совпадает с формулой, приведенной в работе [2], которая была получена из других соображений.
Отметим, что полученные результаты остаются справедливыми только в предположении, что параметр $\varepsilon$ достаточно мал. Если последнее условие не выполнено, границы зоны бареттирования, указанные в неравенстве (16), могут не достигаться.
Заметим также, что режим работы (т. е. величина тока бареттирования $j_0$ и положение точки перехода между высоко- и низкотемпературной зонами) практически не зависит от условий на концах нити. Поэтому замена краевых условий второго рода на условия первого рода приведет к асимптотически малыми изменениям в оценках параметров бареттера.
Методы, развитые в настоящей работе, могут быть применены для асимптотического анализа более сложных систем, содержащих интегральное условие типа (6), например возникающих в дрейфо-диффузионных моделях полупроводников при условии заданного напряжения (см., например, [15], [16]). В частности, интерес для дальнейшего исследования вызывает наличие неоднородности в реактивном слагаемом в системах указанного типа.
Конфликт интересов
Авторы заявляют, что у них нет конфликта интересов.
Список литературы
1.
Н. А. Городецкая, Н. Н. Кралина (ред.), Автоволновые процессы в системах с диффузией, ИПФ АН СССР, Горький, 1981
2.
Ю. Е. Володин, А. И. Вольперт, А. И. Иванова, В. П. Филлипенко, К теории эффекта бареттирования, Препринт, ОИХФ, Черноголовка, 1988
3.
В. В. Барелко, В. М. Бейбутян, Ю. Е. Володин, Я. Б. Зельдович, “Тепловые волны и неоднородные стационарные состояния в системе $\mathrm{Fe} + \mathrm{H}_2$”, Автоволновые процессы в системах с диффузией, ИПФ АН СССР, Горький, 1981, 135–148
4.
D. G. Löffler, L. D. Schmidt, “Steady state multiplicity of an electrically heated iron wire”, Chem. Eng. Sci., 31:12 (1976), 1207–1209
5.
A. Б. Васильева, В. Ф. Бутузов, Асимптотические методы в теории сингулярных возмущений, Высшая школа, М., 1990
6.
А. Б. Васильева, В. Ф. Бутузов, Н. Н. Нефедов, “Контрастные структуры в сингулярно возмущенных задачах”, Фундамент. и прикл. матем., 4:3 (1998), 799–851
7.
Н. Н. Нефедов, “Развитие методов асимптотического анализа переходных слоев в уравнениях реакция-диффузия-адвекция: теория и применение”, Ж. вычисл. матем. и матем. физ., 61:12 (2021), 2074–2094
8.
N. N. Nefedov, E. I. Nikulin, “Existence and stability of periodic contrast structures in the reaction-advection-diffusion problem”, Russian J. Math. Phys., 22:2 (2015), 215–226
9.
N. N. Nefedov, L. Recke, K. R. Schneider, “Existence and asymptotic stability of periodic solutions with an interior layer of reaction-advection-diffusion equations”, J. Math. Anal. Appl., 405:1 (2013), 90–103
10.
Н. Н. Нефедов, Е. И. Никулин, “Существование и асимптотическая устойчивость периодических двумерных контрастных структур в задаче со слабой линейной адвекцией”, Матем. заметки, 106:5 (2019), 708–722
11.
Н. Н. Нефедов, Е. И. Никулин, “Существование и устойчивость периодических контрастных структур в задаче реакция-адвекция-диффузия в случае сбалансированной нелинейности”, Дифференц. уравнения, 53:4 (2017), 524–537
12.
В. Т. Волков, Н. Н. Нефёдов, “Развитие асимптотического метода дифференциальных неравенств для исследования периодических контрастных структур в уравнениях реакция-диффузия”, Ж. вычисл. матем. и матем. физ., 46:4 (2006), 615–623
13.
N. N. Nefedov, K. Sakamoto, “Multi-dimensional stationary internal layers for spatially inhomogeneous reaction-diffusion equations with balanced nonlinearity”, Hiroshima Math. J., 33:3 (2003), 391–432
14.
В. Т. Волков, Асимптотика периодических режимов в системах с малой диффузией и теплопроводностью, Дисс. … канд. физ.-матем. наук, МГУ им. М. В. Ломоносова, М., 1990
15.
Е. И. Никулин, “Контрастные структуры в задаче реакция-адвекция-диффузия, возникающей в дрейфо-диффузионной модели полупроводника, в случае негладкой реакции”, ТМФ, 215:3 (2023), 360–376
16.
C. M. Cuesta, C. Schmeiser, “Stability of solitary waves in a semiconductor drift-diffusion model”, Siam J. Appl. Math., 68:5 (2008), 1423–1438
Образец цитирования:
Е. И. Никулин, В. Т. Волков, А. Г. Никитин, “О контрастных структурах в задаче теории эффекта бареттирования”, ТМФ, 220:1 (2024), 154–163; Theoret. and Math. Phys., 220:1 (2024), 1193–1200