Численные методы решения некоторых краевых задач с обобщенными граничными условиями и их приложения к аэродинамике тема диссертации и автореферата по ВАК РФ 01.01.07, доктор физико-математических наук Сетуха, Алексей Викторович

• Сетуха, Алексей Викторович
• доктор физико-математических наук
• 2003

Оглавление диссертации доктор физико-математических наук Сетуха, Алексей Викторович

ВВЕДЕНИЕ

ГЛАВА 1. ОСНОВНЫЕ И ОБОБЩЕННЫЕ

ФУНКЦИИ. ОБОБЩЕННЫЕ КРАЕВЫЕ ЗНАЧЕНИЯ ФУНКЦИЙ НА КРИВЫХ И ПОВЕРХНОСТЯХ.

1.1 .Используемые обозначения

1.2.Пространства основных функций

1.3.Приближение дифференцируемых функций 45 многочленами

1.4.Пространства обобщенных функций

1.5.Продолжения функционалов и операторов

1.6. Действие функционалов на функции, 58 зависящие от параметра.

1.7.0бобщенные краевые значения функции на поверхности

1.8.0бобщенные краевые значения функций на плоских кривых

1,9.Первообразные обобщенных функций одного аргумента

ГЛАВА 2. ОБОБЩЕННЫЕ КРАЕВЫЕ

ЗНАЧЕНИЯ ГАРМОНИЧЕСКИХ ФУНКЦИЙ.

2.1. Краевые значения поверхностных 72 потенциалов.

2.2.Краевые значения потенциалов в плоском 85 случае

2.3.0 стирании особенностей гармонических функций

Введение диссертации (часть автореферата) на тему «Численные методы решения некоторых краевых задач с обобщенными граничными условиями и их приложения к аэродинамике»

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

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

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

Основы аэродинамической теории несущей поверхности были разработаны в трудах Н.Е. Жуковского и С.А. Чаплыгина [17,60]. В плоском случае задача об обтекании тонкого профиля крыла потоком идеальной несжимаемой жидкости заключается в нахождении векторного поля у(дг) = (у,(д:),у2(д:)),дг = (дг,,д:2) (поля скоростей), определенного на плоскости

Я2 вне разомкнутой кривой Ь, задающей форму профиля. При этом должны выполняться уравнения [23 с.113,с. 130], [36] ОУ, Эу2 . „ ¿У2 ¿У, о2.г (1)

Луу =—!- + —¿- = 0, —1--! = (), хеЯг\Ь, дхх дх2 дх1 дх2 ставится условие на бесконечности

Нт у(*) = Яо > (2) где и^ есть заданный вектор скорости набегающего потока, и в каждой точке контура Ь ставится граничное условие, выражающее отсутствие потока жидкости через поверхность профиля: у±-Я = О, п - вектор нормали на контуре Ь, V* - краевые значения векторной функции у(х) в точках контура со стороны вектора п и с противоположной стороны соответственно.

Рис.1.

В работах Н.Е. Жуковского и С.А. Чаплыгина была выдвинута идея о том, что при решении задачи о стационарном обтекании профиля следует искать решение записанной задачи - векторное поле у(д:) , которое ограничено в окрестности задней кромки профиля (точки В) и может быть неограниченным в окрестности передней кромки (точки А).

В первой половине XX века были разработаны методы аналитического решения плоских задач аэродинамики идеальной жидкости на основе теории функций комплексного переменного [15,41,46]. Однако, эти методы не позволяют выписать поле скоростей в виде, пригодном для количественных расчетов, в случае контура Ь общего вида.

В работах Н.Е.Жуковского, С.А. Чаплыгина и Л. Прандтля (Ь.РгапсШ) была выдвинута идея о том, что несущую поверхность можно моделировать «вихревыми нитями»[ 18,43,61] . Применительно к плоской задаче это означает, что возмущенное поле скоростей й> = V — н^ можно представить в виде суперпозиции особенностей типа вихрь, распределенных вдоль поверхности профиля и-(*) = ¡г^Щх-у^^, хе112\Ь, (3) ь где 5 - естественный параметр на кривой Ь, у^) - точки кривойЬ, - неизвестная плотность вихревого слоя, размещенного на кривой Ь, V есть векторное поле, определяемое законом Био-Савара (поле скоростей, индуцируемое точечным вихрем)

