Научная статья на тему 'ЛИНЕАРИЗОВАННАЯ ДВУМЕРНАЯ ОБРАТНАЯ ЗАДАЧА ОПРЕДЕЛЕНИЯ ЯДРА УРАВНЕНИЯ ВЯЗКОУПРУГОСТИ'

ЛИНЕАРИЗОВАННАЯ ДВУМЕРНАЯ ОБРАТНАЯ ЗАДАЧА ОПРЕДЕЛЕНИЯ ЯДРА УРАВНЕНИЯ ВЯЗКОУПРУГОСТИ Текст научной статьи по специальности «Математика»

CC BY
45
14
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ЛИНЕЙНАЯ ВЯЗКОУПРУГОСТЬ / ОБРАТНАЯ ЗАДАЧА / ДЕЛЬТА-ФУНКЦИЯ / ПРЕОБРАЗОВАНИЕ ФУРЬЕ / ЯДРО / УСТОЙЧИВОСТЬ

Аннотация научной статьи по математике, автор научной работы — Тотиева Жанна Дмитриевна

Представлена линеаризованная обратная задача определения двумерного ядра для системы уравнений линейной динамической вязкоупругости с сосредоточенным источником возмущений на свободной поверхности. Искомой величиной в поставленной задаче является ядро интегрального оператора, моделирующего явление памяти, которое имеет место при распространении волновых процессов в вязкоупругих средах. Прямая начально-краевая задача для вектор-функции смещения содержит нулевые начальные данные и граничное условие Неймана на дневной поверхности специального вида. Для линеаризации искомое ядро разлагается на две составляющие, одна из которых малая по абсолютной величине неизвестная добавка. В качестве дополнительной информации задается отклик линеаризованного поля смещений точек среды на свободной поверхности. В предположении, что коэффициенты системы зависят от одной пространственной переменной, прямая задача сводится к начально-краевой задаче для одного интегро-дифференциального уравнения гиперболического типа второго порядка. Доказывается, что поставленная линеаризованная задача определения сверточного ядра эквивалента некоторой системе линейных интегральных уравнений типа Вольтерра второго рода. К последней применяется обобщенный принцип сжатых отображений. Доказаны теоремы глобальной однозначной разрешимости в пространстве непрерывных функций и устойчивости решения обратной задачи. Приводится теорема о сходимости регуляризованного семейства задач к решению исходной (некорректной) задачи.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Тотиева Жанна Дмитриевна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

LINEARIZED TWO-DIMENSIONAL INVERSE PROBLEM OF DETERMINING THE KERNEL OF THE VISCOELASTICITY EQUATION

A linearized inverse problem of determining the 2D convolutional kernel of the integral term in an integro-differential viscoelasticity equation is considered. The direct problem is represented by a generalized initial-boundary value problem for this equation with zero initial data and the Neumann boundary condition in the form of the Dirac delta-function. The unknown kernel is decomposed into two components, one of which is a small in absolute value unknown additive. For solving the inverse problem, the traces of the solution to the direct problem on the domain boundary are given as an additional condition. It is proved that the linearized problem of determining the convolutional kernel is equivalent to a system of linear Volterra type integral equations. The generalized contraction mapping principle is applied. The main result of the article is the theorem of global unique solvability of the inverse problem in the class of continuous functions. A theorem on the convergence of a regularized family of problems to the~solution of~the~original (ill-posed) problem is presented.

Текст научной работы на тему «ЛИНЕАРИЗОВАННАЯ ДВУМЕРНАЯ ОБРАТНАЯ ЗАДАЧА ОПРЕДЕЛЕНИЯ ЯДРА УРАВНЕНИЯ ВЯЗКОУПРУГОСТИ»

Владикавказский математический журнал 2021, Том 23, Выпуск 2, С. 87-103

УДК 517.958

DOI 10.46698/u2193-3754-6534-u

ЛИНЕАРИЗОВАННАЯ ДВУМЕРНАЯ ОБРАТНАЯ ЗАДАЧА ОПРЕДЕЛЕНИЯ ЯДРА УРАВНЕНИЯ ВЯЗКОУПРУГОСТИ

Ж. Д. Тотиева1

1 Южный математический институт — филиал ВНЦ РАН, Россия, 362027, Владикавказ, ул. Маркуса, 22 E-mail: j annatuaeva@inbox. ru

Аннотация. Представлена линеаризованная обратная задача определения двумерного ядра для системы уравнений линейной динамической вязкоупругости с сосредоточенным источником возмущений на свободной поверхности. Искомой величиной в поставленной задаче является ядро интегрального оператора, моделирующего явление памяти, которое имеет место при распространении волновых процессов в вязкоупругих средах. Прямая начально-краевая задача для вектор-функции смещения содержит нулевые начальные данные и граничное условие Неймана на дневной поверхности специального вида. Для линеаризации искомое ядро разлагается на две составляющие, одна из которых малая по абсолютной величине неизвестная добавка. В качестве дополнительной информации задается отклик линеаризованного поля смещений точек среды на свободной поверхности. В предположении, что коэффициенты системы зависят от одной пространственной переменной, прямая задача сводится к начально-краевой задаче для одного интегро-дифференциального уравнения гиперболического типа второго порядка. Доказывается, что поставленная линеаризованная задача определения сверточного ядра эквивалента некоторой системе линейных интегральных уравнений типа Вольтерра второго рода. К последней применяется обобщенный принцип сжатых отображений. Доказаны теоремы глобальной однозначной разрешимости в пространстве непрерывных функций и устойчивости решения обратной задачи. Приводится теорема о сходимости регуляризованного семейства задач к решению исходной (некорректной) задачи.

Ключевые слова: линейная вязкоупругость, обратная задача, дельта-функция, преобразование Фурье, ядро, устойчивость.

Mathematical Subject Classification (2010): 35L20, 35R30, 35Q99.

