Достаточные условия сходимости положительных решений линеаризованной двумерной задачи транспорта наносов Текст научной статьи по специальности «Математика»

УДК 517.95 10.23947/1992-5980-2017-17-1-5-17

Достаточные условия сходимости положительных решений линеаризованной двумерной

*

задачи транспорта наносов

А. И. Сухинов1, В. В. Сидорякина2, А. А. Сухинов3**

1Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 2Таганрогский институт имени А. П. Чехова (филиал) РГЭУ (РИНХ), г. Таганрог, Российская Федерация 3Южный федеральный университет, г. Ростов-на-Дону-Таганрог, Российская Федерация

Sufficient conditions for convergence of positive solutions to linearized two-dimensional sediment transport problem*** A. I. Sukhinov1, V. V. Sidoryakina2, A. A. Sukhinov3**

1Don State Technical University, Rostov-on-Don, Russian Federation 2 Taganrog Chekhov Institute, RSUE (RINH) branch, Taganrog, Russian Federation 3Southern Federal University, Taganrog - Rostov-on-Don, Russian Federation

Введение. Транспорт наносов является одним из основных процессов, определяющих величины и темпы деформаций донных поверхностей водных объектов. Чаще всего прогностические исследования в этой области строятся на основе математических моделей, которые позволяют сократить, а в ряде случаев исключить дорогостоящие и опасные в экологическом отношении эксперименты. Для прогнозирования изменения рельефа дна в основном используются пространственно-одномерные модели. Для реальных прибрежных систем со сложной формой берега вектор потока наносов в общем случае не ортогонален касательной к береговой линии в каждой из ее точек. Также он может не совпадать с вектором ветровых напряжений. Поэтому для решения многих практически важных задач, связанных с прогнозированием динамики донной поверхности водоемов, необходимо применение пространственно-двумерных моделей транспорта наносов и эффективных численных методов их реализации. Материалы и методы. Авторами (А. И. Сухинов, А. Е. Чистяков, Е. А. Проценко, В. В. Сидорякина) ранее была предложена пространственно-двумерная модель транспорта наносов, удовлетворяющая основным законам сохранения (материального баланса и импульса), которая представляет собой квазилинейное уравнение параболического типа. Были построены и исследованы линейные разностные схемы и решены модельные, а также практические задачи. Однако осталось в тени теоретическое исследование «близости» решений исходной нелинейной начально-краевой и линеаризованной непрерывной задач, на основе которой была построена дискретная модель (разностная схема). Особый интерес представляет исследование корректности линеаризованной задачи и определение достаточных условий положительности решений, т. к. только положительные решения задачи транспорта наносов имеют смысл в рамках рассматриваемых моделей.

Introduction. The sediment transportation is one of the major processes that define the magnitude and back surface changing rate for water bodies. The most used prognostic studies in this field are based on the mathematical models that allow reducing, and in some cases, eliminating the expensive and often environmentally burdensome experiments. Spatially one-dimensional models are mainly used to predict changes in bottom topography. For actual coastal systems with irregular coastal configuration, the flow vector is generally not orthogonal to the tangent line for the coastline at each of its points. It also may not coincide with the wind stress vector. Therefore, to solve lots of practically important problems associated with the prediction of the dynamics of the back surface of water basins, it is necessary to use spatially two-dimensional models of sediment transportation and effective numerical methods of their implementation.

Materials and Methods. The authors (A.I. Sukhinov, A.E. Chisty-akov, E.A. Protsenko, and V.V. Sidoryakina) have earlier proposed a spatially two-dimensional model of sediment transport that satisfies the basic conservation laws (of material balance and momentum) which is a quasilinear equation of parabolic type. The linear difference schemes are constructed and studied; model and some practically important tasks are solved. However, the theoretical study on the proximity of solutions for the original nonlinear initial-boundary value problem and the linearized continuous task - on which basis a discrete model (difference scheme) was built -remained in the shadow. The study of the linearized problem correctness and the determination of sufficient conditions for positivi-ty of solutions are of special interest because only positive solutions to this sediment transport problem have physical value within the framework of the considered models.

c3 «

S X a

<u

'Работа выполнена при финансовой поддержке РФФИ по проектам № 15-01-08619, 16-07-00100, 15-07-08626, 15-07-08408 и по проекту № 00-16-13 в рамках Программы фундаментальных исследований Президиума РАН № 1.33П. **E-mail: sukhinov@gmail.com,cvv9@mail.ru, andreysoukhinov@gmail.com

*** The research is done with the financial support from RFFI (projects nos. 15-01-08619, 16-07-00100, 15-07-08626, 15-07-08408) and on project no. 00-1613 within the frame of the RAS Presidium Program of Fundamental Research no. I.33P.

Результаты исследования. Исследуемая нелинейная двумерная модель транспорта наносов в прибрежной зоне мелководных водоемов учитывает следующие физически значимые факторы и параметры: пористость грунта; критическое значение касательного напряжения, при котором начинается перемещение наносов; турбулентный обмен; динамически изменяемая геометрия дна; ветровые течения и трение о дно. Линеаризация осуществляется на временной сетке — нелинейные коэффициенты параболического уравнения берутся с запаздыванием на один шаг временной сетки. Далее строится цепочка взаимосвязанных по начальным условиям — финальным решениям цепочки линеаризованных смешанных задач Коши на равномерной временной сетке, и таким образом осуществляется линеаризация в целом 2Б нелинейной модели. Ранее авторами были доказаны существование и единственность решения цепочки линеаризованных задач, получена априорная оценка близости решения цепочки линеаризованных задач к решению исходной нелинейной задачи. В данной работе определены условия положительности ее решений и их сходимости к решению нелинейной задачи транспорта наносов в норме Гильбертова пространства Ь1 со скоростью О(т), где т — временной шаг.

Обсуждение и заключения. Полученные результаты исследования пространственно-двумерной нелинейной модели транспорта наносов могут быть использованы при прогнозировании нелинейных гидродинамических процессов, повышения их точности и надежности в силу наличия новых функциональных возможностей учета физически важных факторов, в том числе уточнения граничных условий.