В 50-х годах XX- века С.М.Белоцерковским был предложен численный метод приближенного решения рассматриваемой задачи, названный методом дискретных вихрей [4,5]. В этом методе контур Л, аппроксимируется системой точечных вихрей, которые размещаются в точках л:', / = 1равномерно распределенных по длине контура, и приближенное поле скоростей ищется в виде ад=¿г, ГС*-*'), ы\ где Г, ,/ = 1 \.п - неизвестные циркуляции точечных вихрей. Для нахождения неизвестных Г, ,/ = 1,.« записывается система линейных алгебраических уравнений, аппроксимирующих граничное условие в специальным образом выбранных точках коллокации (контрольных точках). Математическое обоснование данного метода было получено И.К. Лифановым [27]. Он показал, что нахождение неизвестной функции у сводится к решению сингулярного интегрального уравнения с ядром Коши, а метод дискретных вихрей является по существу методом приближенного решения такого уравнения. В его работах метод был обобщен на различные виды краевых задач. Были разработаны подходы, позволяющие сводить плоские краевые задачи Неймана и Дирихле для уравнений Лапласа и Тельмгольца к сингулярным интегральным уравнениям, причем, как для внутренних и внешних задач в области с замкнутой границей, так и для задач в области вне разомкнутой кривой. Для возникающих сингулярных интегральных уравнений разработаны и обоснованы методы численного решения, являющиеся развитием метода дискретных вихрей [29,30,6,32,9].

В 80-90-х годах был разработан метод дискретных вихревых рамок, позволяющий решать трехмерные задачи аэродинамики ([3], см. также [32 с.473-477], [9 с. 439-448]). В этом методе поверхность обтекаемого тела разбивается на ячейки четырехугольной или треугольной формы и по контуру каждой ячейки размещается вихревая нить неизвестной интенсивности. При этом поле скоростей ищется в виде суперпозиции скорости набегающего потока и скоростей, индуцируемых вихревыми рамками в соответствии с законом Био-Савара. Для нахождения неизвестных циркуляций вихревых рамок на каждой рамке определенным образом выбирается точка коллокации и записывается граничное условие равенства нулю нормальной составляющей скорости. Так же, как и в плоском случае, метод сначала был развит на основе эмпирических соображений, а затем под него была подведена теоретическая база.

Трехмерная задача о потенциальном обтекании тела идеальной несжимаемой жидкостью заключается в нахождении векторного поля у(дг), удовлетворяющего всюду вне тела, в области, занимаемой жидкостью, уравнениям [23 с.359], [36] т? = 0, гоЯ = 0, (4) условию (2) на бесконечности и граничному условию уй = 0 на поверхности тела (если тело является тонкой разомкнутой поверхностью, граничное условие должно выполняться для краевых значений с обеих сторон этой поверхности). Удобно перейти к отысканию возмущенного поля скоростей = у-ху^ . Пусть О - область в пространстве /?3, в которой ищется векторное поле й\ В диссертации рассматривается случай, когда обтекаемое тело представляет из себя тонкую разомкнутую поверхность Е, причем, эта поверхность лежит на некоторой плоскости и является открытым множеством на этой плоскости. При этом О = Я3 \Е, где Е есть замыкание множества Е в пространстве Л . Если область О является односвязной, то векторное поле уР является потенциальными для его потенциала и возникает краевая задача Неймана

Аи = 0 в области П, (5) ди , „ (6) = / на поверхностиЕ, дп где / = -н'00Я (граничное условие (6) ставится на обеих сторонах поверхности Е ), причем, функцию и можно искать так, чтобы выполнялось условие

Нш и(х) = 0. (7)

Эффективным методом решения краевых задач для эллиптических уравнений является метод граничных интегральных уравнений, основанный на теории потенциала. В классической теории потенциала обычно рассматриваются краевые задачи в области, границей которой является замкнутая поверхность [12,24,57,58]. Решение краевой задачи Неймана для уравнения Лапласа ищется в виде поверхностного потенциала простого слоя, размещенного на границе области. При этом задача сводится к уравнению Фредгольма второго рода для неизвестной плотности потенциала простого слоя. Однако в случае, когда граница области является разомкнутой (т.е. в случае экранной задачи), решение задачи Неймана нельзя искать в виде потенциала простого слоя, т.к. нормальная производная потенциала простого слоя претерпевает на поверхности, где он размещен, скачок и поэтому граничное условие не может выполняться одновременно на обеих сторонах поверхности.