Образец цитирования: Тотиева Ж. Д. Линеаризованная двумерная обратная задача определения ядра уравнения вязкоупругости // Владикавк. мат. журн.—2021.—Т. 23, вып. 2.—С. 87-103. DOI: 10.46698/u2193-3754-6534-u.

1. Постановка задачи

Рассмотрим при x = (x1,x2,x3) € R+, t € R, R+ = {x € R3 интегро-дифференциальных уравнений

Дад, ¿-12 3 j=1 3

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

«i\t<0 = 0,

: x3 > 0} систему (1.1)

(1.2)

© 2021 Тотиева Ж. Д.

П3\хз=+о = -5-^, ; = 1,2,3, (1.3)

где и(ж,£) = (и^ж, ¿), и2(ж, ¿), и3(ж,£)) — вектор смещений, ¿(-) — дельта-функция Дирака; Ту — тензор напряжений:

г

Ту (ж,£) = ау[и](ж, ^ к(ж2,£ — т)ау [и](ж,т) ^т, (1.4)

= /х ( —— + —— ) + ¿у А Луи. (1.5)

/ дщ ди3 \ дх3 9ж.

Здесь ¿у — символ Кронекера.

Предполагается, что р = р(ж3), ^ = ^(ж3), А = А(ж3) являются функциями одной переменной, удовлетворяющими условиям р(ж3) ^ т > 0, ^(ж3) ^ т > 0, А(ж3) ^ т > 0, причем р'(+0) = ^'(+0) = А'(+0) = 0.

Прямая .задача заключается в отыскании вектор-функции и(ж, ¿) из системы уравнений (1.1) при соответствующих начальных и граничных условиях (1.2), (1.3).

При сделанных предположениях из равенств (1.1)—(1.5) следует, что решение прямой задачи и(ж,£) не зависит от переменной ж1: и(ж,£) = и(ж2,ж3,£) [1].

Предполагаем, что ядро й(ж2,£) можно представить в следующем виде [2]:

к(ж2,^) = ко(ж2^) + (ж2, ¿),

где функция й0(ж2,£) является заданной, а ^(ж^) — неизвестная, малая по абсолютной величине, добавка. Требование малости добавочного ядра понимается как малость по норме, содержащей производные ядра до некоторого порядка. Суть метода линеаризации заключается в следующем. Формально вводится параметр е:

к(ж2,^) = ко(ж2,^) + еЙ1(ж2,^). (1.6)

В рамках данного исследования положим й0(ж2,£) = й0(£).

Решение прямой задачи (1.1)—(1.5) в предположении (1.6) будем искать в виде ряда по степеням е:

те

и(ж2,ж3,£) = ^ иу (ж2,ж3,£). (1.7)

У=о

Подставляя (1.6), (1.7) в (1.1) и приравнивая члены, стоящие при , j = 0,1, получаем две прямые задачи:

(1) Задача определения и0(ж2,ж3,^ = (и°, , и0)(ж2, ж3, ¿) из равенств

М = . = 1,2,3, (1.8)

И сП2 ^ дх^ ' ' '

у=1

и0|г<0 = 0, (1.9)

Т3у |хз=+0 = -¿у ¿(*), ; = 1, 2, 3, (1.10)

где

г

Ту (ж2, ж3, ¿) = ау [и0](ж,£) + У - т)ау [и0](ж,т) ^т.

(л) Задача определения и}(ж2,ж3,£) = (и}, и2, и})(ж2, ж3, £) из равенств

Й = < = ^ (1.П)

у т2 ^ дхт'

и1 ко = 0, (1.12)

Тэ1г|хз=+о = 0, г = 1,2,3, (1.13)

где

(Ж2, Жз, £) = СТ, [и}](ж, £) + / Й}(ж2,^ - [ио](ж,т) + - [и}](ж,т) ^Т. (1.14)

ЧГ

о

Обратная задача заключается в определении ядра ^}(ж2, ¿), £ > 0, входящего в (1.13) посредством формулы (1.14), если относительно решения линеаризованной задачи (1.11)—(1.13) известна дополнительная информация

и

1Ц=+о = Я(ж2,*), (1.15)

где д(ж2,£) — заданная функция.

Определение. Функция й}(ж2,£) из класса непрерывных функций С(М х [0, то)) называется решением обратной задачи (1.11)—(1.15), если соответствующее ей решение прямой задачи (1.11)-(1.14) и(ж,£) из класса обобщенных функций ^'(М+ х М) удовлетворяет (1.15) для д(ж2,£), принадлежащей классу обобщенных функций ^'(М х [0, то)).

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

Задачи определения ядра (зависящего от временной и пространственных переменных) интегрального оператора — направление в теории обратных задач, возникшее в конце прошлого столетия [3-9]. Более подробный анализ источников по данному направлению представлен в монографии [10], которая является одной из последних фундаментальных работ в области исследования обратных задач для сред с памятью (или с последействием). В ней представлены результаты исследования корректности ряда постановок одномерных и многомерных обратных динамических задач для гиперболических интегро-дифференциальных уравнений, возникающих при описании внутренних характеристик сред с последействием по измерениям волнового поля в доступных областях. Доказаны теоремы об однозначной разрешимости поставленных обратных задач, а также получены оценки непрерывной зависимости решений этих задач от входных данных.

Из первых результатов по обратным задачам линейной вязкоупругости, близким к данной, можно отметить [11-13]. Дальнейшее развитие исследований отражено, например, в работах [14-20]. Особый интерес представляют многомерные обратные задачи по определению ядер, когда искомая функция зависит от двух и более переменных. Многомерная обратная задача для (1.1) с начальными и граничными условиями (1.2) и (1.3), дополнительной информацией (1.15) исследована в [21]. В этой работе на основе метода шкал банаховых пространств получена локальная однозначная разрешимость задачи определения ядра &(ж2,£) в классе функций, аналитических по переменной ж2 и гладких по переменной

г

Из результатов по численному исследованию обратных задач для сред с памятью можно отметить работы [22-25].