Research Results. The investigated nonlinear two-dimensional model of sediment transport in the coastal zone of shallow reservoirs takes account of the following physically significant factors and parameters: soil porosity; critical value of the tangent stress at which load transport starts; turbulent interaction; dynamically varying of the bottom geometry; wind currents; and bottom friction. Linearization is carried out on the time grid — nonlinear coefficients of the parabolic equation are taken for the previous time grid step. Next, a chain of tasks connected by the initial data - final solutions of the linearized mixed Cauchy problems chain on a uniform time grid is constructed, and thus, the linearization for the initial 2D nonlinear model is carried out in large. Earlier, the authors have proved the existence and uniqueness of the solution to a linear tasks chain. Prior estimate of the proximity of the linearized problem chain solution to the initial non-linear task solution has been also obtained. The conditions of its solution positivity and their convergence to the nonlinear sediment transport problem are defined in the norm of the Hilbert space Li with the rate O(t) where t is a time step.

Discussion and Conclusions. The obtained research results of the spatially two-dimensional nonlinear sediment transport model can be used for predicting the nonlinear hydrodynamic processes, improving their accuracy and reliability due to the availability of new accounting functionality of physically important factors, including the specification of the boundary conditions.

Ключевые слова: пространственно-двумерная модель транспорта наносов, прибрежная зона, нелинейная задача, линеаризованная задача, положительность решения.

Keywords: spatially 2D model of sediment transportation, coastal zone, nonlinear problem, the linearized problem, solution positivi-

ty.

Введение. В последние десятилетия количество работ, посвященных рассмотрению нелинейных краевых задач, заметно возросло. Этому факту способствовало, в частности, появление новых приложений в физике океана и атмосферы [1-5].

В данной работе рассматривается нелинейная начально-краевая задача для пространственно-двумерной модели транспорта наносов мелководных водоемов и проводится анализ некоторых качественных свойств ее решений.

Исследование указанной задачи осуществляется путем ее линеаризации. С этой целью строится временная сетка, на которой значения нелинейных коэффициентов задаются в левом узле соответствующей временной ячейки. Далее осуществляется решение цепочки связанных по начальным — конечным данным линейных уравнений параболического типа.

Ранее коллектив авторов проводил исследования линеаризованной начально-краевой задачи, которые включали доказательство существования и единственности ее решения с одновременным определением требований к классам гладкости входных данных задачи [6, 7]. Также была получена априорная оценка близости решения цепочки линеаризованных задач к решению исходной нелинейной задачи в норме Гильбертова пространства Ь1 в зависимости от интегральных оценок правой части, граничных условий и нормы начального условия [8]. В предыдущих работах авторов для линеаризованной задачи транспорта наносов была построена и исследована консервативная устойчивая разностная схема, которая численно реализована для модельных и реальных задач прибрежной зоны [9-11]. В данной работе определены достаточные условия положительности решений линеаризованной задачи и их сходимости к ре> шению исходной нелинейной начально-краевой задачи транспорта наносов в норме Гильбертова пространства Ь1 при ^ стремлении к нулю параметра — шага временной сетки, на которой осуществлялась линеаризация.

й

л

1. Постановка нелинейной начально-краевой задачи для двумерной модели транспорта наносов. Следуя [4], рассмотрим уравнение транспорта наносов:

3H

sin ф0