В работах И.К.Лифанова предложен подход к решению трехмерных краевых задач Неймана, в котором решение ищется в виде потенциала двойного слоя [28,32]. В этом случае для неизвестной плотности потенциала двойного слоя возникают гиперсингулярные интегральные уравнения, причем, такие уравнения были выписаны как для внутренних и внешних задач в области с замкнутой границей, так и для экранной задачи. И.КЛифанов и Л.Н.Полтавский показали, что метод дискретных вихревых рамок по существу является методом решения гиперсингулярных интегральных уравнений [42,9]. В этих работах была исследована сходимость квадратурных сумм типа метода дискретных вихревых рамок для вычисления гиперсингулярных интегралов в случае, когда область интегрирования есть гладкая разомкнутая поверхность. В случае, когда поверхность Е является частью плоскости, Л.Н. Полтавский доказал сходимость численных решений гиперсингулярного интегрального уравнения, получаемых с использованием метода дискретных вихревых рамок, к точному решению. Этим исследованиям посвящен цикл работ [42,31,33,34]. Систематическое изложение полученных в этих работах результатов можно найти в [9]. Опишем основные результаты указанных работ.

Рассмотрим краевую задачу Неймана (5)-(7) в случае, когда О есть множество точек трехмерного арифметического пространства Л3, лежащих вне разомкнутой поверхности £, которая является частью координатной плоскости Оххх2. При этом предположим, что £ есть выпуклое открытое ограниченное множество на плоскости, границей которого является бесконечно гладкая кривая дИ. Решение такой задачи ищем в виде потенциала двойного слоя л:еО, где д!дпу есть производная по направлению вектора Я = (0,0,1), перпендикулярного к поверхности вычисляемая по координатам точки у, интеграл понимается как поверхностный интеграл первого рода (по площади поверхности), V - неизвестная плотность потенциала двойного слоя. Запись граничного условия (6) приводит к интегральному уравнению для функции V (в работе [28] это уравнение названо уравнением Прандтля)

Т~ Ку), 1 ,3 = Ах), х е I, ( ^ где интеграл понимается в смысле конечного значения по Адамару [2] пг ,, , 1 Му)с1у 1 . | Г у(у)4у 2тСУ(х)\ ^

5(е,х)- открытый круг на плоскости радиуса б с центром в точке х (см. глава 2, формула (2.10)).

Вычислительная схема для приближенного решения уравнения (8) строится следующим образом. Плоскость, на которой лежит поверхность £ разбивается на квадратные ячейки со стороной к. Пусть ак, к = 1,.,«, есть ячейки, целиком лежащие в множестве Приближенное решение ищется в виде функции И, заданной на множестве £л = ак и принимающей на каждой ячейке постоянное значение ук. В центре каждой ячейки разбиения ак берется точка коллокации х* =(^,^2), к-\,.,п и записывается система линейных алгебраических уравнений.

V /■ • 1 (10) диа[е\ 1 . 1 где-/=/(*'), —(*) = —]:-¡г^сг / = = 1,[е] дп есть потенциал двойного слоя с плотностью равной 1, размещенный на ячейке сгу. При 7 = у интеграл понимается в смысле (9), причем, в п.5.2.2 показано, как выражение для коэффициентов аи, / = 1,.,«, преобразовать к обычному несобственному интегралу.) При этом приближенное решение краевой задачи Неймана (5)-(7) представляется в виде

ПП л,

Отметим, также, что градиент функции и (х) записывается в виде [23] 1 причем, для каждого у = 1,.,« функция и^Дл:) может быть представлена как векторное поле, индуцируемое вихревой нитью, размещенной по контуру ячейки сту, которое определяется законом Био-Савара. В случае рамки прямоугольной формы векторное поле йДл:) может быть вычислено аналитически [32].

Л.Н. Полтавский доказал (см. [42,9]), что если V есть заданная функция на множестве £, представляющаяся в виде г(г) у\Х) (П)