В работах [26, 27] рассмотрены задачи определения двумерного ядра интегро-диф-ференциального уравнения со слабо горизонтальной однородностью.

Прямая задача (1.8)-(1.10) изучена в работе [1]. Показано, что она распадается на две независимо решаемые задачи по определению и°(ж2,жз,*) и (и°,■и°)(ж2,ж3,*) соответственно. Причем для второй задачи получаем однородную систему уравнений с однородными граничными условиями и нулевыми начальными данными. Поэтому и2 = и° = 0. Кроме того, -и°(ж2,ж3,*) = -и°(ж3,*). Его структура имеет вид

и?(гЧу),*) = Ку)

ъ

V°(у,*) + У г°(* - тЖ°(у,т) ^т

V°(у,*) = ^ -

+ «(у,*)

где

жз I-

У = ФЫ) := / ^у, 8(у) :=

с(+0)р(+0)

с (ф-(у)) р (ф-?(у))

гоСО = - / - т)го(г) (1т, а = И+0)р(+0)]" ,

6>(*) = 1, * ^ 0, 6>(*) = 0, * < 0.

Функция г>(у,*) — решение интегрального уравнения Вольтерра второго рода в области Бт = {(у,*) : |у| ^ * ^ Т — |у|} для фиксированного Т > 0

УШ)= ] I (1.

16)

.= *"(у) 2

*>'(У)

«(у)

При предположении к°(*) € С3[0,Т] с помощью дифференцирования по параметру можно показать, что г>(у,*) € С2 (* ^ |у|). Из (1.16) следует, что ^|ъ=|у| = 0.

В работе [28] для системы уравнений (1.8)-(1.10) изучена обратная задача определения к°(*). В дальнейшем будем считать функции к°(£) и и°(жз,*) известными величинами и допустим, что к°(0) = 0. Это условие существенно облегчает выкладки и не меняет сути исследования.

Основная цель представленной работы — построение метода нахождения к? (ж2, *).

Результатами исследования являются теоремы глобальной однозначной разрешимости и устойчивости решения обратной задачи.

2. Задача определения функций к? и и?

Определим билинейный интегральный оператор Ь по формуле

ъ

Ь[к°(*),и(ж2,ж3,*)] = и(ж2, ж3, *) + У к°(* — т)и(ж2, ж3, т) ^т,

а

2

ъ

где и — скалярная функция. В дальнейшем для сокращения записи иногда не будем в операторе Ь указывать зависимость функций от переменных, подразумевая зависимость первой функции от переменной а второй — от ж2, жз,

Из (1.11)-(1.15) следует, что компонента вектор-функции и} = и1(ж2,ж3,£) при ж3 > 0, (ж2,£) € М2, удовлетворяет следующим равенствам:

д2и1(ж2,жз,£) Р&*)-Ы2- =Ь

д ( ди1 \ д2и1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

о

и11г<о = 0,

Ь

ди1

джз

г

+ / ^1(ж2,£ - Т)^(жз)

ди1(ж3, т ) дх3

0,

хз=+о

(2.2)

(2.3)

(2.4)

и1(ж2, 0, £) = #(ж2,£).

Перейдем от функции и1(ж2,жз,£) к ее образу Фурье и!(^, жз,£) := [и1](^, жз, £), причем вирри1(^, жз,£) С [—А, А], где Л — фиксированное положительное число.

Будем считать, что £) = [й1](^, £) € Л(А, Т) тогда и только тогда, когда А;1 (V, £) € С1 (М х М+) и для любого фиксированного £ € [0, Т] вирр А;1 (V, £) С [—Л, Л]. Соответственно, й1(ж2,£) € Л(А, Т) тогда и только тогда, когда А;^, £) = [А^]^, £) € Л(А,Т). Тогда обратная задача (2.1)-(2.4) в терминах функции и1 перепишется в виде

р(жз)

д2и!^, жз, £)

дР

Ь

+ Мг, (2.5)

о

и 11 г< о = 0,

М+0) ь

ди1

джз

+ / ЫМ - т)

ди\{ь>, х3,т) дх3

0,

хз=+о

Ul(v, 0,£) = ^ ЫМ) := ¿г(^*).

(2.6)

(2.7)

(2.8)

Пусть

У1(^у,*) := Ь

Ло,

Тогда (2.5)-(2.8) для у > 0, £ € М примут вид

д£2

ду2

г

+ Я(V, у)У1 -I <(* - т)^,у,т) ^т

+ у А;^ - т) о

Э2Ь[го,Уо(у,т)]

ду2

+ 9(у)Ь[го, ^о(у,т)]

^Т,

(2.9)

г

г

г

у=+°

где

V |у=°= Ь [к°,д^, *)] ,

Н(у^) := д(у) — V2с2(^-!(у)).

(2.10) (2.11) (2.12)

Имеем

Ь[г°,^(у,*)] = 0(* — у)

+ и(у,*)

ъ

+ I г°(* — т)0(т — у)

+ и(у,т)

^т,

9Ь[го,Уо(у,г)] Л к^и^- т)—-——-- йт

д ъ ( Т

— У кг(и, "(у,т) + / г°(т ~ ^

у у

+ и(у,п)

¿П ? ¿т (2.13)

= + / ЫМ-тМ ¿о

ди

~ 2Го(т - у) > ¿т.

Поступая аналогично,

— т)

д2Ь[г°, М°(у,т)] _а

ду

2 (¿Г = -к^г - у) + - у)Ь0

Го,^(у,у) ду

+ у Л!^,* — т К Ь°

у

Г°' Эу2 Т

+ 2Го(^-у) >

(2.14)

Здесь и далее, к примеру, А;?, д" означают операции однократного и двукратного дифференцирования по переменной * соответствующих функций.

С учетом (2.13), (2.14) задача (2.9)-(2.12) перепишется в следующем виде для у > 0, * € М:

д2^ д2У? тт. ,т^ а. . ~ , ,т

го,^(У,У)

ъ ъ

/

(2.15)

— у Г° (* — т)У!^,у,т) ^т + ^ * — т)р(у, т) ¿т,