(1 "s)'— = div k-—ф--gradH -div(k-xb). (1)

6

Здесь H = H (x, y, t) — глубина водоема; е — пористость донных материалов; xb — вектор касательного тангенциального напряжения на дне водоема; xbc — критическое значение тангенциального напряжения; xbc = a sin ф0, ф0 — угол естественного откоса грунта в водоеме; k = k (H, x, y, t) — нелинейный коэффициент, определяемый соотношением:

М

Aro d k =-

\P

4 -—Г" • gradH sin ф0

(2)

(p( -P0 )gd)

где p1, p0 — плотности частиц донного материала и водной среды соответственно; g — ускорение силы тяжести; со — частота волны; A и Р — безразмерные постоянные; d — характерные размеры частиц грунта.

Пусть D с Rn — область, где происходит процесс, и S — ее граница, которая является кусочно-гладкой линией. Областью задания уравнения (1) считаем цилиндр ЦТ = D х(0, T) высоты T с основанием D . Его граница состоит из боковой поверхности S х[0, T] и двух оснований: нижнего D x{0j и верхнего D x{Tj. Далее для простоты

уравнение (1) рассматривается в прямоугольной области D (x, y) = {0 < x < Lx ,0 < y < Ly j .

Дополним уравнение (1) начальным условием, предполагая, что функция начальных условий принадлежит соответствующему классу гладкости:

H (x, y,0) = H0 (x, y), H0 (x, y) e С2 (D)n С (D), gmd{x y }H0 e С (D), (x, y) e D . (3)

Сформулируем условия на границе области D :

И = 0, (4)

I I y=0

H (Lx, y, t) = H2 (y, t), 0 < y < Ly, (5)

H (0, y, t) = H (y, t), 0 < y < Ly, (6)

H (x, 0, t) = H3 (x), 0 < x < Lx, (7)

H (x,L'y, t) = H4 (x, t), 0 < x < Lx, L'y < Ly. (8) Дополнительно к граничным условиям (5)-(8) предполагаем выполнение условий их гладкости — существование непрерывных производных на границе области D :

grad{x y)H e С (Цт )п С1 (Ц ). (9)

Считаем, что всегда есть слой жидкости конечной толщины в рассматриваемой области и для указанного временного промежутка не происходит осушения области, т. е.

H (x, y, t)> c0 = const > 0, 0 < x < Lx, 0 < y < L'y, 0 < t < T. (10)

Условие невырожденности оператора задачи имеет вид:

k > k0 = const > 0, V(x,y) e D, 0 < t < T. (11)

Вектор тангенциального напряжения на дне выражается с использованием единичных ортов системы координат естественным образом:

tb = itbx + j4y , tbx=tbx (x,y,t) , tby=tby (x,y,t). (12)

2. Линеаризация нелинейной начально-краевой задачи транспорта наносов. Для получения линеаризованной модели на временном отрезке 0 < t < T построим равномерную сетку ют с шагом t, т. е. множество точек

ют = {tn = n-t, n = 0,1,...,N, N-t = Tj.

Осуществим линеаризацию указанной начально-краевой задачи, используя равномерную временную сетку

f ^ \ Й

= {tn = n-x, n = 0,1,..., N, N •x = T}. Линеаризацию члена div k ——— • gradH

sin ф0

s

(коэффициента k ) выполним пу- К

ей

тем выбора его значений в момент времени t = tn, п = 0,1,..., N и рассмотрения уравнения (1) на временном промежутке ^ < t < tп+1, п = 0,1,2,..., N -1.

го

x

тз М й

При этом предполагается, что известны функция H (п'(х, у, tn_1 ) = H (" ^(х, у, tи_1) и ее частные производные по пространственным переменным.

Если п = 1 , то в качестве H (1'(х, у, t0) достаточно взять функцию начального условия, т. е.

H(1' (x,у,t0) H0 (х,у). Если же п = 2,...,N , то функция H(п' (x,у,tn_1) = H(п_1' (x,у,tn_1) предполагается известной, поскольку предполагается решенной задача (1)-(14) для предыдущего временного промежутка tn_2 < t < tn_1 . Введем обозначения:

Р_1

, п = 1,2,...,N. (13)

k

и_1) = A& d

((Pi -Ро )d)

lb • gradH(n-1) (x,y,tn-1)

sin ф0

Тогда уравнение (1) после линеаризации запишем в виде:

~,тт{п) / Л _

(1 = div k.gradH(n) -div (k(n-1Ub I, tn 1 < t < tn, n = 1,..., N (14)

dt sin ó„ V '

• a

и дополним его начальными условиями:

H(x,у,t0 ) = Hо (x,у),H(n) (x,у,/„_!) = H(n_1) (x,у,tn_i), (x,у) 6 D, n = 2,..., N. (15)

Член вида div(k(n 1)-хг,) является при такой линеаризации известной функцией правой части. Граничные условия (4)-(8) предполагаются выполненными для всех промежутков времени tn_1 < t < tn, n = 1,2,..., N.

Отметим, что коэффициенты k(n 1, n = 1,2,...,N зависят от пространственных переменных x,у и временной переменной tn_, n = 1,2,...,N , определяемой выбором шага х сетки юх, т. е. k(n_l) = k(n_1) (x,у,tn_l), n = 1,2,...,N.

3. Положительность решений линеаризованной начально-краевой задачи транспорта наносов. Существование и единственность решения линеаризованной задачи (14) с начальными условиями (15) и с условиями на

границе (4)-(8) обосновано в работах [6-8]. Обозначим решение линеаризованной задачи через H (x, у, t), (x, у) 6 D . Решение H(x,у,t), (x,у) 6 D зависит от шага х сетки ют = jtn = nx, n = 0,1,...,N, N-х = Г} . По существу имеем семейство решений jH(n) (x,у,t)|, n = 1,2,...,N, зависящее от параметра х . Убедимся в положительности решений указанной линеаризованной задачи.

Под положительным решением рассматриваемой начально-краевой задачи понимается функция

H(n)(x,у,t), tn_ < t < tn, n = 1,2,...,N классаC2 (цт )nC(цт), grad(xy6 C(цт ), удовлетворяющая уравнению

(14) с начальными и граничными условиями (4)-(8), (15), для которой H(n)(x,у,t)> с0 > 0, с0 = const, n = 1,2,..., N при любом n = 1,2,...,N.

Относительно рассматриваемой задачи будем предполагать выполнение условия

div (k (n_1)-Xb )< 0 (16)

для слагаемого левой части уравнения (14), а также положительность начального условия (15):

H0 (x, у), > c0 = const > 0, (x, у) 6 D (17)

g и граничных условий (5)-(8):

3 Hl (у, t), H2 (у, t), H3 (у, t), H4 (у, t)> C0 - const > 0. (18)

Д

о Рассмотрим вначале случай n = 1, для которого существует единственное решение уравнения (14). Обозна-

чим данное решение H(x, у, t). Указанное решение принадлежит классу C2 (Ц) n C (ц^), grad(x у)H6 C (Цt1).

¡> Предположим, что функция H(x, у, t) достигает своего минимального значения c0 или значения меньшего

ji его в некоторой внутренней точке параллелепипеда р = D х [0, t1 ] или в точке его «верхней крышки». Тогда это противоречит принципу максимума (минимума): решение начально-краевой задачи не может достигать экстремума, в том числе минимума во внутренних точках области P1 = D х[0, t1 ] и на ее «верхней крышке», где во всех точках, кроме

относящихся к боковой поверхности £ х[0,] параллелепипеда Р1 = Б х[0,г1 ] (рис. 1), ищется решение начально-краевой задачи для уравнения (14) при п = 1 [12].

T

ti-i

*2 h

О

и-

/ У

D

Рис. 1. Расчетная область Fig.1. Computational domain

Следовательно, функция Й (х, y, t) достигает экстремума, в том числе минимума, на нижнем основании Dх{0} параллелепипеда P1 = Dх[0,t1 ] или на его боковых границах Sх [0,t1 ], т. е.:

H (1)(х, y, t )> min В силу условий (17), (18), имеем

min Й (1)(х, y, t), min H (1)(х, y, t)

(х, У )eD,t=0 (х, y )e5,0<t<t1

К1)!

min H(1) (х, y, t) = Со > 0, min H(1)(х, y, t) = Co > 0,

(х, y )eD,t=0 (х, y )eS,0<t<t1

(19)

и потому

Й(1) (х, y,t) > c0 = const > 0.

К1)

(20)

Таким образом, финальное решение Н(' (х, у, г) для 0 < г < г1 будет положительным и большим

При п = 2 имеем начально-краевую задачу с начальным условием Н(2)(х, у, г1 ) = Н(1)(х, у, г1), имеющим гладкость такую же, как и начальное условие для уравнения (14) номера п = 1. Рассуждая аналогично, решение уравнения (14) г1 < г < г2 существует и единственно. Обозначим данное решение Н(2)(х, у, г) . Известно, что функция

Н(2) (х,у,г) принадлежит классу С2 (Ц^ )п С(ц^ ), )Н(2) е С(Ц^ ).

Очевидно, что решение Н(2)(х, у, г) второй задачи цепочки начально-краевых задач вида (4)-(8), (14), (15) достигает своего минимального значения с0 или значения меньшего с0 в некоторой внутренней точке параллелепипеда Р2 = Б х[г1, г2 ] или в точке его «верхней крышки», что противоречит принципу максимума (минимума). Тогда решение Н(2) (х,у,г) начально-краевой задачи достигает экстремума на нижнем основании указанного параллелепипеда или его боковых границах £ х г2 ] , и потому выполняется неравенство:

i= Й 0)<

Й(2) (х, y,t) > c0 = const > 0.

(21)

Отсюда следует, что решение Н(1) (х, у, г) для 0 < г < г1 будет положительным и большим величины

Нетрудно видеть, что вообще при любом п = к, к = 1,2,...,N получаем смешанную задачу для линейного уравнения параболического типа, имеющего достаточную гладкость начальных и граничных условий. Для любого

eö И К X а X к

t

c

0

c

0

номера k, к = 1,2,..., N функция решения Й (k '(х, y, t), tk-1 < t < tk, к = 1,2,..., N класса C2 (tk )n C (ц^) ,

grad(xy)Й(k) e C (цtt) начально-краевой задачи (4)-(8), (14), (15) в параллелепипеде Pk = D x[tk_1, tk ] достигает экстремума на нижнем основании указанного параллелепипеда или на его боковых границах, где в силу условия (17)-(18) будет выполняться

Й(k'(х,y,t)> c0 = const > 0, k, k = 1,2,...,N. (22)

Учитывая полученные соотношения (20)-(22), приходим к теореме. Теорема 1. Пусть даны уравнения (14):

(1 -е)

в ЦТ = Dx (0,T), D с Rn .

дЙ

(n)

dt

• = div

bc

Здесь k (n-1U-

Aco d

Й (1)(x, y, t0 )= Й0 (x, y), Й (n)(x, y, tn-1 ) = Й(n :)(x, y, tn-1), (x, y)e D, n = 2,..., N и c граничными условиями (4)-(8),

((P1 _P0 )gd)

r(

sin ф0

sin ф0

gmrn(n 1)(x, y, tn _j)

• gradЙ (n)j-div (k(n_1)-Tb ), tn_j < t < tn, n = 1,2,...

N,

P-1

div

(k (n-1)Xb )

< 0 с начальными условиями

n = 1,2,..., N

удовлетворяющими условию (17)-(18).

Тогда, если k(n-1) > k0 > 0, k(n-1) e C1 (d) , то для любого n = 1,2,...,N решение

Й(n)(x,y,t), tn-1 < t < tn, n = 1,2,...,N уравнения номера n в цилиндре ЦТ = Dx(0,T) является положительным и выполняется оценка

Й(n)(x,y,t)> с0 = const > 0, tn-1 < t < tn, n = 1,2,...,N.

(23)

4. Сходимость решений линеаризованной задачи к решению нелинейной начально-краевой задачи.

Убедимся в том, что решение линеаризованной задачи стремится к решению нелинейной задачи в норме Гильбертова пространства Ь1 (Б) при т^0, N-т = Т . С этой целью обозначим решение нелинейной задачи через

Н (х, у, t), (х, у)е Б , а решение линеаризованной задачи — Н (х,у,t), (х, у)е Б . Будем считать, что функция

Н (х, у, t) ограничена на промежутке 0 < t < Т , а также существуют и ограничены смешанные частные производные

до третьего порядка включительно по пространственной и временной переменной t от функций Н (х, у, t), Н (х, у, t)

для всех (х, у) е Б, 0 < t < Т :

dx (grad ( dt ]), dy (grad ( dt ]), dx

grad\dL

д_

dy

( Jm)) grad\ -df

Заметим, что функция Н(х,у,t) ограничена на каждом временном промежутке tn< t < tn, п = 1,2,..., N, поскольку является решением линейной начально-краевой задачи для уравнения параболического типа с «гладкими» входными данными.

Обратимся к нелинейной задаче (1)-(8). При подстановке функции решения Н (х, у, t) в уравнение (1), получим:

(1 -е)-дЙ = div v ' dt

k -4c

sin ф0

•gradЙ

- div (k • ть ),

(24)

й о тз

"¡3

и (U

o, £

где k =-

Aco d

(((51 -P0 )gd)

Tb --

4c

sin ф0

•gradЙ

P-1

Умножим обе части уравнения (24) на функцию Н (х, у, t), (х, у)е Б и проинтегрируем по переменной ^ 0 < t < Т и по переменным (х, у) в области Б . Получим

= J

о

t

J

о

íbc

J| (1 -e))JH—dxdy Idt =

d dt

JJ H •div | k--—— • gradH I dxdy

d \ sin фо

dt-J| JJ H •div (k-Tb )dxdy dt.

о VD V )

(25)

Выполняя ряд преобразований в равенстве (25), получим:

2(1 - s) А (я2 (x, y, T) - H2 (x, y, 0)) =

T( f - T f . _ > Л

• gradH dxdy dt - M JJ H • div я • -ь ) dxdy Idt.

= J

JJ H •div

k -

%c

sin ф0

(26)

Далее рассмотрим линеаризованную задачу. При подстановке Н (х, у, t) в уравнение (14) получим:

где к

(1 -в)

A(Ьd

dH(n) J (

- = div

dt

n-1) Tbc

БШ ф0

gradH(n) I - div (k(n-1) • ть ), tn-1 < t < tn

n = 1,..., N,

(27)

((Р -Ро )gd)

Tb —

Lbc

бШ ф(

■• gradH

P-1

Умножим обе части уравнения (27) на функцию Н(п) = Н(п) (х, у, t), (х,у) е Б и проинтегрируем обе части этого уравнения по области Б , а после этого — по переменной t, tn_l < t < tn. Получим:

J

(1 -B))J h (

dH

(

= J

tn—1

(

JJ H(n) •div

D

k n• -^ •gradH(

dt

dxdy

) )

dxdy

dt =

dt - J |JJ H(n) -div (k(

tn-1 V D V

-1

sin ф0

Просуммировав обе части соотношения (28) по n, n = 1,...,Ж , находим:

А „ , „Лín) Л

tb

)drdy

(28)

dt.

JJ(-в)

D

E J H

n=1 t ,

ln-1

N

E J jj h(n)• div

n=1 _tn-1 VD V

^•J^ •gradH(

3

sin ф0

dt

Л I

dt

dxdy =

dxdy

) )

'(((n)

dt - J | JJ Hní •div \tn-4 •Tb )dxdy | dt

tn-1 V D

(29)

Преобразуя левую часть равенства (29), запишем:

-2(1-е)Я(я2 (х,у,Т)-Н2 (х,y,0))dxdy ■■

N tn

E J jj h(n)• div

n=1 tn-1 lD V

k"1 •gradH(

sin ф0

Л Л dxdy

) )

dt - J ^JJ H(n) •div (k(n-1)-Tb )xdy ) dt

(30)

Вычитая из равенства (30) выражение (26) и учитывая, что Н (х, у,0) = Н (х, у,0) , получаем:

2(1 -8)Я(Я2 (х,у,Т)-Н2 (х,у,Т)) =

= J JJ H •div

' T Л i

k--• gradH dxdy

sin Ф0 )

dt - J |jjH • div (k • Tb ) dxdy |dt -

(31)

N

-E

n=1

t |

J JJ H (n)•div|k (

Б1П ф0

•gradH(n) Idxdy

Л t„ /

dt - j" IJJ H(n '•div (k (n-:1 )dxdy dt

V-1 V D

Выполним преобразования правой части равенства (31). С этой целью добавим, а затем вычтем выражения

(

Hn> • div

k •-

sin ф0

• gradH

Н(п) • div(к • Хб ) под знаком интеграла. Комбинируя слагаемые, получим:

ев «

S X ев

X

К

1 (1 -в)Я(#2 (x,y,T)-H2 (x,y,T))xdy --

N

- z

n=1

tf(( - H (n))div

(

k -4c

sin ф0

• grad H

+tf H(n) • div I - • (k • grad H - k(n-1) • gradH(n)) D У sin ф0 V )

Л Л dxdy +

J J

Л

dxdy +

(32)

+я(я("У -HУ div(k• ть )dxdy + JJH(nУ •div((k(n-1)-kjdxdyjdt Оценим каждый из интегралов, стоящих в правой части равенства (32) под знаком суммы. Введем обозначе-

ния:

/in = 1

n ( , (

k •- bc

л (я - H(n)) div

• gradH

Л Л dxdy J J

in - í

¿и-1

л H(n)• div

D

lbc

sin ф0

(k • gradH - k(n-1)- gradH(n))

dt.

1n -14 -

sin ф0

I (я(H(n) -H)• div(k-Tb)dxdyJjdt

J (flH(n) .div((k(n-1) -k)• Tb)dxdyЛdt

Л | dxdy J J

dt,

1 n -13 -

tn-1 y D

n-1,..., N.

(33)

(34)

(35)

(36)

Полагая п = 1 в равенствах (33)-(36), рассмотрим интегралы по первому временному интервалу /0 < t < ^ . Для начала выполним оценку интеграла //. Используя разложения функций Н (х, у, t) и Н (1)(х, у, t) по фор муле Тейлора относительно ^ -10) с остаточным членом в форме Лагранжа, с учетом условия (3), можем записать:

'ЭН (х,у, ^) дН(1) (х, у, %2)

H (x, y, t)-H (1)(x, y, t )-

dt

dt

•(t - to),

(37)

где to <5i <ti, to <^2 <A, to <t<ti.

Используя равенство (37), получаем выражение для интеграла //:

/} - j dt H(( - H(1)) div

k • -Tbc

( (

-1 (-to )

dH (x, y, ^) dH (1)(x, y, ^ )) (

У D У

dt

dt

sin ф0

div

gradH

dxdy

J J

(38)

k -4c

sin ф0

• gradH

dxdy

dt.

С учетом равенства t -10 - t из (38) можно перейти к оценке:

i1 < L T2 • v v м;,

(39)

где

й о тз

"¡3

и

<и

Ü £ л

М; = max ^ max_

to <t<t; I (x,y)eD

dH (x, y, ) dH (1)(x, y,

dt

dt

• div

k -4c

sin ф0

• gradH

Далее получим оценку для интеграла /^. Заметим, что справедливо равенство:

div

k •-

Lbc

•gradH

- div

- div

/ \ | (

+ div

k(°) •^^•gradH(1) sin ф0

sin ф0

г

(k - k gradH

O ' sin фо

Нетрудно видеть, что интеграл 11 с учетом выражения (40) можно представить в виде:

k(0) Tbc

sin ф0

• (gradH - gradH(1))

¡2 = 1 я И (1) • ¿¡у

Ьс

sin ф(

-• ((к - к(0) ) • gradH + &(0) • {^тёИ - gradH(1))) (Ыу) ¿й\

(41)

Разложим функции к (х, у, t) , И (х, у, t) по формуле Тейлора и, воспользовавшись равенствами

к (х, у, ^ ) = £(0) (х, у, ^ ) , И (х, у, ^ ) = ИИ(1) (х, у, t0 ) , находим:

к - к (0)=^1у!3) - to ),

д/

grad И - grad И(1) = grad

дИ (х, у, )

д/

(t - to )

(42)

(43)

где to <§э <4, to <^4 <А, ^ <t <tl.

Используя выражения (42) и (43) для интеграла вида (41), получаем:

12 = 1(-ч)

11И(1) -¿¡у

в

Ьс

sin ф0

дк (x, y, ^3)

дt

gradH + к((0' • grad

дИ (х, у,

ЛЛ

дt

dxdy

уу

(44)

Предполагая ограниченность соответствующих производных, введем постоянную величину:

М 2 = тах \ тах

to^<¿1 I (х,у)еВ

И (1)Лу

А

ть ( дк (х, у, Е3) ~(0 \ (дИ (х, у, Е4 )

Ьс ^ 1 ^ 37 • gradH + кИ(0gradl 1 У

sin ф

дt

дt

ЛЛ

dxdy

Из равенства (44), используя обозначение для величины М 2, получим

4 < т2 • м1 • ьх • Ьу. Для оценки интеграла ¡\ воспользуемся равенством (37), с учетом которого имеем

(45)

¡1 = 1 (-'0 )

Г Г дИ(!)

в

V V

дИ(1)(х, у, ^ ) дИ (х, у, ^ )

дг

дг

• ¿¡у (к • тЬ ) dxdy

dt,

где ?0 <^5 <?1, ^ <^6 ?0 <1 <1. Вводя обозначение

М1 = тах <; тах

to^^ I (х,у)еВ

дИ(1)(х, у, ^ ) дИ (х, у, )

д(

д(

¿¡у (к • хЬ |

можем записать:

¡1 < ^ т2 • Ьх- Ьу- м1.

(46)

Далее перейдем к оценке интеграла ¡4 . С этой целью выполним несложные преобразования ¡4 с учетом вы-

ражения вида (42). Получим равенство

¡4 = 1 ( -0 )•

11И(1) • ¿¡У-

(

VV

дк (х, у, ^7 )

дt

Л _ Л ТЬ

у у

dxdy

dt.

(47)

В силу ограниченности функции И(1)(х, у, t), /0 < / < ^ и смешанных производных функции И (х, у, t) по пространственной и временной переменной из (47) получаем неравенство:

¡1 < 2 • т2 • М4 • Ьх • Ьу,

(48)

где

М4 = тах <! тах <

to ^ I (х,у )еВ

И(1) •¿¡у

Гг

VV

дк ( y, )

д/

Л _ л

ТЬ

У у

ев И К

к

ев

X

К

Собирая полученные оценки из (39), (45), (46), (48) и пользуясь неравенством многоугольника для модулей величин, получаем неравенство вида

Я(Я2 (х, у, ) - Н(1)2 (х, у, ^у < т2 -Ьх-Ьу- М;

(49)

где

м ; = —-((, + м 1 + м, + м4). 1 -8 ^ '

Меняя местами функции Н2 (х,у,), Н(1)2 (х,у,) и проводя рассуждения, аналогичные предыдущему, начиная с равенства (24), можно получить оценку:

Ц (((1)2 (х, у, /1) - Н2 (х, у,/1 ))у < т2 - Ьх - Ьу -м;. (50)

Из неравенств (49) и (50) следует неравенство

я Н(1)2 (х,у,/1)-Н2 (х,у,/1 )с1хёу <т2 -Ьх -Ьу -м;.

5

Преобразуем левую часть неравенства (51):

я я(1)2 (х,у,/1)-Н2 (х,у,/1 )dxdy =

= Я Я(1)(х,у,/1)-Н(х,у,/1) - Я(1)(х,у,/1) + Н(х,у,/1 )dxdy.

и

(51)

(52)

Далее будем предполагать, что выполняются условия теоремы 1 и неравенство вида (23) выполняется для функций Н (х, у, /) и Н(п) (х,у,/), п = 1,2,..., N, т. е.

Н (х, у, /) > с0 > 0, (х, у) е 5, 0 < / < Т ,

Н(п) (х,у,/) > с0 > 0, (х,у) е 5, 0 < / < Т . С учетом выражений (53), (54) находим

Я Н(1)2 (х, у, /1) - Н2 (х, у,/1) ёхёу > 2с0 - Ц Н) (х,у,/1) - Н (х,у,/1)ёхёу. и и

Из соотношений (51), (52) и (55) следует оценка

я Н(1)(х,у,/1)-Н(х,у,/1 )скёу <-1--т2-Ьх-Ьу -м;.

5 2с0

Требуемая оценка для п = 1 получена, т. к. неравенство (56) эквивалентно

К1)!

Н (1)(х, у, /1)-Н (х, у, /1)

( 2--т •Ьх-Ьу -м;.

1,(5) 2С0

Фактически сделан первый шаг индукции при п = 1; к этому относятся выкладки (37)-(57). Пусть п = 5, 5 = 2,..,N . Тогда, рассуждая аналогичным предыдущему образом, получаем оценку:

Н (Т )-Н <Т >1М» )< - '2- ^ Ьу С' N• (<1 -') + (1 -т2 )+.+(1 )).

Здесь , — некоторые постоянные функции, определяемые для временного промежутка соответственно.

Т

Поскольку N = —, то из неравенства (58) получаем: т

||я (Т )-н (Т )|| ь

- Ьх- Ьу-С - Т - ((1 -т) + (1 -т2 ) + ... + (1 -т5+1)).

й о тз

"¡3

и

и

О, £ -Й

|1ь1 (») 2с0 - (1 -т) При т —у 0 из (59) находим

\\Н (Т)-Н (Т ) А(»)— 0, что и завершает доказательство сходимости линеаризованной задачи к решению исходной задачи. Теорема 2. Пусть дана начально-краевая задача:

(1 - 8)—Н- = йп>

к-gradH (кть ) в ЦТ = 5 х( 0, Т), 5 с Яп . I 8Ш Ф0 ) V >

(53)

(54)

(55)

(56)

(57)

(58)

(59)

5

1

Здесь k = -

A со d

((Pi -Po )gd)

ть --

sin ф0

• gradH

P-i

div I к

(k ("-1)Tb)

< 0 с начальными условиями:

H (х,у,0) = Hо (х,у), Hо (х,у) е С2 (Р)п С (D) , grad(x у}Но е C(D), (х,у) е D , условиями (4)-(8) на границе области Р . Функция Н (х,у,t) — решение задачи класса С2 (цт )пС (цт ) ,

grad{x,у)Н е С (цт ) .

Пусть далее на временной сетке ют = {/и = п -т, п = 0,1,...,N, N -т = Т} построена линеаризованная начально-краевая задача:

(1 )

дИ

(n)

dt

■ = div

(n-i)

Lbc

Здесь k

n-1)

A со d

((Pi -Po )gd)

sin ф0

т>

• gradH(

- div

(к (n-14),

n = 1,2,..., N.

Tb --

bc

И(1) (x,y,to) = Иo (x,y), И(n> (x,y,tn-i) = И(n~1) (x,y,tn-i) , (x,y) e D, n = 2

sin ф0

J ( n

grnM(n )x, y, tn-i)

P-i

с начальными условиями:

, N.

граничными условиями (4)-(8), удовлетворяющими условию (17)-(18). Функция Н (х, у, /) — положительное решение задачи класса С2 (цт )п С (цт), grad( x у)Н е С (ЦТ ) . Также предположим, что существуют и непрерывны смешанные частные производные до третьего порядка включительно по пространственной и временной переменной t от функций Н (х,у, t) , Н (х,у, {) для (х,у) е Р, 0 < t < Т :

dx (grad ( dt ]], ddy [grad ( dt )), dx

fgrad (?]], dy [grad (fT

Тогда, если к(n 1) > к0 > 0, к0 = const, к(n 1)

||И (T)-И (T)||

е С1 (D), n = 1,2,..., N , то имеет место оценка 1

/ • Lv •С •Т

lLi(D) 2c0 •(i -т) x y

(60)

где С — постоянная величина, большая 0, т. е. положительное решение линеаризованной задачи И (x, y, t) стремится к решению И (x,y,t) нелинейной задачи в норме Гильбертова пространства L1 (D) при т ^ 0, N•т = T .

Заключение. Исследована линеаризованная двумерная математическая модель транспорта донных наносов в прибрежных системах, удовлетворяющая основным законам сохранения и учитывающая следующие физические процессы и параметры: пористость донного материала; критическое значение тангенциального напряжения, при котором начинается транспорт донных материалов; динамически изменяемая геометрия донной поверхности за счет движения водной среды.

Получены условия положительности решений линеаризованной начально-краевой задачи и их сходимости к решению исходной нелинейной начально-краевой задачи в норме Гильбертова пространства L1, сформулированные в виде теорем.

Библиографический список

1. Марчук, Г. И. Математические модели в геофизической гидродинамике и численные методы их реализации / Г. И. Марчук, В. П. Дымников, В. Б. Залесный. — Ленинград : Гидрометеоиздат, 1987. — 296 с.

2. Sukhinov, A. I. Numerical realization of the three-dimensional model of hydrodynamics for shallow water basins on a high-performance system / A. I. Sukhinov, A. E. Chistyakov, E. V. Alekseenko // Mathematical Models and Computer Simulations. — 2011. — Vol. 3, is. 5. — P. 562-574.

3. Леонтьев, И. О. Прибрежная динамика: волны, течения потоки наносов / И. О. Леонтьев. — Москва : ГЕОС, 2001. — 272 с.

4. Predictive modeling in sediment transportation across multiple spatial scales in the Jialing River Basin of China / Xiaoying Liu [et al.] // International Journal of Sediment Research. — 2015. — Vol. 30, is. 3. — P. 250-255.

5. Sukhinov, A. I. Mathematical modeling of sediment transport in the coastal zone of shallow reservoirs / A. I. Sukhinov, A. E. Chistyakov, E. A. Protsenko // Mathematical Models and Computer Simulations. — 2014. — Vol. 6, is. 4. — P. 351-363.

ce

и к

X а X <u

6. Сухинов, А. И. Существование и единственность решения линеаризованной двумерной задачи транспорта наносов / А. И. Сухинов, В. В. Сидорякина // Теория операторов, комплексный анализ и математическое моделирование : тезисы докладов XIII междунар. науч. конф. — Владикавказ : ЮМИ ВНЦ РАН, 2016. — С. 184-185.

7. Сухинов, А. И. О единственности решения линеаризованной двумерной начально-краевой задачи транспорта наносов / А. И. Сухинов, В. В. Сидорякина // Вестник Таганрог. ин-та им. А. П. Чехова. — 2016. — № 2. — С. 270-274.

8. Сухинов, А. И. Априорная оценка решения двумерной задачи транспорта наносов / А. И. Сухинов, В. В. Сидорякина, А. А. Сухинов // Актуальные проблемы прикладной математики и автоматизации : мат-лы между -нар. конф. ; Нелокальные краевые задачи и современные проблемы анализа и информатики : мат-лы XIV школы молодых ученых. — Терскол, 2016. — 363 с.

9. Sukhinov, A. I. Adaptive modified alternating triangular iterative method for solving grid equations with a non-self-adjoint operator / A. I. Sukhinov, A. E. Chistyakov // Mathematical Models and Computer Simulations. — 2012. — Vol. 4, is. 4. — P. 398-409.

10. Сухинов, А. И. Математическое моделирование транспорта наносов в прибрежных водных системах на многопроцессорной вычислительной системе / А. И. Сухинов, А. Е. Чистяков, Е. А. Проценко // Вычислительные методы и программирование. — 2014 — Т. 15, вып. 4. — С. 610-620.

11. Сравнение вычислительных эффективностей явной и неявной схем для задачи транспорта наносов в прибрежных водных системах / А. И. Сухинов [и др.] // Вычислительные методы и программирование. — 2015. — Т. 16, вып. 3. — С. 328-338.

12. Годунов, С. К. Уравнения математической физики / С. К. Годунов. — 2-е изд., исправл. и дополн. — Москва : Наука, 1979. — 392 с.

References

1. Marchuk, G.I., Dymnikov, V.P., Zalesny, V.B. Matematicheskie modeli v geofizicheskoy gidrodinamike i chislen-nye metody ikh realizatsii. [Mathematical models in geophysical hydrodynamics and numerical methods for their implementation.] Leningrad: Gidrometeoizdat, 1987, 296 p. (in Russian).

2. Sukhinov, A.I., Chistyakov, A.E., Alekseenko, E.V. Numerical realization of the three-dimensional model of hydrodynamics for shallow water basins on a high-performance system. Mathematical Models and Computer Simulations, 2011, vol. 3, iss. 5, pp. 562-574.

3. Leontyev, I.O. Pribrezhnaya dinamika: volny, techeniya potoki nanosov. [Coastal dynamics: waves, moving streams, deposits drifts.] Moscow: GEOS, 2001, 272 p. (in Russian).

4. Xiaoying Liu, et al. Predictive modeling in sediment transportation across multiple spatial scales in the Jialing River Basin of China. International Journal of Sediment Research, 2015, vol. 30, iss. 3, pp. 250-255.

5. Sukhinov, A.I., Chistyakov, A.E., Protsenko, E.A. Mathematical modeling of sediment transport in the coastal zone of shallow reservoirs. Mathematical Models and Computer Simulations, 2014, vol. 6, iss. 4, pp. 351-363.

6. Sukhinov, A.I., Sidoryakina, V.V. Sushchestvovanie i edinstvennost' resheniya linearizovannoy dvumernoy zadachi transporta nanosov. [Availability and uniqueness of solution to linearized two-dimensional sediment transportation problem.] Teoriya operatorov, kompleksnyy analiz i matematicheskoe modelirovanie: tezisy dokladov XIII mezhdunar. nauch. konf. [Theory of operators, complex analysis and mathematical modeling: Abstracts of XIIIth Int. Sci.Conf.] Vladikavkaz: YuMI VNTs RAN, 2016, pp. 184-185 (in Russian).

7. Sukhinov, A.I., Sidoryakina, V.V. O edinstvennosti resheniya linearizovannoy dvumernoy nachal'no-kraevoy zadachi transporta nanosov. [On the uniqueness of solutions of the linearized two-dimensional bottom deposit transportation problem.] Taganrog Institute Journal, 2016, no. 2, pp. 270-274 (in Russian).

g 8. Sukhinov, A.I., Sidoryakina, V.V., Sukhinov, A.A. Apriornaya otsenka resheniya dvumernoy zadachi transporta

ja nanosov. [A priori estimate of solution to two-dimensional sediment transportation problem.] Aktual'nye problemy prikladnoy д matematiki i avtomatizatsii: mat-ly mezhdunar. konf.; Nelokal'nye kraevye zadachi i sovremennye problemy analiza i in-Td formatiki: mat-ly XIV shkoly molodykh uchenykh. [Challenging issues of applied mathematics and automation: Proc. Int. rg Conf.; Nonlocal boundary-value problems and current problems of analysis and informatics: Proc. XIV School of young и scholars.] Terskol, 2016, 363 p. (in Russian).

.> 9. Sukhinov, A.I., Chistyakov, A.E. Adaptive modified alternating triangular iterative method for solving grid equa-

jd, tions with a non-self-adjoint operator. Mathematical Models and Computer Simulations, 2012, vol. 4, iss. 4, pp. 398-409. ja 10. Sukhinov, A.I., Chistyakov, A.E., Protsenko, E.A. Matematicheskoe modelirovanie transporta nanosov v

pribrezhnykh vodnykh sistemakh na mnogoprotsessornoy vychislitel'noy sisteme. [Sediment transport mathematical modeling in a coastal zone using multiprocessor computing systems.] Numerical Methods and Programming, 2014, vol. 15, iss. 4, 16 pp. 610-620 (in Russian).

11. Sukhinov, A.I., et al. Sravnenie vychislitel'nykh effektivnostey yavnoy i neyavnoy skhem dlya zadachi transporta nanosov v pribrezhnykh vodnykh sistemakh. [Comparison of computational efficiency of explicit and implicit schemes.] Numerical Methods and Programming, 2015, vol. 16, iss. 3, pp. 328-338 (in Russian).

12. Godunov, S.K. Uravneniya matematicheskoy fiziki. [Equations of mathematical physics.] 2nd revised and enlarged ed. Moscow: Nauka, 1979, 392 p. (in Russian).

Поступила в редакцию 25.12.2016 Сдана в редакцию 25.12.2016 Запланирована в номер 11.01.2017

Об авторах:

Сухинов Александр Иванович, проректор по научно-исследовательской работе и инновационной деятельности Донского государственного технического университета (РФ, 344000, Ростов-на-Дону, пл. Гагарина, 1), доктор физико-математических наук, профессор, ORCID: http://orcid.org/0000-0002-5825-1523 , Scopus Author ID: 8573972700,

WoS Researcher ID: I-1091-2016, sukhinov@gmail.com

Сидорякина Валентина Владимировна, доцент кафедры «Математика» Таганрогского института имени А. П. Чехова (филиала) РГЭУ (РИНХ), (РФ, 347936, г. Таганрог Ростовской области, ул. Инициативная, д. 48), кандидат физико-математических наук, доцент, cvv9@mail.ru

Сухинов Андрей Александрович, аспирант кафедры «Информационно-аналитические системы безопасности» Южного федерального университета (РФ, 347928, г. Таганрог, ул. Чехова, 2), honey.lusine@mail.ru

Received 25.12.2016 Submitted 25.12.2016 Scheduled in the issue 11.01.2017

Authors:

Sukhinov, Alexander I., pro-rector for Research and Innovation Activities, Don State Technical University (RF, 344000, Rostov-on-Don, Gagarin sq., 1), Dr.Sci. (Phys.-Math.), professor, ORCID: http://orcid.org/0000-0002-5825-1523 , Scopus Author ID: 8573972700,

WoS Researcher ID: I-1091-2016, sukhinov@gmail.com

Sidoryakina, Valentina V., associate professor of the Mathematics Department, Taganrog Chekhov Institute, RSUE (RINH) branch ( RF, Taganrog, Rostov Region, Initsiativnaya St., 48), Cand.Sci. (Phys.-Math.), associate professor, cvv9@mail.ru

Sukhinov, Andrey A., postgraduate student of the Information Analysis Safety Systems Department, Southern Federal University (RF, Taganrog, Chekhov St., 2), honey.lusine@mail.ru

CÖ «

S X a X <u