ТрС^аЁ)' где V* имеет на множестве 2 первые производные, непрерывные по Гельдеру, р(х,д1.) - расстояние от точки х до края поверхности £, то выполнена оценка

Е1-Г 7=1 л/р(*у,д£) где С и ¡л есть некоторые положительные константы.

Обоснование сходимости численных решений уравнения (8), получаемых на основании уравнений (10) впервые было получено в докторской диссертации Л.Н. Полтавского [42], в предположении, что уравнение (8) имеет решение, представляющееся в виде (13). При этом была доказана сходимость приближенных решений к точному в среднем, а также равномерная сходимость приближенных решений к точному и разностных отношений, получаемых при использовании приближенного решения, к производным точного решения первого порядка на любом замкнутом подмножестве внутренности множества £ (здесь множество £ рассматривается как множество на плоскости). Одним из результатов, полученных автором настоящей диссертации, явилось доказательство равномерной сходимости приближенных решений к точному на всем множестве £ (см. п.5.2.2 настоящей диссертации). Опираясь на этот результат, Л.Н. Полтавский получил оценки для численных значений производной неизвестной функции, справедливые вплоть до границы множества £ [9]. (Нахождение приближенных значений производных функции у необходимо в задачах аэродинамики для расчета сил давления, действующих на поверхность).

Новый класс постановок краевых задач и граничных интегральных уравнений возник при моделировании обтекания тел потоком идеальной несжимаемой жидкости при наличии отсоса внешнего потока [7]. Так при моделировании плоского обтекания тонкого профиля с отсосом потока, осуществляемым на одной стороне профиля в заданной точке <? = (<?1,<?2)е^' на профиле, возникает задача о нахождении векторного поля скоростей, имеющего в окрестности точки д со стороны, где осуществляется отсос потока, особенность типа сток. С другой стороны профиля поле скоростей должно быть непрерывным (рис. 2).

В работе В.И. Бушуева и И.К. Лифанова [7] предложено искать такое поле скоростей в виде V = + й^ + й», где й>да есть скорость потока на бесконечности, - векторное поле, индуцируемое особенностью типа сток: *2> x = {xvx2)eR2,xФq, - заданная константа, выражающая объемный расход отсасываемой жидкости за единицу времени, й есть векторное поле, индуцируемое вихревым слоем неизвестной интенсивности у(з), которое определяется выражением (3), причем, функция у должна представляться в виде

2 1 . О*)

К 5 где - значение естественного параметра на контуре Ь, соответствующее точке <7, у* - функция, непрерывная по Гельдеру в окрестности точки(предполагается, что отсос потока осуществляется со стороны вектора нормали Я и пара векторов (г,Я), где т = ¿у^я)/ей, является правой). При этом доказано [35], что в окрестности точки q со стороны, противоположной вектору Я, поле V является непрерывным, а со стороны вектора Я представляется в виде

-/ \ 0. -*, х где \>*(л:) есть непрерывная функция в окрестности точки ц со стороны вектора Я . Задача сводится к сингулярному интегральному уравнению с ядром Коши на отрезке для функции у ГК(8,50)у{8)аз = /(50 ), (15) которое должно выполняться во всех точках е Ь, кроме концов кривой Ь и точки <7, правая часть уравнения и ядро есть некоторые заданные гладкие функции, причем ищется решение записанного уравнения, имеющее вид (14). Такие решения были названы сингулярными. Для приближенного нахождения сингулярных решений была предложена численная схема типа метода дискретных вихрей и доказана ее сходимость (обоснование сходимости см.[32 с.350-354]). При этом нахождение приближенных значений функции у в узлах сводится к решению алгебраической системы линейных уравнений, в которых заложена ассимптотика функции у в окрестности точки ^.

В дальнейших работах, выполненных Лифановым И.К. с рядом соавторов [65-69,35], были рассмотрены различные варианты постановок краевых задач и сингулярных уравнений с сингулярными решениями. Так оказалось, что если интенсивность стока считать неизвестной, то можно поставить краевую задачу об отыскании векторного поля V, ограниченного в окрестностях обоих концов кривой £. С точки зрения аэродинамики организация такого течения с безударным обтеканием передней кромки представляет значительный интерес, поскольку за счет этого имеется возможность предотвратить отрыв потока и потерю несущих свойств профиля на больших углах атаки. Систематическое изложение различных постановок краевых задач и численные схемы их решения приведено в монографии [32].

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