уу

<9у

а

у=+°

Г°, А! (V, *)

(2.16)

0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

а

а

2

2

ъ

а

2

ъ

ъ

ъ

У1|у=О = Ь [Л£)]

(2.17)

где

Р(у,£) := Ьо

=0, 1 \г=у '

д2-и

Г°> ¿^2 (У> Т) + ^ | + Т)

(2.18)

В равенстве для р(у,£) оператор Ьо отличается от оператора Ь тем, что нижний индекс в интеграле оператора заменен на у. С помощью формулы Даламбера получаем

г+у г—у

+ Я (^СЖ(^С,Т ) (2.19)

у ¿-Ну—^ ✓

+ \ [ [ ¡^(чт-О+Ычт-ОЬо

о г—у+?

ду

+ У - п)р(С,п) - Го'(т - (V, С, п))

Переходя в равенстве (2.19) к пределу £ ^ у + 0 с учетом (2.18) и (/(V, 0) = 0, имеем

у 2у—

+ / Мо^'г-О + Ы^т-О^о

о ?

ду

+ Я (^сЖ^т) (2.20)

+ у [- п)р(С,п) - Го'(т - п)^1 (^С,п)) ¿пМ^С Заменяя 2у на £ и дифференцируя (2.20) по получаем

г—?

ду

+ Я- С)

+ у (^(М - С - Т )р(С,т) - го'(£ - С - Т )^,С,Т)) ¿ТК.

(2.21)

Дифференцируя по * (2.21) (предварительно сделав замену переменной во втором интеграле * — 2{ = т), получаем

ъ

2

к[(и,1) = [**),£>,*)] +*о(*№0)) -Ц Го(Ь-т)к1(и, т)йт

°

°

дМ

V ^ V ^ ^

+а?^, 0)р(е, * — о — г°'(0) V? (V, е,*—е)

— ъ-2? .

+ У У^-а-тШ^йт- I ? ° ^

Уравнение для получается дифференцированием по £ уравнения (2.19)

а / ~ ~ \ а -

+ + 2/)] - ¿[г0,А;1(1/,4 -у)]] + ^уА^М -у)

у (

+я (V, е) (V? +у—е) — V? —у+о) (2.23)

+ ( +у — е — п)р(е,п) — г°'(*+у—е — пж? ^,е,п)) ¿п

ъ+у-?

ъ-у+?

( (V, *—у+е — п)р(е, п) — г°'(* — у+е — пж? (V, е, п)) ¿п ? ¿е.

Добавляя к уравнениям нижеприведенное

ъ

й?^,*) = й?^, 0) + у"(г — т)а! (V, т)^т, (2.24)

°

получаем замкнутую линейную систему интегральных уравнений Вольтерра второго рода (2.19), (2.22)-(2.24) в области £>т относительно ^1(1/,у,, ЫМ)-Так как в данную систему входит неизвестная величина А;?^, 0), то ее можно найти из уравнения (2.21) при * = 0:

2

^1(^,0) = ——д(и,0).

Основными результатами этого параграфа являются следующие теоремы однозначной глобальной разрешимости и устойчивости обратной задачи определения А;?^, *).

Теорема 2.1. Пусть Л, Т — фиксированные положительные числа. Для существования и единственности решения обратной задачи (1.11)—(1.15) /^(ж^) € Л(Л,Т), необходимо и достаточно, чтобы А0(£) € С3[0,Т], € С42(М х [0, Т]), 0) = 0, и для любого фиксированного £ € [0, Т] вирр£) С [—Л, Л].

Теорема 2.2. Пусть /(1)(ж2, £), /(2) € Л(Л, Т) — решения обратной зада-

чи (1.11)—(1.15), отвечающие информациям $(1)(ж2,£), $(2)(ж2,£) соответственно. Тогда при выполнении условий теоремы 2.1 имеет оценка устойчивости

^ — А^П2

С[0,т ]

^ж2 ^ С

|д(1) — й(2) П2 ^ I9 9 II с2[о,т]

(2.25)

где С — некоторая константа, зависящая от величин Л, Т и значений функций ^(жз), р(жз), Ао(£).

< Доказательство теоремы 2.1. Обратная задача (1.11)-(1.15) эквивалентна системе интегральных уравнений (2.19), (2.22)-(2.24). Данная система является замкнутой линейной системой интегральных уравнений Вольтерра второго рода с непрерывными свободными членами и ядрами относительно неизвестных функций в области при V € М. Идея доказательства существования единственного решения данной системы состоит в применении обобщенного принципа сжатых отображений. Запишем систему (2.19), (2.22)-(2.24) в виде операторного уравнения

ф = Вф, (2.26)

Ф = [Ф1^,у,£), ф2(М), фз(^у,£), ф4(М)] :=

VI(V, у, £), (V, £),

^ У, - | (^1 * + 2/) - * - 2/)) - | ("> * - У)> (">

Оператор В = (В^ В2, В3, В4) определен на множестве функций ф € С(М х ), где

у

у

22

о

+ Н (^)ф1(^,т)

+ ^ ^Ф4(V, т — п)р(С,п) — Го'(т — п)Ф1^,£,п)) ^Мт^,

в2ф = фо2(», I) - ^ I го- 2£)ф420 ^

