MSC 47J05 DOI: 10.14529/mmp140409
ON THE UNIQUENESS OF A NONLOCAL SOLUTION IN THE BARENBLATT - GILMAN MODEL
E.A. Bogatyreva, South Ural State University, Chelyabinsk, Russian Federation, [email protected],
I.N. Semenova, Ural State Pedagogical University, Yekaterinburg, Russian Federation, semenova i [email protected]
This article deals with the question of uniqueness of a generalized solution to the Dirichlet-Cauchy problem for the Barenblatt-Gilman equation, which describes nonequilibrium countercurrent capillary impregnation. The unknown function corresponds to effective saturation. The main equation of this model is nonlinear and implicit with respect to the time derivative, which makes it quite hard to study. In a suitable functional space, the Dirichlet-Cauchy problem for the Barenblatt-Gilman equation reduces to the Cauchy problem for a quasilinear Sobolev-type equation. Sobolev-type equations constitute a large area of nonclassical equations of mathematical physics. The techniques used in this article originated in the theory of semilinear Sobolev-type equations. For the Cauchy problem we obtain a sufficient condition for the existence of a unique generalized solution. We establish the existence of a unique nonlocal generalized solution to the Dirichlet-Cauchy problem for the Barenblatt-Gilman equation.
Keywords: Barenblatt-Gilman equation; quasilinear Sobolev-type equation; generalized solution.
Introduction
The Barenblatt-Gilman mathematical model [1] amounts to the equation
ut - Aa(A$(u))t = aA$(u) (0.1)
with the Dirichlet-Cauchy conditions
u(x, t) = 0, (x, t) e dQ x (0, t), (0.2)
u(x, 0) = u0(x),x e Q, (0.3)
where Q C is a bounded region with boundary of class C™, the function $(u) = \u\p-2u with p > 2 is monotonely increasing and smooth, while t, a, and A are positive real parameters.
This equation models nonequilibrium countercurrent capillary impregnation. In a suitable functional space we reduce (0.1) - (0.3) to the Cauchy problem
u(0) = uo (0.4)
for the abstract operator differential equation
d
— (L(u)) + M (u) = 0, (0.5)
dt
where L(u) = Au + XM(u) with A G R+. This is a quasilinear Sobolev-type equation. Implicit equations attract many researchers [2-5]. This article continues [7], which dealt with the question of the existence of a generalized solution to (0.1)-(0.3).
1. Abstract Cauchy Problem
Consider a real Hilbert space H = (H, (•, •}) identified with its dual and equipped with a dual pair of reflexive Banach spaces U = (U, || • ||) and U* = (U, || • ||*) such that we have a continuous dense embedding U ^ H ^ U*. Take an operator M G Cr+1(U; U*) with r G N and an operator A G L(U; U*).
Definition 1. A vector function u G L^(0,r; U) such that G L2(0,t; H) for t G R+ is called a generalized solution to problem (0.4), (0.5) whenever it satisfies
dt
(JtL(u(t)),v^ + (M(u(t)),v) = 0 for almost all t G (0,т)
and
(u(0) — u0, v) = 0 for a 11 v G U.
Theorem 1. [6] Take two operators A G L(U; U*) and M G Cr+1(U; U*) with r G N and suppose that at u0 G U the operator N(u0) : U ^ U* is a toplinear isomorphism. Then there exists a unique local solution to problem (0.4), (0.5) passing through the point u0.
Theorem 2. [7] Take an operator M G Cr+1 (U; U*) with r G N which is s-monotone, p-coercive, and homogeneous of degree k G R+, and, furthermore, its Frechet derivative is symmetric, and a positive definite symmetric operator A G L(U; U*). If problem (0.4), (0.5) has a unique local solution on some interval (—т0, т0) with т0 G R+ then it has a generalized solution.
Condition 1. If there exists a function F G C[0, то), possibly after a change on a negligible set, such that F(s) > 0 for almost all s G [0, то) and
\\M(u) — M(v)||* < F(so)||u — v\\
for almost all s0 G [0, то) and all u = u(s0),v = v(s0) G U then we say that the operator M : U ^ U* satisfies condition 1.
Definition 1 of a generalized nonlocal solution to problem (0.4), (0.5) is equivalent [8] to
Definition 2. Refer as a weak generalized solution to the Cauchy problem (0.4), (0.5) to a function u(t) G Ь^(0,т,U) with du G Ь2(0,т,U) satisfying
T
f d
I (- (L(u),w) + (M(u),w))<p(t)dt = 0,
u(0) = uo, Vw G U, Уф G L2(0,т). 114 Вестник ЮУрГУ. Серия «Математическое моделирование и программирование»
o
Theorem 3. Suppose that for u0 G U a generalized solution to problem, (0.4), (0.5) exists. If the operator M G Cr+1(U;U*), with r G N, ¿s s-monotone and satisfies condition 1, while the operator A G L(U; U*) is positive definite, then this solution is unique.
Proof. Consider two distinct weak generalized solutions ui and u2 to problem (0.4), (0.5) in the sense of Definition 2.
Take p(t) = 1 E L2(0,t). Then for almost all t E (0,t) we infer from (1.1) that
t
f d
(— (L(ui),w) + (M(ui),w))ds = 0, Vw E U,i = 1, 2; (1.2)
J ds 0
where L(u) = Au + XM(u). Integrating (1.2) yields
t
(Aui(t) + XM(ui(t)),w) + j (M(ui(s)),w) ds = (Au0 + XM(u0),w), Vw E U,i = 1, 2.
0
(1.3)
Subtract from (1.3) for u the same equation for u2 and put v = u\ — u2- Then
t
(Av,w) + X (M(ui(t)) — M(u2(t)),w) + j (M(ui(s)) — M(u2(s)),w) ds = 0,Vw E U.
0
(14)
Since ui E L^(0,t,U), it follows that v(a) E U for almost all a E (0,t). Therefore, (1.4) yields
a
(Av,v) + X (M(ui(a)) — M(u2(a)),v) = j (MMs)) — M(ui(s)),v(a)) ds. (1.5)
0
Since A is a positive definite operator for all u E U, we have
(Au,u) > m\\u\\2H,m> 0. (1.6)
Ms
(M(ui(a)) — M(u2(a)),v)> 0. (1.7)
By (1.5)—(1.7) we arrive at
a
\\v(a)\\2 < m- j (M(u2(s)) — M(ui(s)),v(a)) ds. 0
Now Holder's inequality implies that
a
\\v(a)\\2 < m-i\\v(a)\\ [ \\M(^(s)) — M(ui(s))\\ds.
Since M satisfies condition 1, we infer that
\\v(a)\\2 <\\v(a)\\ F(S)\\v(S)\\dS,
and then the Gronwall-Bellman lemma [9] shows that u1 = u2 almost everywhere.
□
2. The Dirichlet-Cauchy Problem
for the Barenblatt-Gilman Equation
Let us reduce problem (0.1) - (0.3) to the abstract problem (0.4), (0.5). Take H = W-l(Q) with the inner product
(u,v) = u(—A) vdx,u,v E H
(2.1)
where (—A)"1 is the Green operator of the homogeneous Dirichlet problem for the Poisson equation Au = / in Q. Choose the space Lp(Q) as U sinee Lp(Q) — L2(Q) for p E [2, in the case of a bounded region Q, which we have, while L2(Q) — W"l(Q). Then we can take as U* the dual space to U with respect to the duality (2.1). In this case the embedding U c—y H ^ U* for p E [2, is dense and continuous. Define the operator A as
(Au,v) = (u,v) = u(-A) vdx,u,v E U
M
(M(u),v) = a \u\p 2uvdx, u, v E U.
M
Proof. There exists a nonnegative continuous function g : R2 ^ R such that
\\u\p-2u - \v\p-2v\ < g(u,v)\u - v\
u, v E R
g(u,v)
{
for u = v,
| u—v \ ' '
(p — 2)\u\p—l fo ru = v
sup a \\u\p 2u — \v\p 2v\\w\dx < sup a f (x, s)\u — v\\w\dx.
IMI=i
(2.2)
satisfies these conditions. By (2.2), for all real-valued functions u = u(x,s), v = v(x, s) E U
\\u\p"2u — \v\p-2v\ < /(x, s)\u — v\ (23) almost everywhere, where /(x,s) = g(u(x, s),v(x, s)). It follows from (2.3) that
(2.4)
a
116 BecTHHK lOypry. Cepna «MaTeMaTHHecKoe MOflenHposaHHe h nporpaMMHposaHHe »
In the left-hand side we have
sup a \\u\p 2u — \v\p 2v\\w\dx > sup a
MM
\H\ = 1
\u\p 2uwdx — \v\p 2vwdx
sup \ (M(u) — M(v),w) \ = \\M(u) — M(v)
\\w\\ = 1
In the right-hand side we have
sup a f (x, s)\u-v\\w\dx = sup a (f (x,s)\u — v\, \w\) < sup a\\f (x, s)\u—v\\l2 \\w\\l2 <
\N\=i
H\ = 1
\\w\\ = 1
Hence,
< aa \\f (x, s)\u — v\\\b2 < aC\\f (x, S)\\\\u — v\\ = F (S)\\u — v\\. \\M(u) — M(v\ < F(s)\u — v\\.
□
In the case p > 2 the existence of a generalized solution to problem (0.1) - (0.3) for arbitrary initial conditions is established in [7]. Thus, Lemma 1 implies
Theorem 4. Assume that p > 2 and take a,X E R+. Then for every u0 E U and every t E R+ there exists a unique generalized solution to problem (0.1)-(0.3).
References
1. Barenblatt G.I., Gilman A.A. [Mathematical Model of the Countercurrent Capillary Impregnation]. Inzhenerno-fizieheskiy zhurnal [Journal of Engineering Physics and Thermophysics], 1987, vol. 52, no. 3, pp. 456-461. (in Russian)
2. Zagrebina S.A., Sagadeeva M.A. [The Generalized Showalter-Sidorov Problem for the Sobolev Type Equations with strongly (L,p)-radial operator] Vestnik MaGU. Matematika [Bulletin of the Magnitogorsk State University], 2006, no 9, pp. 17-27. (in Russian)
3. Keller A.V. [The Leontief Type Systems: Classes of Problems with the Showalter-Sidorov Intial Condition and Numerical Solving] Izvestiya Irkuiskogo gosudarstvennogo universitita. Seriya: Matematika [The Bulletin of Irkutsk State University. Series: Mathematics], 2010, no 2, pp. 30-43. (in Russian)
4. Sviridyuk G.A., Zagrebina S.A. Verigin's Problem for Linear Equations of the Sobolev Type with Relatively P-Sectorial Operators. Differential Equations, 2002, vol. 38, no. 12, pp. 1745-1752. DOI: 10.1023 A: 1023812213901
5. Sviridyuk G.A., Zamyshlyaeva A.A. The Phase Spaces of a Class of Linear Higher-Order Sobolev Type Equations. Differential Equations, 2006, vol. 42, no. 2, pp. 269-278. DOI: 10.1134/S0012266106020145
6. Manakova N. A., Bogatyreva E.A. [Numerical Research of Processes in the Barenblatt - Gilman Model]. Vestnik MaGU. Matematika [Bulletin of the Magnitogorsk State University], 2013, issue 15, pp. 58-67. (in Russian)
*
7. Manakova N. A., Bogatyreva E.A. [On a Solution of the Dirichlet - Cauchy Problem for the Barenblatt - Gilman Equation]. Izvestiya Irkutskogo gosudarstvennogo universitita. Seriya: Matematika [The Bulletin of Irkutsk State University. Series: Mathematics], 2014, vol. 7, pp. 52-60. (in Russian)
8. Al'shin A.B., Korpusov M.O., Sveshnikov A.G. Blow-up in Nonlinear Sobolev-Type Equations. Berlin, N.-Y., Walter de Gruyter GmbH & Co. KG, 2011. DOI: 10.1515/9783110255294
9. Demidovich B.P. Lektsii po matematicheskoy teorii ustoychivosti [Lectures on the mathematical theory of stability]. Moscow, Nauka, 1967. (in Russian)
Received May 16, 2014
УДК 517.9 DOI: 10.14529/mmpl40409
О ЕДИНСТВЕННОСТИ НЕЛОКАЛЬНОГО РЕШЕНИЯ МОДЕЛИ БАРЕНБЛАТТА - ГИЛЬМАНА
Е.А. Богатырева, И.Н. Семенова
В работе рассматривается вопрос единственности обобщенного решения задачи Дирихле - Коши для уравнения Баренблатта - Гильмана. Это уравнение описывает неравновесную противоточную капиллярную пропитку. Неизвестная функция соответствует функции эффективной насыщенности. Основное уравнение модели является нелинейным и не разрешимо относительно производной по времени, что создает значительные трудности при его рассмотрении. Задача Дирихле - Коши для уравнения Баренблатта - Гильмана в подходящих функциональных пространствах редуцируется к задаче Коши для квазилинейного уравнения соболевского типа. Уравнения соболевского типа составляют обширную область неклассических уравнений математической физики. Методы исследования, которые используются в работе, первоначально возникли в теории полулинейных уравнений соболевского типа. Для задачи Коши получено достаточное условие существования единственного обобщенного решения. Показано существование единственного нелокального обобщенного решения задачи Дирихле -Коши для уравнения Баренблатта - Гильмана.
Ключевые слова: уравнение Баренблатта-Гилъмана; квазилинейное уравнение соболевского типа; обобщенное решение.
Литература
1. Баренблатт, Г.И. Математическая модель неравновесной протнвоточной капиллярной пропитки / Г.И. Баренблатт, A.A. Гильман // Инженерно-физический журнал. - 1987. - Т. 52, № 3. - С. 456-461.
2. Загребина, С.А. Обобщенная задача Шоуолтера - Сидорова для уравнений соболевского типа с сильно (L,р)-радиальным оператором / С.А. Загребина, М.А. Ca-гадббва // Вестник Магнитогорского государственного университета. - 2006. -до д. _ с. 17-27.
118 Вестник ЮУрГУ. Серия «Математическое моделирование и программирование»
3. Келлер, A.B. Системы леонтьевского типа: классы задач с начальным условием I Поуо. пера Сидорова и численные решения / A.B. Келлер // Известия Иркутского государственного университета. Серия: Математика. - 2010. - № 2. - С. 30-43.
4. Свиридюк, P.A. Задача Веригина для линбиных уравнении соболевского типа с относительно р-секториальными операторами / P.A. Свиридюк, С.А. Загребина // Дифференциальные уравнения. - 2002. - Т. 38, № 12. - С. 1646.
5. Свиридюк, P.A. Фазовые пространства одного класса линейных уравнений соболевского типа высокого порядка / P.A. Свиридюк, A.A. Замышляева // Дифференциальные уравнения. - 2006. - Т. 42, № 2. - С. 252-260.
6. Манакова, H.A. Численное исследование процессов в модели Баренблатта-Гильмана / H.A. Манакова, Е.А. Б огатырбва / / Вестник МаГУ. Математика. - 2013. - Вып. 15. - С. 58-67.
7. Манакова, H.A. О решении задачи Дирихле-Коши для уравнения Баренблатта-Гильмана / H.A. Манакова, Е.А. Б огатырбва // Известия Иркутского государ-ствбнного университета. Серия: Математика. - 2014. - Т. 7. - С. 52-60.
8. APshin, A.B. Blow-up in Nonlinear Sobolev-Type Equations / A.B. APshin, M.O. Korpusov, A.G. Sveshnikov. - Berlin, N.-Y.: Walter de Gruyter GmbH & Co. KG, 2011.
9. Демидович, Б.И. Лекции по математической теории устойчивости / Б.П. Деми-дович. - М: Наука, 1967.
Статья подготовлена в рамках выполнения работ по госзаданию МОиН РФ 200Ц-392, проект 1942.
Екатерина Александровна Б огатыревсЦ сюсистбнт^ кафедра «Уравнения математической физики», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Ирина Николаевна Семенова, кандидат педагогических наук, кафедра «Теория и методика обучения математике:», Уральский государственный педагогический уни-вврситбт (г. Екатеринбург, Российская Федерация), [email protected].
Поступила в редакцию 16 мая 2014 г.