В предлагаемой диссертации развит новый подход к постановке и решению задач указанного типа. Основная идея этого подхода заключается в том, чтобы рассматривать краевые значения неизвестных функций и их нормальных производных как обобщенные функции. Так в задаче об обтекании профиля с отсосом потока в точке нормальная составляющая вектора скорости в каждой точке профиля, отличной от точки <7, должна быть равна нулю. В то же время суммарный поток вектора скорости через контур Ь со стороны вектора п должен быть равен -(). Это наводит на мысль поставить со стороны вектора Я краевое условие vn = где 8 есть дельта функция Дирака. Для того, чтобы придать записанному равенству строгий математический смысл, в диссертации вводится понятие обобщенных краевых значений и обобщенных нормальных производных. Дается постановка плоской краевой задачи для векторной функции, удовлетворяющей уравнениям (1), и краевой задачи Неймана (5)-(7) в пространственном случае при условии, что правая часть в граничном условии есть обобщенная функция. Поставленные краевые задачи сводятся к граничным интегральным уравнениям, причем, в плоском случае возникает сингулярное интегральное уравнение (15), а в пространственном случае гиперсингулярное интегральное уравнение (8). Для записанных уравнений вводится понятие обобщенного решения. Рассмотрены вопросы разрешимости поставленных краевых задач и сингулярных интегральных уравнений. Предложены численные схемы решения записанных интегральных уравнений в случае, когда правая часть есть обобщенная функция некоторого частного вида и доказана сходимость приближенных решений к точному в некоторой слабой топологии. При этом для решения соответствующих краевых задач доказана сходимость приближенных решений и их частных производных всех порядков к точным значениям во всех внутренних точках области определения.

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

Способ постановки краевых задач с обобщенными граничными условиями и методы теоретического исследования разрешимости этих постановок отличаются от известных ранее. Возможный способ такой постановки краевых задач для эллиптических уравнений был предложен в работе М.И.Вишика и С.Л.Соболева [11]. В этой работе рассматривалась внутренняя краевая задача для уравнения Пуассона. При этом неизвестная функция рассматривалась как обобщенная функция, заданная во всем пространстве и равная нулю вне замыкания области, в которой ищется решение. Левые части уравнения и граничного условия рассматривались как единый оператор над такой обобщенной функцией. Однако такой подход не применим к экранным задачам.

В настоящее время достаточно широко, также, развиты постановки и теория краевых задач, в которых решение является обобщенной функцией в области, а на границе краевые значения неизвестной функции и ее производных лежат в пространствах Соболева интегрируемых функций Wrp, р e[l,oo),r е Z+

55,38]. Но эти пространства не являются достаточно широкими для применения в задачах аэродинамики с отсосом потока ( так эти пространства не содержат 8 -функцию).

Другой подход к введению обобщенных краевых значений функций и исследованию поведения функции у границы области основан на представлении функции через интегралы по границе. Плоские краевые задачи для эллиптических уравнений тесно связаны с краевыми задачами для аналитических функций комплексного переменного, которые, в свою очередь, сводятся к сингулярным интегральным уравнениям. Обычно сингулярные интегральные уравнения рассматриваются в классах функций, непрерывных по Гельдеру [40,13]. Интегрируемые по Лебегу краевые значения аналитических функций комплексного переменного и плоских гармонических функций, возникли при исследовании свойств интеграла Коши с плотностью, интегрируемой по Лебегу [70,19,10,44,59]. В работах [53,54,16] рассматривались краевые задачи Римана и Гильберта для аналитических функций с измеримыми коэффициентами. К этим вопросам примыкают исследования по теории сингулярных интегралов и сингулярных интегральных уравнений в классах измеримых функций [39, 37]. Эти результаты могут быть использованы при постановке и исследовании двумерных краевых задач.