+ Н(и, О^(г/, + Ф4(*Л 0)р(£, 4 - О

■Л^Олт-ОЬо

0

——? — ч

л

где

ду

Взф = 0Оз(*Л y,t) + j j r0(t - т) [04(z/, Т + у) - 04(г/, г - у)] (1т

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

о

У (

+ ^ J<%Mv,t + V-4) + (Mv,t + V-4)-Mt-y))Lo

о ^

+ H(v, С) (01 (v, + у - О - 0i(v,£,i - у + о)

t+y-?

+ J (04(v,t + у - С - п)р(С,п) - Го'(t + у - С - n)0i(v,£,n)) dn - / (04(v,t - у + С - n)p(C,n) - Го'(t - у + С - n)0i(v,C,n))dn>d£,

t 2

B40 = 004(v) + y"(t - 2C)02(v, 2C) dC,

dvi a

-QffaV't) =fo(v,y,t) + j('il)4(v,t + y) -Ip4(v,t -y) +yip2(v,t -y)), 0oi(^,y,i) = ^L[Äo,5f(i/,i -y)] +L[fc0,5(i/,i+ 2/)]),

002(i) = "¿(b [WM)]

0оз(гЛУ,£) = ^(¿[fo^'M-у)] + y)]),

004 (v, у, t) = fci(v, 0).

Покажем теперь, что некоторая степень n (n — натуральное число) линейного отображения B0 является сжатием. Положим

V) = та^ тах |0? у, *))|, тах |^2 4(V, |

I (у , ъ)еДТ (ъ)е[°,т]

Пусть 0(?), 0(2) — две непрерывные вектор-функции в М х Бт, удовлетворяющие линейной системе интегральных уравнений (2.26). Обозначим

ЛЫ) = |(С,т) : 0 < С < у, t - у + С < т < t + у - С^, ЦуЛС) = {т : (С,т) € Л(у,t)}.

Тогда для (V, у, £) € М х имеем (в оценках используем факт, что в уравнениях (2.22), (2.24) £ = 2у)

—у

^ дЛ тах ^ тах

' 1 Т е£(у,а)

0

В1ф(1) — В1ф(2) (V, у, £)

ф(1)^,е,т)—ф(2)^,е,т), ф214(2О — ф224(2О

¿е

—у

^ д 1 там тах

7 (?,т)€дт

0

ф(1)^,е,т) — фГ ^,е,т)

(2) 1

, т ах

ф$(20 — ф224 (20

¿е

< Д1У

ф(1) — ф(2)

далее

В2ф(1) — В2ф(2) (2у)

^ Д2 I та^ тах

I Ьф,|]

ф31) (е,т — е^)—ф32)(е,т—е^)

т ах

^ [о,?]

ф214(2е) — ф224 (2е)

¿е

< Д2у

ф(1) — ф(2)

—у

^ д^ тах ^ тах

т €Е(у,*,0

0

В3ф(1) — В3ф(2)|(^у,£)

ф(1) (V, е,т)—ф(2) (^е,т) |, | ф214 (2е) — ф224 (2е) |} ¿е

< Д3у ф(1) — ф(2)

| —у | |

В4ф(1) — В4ф(2) (2у) < дЛ 1Пах п ф21)(2е) — ф22)(2е) ^ < Д4у ф(1) — ф(2)

{

где д^- — константы, зависящие от величин, входящих в С (теорема 2.2). Полагая М = тах{д1, д2, д3, д4}, получаем, что

тах 1<К4

В,ф(1) — В,ф(2) (V, у, £) < Му ф(1) — ф(2) , (V, у,£) € М х

Далее,

^ ди та^ тах

7 I ?€£(у,4,т) 0

В2ф(1) — В2ф(2) (^у,£)

В1ф(1) (v,е,т) — В1ф(2) (v,е,т), В2,4ф(1) (2е) — В2,4ф(2)(2е)

¿е

у

< д1^у е

0

ф(1) — ф(2)

ф(1) — ф(2)

у

у

—у

^ д2 та^< тах

У I Т€[£,2у-£]

0

В|ф(1) — В|ф(2) (2у)

Bзф(l)(v,е,т) — В3ф(2)(v,е,т), В2,4ф(1)(2е) — В2,4ф(2)(2е)

у

< д2М у е

0

ф(1) — ф(2)

¿е < д2М

ф(1) — ф(2)

^ дн та^< тах

Т еЕ(у,а)

0

В2ф(1) — В2ф(2) (V, у, £) Blф(l)(v,е,т) — В1ф(2)(v,е,т), В2,4ф(1)(2е) — В2,4ф(2)(2е)

< дзму е 0

ф(1) — ф(2)

ф(1) — ф(2)

| —у | В42ф(1) — В42ф(2) (2у) < д^ тах В2ф(1)(2е) — В2ф(2)(2е)

¿е

—у

< д4^ е

0

ф(1) — ф(2)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Отсюда

тах

ККз

¿е < д4М

ф(1) — ф(2)

и, вообще,

тах

КК4

б2ф(1) - В2ф(2) (г/, у, *) < М2|у ф(1) - Ф{2) , (V, г/, *) € М х Дт

(V, у,£) € М х ,

п!

В пф(1) — В пф(2)

^ М

т\п

га v 2 /

П!

ф(1) — ф(2)

¿е

¿е

При любом фиксированном Т число п можно выбрать настолько большим, что

/у\га

МгаЩ- := а < 1. п!

Тогда отображение Вп является сжатием. Согласно обобщению принципа сжимающих отображений уравнение Вф = ф имеет одно и только одно решение, принадлежащее С (М х ). Данное решение может быть найдено методом последовательных приближений. >

< Доказательство теоремы 2.2. Пусть ф(:,) — вектор-функции, которые являются решениями (2.26), т. е. справедливы уравнения ф(:,) = Вф(:,).

Переходя в этом выражении к разностям ф(1) — ф(2), 9(1) — 9(2), из рассуждений, проведенных при доказательстве теоремы 3.1, получим оценку

ф(1) — ф(2) (V) ^

а

9(1) — 9(2)

с2[о,т ]

—у

(V)+^ ||ф(1) — ф(2)| (v,е) ¿е, (2.27)

2

у

у

2

2

где постоянные а, в зависят от величин, входящих в С. Из неравенства (2.27) следует, что

^(1) - ^(2) ) < С

^(1) - ^(2)

С2 [0,Т ]

(V).

Неравенство (2.25) вытекает из следующего свойства изометричности оператора преобразования Фурье .Х2:

(V,*) — (V, ¿)|2 ^ = I |к(1)(Х2,^) — к(2)(Ж2,^)|2 ^2. >

Теперь рассмотрим вопрос о сходимости последовательности решений некоторого семейства задач к искомому решению. Пусть существует решение ^1(^2,*) € ^(К х К+), отвечающее информации (/(V, *).

Определим множество функций *) по правилу *) := 0(Л — IV*).

Выделим семейство обратных задач: определить функцию кЛ(х2,*) = .Р—1^^, *)] по информации *).

Теорема 2.3. Данное семейство является регуляризованным, т. е.

1) для каждого Л > 0 обратная задача корректна;

2) если данные таковы, что решение исходной (корректной) задачи существует, то при Л ^ то последовательность решений семейства задач с этими данными стремится к решению исходной (некорректной) задачи.

< Теоремы 2.1 и 2.2, доказанные ранее, утверждают о корректности обратной задачи. Теперь покажем, что

Иш вир / IкЛ(х2,*) — к1(ж2,*)|2= 0.

Л^те 4€[0,ТЫ

ге[о,т ]

Действительно,

Л

Иш вир I |кЛ(х2,*) — к1 (х2, *)|2 = Иш вир I / *) — (V, *)|2 ^

Л^те 4е [о,т ] К Л^те [о,т ] Л

+те

+ / |fc1(v,í)|2dv + / [й^,*)!2^ | = Иш вир

Л^те ¿е[о,т ]

— те Л / \—те

—Л +те