Дальнейшее развитие обобщенных постановок граничных условий оказалось возможным путем рассмотрения граничных интегральных уравнений в классах обобщенных функций. При этом широкое распространение получила теория псевдо-дифференциапьных операторов, основанная на преобразовании Фурье [1,62,64,63]. В работах Л.Н. Полтавского (см.[42,31,33,34,9]) теория псвдодифференциапьных операторов была применена к исследованию разрешимости гиперсингулярного интегрального уравнения (8), в котором правая часть и неизвестная функция рассматривались в пространствах Соболева-Слободецкого. Было получено и обоснование численного метода типа метода дискретных вихревых рамок для отыскания таких решений. Однако порядок обобщенных функций в правой части уравнения (8), для которых были получены указанные результаты, не достаточен для приложения к описанным выше задачам аэродинамики. В частности случай, когда в правой части стоит дельта-функция, этими работами не охватывается.

В работах Лифанова И.К. и Вайникко Г.М. рассмотрено сингулярное интегральное уравнение с ядром Гильберта в случае, когда правая часть в этом уравнении есть дельта функция см. [8,9]. Для исследования такого уравнения применена теория псевдо-дифференциальных операторов и показано, что решение уравнения имеет неинтегрируемую особенность вида (14). Такое уравнение рассматривалось в связи с приложением к плоским краевым задачам аэродинамики, описанным выше. Однако трехмерные краевые задачи при этом не рассматривались.

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

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

В плоском случае при исследовании обобщенных решений сингулярного интегрального уравнения с ядром Коши на отрезке использован подход, основанный на сведении этого уравнения к краевой задаче Римана для аналитических функций. Для уравнений, рассматриваемых в классах функций, непрерывных по Гельдеру, этот подход описан в монографиях Н.И. Мусхелишвили и Ф.Д. Гахова [40,13]. В настоящей диссертации введено понятие краевой задачи Римана с обобщенными граничными условиями и получены аналитические представления для решений такой задачи. Затем на основе этих представлений получены формулы для обращения сингулярного интегрального уравнения на отрезке с обобщенной правой частью.

Содержание диссертации.

По структуре диссертация состоит из введения, 5 глав и заключения.

Заключение диссертации по теме «Вычислительная математика», Сетуха, Алексей Викторович

ЗАКЛЮЧЕНИЕ

В диссертации получены следующие основные результаты, которые выносятся на защиту:

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

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

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

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

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

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

Список литературы диссертационного исследования доктор физико-математических наук Сетуха, Алексей Викторович, 2003 год

1. Агранович М.С. Эллиптические сингулярные интегро-дифффернциалъные операторы. И Успехи математических наук, т.20, №5. 1965г.- с.3-120.

2. Адамар Ж. Задача Коши для; линейных уравнений с; частными производными гиперболического типа. М.Наука, 1978. - 353с.

3. Апаринов В.А., Белоцерковский С.М., Лифанов И.К., Михайлов A.A. -Расчет нестационарных аэродинамических характеристик тел при отрывном обтекании. -ЖВМиМФ, 1988. -т.2, №16. -с.1558-1566.

4. Белоцерковский С.М. Исследования по аэродинамике современных несущих поверхностей. Дис.д-ра.тех.наук. М.1955.

5. Белоцерковсеий С.М; Тонкая несущая поверхность в дозвуковом потоке газа. М.Наука. 1965. 244с.

6. Белоцерковский С.М., Лифанов ПК. Численные методы в сингулярных инетегральных уравнениях. М.Наука.- 1985. — 256с.

7. Бушуев.В.И., Лифанов И.К. Численное решение сингулярных интегральных уравнений в классе сингулярных функций и задача отсоса потока в аэродинамике. //ЖВМ иМФ. 1986. -№10. 1572-1577.

8. Вайникко Г.Е., Лифанов И.К. Моделирование задач аэродинамики и дифракции волн и расширение интегральных операторов типа Коши на замкнутых и открытых кривых. II Дифференциальные уравнения. 2000г. Т.36, №9.- с. 1184-1195.

9. Вайникко Г.М., Лифанов И.К., Полтавский Л.Н. Численные методы в гиперсингулярных интегральных уравнениях и их приложения. М.: "Янус-К", 2001,-508с.

10. Ю.Векуа И.Н. Обобщенные аналитические функции. М.: Физматгиз, 1959.