[ |Й1(М)|2^ + / I =0.

Л

Тем самым теорема 2.3 доказана. >

Литература

1. Туаева Ж. Д. Многомерная математическая модель сейсмики с памятью // Исслед. по диф. ур-ям и мат. моделированию.—Владикавказ: ВНЦ РАН, 2008.—C. 294-303.

2. Романов В. Г. Обратные задачи математической физики.—М.: Наука, 1984.

3. Lorenzi A., Sinestrari E. An inverse problem in the theory of materials with memory // Nonlinear Anal. TMA.—1988.—Vol. 12, № 12.—P. 1317-1335. DOI: 10.1016/0362-546X(88)90080-6.

4. Lorenzi А. An inverse problem in the theory of materials with memory II // J. Semigroup Theory and Applications. Ser. Pure and Appl. Math.—1989.—Vol. 116.—P. 261-290.

5. Дурдиев Д. K. Обратная задача для трехмерного волнового уравнения в среде с памятью // Мат. анализ и дискретная математика.—Новосибирск: Изд-во Новосибирского ун-та, 1989.—C. 19-27.

6. Lorenzi A., Paparoni E. Direct and inverse problems in the theory of materials with memory // Rend. Sem. Mat. Univ. Padova.—1992.—Vol. 87.—P. 105-138.

7. Bukhgeym A. L. Inverse problems of memory reconstruction // J. of Inverse and Ill-posed Problems.— 1993.—Vol. 1, № 3.—P. 193-205. DOI: 10.1515/jiip.1993.1.3.193.

8. Дурдиев Д. К. Многомерная обратная задача для уравнения с памятью // Сиб. мат. журн.—1994.— Т. 35, № 3.—С. 574-582.

9. Bukhgeim A. L., Dyatlov G. V. Inverse problems for equations with memory // SIAM J. Math. Fool.— 1998.—Vol. 1, № 2.—P. 1-17.

10. Дурдиев Д. К. Обратные задачи для сред с последействием.—Ташкент: Турон-Икбол, 2014.

11. Lorenzi A., Ulekova J. Sh., Yakhno V. G. An inverse problem in viscoelasticity // J. of Inverse and Ill-posed Problems.—1994.—Vol. 2, № 2.—P. 131-164. DOI: doi.org/10.1515/jiip.1994.2.2.131.

12. Janno J., Von Wolfersdorf L. Inverse problems for identification of memory kernels in viscoelasticity // Math. Methods in Appl. Sciences.—1997.—Vol. 20, № 4.—P. 291-314. DOI: 10.1002/(SICI)1099-1476(19970310)20:4<291::AID-MMA860>3.0.C0;2-W.

13. Janno J., Von Wolfersdorf L. An inverse problem for identification of a time- and space-dependent memory kernel in viscoelasticity // Inverse Problems.—2001.—Vol. 17, № 1.—P. 13-24. DOI: 10.1088/0266-5611/17/1/302.

14. Lorenzi A., Messina F., Romanov V. G. Recovering a Lame kernel in a viscoelastic system // Applicable Analysis.—2007.—Vol. 86, № 11.—P. 1375-1395. DOI: 10.1080/00036810701675183.

15. Romanov V. G., Yamamoto M. Recovering a Lame kernel in a viscoelastic equation by a single boundary measurement // Applicable Analysis.—2010.—Vol. 89, № 3.—P. 377-390. DOI: 10.1080/00036810903518975.

16. Lorenzi A., Romanov V. G. Recovering two Lame kernels in a viscoelastic system // Inverse Probl. Imaging.—2011.—Vol. 5, № 2.—P. 431-464. DOI: 10.3934/ipi.2011.5.431.

17. Романов В. Г. Двумерная обратная задача для уравнения вязкоупругости // Сиб. мат. журн.— 2012.—T. 53, № 6.—C. 1401-1412.

18. Romanov V. G. Inverse problems for differential equations with memory // Eurasian J. of Mathematical and Computer Applications.—2014.—Vol. 2, № 4.—P. 51-80.

19. Дурдиев Д. К., Тотиева Ж. Д. Задача об определении одномерного ядра уравнения электровяз-коупругости // Сиб. мат. журн.—2017.—Т. 58, № 3.—С. 553-572. DOI: 10.17377/smzh.2017.58.307.

20. Дурдиев Д. К., Рахмонов А. А. Обратная задача для системы интегро-дифференциальных уравнений SH-волн в вязкоупругой пористой среде: глобальная разрешимость // Теор. и мат. физика.— 2018.—T. 195, № 3.—C. 491-506. DOI: 10.4213/tmf9480.

21. Дурдиев Д. К., Тотиева Ж. Д. Задача об определении многомерного ядра уравнения вязкоупругости // Владикавк. мат. журн.—2015.—T. 17, № 4.—С. 18-43. DOI: 10.23671/VNC.2015.4.5969.

22. Карчевский А. Л., Фатьянов А. Г. Численное решение обратной задачи для системы упругости с последействием для вертикально неоднородной среды // Сиб. журн. вычисл. матем.—2001.—Т. 4, № 3.—С. 259-268.

23. Дурдиев У. Д. Численное определение зависимости диэлектрической проницаемости слоистой среды от временной частоты // Сиб. электрон. мат. изв.—2020.—Т. 17.—С. 179-189. DOI: 10.33048/semi.2020.17.013.

24. Bozorov Z. R. Numerical determining a memory function of a horizontally-stratified elastic medium with aftereffect // Eurasian J. of Mathematical and Computer Applications.—2020.—Vol. 8, № 2.—P. 4-16.

25. Kabanikhin S. I., Karchevsky A. L., Lorenzi A. Lavrent'ev regularization of solutions to linear integro-differential inverse problems // J. of Inverse and Ill-posed Problems.—1993.—Vol. 1, № 2.—P. 115-140. DOI: 10.1515/jiip.1993.1.3.193.

26. Дурдиев Д. К., Бозоров З. Р. Задача определения ядра интегро-дифференциального волнового уравнения со слабо горизонтальной однородностью // Дальневост. мат. журн.—2013.—Т. 13, № 2.— C. 209-221.

27. Тотиева Ж. Д. Определение ядра уравнения вязкоупругости в слабо горизонтально-неоднородной среде // Сиб. мат. журн.—2020.—Т. 61, № 2.—С. 453-475. DOI: 10.33048/smzh.2020.61.217.

28. Дурдиев Д. К., Тотиева Ж. Д. Задача об определении одномерного ядра уравнения вязкоупругости // Сиб. журн. индустр. матем.—2013.—Т. 16, № 2.—С. 72-82.

Статья поступила 11 января 2021 г. Тотиева Жанна Дмитриевна

Южный математический институт — филиал ВНЦ РАН, старший научный сотрудник отдела математического моделирования РОССИЯ, 362027, Владикавказ, ул. Маркуса, 22; E-mail: j annatuaeva@inbox. ru

Vladikavkaz Mathematical Journal 2021, Volume 23, Issue 2, P. 87-103

LINEARIZED TWO-DIMENSIONAL INVERSE PROBLEM OF DETERMINING THE KERNEL OF THE VISCOELASTICITY EQUATION

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Totieva, Zh. D.1

1 Southern Mathematical Institute VSC RAS, 22 Marcus St., Vladikavkaz 362027, Russia E-mail: j annatuaeva@inbox. ru

Abstract. A linearized inverse problem of determining the 2D convolutional kernel of the integral term in an integro-differential viscoelasticity equation is considered. The direct problem is represented by a generalized initial-boundary value problem for this equation with zero initial data and the Neumann boundary condition in the form of the Dirac delta-function. The unknown kernel is decomposed into two components, one of which is a small in absolute value unknown additive. For solving the inverse problem, the traces of the solution to the direct problem on the domain boundary are given as an additional condition. It is proved that the linearized problem of determining the convolutional kernel is equivalent to a system of linear Volterra type integral equations. The generalized contraction mapping principle is applied. The main result of the article is the theorem of global unique solvability of the inverse problem in the class of continuous functions. A theorem on the convergence of a regularized family of problems to the solution of the original (ill-posed) problem is presented.

Key words: linear viscoelasticity, inverse problem, delta function, Fourier transform, kernel, stability.

Mathematical Subject Classification (2010): 35L20, 35R30, 35Q99.

For citation: Totieva, Zh. D. Linearized Two-Dimensional Inverse Problem of Determining the Kernel of the Viscoelasticity Equation, Vladikavkaz Math. J., 2021, vol. 23, no. 2, pp. 87-103 (in Russian). DOI: 10.46698/u2193-3754-6534-u.

References

1. Tuaeva, Zh. D. Multidimensional Mathematical Model of Seismic Memory, Issledovaniya po differen-cial'nym uravneniyam i matematicheskomu modelirovaniyu [Research on Differential Equations and Mathematical Modeling], Vladikavkaz, VNC RAN, 2008, p. 297-306 (in Russian).

2. Romanov, V. G. Obratnye zadachi matematicheskoj fiziki, Moskva, Nauka, 1984 (in Russian).

3. Lorenzi, A. and Sinestrari, E. An Inverse Problem in the Theory of Materials with Memory I, Nonlinear Analysis: Theory, Methods and Applications, 1988, vol. 12, no. 12, pp. 1217-1335. DOI: 10.1016/0362-546X(88)90080-6.

4. Lorenzi, А. An inverse Problem in the Theory of Materials with Memory II, Journal of Semigroup Theory and Applications, Series on Pure and Applied Mathematics, 1989, vol. 116, pp. 261-290.

5. Durdiev, D. K. The Inverse Problem for a Three-Dimensional Wave Equation in a Memory Environmentu, Matematicheskij analiz i diskretnaya matematika, Novosibirsk, Izd-vo Novosibirskogo Universiteta, 1989, pp. 19-27 (in Russian).

6. Lorenzi, A. and Paparoni, E. Direct and Inverse Problems in the Theory of Materials with Memory, The Mathematical Journal of the University of Padua, 1992, vol. 87, pp. 105-138.

7. Bukhgeym, A. L. Inverse Problems of Memory Reconstruction, Journal Inverse Ill-Posed Problems, 1993, vol. 1, no. 3, pp. 193-206.

8. Durdiev, D. K. A Multidimensional Inverse Problem for an Equation with Memory, Siberian Mathematical Journal, 1994, vol. 35, pp. 514-521. DOI: 10.1007/BF02104815.

9. Bukhgeim, A. L. and Dyatlov, G. V. Inverse Problems for Equations with Memory, SIAM J. Math. Fool., 1998, vol. 1, no. 2, pp. 1-17.

10. Durdiev, D. K. Obratnye zadachi dlya sred s posledejstviem, Tashkent, Turon-Ikbol, 2014.

11. Lorenzi, A., Ulekova, J. Sh. and Yakhno, V. G. An Inverse Problem in Viscoelasticity, Journal Inverse Ill-Posed Problems, 1994, vol. 2, no. 2, pp. 131-165. DOI: doi.org/10.1515/jiip.1994.2.2.131.

12. Janno, J. and Von Wolfersdorf, L. Inverse Problems for Identification of Memory Kernels in Viscoelasticity, Mathematical Methods in the Applied Sciences, 1997, vol. 20, no. 4, pp. 291-314. DOI: 10.1002/(SICI)1099-1476(19970310)20:4<291::AID-MMA860>3.0.C0;2-W.

13. Janno, J. and Von Wolfersdorf, L. An Inverse Problem for Identification of a Time- and Space-Dependent Memory Kernel in Viscoelasticity, Inverse Problems, 2001, vol. 17, no. 1, pp. 13-24. DOI: 10.1088/02665611/17/1/302.

14. Lorenzi, A., Messina, F. and Romanov, V. G. Recovering a Lame Kernel in a Viscoelastic System, Applicable Analysis, 2007, vol. 86, no. 11, pp. 1375-1395. DOI: 10.1080/00036810701675183.

15. Romanov, V. G. and Yamamoto M. Recovering a Lame Kernel in a Viscoelastic Equation by a Single Boundary Measurement, Applicable Analysis, 2010, vol. 89, no. 3, pp. 377-390. DOI: 10.1080/00036810701675183.

16. Lorenzi, A. and Romanov, V. G. Recovering Two Lame Kernels in a Viscoelastic System, Inverse Problems and Imaging, 2011, vol. 5, no. 2, pp. 431-464. DOI: 10.3934/ipi.2011.5.431.

17. Romanov, V. G. A Two-Dimensional Inverse Problem for the Viscoelasticity Equationi, Siberian Mathematical Journal, 2012, vol. 53, pp. 1128-1138. DOI: 10.1134/S0037446612060171.

18. Romanov, V. G. Inverse Problems for Differential Equations with Memory, Eurasian Journal of Mathematical and Computer Applications, 2014, vol. 2, no. 4, pp. 51-80

19. Durdiev, D. K. and Totieva, Z. D. The Problem Of Determining The One-Dimensional Kernel of the Electroviscoelasticity Equation, Siberian Mathematical Journal, 2017, vol. 58, no. 3, pp. 427-444. DOI: 10.1134/S0037446617030077.

20. Durdiev, D. K. and Rahmonov, A. A. Inverse Problem for a System of Integro-Differential Equations for SH Waves in a Visco-Elastic Porous Medium: Global Solvability, Theoretical and Mathematical Physics, 2018, vol. 195, pp. 923-937. DOI: 10.1134/S0040577918060090.

21. Durdiev, D. K. and Totieva, Zh. D. The Problem of Determining the Multidimensional Kernel of Viscoelasticity Equation, Vladikavkaz Mathematical Journal, vol. 17, no. 4, pp. 18-43 (in Russian). DOI: 10.23671/VNC.2015.4.5969.

22. Karchevsky, A. L. and Fatianov, A. G. Numerical Solution of the Inverse Problem for a System of Elasticity with the Aftereffect for a Vertically Inhomogeneous Medium, Sibirskii Zhurnal Vychislitel'noi Matematiki [Numerical Analysis and Applications], 2001, vol. 4, no. 3, pp. 259-268 (in Russian).

23. Durdiev, U. D. Numerical Method for Determining the Dependence of the Dielectric Permittivity on the Frequency in the Equation of Electrodynamics with Memory, Sibirskie Elektronnye Matematicheskie Izvestiya [Siberian Electronic Mathematical Reports], 2020, vol. 17, pp. 179-189 (in Russian). DOI: 10.33048/semi.2020.17.013.

24. Bozorov, Z. R. Numerical Determining a Memory Function of a Horizontally-Stratified Elastic Medium with Aftereffect, Eurasian Journal of Mathematical and Computer Applications, 2020, vol. 8, no. 2, pp. 4-16.

25. Kabanikhin, S. I., Karchevsky, A. L. and Lorenzi, A. Lavrent'ev Regularization of Solutions to Linear Integra-Differential Inverse Problems, Journal Inverse Ill-Posed Problems, 1993, vol. 1, no. 2, pp. 115140. DOI: 10.1515/jiip.1993.1.3.193.

26. Durdiev, D. K. and Bozorov, Z. R. A Problem of Determining the Kernel of Integrodifferential Wave Equation with Weak Horizontal Properties, Dal'nevostochnyi Matematicheskii Zhurnal, 2013, vol. 13, no. 2, pp. 209-221 (in Russian).

27. Totieva, Z. D. Determining The Kernel of the Viscoelasticity Equation in a Medium with Slightly Horizontal Homogeneity, Sibirskii Matematicheskii Zhurnal [Siberian Mathematical Journal], 2020, vol. 61, no. 2, pp. 359-378 (in Russian). DOI: 10.33048/smzh.2020.61.217.

28. Durdiev, D. K. and Totieva, Zh. D. The Problem of Determining the One-Dimensional Kernel of the Viscoelasticity Equation, Sibirskii Zhurnal Industrial'noi Matematiki [Journal of Applied and Industrial Mathematics], 2013, vol. 16, no. 2, pp. 72-82 (in Russian).

Received January 11, 2021 ZHANNA D. ^HEVA

Southern Mathematical Institute VSC RAS, 22 Marcus St., Vladikavkaz 362027, Russia, Senior Researcher of the Department of Math. Modeling E-mail: j annatuaeva@inbox. ru

i Надоели баннеры? Вы всегда можете отключить рекламу.