11. И.Вишик М.И., Соболев СЛ. Общая постановка некоторых краевых задач для эллиптических дифференциальных уравнений в частных производных. II ДАН СССР, 1956., т. 111, №3, - с. 521-523.

12. Владимиров B.C. Уравнения математической физики. — М.Наука. 1976 -528с.

13. Гахов Ф.Д. Краевые задачи. М.: Наука, 1977. - 640с.

14. М.Гельфанд И.М., Шилов Г.Е. Обобщенные функции и действия над ними. М.: Гостехиздат. 1958. - 440с.

15. Голубев В.В. Теория крыла аэроплана в плоско параллельном потоке. -М.: ОНТИ, 1936.

16. Данилюк W.W. О задаче Гильберта с измеримыми коэффициентами. II Сиб. мат. журнал, т.Г, №2, 1960, с. 171-197.

17. Жуковский Н.Е. Теоретические основы воздухоплавания. Собрание сочинений. Т.6, Гостехиздат, Москва , 1948.

18. Колтон Д., Кресс Р. Методы интегральных уравнений в теории рассеяния. М.:Мир, 1987. - 311с.

19. Кочин Н.Е. Кибель И.А. Розе Н.В. Теоретическая гидромеханика. 4il. -М:Физматгиз. 1963г. 584 с.

20. Курант Р. Уравнения с частными производными. М.Мир, 1964. - 882с.

21. Куратовский К. Топология. Т. 1. М.Мир, 1966, 596 с.

22. Лаврентьев М.А., Шабат B.BL Методы теории функций комплексного переменного. М.Наука. 1987.

23. Лифанов И.К., Полонский Я.Е. Обоснование численного метода дискретных вихрей решения сингулярных интегральных уравнений II ПММ. 1975, т39, №4 - с.742-746.

24. Лифанов И.К. О методе дискретных вихрей для крыла бесконечного размаха и уравнении Прандтля для крыла конечного размаха. II Известия вузов. Математика. -1980. №6. - с.44-51.

25. Лифанов И.К. О численном решении сингулярных интегральных уравнений // Дифференц. уравнения. 1981., т17, №12.

26. Лифанов И.К. Сингулярное интегральное уравнение первого рода задачи Неймана. II Дифференц. уравнения. 1988., т24, №1. -с.110-115.

27. Лифанов И.К. Полтавский Л.Н. Обобщенный оператор Фурье и его применение в обосновании метода дискретных вихрей. II Математический сборник. 1992., т 183., №5 -с.79-114.

28. Лифанов И.К. Метод сингулярных интегральных уравнений и численный эксперимент. М.:ТОО "Янус", - 1995. - 520с.

29. Лифанов И.К. Полтавский Л.Н. Пространства дробных отношений, дискретные операторы и их приложения. I. II Математический сборник. 1999., т190., №9 -с.41-98.

30. Лифанов И.К. Полтавский Л.Н. Пространства дробных отношений, дискретные операторы и их приложения: II. II Математический сборник. 1999., т190., №11 -с.67-134.

31. Лифанов И.К., Сетуха A.B. О сингулярных решениях некоторых краевых задач и сингулярных уравнений. // Дифференциальные уравнения т.35, №9, 1999г с. 1227-1241.

32. Лойцянский Л.Г. Механика жидкости и газа. М;: Наука, 1978 г. - 736с.

33. Магнарадзе Л.Г. Об одном обобщении теоремы Племели-Привалова. Н Сообщ. АН Груз. ССР. т.8, 1947г. с.509-516.

34. Михайлов. Дифференциальные уравнения в частных производных. -М. Наука. 1983.

35. Михлин С.Г. Сингулярные интегральные уравнения. И УМН III, вып.3(25), l948.-c.19-l 12.

36. Мусхелишвили Н.И. Сингулярные интегральные уравнения. М.: Наука. 1968., 512 с.

37. Некрасов А.И. Теория крыла в нестационарном потоке газа. М.: изд-во АН СССР, 1947.

38. Полтавский Л.Н. Математическое обоснование некоторых численных схем в аэродинамике. Диссертация д-ра физ. мат. наук. М. 1991.

39. Прандтль Л. Теория несущего крыла. Перевод, с немецкого. М.: 1931.

40. Привалов И.И. Граничные свойства аналитических функций. -Гостехиздат. 1950.

41. Привалов И.И. Введение в теорию функций комплексного переменного. М.: Наука. 1984ю - 432с.

42. Седов Л.И. Плоские задачи гидродинамики и аэродинамики. Изд.2-е, М.: Гостехиздат, 1950.

43. Сетуха A.B. Краевая задача Неймана в полупространстве. // Дифференциальные уравнения, т.36, №9, 2000 г. с. 1172-1185.

44. Сетуха A.B. Краевая задача Неймана в с граничным условием на плоской разомкнутой поверхности. // Дифференциальные уравнения. т.37, №10, 2001г. с.1311-1329.

45. Сетуха А.В: О построении фундаментальных решений краевой задачи Неймана в области вне разомкнутой плоской поверхности. // Дифференциальные уравнения т.З 8, №4, 2002г. с505-515.

46. Сетуха A.B. О плоской краевой задаче Неймана с обобщенными граничными условиями. // Дифференциальные уравнения т.З 8, №9, 2002г . с. 1172-1182.

47. Сетуха A.B. Фундаментальные решения плоской краевой задачи Неймана для уравнения Лапласа. // Дифференциальные уравнения т.39, 2003 г. (сдана в печать).

48. Сетуха A.B. Трехмерная краевая задача Неймана с обобщенными граничными условиями и уравнение Прандтля. // Дифференциальные уравнения т.39, №9, 2003 г.

49. Симоненко И.Б. Краевая задача Римана с измеримым коэффициентом. II ДАН СССР, т.135, №3, 1960. с.538-541.

50. Симоненко И.Б. Краевая задача Римана для п пар функций с измеримыми коэффициентами и ее применение к исследованию сингулярных интегралов в пространстве Lp с весами. II Изв. АН СССР, сер. Матем. т.28, №2, 1964. с.277-306.

51. Соболев С.JI. Methode nouvelle a resoudre le probleme de Cauchy pour les equations linears hyperboliques normales. II Мат. сб. т. №1(43), 1936г. -с. 39-72.

52. Соболев С.Л. Некоторые применения математического анализа в математической физике. М.: Наука. 1988г.

53. Соболев С.Л. Уравнения математической физики. М.:Наука. 1992, -432с.

54. Тихонов А.Н. Самарский А.А. Уравнения математической физики. -Учебное пособие. 6-е изд., М.: Изд-во МГУ, 1999г. 798 с.

55. Чаплыгин. С.А. О давлении плоско параллельного потока на преграждающие тела. Собрание сочинений. Т.2, Гостехиздат, Москва , 1948.

56. Шилов Г.Е. Математический анализ. Второй специальный курс. М., Наука. 1965.

57. Шубин M.A. Псевдодифференциальные операторы и спектральная теория. М.Наука., 1978.

58. Эскин Г.И. Краевые задачи для эллиптических псевдодифференциальных уравнений. М., Наука. 1973,232с.

59. Gulyaev V.V., Lifanov I.K., Misko V.A. Mathematical model of the flow around airfoil with energy hight-lift devicis. II Russ J. Numer. Anal. Math. Modelling. 1996. V.9. №3. - p237-251.

60. Lifanov I.KSingular solutions of singular integral equations end flow ejecting for an arbitrary contur. И Sov. J. Numer. Anal. Math. Modelling. -1989. V.4. №3. — p239-252.

61. Lifanov I.K., Mikhailov A.A. Titsky S.V. Mathematical modelling of airfoil flowing control by ejection. II Sov. J. Numer. Anal. Math. Modelling. 1990. V.5. №3. - p209-220.

62. Lifanov I.K., Mikhailov A.A. Titsky S.V. External flow entrainment and boundary layer separation control for thick airfoil I I Sov. J. Numer. Anal. Math. Modelling. 1991. V.5. №4. -p325-334.

63. Lifanov I.K., Panarin L.N., Setukha A.V. Mathematical modeling of the rotor of a helicopter with ejection devices. II Russ. J. Numer.Anal. Math. Modelling. Vol.14, No 3, p.237-264.

64. Reisz M. Sur les fonctions conjugates II Math. Zeit. т.27., c.218-244. 1927.

65. Schwartz L. Theorie des distributions I,II. Paris., 1950-51.