автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическое моделирование обтекания профилей с отсосом и численное решение сингулярных интегральных уравнений в классе обобщенных функций

ВОЕННО-ВОЗДУШНАЯ ИНЖЕНЕРНАЯ АКАДЕМИЯ имени профессора Н.Е. Жуковского

На правах рукописи

Лебедева Нина Васильевна

УДК 517

Математическое моделирование обтекания профилей с отсосом и численное решение сингулярных интегральных уравнений в классе обобщенных функций

05.13.18 - математическое моделирование, численные методы и комплексы программ

Автореферат диссертации на соискание учёной степени кандидата физико-математических наук

Москва - 2005

амо

ВОЕННО-ВОЗДУШНАЯ ИНЖЕНЕРНАЯ АКАДЕМИЯ имени профессора Н.Е. Жуковского

На правах рукописи

Лебедева Нина Васильевна

УДК 517

Математическое моделирование обтекания профилей с отсосом и численное решение сингулярных интегральных уравнений в классе обобщенных функций

05.13.18 - математическое моделирование, численные методы и комплексы программ

Автореферат диссертации на соискание учёной степени кандидата физико-математических наук

Москва , ^,

2101873

Работа выполнена на кафедре высшей математики Военно-технического университета при Федеральном агентстве специального строительства

Научный руководитель:

доктор физико-математических наук, профессор И.К. Лифанов Официальные оппоненты:

доктор физико-математических наук,

профессор

Е.В. Захаров

кандидат физико-математических наук, доцент

Крутицкий П.В. Ведущая организация:

ЦАГИ им. Н.Е. Жуковского.

Защита состоится ¿^/Я>Р 2005 года в

_часов на заседании диссертационного совета Д

215.001.01 Военно-воздушной инженерной академии имени профессора Н.Е. Жуковского по'адресу: 125190, г. Москва, ул. Планетная, д.З

С диссертацией можно ознакомиться в библиотеке Военно-воздушной инженерной академии имени профессора Н.Е. Жуковского.

Автореферат разослан «__»

2005

года.

Учёный секретарь диссертационного с кандидат физико-математических нау

С. Ненашев

«лыкжяии

$ А вСР й МДОНЛЛЫКЯ»|1 е р о в л н о 22 страницы БИБЛИОТЕКА СП«

09 "V---'М

Обицая характеристика работы.

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

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

Целью работы является :

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

2. Доказательство сходимости приближенных решений, получаемых по методу дискретных вихрей к точным а) для сингулярного и гиперсингулярного интегральных уравнений на отрезке с ядром Коши; б) для сингулярного и гиперсигулярного интегральных уравнений в классе периодических функций.

Научная новизна и теоретическая ценность работы состоят в следующем :

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

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

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

Результаты, полученные в этой работе, показывают, что математические модели, исследуемые в данной работе, позволяют

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

Практическая значимость работы.

Построенные модели возникли из актуальных задач практики. С помощью математических доказательств

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

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

Апробация работы.

Работа в целом докладывалась на заседании научного семинара по интегральным уравнениям на ВМиК МГУ (рук. Е.В. Захаров, И.К. Лифанов). Частичные результаты диссертационной работы в виде докладов были представлены на X и XI Международных симпозиумах «Методы дискретных

особенностей в задачах математической физики» (Херсон, МДОЗМФ-2001, МДОЗМФ-200Э), а также в сборнике «Труды Международных школ - семинаров МДОЗМФ» (0рел-2002).

Основные результаты диссертационной работы были опубликованы в статье Вайникко Г.М., Лебедева Н.В., Лифанов И.К. «Численное решение сингулярных и гиперсингулярных интегральных уравнений на отрезке и дельта - функция », Матсборник, 2002, т.193, №10, с.3-16, а также в статве Lebedeva N.V. and Lifanov I.K. "Singular and hipervsingular integral equations with the Hilbert kernel, delta - function, and metod of discrete vortices" Proc. of the Estonian Acad, of sciences physics, mathematics, - 2004, т.53, №2, с.84-91.

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

Структура и объем работы. Диссертационная работа состоит из введения, четырех глав, заключения, списка использованной литературы. Общий объем работы составляет 77 страниц. Библиография составляет 56 наименований.

Основное содержание работы.

Во введении обоснована актуальность темы диссертации, приведен обзор литературы по теме исследования, указана цель работы и её новизна, показана теоретическая значимость результатов исследований. Приведено краткое изложение основного содержания диссертационной работы.

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

Рассматривается плоскопараллельное обтекание изолированного непротекаемого профиля установившимся потоком идеальной несжимаемой жидкости, скорость которого й0 =и0„1+и0у]. Сам профиль является неподвижным. Под профилем понимаем цилиндрическую поверхность с образующей, параллельной оси ОХ, и направляющей кривой I в плоскости ОХУ (рис.1). Так как все параметры набегающего потока и возмущенного не зависят от координаты г, то в дальнейшем профиль изображается только кривой I в плоскости ОХУ. Если кривая ¿ является простой разомкнутой (рис. 1а), то такой профиль называется тонким, а если кривая £ является простой замкнутой кусочно - гладкой (рис. 16), то такой профиль называется телесным. Поскольку рассматривается стационарная задача, то следа за профилем нет, и возмущенное течение будет потенциальным везде вне профиле

а)

б)

Рис.1 Разомкнутый (а) и замкнутый (б) (телесный) профили

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

± р, (Л ' У) + <■(*„ - *) = ДМ х (1)

2яг I гОш,

где /(М0) = -й„пк, г(м)- интенсивность распределенного вихревого слоя в точке Л/е£ и /(м0) гладкая функция, Дл/0)еЯ(а) на I.

Далее рассматриваются сингулярные интегральные уравнения каждого типа, к которым сводится уравнение (1).

Если Ь - гладкая разомкнутая кривая, заданная параметрически х = *(/), у = <е[-1Д],то уравнение (1) может быть записано в виде

е[-1,1], (2)

'о~' -I

где /(/)е Я(а) на [-1,1], *,(/„,/)е Я(а) на [-иМ-М]. Это уравнение имеет решения трех типов:

= , где Ьдод-о, (4)

VI-/ £

(5)

где функция (^(0 не имеет особенностей на Ь, соответствующие циркуляционному, бесциркуляционному и безударному обтеканию профиля I соответственно.

Если профиль Ь является телесным без острой кромки, т.е. кривая £ является простой замкнутой гладкой, и параметр / изменяется на отрезке [0,2л-], то уравнение (1) можно записать в виде

+- до, /„€[0,2*] (6)

™ о * о

Так как набегающий поток является потенциальным во всей плоскости, то для любого телесного профиля выполняется равенство

]/(ЛОА = 0 (7)

Для телесного профиля решение уравнений (1) определенно с точностью до константы, т.е. в этом уравнении константа является собственным решением.

Далее рассматривается задача с отсосом внешнего потока. Пусть в точке Мч(хч,уч) профиля Ь, =*(',). у, =Я',)> е(°Д происходит всасывание внешнего потока внутрь оболочки профиля (рис 2). Устройство отсоса внешнего потока будем моделировать стоком интенсивности (?. Поле скоростей порождённое стоком интенсивности (?, имеет особенность в точке Мч и задаётся формулой

Л* Ху

V*

Рис 2. Стационарное обтекание разомкнутого гладкого профиля (а) и замкнутого профиля (б) с отсосом внешнего потока.

Традиционное условие непротекания на V даёт уравнение

где

-у,

М0еЛ, М0*М„,

\ - £? {? I . \

г2 2л

(9)

(10)

Так как из физики следует, что на V поле скоростей должно иметь особенность типа стока, а на 17- оно является гладким, то можно получить, что решение уравнения (9) имеет вид

если Ь является разомкнутой кривой и /7,(/)е#*на I, и

если £ является замкнутой кривой, где ^(/)€#(«) на [о,/]

Если £ является единичной окружностью, то решением уравнения (6), где в этом случае должна быть функция,

т.е. в смысле главного значения по Коши должно быть верно равенство

1 » « г\ 1-1 /п

(П)

1 г 'о-'б Я

2 2 л

4 л

Если теперь взять нетрадиционное условие непротекания, то есть на I:, то те же рассуждения приводят к уравнению

М0 6/. (12)

Оно отличается от уравнения (9) тем, что оно написано для любой точки М0е£, в частности, и для точки Л/,, но теперь в

правой части появилась S- функция (обобщённая функция). Поэтому теперь интегральный оператор в равенстве (12) надо понимать в смысле теории п.севдодиференциальных операторов. В частности, если L единичная окружность, то уравнение (12) приобретает вид (г70(л/0)=0)

Поэтому из теории псевдодифферециальных операторов без исследования поведения касательных скоростей на Г сразу следует, что решением уравнения (13) будет функция

Вторая глава посвящена квадратурным формулам, которые имеют более физический смысл относительно тех же формул, рассмотренных в книге С.М. Белоцерковского и И.К. Лифанова «Численные методы в сингулярных интегральных уравнениях», 1985г.

Классически в работах С.М. Белоцерковского метод дискретных вихрей формулируется для циркуляционного обтекания в виде

И р

£—'-— = /(*„) j=l,. n , где Г,- интенсивность

i«l хоt ~~

/ дискретного вихря. Затем в работах И.К.Лифанова вместо Г, используется представление в виде y(x,) h. Поскольку в работах Белоцерковского Г, это есть дискретный вихрь равный завихренности на отрезке [х01_,,*„,]» то возникает целесообразность исследования квадратурной формулы для сингулярного интеграла следующего вида:

-1 X0j Xl '-I X0j Xi Ко, ,

В §1 рассматривается сингулярный интеграл:

/(/,)= J^ (15)

L '"'о

по окружности L единичного радиуса с центром в начале координат;

<p(t) - функция класса Н(а), удовлетворяет условию Гёльдера степени а на L, 0<а<1.

На L выбирается два множества узлов (см. рис. 3)

E={tk, к=1,...,п} и Eo={tok, к=1,...,п}, таких, что 4, к=1.....и, разбивают

окружность (рис. 3) на п равных частей, а точка tot является серединой дуги

где полагаем /„,/"'/• Выбранные таким образом множества Е и Е0 называются каноническим разбиением окружности Ь (рис. 3 ).

Рис 3. Каноническое разбиение окружности на части.

Теперь для сингулярного интеграла на окружности берется квадратурная формула:

)=Е—У=1 (16)

В этом случае справедлива теорема:

Теорема 2.1. Пусть <р(1) удовлетворяет условию Н(а) на I. Тогда выполняется неравенство

1пя I I 1 7 = !»••.«.

В §2 рассматривается сингулярный интеграл с ядром Гильберта

1(ва)~- Й?^1 ъФХЮ, (17)

о ^

где <рх{в) - функция с периодом 2я.

Квадратурная формула, аналогичная квадратурной формуле (16), для сингулярного интеграла с ядром Гильберта имеет вид:

) = 7 = 1, п (18)

ы * о,

Доказывается теорема:

Теорема 2.2. Пусть функция на отрезке [0,2я] и

имеет период 2я. Пусть точки Е={вц, к = 1,...,п} и Е0 = { вок, к = 1,...,п} выбраны на [0,2л] следующим образом: -вк = * = 1,...,л-1, вы =вк = . ,п, т.е. точки 1к = е'"' и /01=е'в°*, к=1,...,п, образуют каноническое разбиение I. Тогда справедливо следующее неравенство:

) - (во, \ * о(п а 1п4 ) = 1,.., п

В §3 рассматривается сингулярный интеграл с ядром Коши:

/(<„)= (19)

I 1 'о

где L=[a,b] на действительной оси, а функция q>(t)eЯ* на L, т.е. A'h (t_ayfb_ty ' Г-Р*1- н(а) на [a,b).

Пусть точки a = t0,t„ =Ь разбивают отрезок [а,ь] на л + 1 равных

частей длины л = ~~> а точка /„^является серединой отрезка l^.ViJ» 7 = 0,1, ...п.

Точки множеств E=ftk, к=1. ,п} и E0^{t0j, j=0, ,п} образуют каноническое разбиение отрезка [а,ь\ с шагом h (см. рис 4)

too toi t()2 •on . .

V=a i - ■■ • - f - « ■■ b=t„+1

Ц 4Z И

Рис.4. Каноническое разбиение отрезка [а,б]. Квадратурная сумма для интеграла /(г0) в точке tai следующая:

'Xhtr1!- (20)

к. 1 '» - 'о, ,

На основании теоремы, доказанной А.Ю. Евсеевой*: Теорема 2.4. Пусть Я" на [а,ь\ и множества Е и Е0 образуют каноническое разбиение этого отрезка. Тогда для любой точки t0j е Е0 справедливо неравенство:

ь

~ * »il / А=1 *П > •

^('оД у = 0,1,.. ,п,

•оу >0, ,

где величина #(/„,) удовлетворяет неравенствам:

1) Д, 51, для всех точек /0/ е[а+<У,&-<?], где ¿¡-0;

2) Дг <1, для всех точек /0,

¡-о

далее рассматривается характеристическое уравнение

*#7 = /('о), ^М)- (21)

а 0

Известно, что решение индекса /с=0,1,-1 этого уравнения задается формулой

где у, =l>v0 = v_, =о,

* Евсеева А Ю Об одной модификации истода дискретных вихрей численного решения сингулярных интегральных уравнений на отрезке //Тр IX МДОЗМФ,-2(Х>Ог, Орел, с 209-219

Теорема 2.5. Пусть на [а,ь\ и множества Е и £0 образуют

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

£ с

*-1 0/

7 = 1 ,...п,

) = 1,

7 = 0,1,. .п,

и решениями индекса *■ = 0,1,-1 уравнения (21) выполняется

неравенство

в котором величина 0„('*) удовлетворяет неравенствам: О^О^оДл'"), Я, >-0, для всех точек е[а + 8,Ь-б\, где 3>0\

2)£), Л^О, 'для всех точек ¡к е[а,ь].

4-1

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

Вначале рассматривается сингулярное интегральное уравнение

= /(*„), *ое(-1,1), (22)

где /(х0)может быть гельдеровской функцией степени а т.е. /(х0)е н(а) на [-1,1] или быть дельта - функцией.

Случай, когда /(*0)еЯ(а) был подробно рассмотрен в книге С.М. Белоцерковского и И.К. Лифанова «Численные методы в сингулярных интегральных уравнениях», 1985г., поэтому в данной работе рассматривается случай, когда в правой части стоит функция следующего вида

/(*о) = /■,.,(*„) = в ■ -Ч), *0>9 е (-1,1), где ¿(jCj-^lim^O^-i?),

1 h h

/Г х0 е q—,q + ~ 2 2.

0, h h

*о * q—,q + — . 2 2.

*0e(-U).

(24)

Далее рассматривается метод численного решения уравнения (22) для различных индексов i<-0,1,-1 На отрезке [-1,1] задаются два множества точек Е = {хк,к = 0Х ,п) и Е0 = {*Ojj = 0,l,..,n} удовлетворяющих условиям:

2

•*»!=*' k=l> h=—v

n + Z

Л Л , ,

где pt и р2 - заданные числа и точка при любом п лежит в множестве Е0, д = хь<.

Вначале рассмотрим решения индекса к =1. При моделировании-тонкого разомкнутого профиля вихревым слоем в случае бесциркуляционного обтекания с отсосом внешнего потока, используя нетрадиционное условие непротекания, приходим к сингулярному интегральному уравнению (22), которое рассматриваем для решений класса индекса к =/, т.е. обращающиеся на обоих концах отрезка [-1,1] в бесконечность, а правая часть представляет собой дельта - функцию (см.(23)). Справедлива следующая

Лемма 3.1. Пусть в правой части уравнения (22) стоит функция fs,4(x„), (-1,1). Тогда решение системы линейных

алгебраических уравнений

1 h г, X / \ I 1

я 1.1 JC0/ - xt

X>„(*t)A = 0, (25)

«•I

и значения функции

Ю--1Щ.-0- (26)

* Ь-хг х-Ч связаны соотношением

* = (27)

в котором величина в^х^ удовлетворяет неравенствам:

1) для всех точек х^[-\ + 8,д-8]и[д + 8,\-5\ где 8>0 сколь угодно мало,

<?„<**)£С, >>\ Л, >-0, (28)

2) для всех точек хк е [- ц]

Л^о, (29)

где Сд, С - некоторые константы, не зависящие от п.

Чтобы показать, что функция у(х) в (26) является решением уравнения (22), сингулярный интеграл в нем рассмотрим как оператор в пространствах Н* (типа Соболевских с весом).

Напомним вначале понятие пространства ь1рдействительных функций, квадрат модуля которых интегрируем на [-1,1] с весом

I

Р - р{х), т.е. таких функций /(*), х е [-1,1] , что [р(х) \/(х\'сЫ-<, где

-1

й(х)>0 почти всюду на [-1,1]. Естественным образом вводится скалярное произведение функций в этом пространстве

ШМ= '¡рШШАЬ-

После этого пространство Ь2,р становится гильбертовым, и поэтому можно взять полную систему ортонормальных функций (Рп(х), п=0,1,.. }, являющихся базисом в т.е.

где 8"=0,п*т, и 8" =1,да,я = 0,1, > - символ Кронекера. Любая функция /(х)ьЬ1р представляется рядом Фурье по системе этих функций, сходящимся по норме Ь2,р:

/(«)= \р{х)Ах)рАх\1х

Теперь можно дать следующее определение пространства (типа соболевского) с весом р, связанного с пространством Ь2,р.

Определение. Пусть Лей - произвольное действительное число. Весовым пространством Н1Р = Н*{{Р„(х)}) на [-1,1] называется множество таких функции (или обобщенных функций дляЛчО) и(х), что функция

п-0

и{п) = ]р{х)и{х)Р„ {¿уь = (и(4 Р. М),

-I

принадлежит пространству где п=тах{],п}. Нормой функции и(х) в пространстве Нхр назовем число =^.0??и|й(л)|!.Скалярным произведением функций и(х) и \(х) из пространства Я* назовем число (и, = Относительно этого скалярного

произведения пространство Н* становится гильбертовым.

Рассматривая интеграл в (22) как оператор в пространствах

нп гДе Р\ ~ I—-—- и Рг = л/1-*2, показывается, что функция

VI-*2

(26) является единственным решением (22).

Отметим, что справедливо следующее замечание: Замечание 3.1. Если рассмотреть гладкий слабоизогнутый тонкий профиль, то нетрадиционное выполнение условия непротекания приводит к тому, что правая часть в уравнении (22) будет представлять сумму гёльдеровской и дельта -функции: /(х0)+/^(х0), где /(*0)б£2й. В этом случае все решения г(х) индекса к= 1 в паре пространств (Я* ;//£),Л-¡-1/2, даются формулой.

1 'гУ'-*ог ЖК | с 11 УГу

Я-4Х-Х1 }, х-х0 я-л/1-*2 * VI-*2 9-*'

хе(-М),

где функция у(х) удовлетворяет равенству

При моделировании тонкого, разомкнутого профиля вихревым слоем в случае безударного обтекания с отсосом внешнего потока, используя нетрадиционное условие непротекания, приходим вновь к уравнению (22), которое рассматривается для решений индекса к =-/, при наличии в правой части функции }д,ч(хо), в паре пространств Н"п и Я* и

показывается, что при этом должно выполняться условие разрешимости для индекса к ='-7 .

j/^y^o, (30)

-i VI -хг

где /(*)- гельдеровская функция, то из (30) следует Далее рассматривается уравнение

Ií*2£ = /(Xe) + / х0 е(-1,1), (31)

л дс0 - X

Для численного решения уравнения (31) на основе идеи метода дискретных вихрей воспользуемся множествами Е = {хк, к = 1,. ,.,п} и Ео = {x0¡, j = 0,1,.. , п}, построенными выше. Справедлива следующая

Лемма 3.3. Пусть функция f(x) принадлежит классу Н(а) на [-1,1], а функция fi,q(x) имеет вид (23). Тогда решение системы линейных алгебраических уравнений

= J = 0X...,n, (32)

я Л ».I *оJ ~хк

где yon - регуляризирующая переменная , и значения функции

„(_и)> (33)

я u4\-xl(x-x¡¡) nJxTqix-q

связаны соотношением (27), в котором величина в„(Хк) удовлетворяет неравенствам (28) и (29).

Наконец, моделируя тонкий, разомкнутый профиль вихревым слоем в случае циркуляционного обтекания с отсосом внешнего потока, и используя нетрадиционное условие непротекания, рассмотрим уравнение (22) для решений индекса к =0 в паре пространств Я* и Яр ,где

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

Е = {хк, к = 1, ,п}, Е0 = {x0j, j = 0,1,. , п},

построенные выше, аналогично леммам 3.1 и 3.3 с использованием теоремы 5.1.1 из книги С.М. Белоцерковского и И.К. Лифанова «Численные методы в сингулярных интегральных уравнениях», 198 5 г. доказывается

Лемма 3.5. Пусть в правой части уравнения (22) стоит функция /д,ч(х0), х0 е (-1,1). Тогда решение системы линейных алгебраических уравнений

= 7 = 1...,", (34)

Я к-1 хи, ~ Хк

и значение функции

ДГ0€(-Ц), (35)

связаны соотношением (27), в котором величина ва(хк) удовлетворяет неравенствам (28) и (29).

Если интеграл в уравнении (22) рассматривать как оператор в

паре пространств то доказывается, что уравнение

(_Ц) (36)

я м х0-х

имеет единственное решение у(х), даваемое формулой (35).

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

1(Нх)Ж + | К(х^х)г{х)ск = {(Хв) + 1 (д:д х0 е (-1,1), (37)

я м х0 - X "

где функции /(х) и К(х0, х) принадлежат классу Н(а) на соответствующих множествах. Предполагается, что уравнение (37) имеет единственное решение в паре пространств ИД).

Пусть теперь оператор (37) есть оператор из Я* в Я^ для решений индекса к =1 Решение у(х) ищется в виде

= + хое{-Ц), (38)

я^-х^х-Ч я< 1-х2

Тогда для функции получаем интегральное уравнение

Фредгольма второго рода

4>(х)-[к(х,тМг)с/г = Д(дг,с) + /<1(дг), *„ € (-1,1),

«м г — v.

а-с

Я" л/1-Г2 *-*»

которое, в силу сделанного предположения, имеет единственное решение.

Аналогичные рассуждения сведения сингулярного интегрального уравнения к уравнению Фредгольма можно провести и для решений индексов к=-1, к=0. Только теперь при к=-1 условие разрешимости несколько усложняется для поиска числа Q.

Для численного решения уравнения (37) можно опять применить идеи метода дискретных вихрей и рассмотреть для нахождения решений индексов к = 1,-1,0 системы линейных алгебраических уравнений вида (25), (32) и (34) соответственно. Справедлива следующая

Теорема 3.1. Пусть в уравнении (37) функция /(*„) и ядро к(хй, ^принадлежат классу Н(а) на [-1,1] и [-1,1]х [-1,1] соответственно, а функция /¿гЧ(хо) имеет вид (23). Тогда решения систем линейных алгебраических уравнений

% к-1 Х0/~~Хк кш 1

= Д*о,) + /.м(*оД У = (40)

я к=I х0) ~Хк »-I

= /(*о,) + /.,,(*оу)>

|>„(х,)/, = С, (41)

Уо^-ЪЪ^-^Кх,,,^) гЛ*>) к-

я я 75 х,

4-1 |

= /(*») + /,„(*.,). 7 = 1.....и. (42)

и соответствующие решения индекса к =0,1,-1 уравнения (37) связаны соотношением (27), в котором величина удовлетворяет неравенствам (28) и (29).

Если в задаче обтекания профиля моделировать профиль слоем диполей, то задача нахождения интенсивности слоя диполей

сводится к следующему гиперсингулярному интегральному уравнению

11 л

-/(*„), хае(-Ц), (43)

4*0 -■*)'

которое эквивалентно уравнению

е'(х)сЬс

Ч!

я «

/(х„), х0 € (-1,1), (44)

к

при условии

£*'(*)<& = 0. (45)

Если в уравнении (43) справа стоит функция /¿,ч(х0) (см.(37)), тогда уравнение (44) имеет своим решением функцию

*'(*) = ,ое(-щ (46)

я4\-хг х-Ч

а, следовательно, в силу равенства (45) уравнение (43) имеет своим решением функцию

в(,)я1|п-^- *0е(-1,1) (47)

Для численного нахождения решения (47) уравнения (43) в случае, когда в правой части стоит функция /».ч(хо), опять воспользуемся идеей метода дискретных вихрей с использованием множеств Е ~ {хк, к = /,. ..,п} и Ед = {хц. ] ~ 0,1, ., п}, построенных в § 1,гл. 3. Справедлива

Теорема 3.2. Пусть в правой части уравнения (43) стоит функция /¿^(хо), дг0,?б(-1Д) Тогда решение системы линейных алгебраических уравнений

1 "

П 4.1

1 1

лоу

>-9). ) = <д • (48)

и значения решения %(х), даваемое формулой (47), уравнения (43) связаны соотношением

(*<*«)- 8.С*»)]* С&\ Х5 > 0, (49)

где С не зависит от п.

В главе 4 рассматриваются схемы метода дискретных вихрей для сингулярных интегральных уравнений с ядром Гильберта первого рода в классе обобщенных функций и их обоснование.

В §1 рассматривается сингулярное интегральное уравнение с ядром Гильберта, которое мы получаем, моделируя телесный профиль вихревым слоем с отсосом внешнего потока, используя нетрадиционное условие непротекания

где в правой части есть дельта - функция

0- произвольное действительное число, а дельта - функция <У(0о-<1) с носителем в точке q определяется условиями:

^-•ЧЛ V; (51)

[+0О, в0=ч,

р(0о-*>Юо=1. (52)

о

В главе 1 показано, что уравнение (50), рассмотренное в обобщенных пространствах Нх, имеет решение только при выполнении условия:

(53)

о

При выполнении условия (53) решение уравнения (50) дабтся формулой

(54)

1 2'

где — Ш)ав = С, б»е[0;2я-]. (55)

2я о

Используя идеи метода дискретных вихрей, построим метод численного решения уравнения (50). Пусть отрезок [о;2л-] разбит двумя множествами точек £ = = 1,2,...,«} и £0=ф ,у = 1,2,...,л)

1п

таким образом, чтобы вул-вк=Ь = —,

п

к-1,2,...,п; 0^=0,; 0Ь = <?,+- и ?е£0, при этом я = в0^.

Теорема 4.1. Пусть в правой части уравнения (50) стоит сумма функций /(ва)<= н(а) и /¡ч(вй), 0о>де[о,2я]. Тогда между решением системы линейных алгебраических уравнений

2гг*-> 2 (56)

.¿я кт1 п где

(0, 1*], .

1 в л л \ . * = (57)

и решением уравнения (50), задаваемым формулой (54) при условии (55), выполняется соотношение

п

В §2 исследуется полное сингулярное интегральное уравнение, которое получается при моделировании обтекания произвольного телесного профиля вихревым слоем с отсосом внешнего потока с использованием нетрадиционного условия непротекания. Задача нахождения интенсивности вихревого слоя сводится к решению сингулярного интегрального уравнения

~' ¡СЛ?Зс*у(в)ао+г]фп,в)у{в)ав=М) + /,,{в0\ 0о,9е[о,2ж]. (58)

¿71 0 X „

Его решение ищется в виде

0,де[о,2г] (59)

¿к I

Тогда для функции ц/(в) уравнение (58) перейдет в уравнение

± (к1±Ч1{о)ав + |А- (ва, =

о 2 о

= /И+ £ + ^'¡Фо^Ц^-М, 00 <= [0,2ж\ (60)

2.Л 2я * 2

При этом условие разрешимости уравнения (58) в классе функций вида (59) получит вид

|/(Ф - = О. (61)

Для уравнения (58) предлагается следующий численный метод решения

4-1 * 4-1

= 7 = 1,...«.

¿Я 4»,

Доказывается, что /„(^)сходится в точках Ок к точному решению.

Далее в §3 рассматривается уравнение с гиперсингулярным периодическим ядром. Если в задаче обтекания профиля моделировать профиль слоем диполей, то задача нахождения интенсивности слоя диполей сводится к следующему гиперсингулярному интегральному уравнению

¿1—5= г>0е [0,24 (62)

2

Уравнение (62) эквивалентно (в силу определения гиперсингулярного интеграла)уравнению

= /(0О)+Л>о), 0О 6 [о,2/г].

где0о е[0,2ж]. Решение уравнения (62 )получаем в виде

$т((в~д)/2)

5ГП(9/2)

+ С,

(63)

(64)

Для численного решения уравнения (62) вновь рассматриваем два множества точек Е к Е0, полагая д = 0„и е£0. Используя идеи метода дискретных вихрей, построим метод численного решения уравнения (62). Заменим (62) следующей системой линейных алгебраических уравнений

2 лы

---

(65)

Систему (65) можно записать в виде системы (56) относительно функции

е,

¿Л

±8'М)Ь = О

Для однозначного нахождения служат равенства

(66)

*-]

Разрешая систему (66) и используя квадратурные формулы для интеграла с ядром Гильберта имеем:

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

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

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

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

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

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

Автор выражает благодарность профессору И.К. Лифанову за постоянную поддержку, внимание и помощь ' при подготовке диссертационной работы.

Основные результаты диссертационной работы отражены в следующих публикациях:

1. Лебедева Н.В., Сингулярный интеграл по замкнутому контуру // Тр. X Междунар.симп. «МДОЗМФ-2001». Харьков-Херсон, 2001г., с.192-195

2. Вайникко Г.М., Лебедева Н.В., Лифанов И.К. Численное решение сингулярного и гиперсингулярного интегральных уравнений на отрезке и дельта - функция. // Математический сборник , 2002г., т. 193, №10, с.3-16

3. Лебедева Н.В., Вопросы численной реализации решений сингулярного интегрального уравнений на отрезке в классе обобщенных функций. // Тр. Междунар. Школ-семинаров «МДОЗМФ». Орел, 2002. с.49-53

4. Лебедева Н.В. .Численное решение сингулярного и гиперсингулярного интегрального уравнений с ядром Гильберта первого рода и дельта -

функция. // Тр. XI Междунар. симп. «МДОЗМФ-2003». Харьков-Херсон, 2003г., с. 165-169 5. Lebedeva N.V. and Lifanov I.К. Singular and Hipersingular integral

equations with the Hilbert kernel, delta-function, method of discrete vortices // Proc. Estonian Acad. Sci. Phys. Math., 2004, 53/2, p.84-91

Лебедева Нина Васильевна Математическое моделирование Автореферат Под авторской редакцией Подписано в печать 8.07.2005. Формат 60х84\16. Бумага для офисной техники. Гарнитура Тайме. Уч.-изд. 1,1. Тираж 100 экз. Заказ № 199. Цена договорная. РИО. ВТУ. Балашиха

»14234

РНБ Русский фонд

2006-4 21180

Введение.

Глава 1. Математическая модель обтекания профиля с отсосом внешнего потока и сингулярное интегральное уравнение.

§ 1 .Моделирование профиля вихревым слоем

§2. Сингулярное интегральное уравнение для разомкнутого контура.

§3. Сингулярное интегральное уравнение для замкнутого контура.

§4.0тсос внешнего потока и требуемый класс решений для сингулярного интегрального уравнения.

§5. Отсос внешнего потока и дельта - функция в правой части сингулярного интегрального уравнения.

Глава 2. О некотором обобщении квадратурных формул типа метода дискретных вихрей для сингулярных интегралов.

§1. Сингулярный интеграл на окружности.

§2. Сингулярный интеграл с ядром Гильберта.

§3. Сингулярный интеграл с ядром Коши.

Глава 3. Схемы метода дискретных вихрей для сингулярных интегральных уравнений первого рода на отрезке в классе обобщенных функций и их обоснование.

§ 1. Характеристическое сингулярное интегральное уравнение на отрезке.

§ 2.Сингулярное интегральное уравнение первого рода на отрезке.

§ 3. Характеристическое гиперсингулярное интегральное уравнение на отрезке.

§ 4. Результаты численного эксперимента.

Глава 4. Схемы метода дискретных вихрей для сингулярных интегральных уравнений с ядром Гильберта первого рода в классе обобщенных функций и их обоснование.

§ 1. Характеристическое сингулярное интегральное уравнение с ядром Гильберта.

§2 Полное сингулярное интегральное уравнение.

§3 Характеристическое гиперсингулярное интегральное уравнение.

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

Следует отметить, что при - аналитических исследованиях в приложениях уже давно некоторые задачи стали сводить к сингулярным интегральным уравнениям, т.к. для них в одномерном случае получена хорошая теоретическая база, которая довольно полно изложена в монографиях [30] Гахова Ф.Д. «Краевые задачи» - Москва, изд. Наука, 1977г. и [48] Мусхвелишвили Н.И. «Сингулярные интегральные уравнения » - Москва, изд. Наука, 1968г.

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

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

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

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

Идея метода дискретных вихрей состоит в следующем. Непрерывный вихревой слой, моделирующий несущую поверхность и след за нею, заменяется системой дискретных вихрей. На несущей поверхности выбираются точки, называемые расчетными, в которых выполняется условие непротекания (сумма нормальных составляющих скоростей, индуцируемых вихрями, и набегающего потока равна нулю). Задача нахождения неизвестных циркуляций дискретных вихрей сводится к системе алгебраических уравнений. Решение задачи не единственно и может иметь особенности на кромках и изломах несущей поверхности. Нужный класс решения определяется физическим содержанием задачи и выделяется выбором взаимного расположения множества дискретных вихрей и расчетных точек. К тем кромкам, где решение должно быть неограниченным, ближайшими располагаются дискретные вихри, а к тем, где оно должно быть ограниченным - расчетные точки. Кроме того, суммы, которыми заменяются сингулярные интегралы в теории несущей поверхности, должны соответствовать главным значениям интегралов в смысле Коши. Для этого внутренние расчетные точки должны лежать посередине между вихрями на поверхности (или стремиться к этим положениям). Именно в таком виде впервые был сформулирован метод дискретных вихрей в 1955г. в докторской диссертации С.М.Белоцерковского, после чего началась систематическая реализация его в аэродинамике [1,2,3,6,8].

Так как по существу метод дискретных вихрей использует вычисление сингулярного интеграла с помощью специальных квадратурных сумм типа прямоугольников, а при выборе класса решения не используется явное выделение особенности на кромках, то он начал подвергаться критике [50].

Строгое математическое обоснование метода дискретных вихрей было выполнено И.К. Лифановым в ряде работ [36-46], завершившихся защитой в 1981г. докторской диссертации. После этого удалось перенести идеи метода дискретных вихрей в теорию упругости и электродинамику [4, 5, 7, 19-29]. Эти первые результаты бьши изложены в [6, 52].

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

В последнее время в аэродинамике начали исследовать механизацию крыльев с использованием устройств отсоса внешнего потока. Было показано для тонкого профиля, что эта задача сводится к решению сингулярного интегрального уравнения первого рода на отрезке в классе функций, которые в точке отсоса имеют неинтегрируемую особенность. В некоторых частных случаях было дано численное решение таких задач с помощью метода дискретных вихрей. Однако каждый частный случай требовал особого подхода. При численном решении задачи с отсосом внешнего потока методом дискретных вихрей расчетные точки выбирают так, чтобы точка отсоса являлась одной из них, и при написании системы линейных алгебраических уравнений то уравнение, что соответствует точке отсоса, пропускается. Поэтому, например, при решении методом дискретных вихрей циркуляционной задачи для тонкого профиля получается система с п неизвестными и п-1 уравнениями. Затем из тех или иных физических соображений ищется дополнительное уравнение, чтобы получить разрешаемую однозначно систему п уравнений с п неизвестными [18], [35]. Это показывает, что для каждой новой задачи с отсосом (несколько новых точек отсоса внешнего потока) надо думать, как замкнуть полученную систему линейных алгебраических уравнений. Однако, в задачах обтекания телесного профиля до сих пор для некоторых случаев непонятно было как это сделать. Взяв несколько иную физическую трактовку задачи с отсосом [16], приходим к тому же самому сингулярному интегральному уравнению первого рода, но в правой части этого уравнения появляется дельта — функция с носителем, расположенным в точке отсоса. Оказалось, что наличие дельта - функции в правой части позволяет строить более удобные методы численного решения соответствующих задач аэродинамики. Действительно, в старом подходе к решению задач аэродинамики [15,52,35,53] каждая новая точка отсоса на профиле приводила к необходимости изменения системы линейных алгебраических уравнений . При новом же подходе можно построить вычислительный алгоритм решения этой задачи, в котором с.л.а.у. остается всегда неизменной и меняется только правая часть. Более того, доказательство сходимости численного решения к точному в этом случае тоже стало ближе к классическому доказательству при хороших правых частях этого интегрального уравнения. Именно таким доказательствам и посвящена представленная работа. Таким образом, актуальность выбранной темы не вызывает сомнений.

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

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

Рассматривается плоскопараллельное обтекание изолированного непротекаемого профиля установившимся потоком идеальной невязкой несжимаемой жидкости, скорость которого и0 =и0х1 + и0у]. Сам профиль является неподвижным. Под профилем понимаем цилиндрическую поверхность с образующей параллельной оси 02 и направляющей кривой Ь в плоскости ОХУ (рис.1). Так как все параметры набегающего потока и возмущенного не зависят от координаты г, то в дальнейшем профиль изображается только кривой Ь в плоскости ОХУ. Если кривая Ь является простой разомкнутой (рис. 1а), то такой профиль называется тонким, а если кривая Ь является простой замкнутой кусочно - гладкой (рис. 16), то такой профиль называется телесным. Поскольку рассматривается стационарная задача, то следа за профилем нет, и возмущенное течение будет потенциальным везде вне профиля Ь.

Рис. 1 Разомкнутый (а) и замкнутый (б) (телесный) профили

При моделировании профиля вихревым слоем задача сводится к решению следующего сингулярного интегрального уравнения где /(М0) = -и0пм<, у{м)~ интенсивность распределенного вихревого слоя в точке М е I и /(М0) гладкая функция, /(М0)е Н(а) на Ь.

Далее рассматриваются сингулярные интегральные уравнения каждого типа, к которым сводится уравнение (1).

Если Ь - гладкая разомкнутая кривая, заданная параметрически х = х(0, у = ХО, ' е [-1,1], то уравнение (1) может быть записано в виде х а) б) О где/(0еЯ(ог) на [-1,1], Н(а) на [-1,1]х[-1,1].

Это уравнение имеет решения трех типов: т= где |/(А/)й& = 0.

4) г«)7><о (5) где функция не имеет особенностей на Ь, соответствующие циркуляционному, бесциркуляционному и безударному обтеканию профиля I соответственно.

Если профиль Ь является телесным без острой кромки, т.е. кривая Ь является простой замкнутой гладкой, и параметр ( изменяется на отрезке [0,2я], то уравнение (1) можно записать в виде

4я- о о

Так как набегающий поток является потенциальным во всей плоскости, то для любого телесного профиля выполняется равенство

6)

М)Л = 0

7)

Для телесного профиля решение уравнений (1) определенно с точностью до константы, т.е. в этом уравнении константа является собственным решением.

Далее рассматривается задача с отсосом внешнего потока. Пусть в точке Мч(хч,уд) профиля Д хд = х(гд), уч = у(гд), г,е(о,/), происходит всасывание внешнего потока внутрь оболочки профиля (рис 2).

Устройство отсоса внешнего потока будем моделировать стоком интенсивности Поле скоростей ' Уч(м0), порождённое стоком интенсивности Q , имеет особенность в точке Мч и задаётся формулой

1п

8) б)

Рис 2. Стационарное обтекание разомкнутого гладкого профиля (а) и замкнутого профиля (б) с отсосом внешнего потока.

Традиционное условие непротекания на V даёт уравнение К{М0)пМп =-и0 {М0)пМа -уч (М0)пМа, М0<ЕЬ, М0*МЧ, (9) где

7<M0) =f iMbJ.fi. (10) l7t rMoM, 2;r м0 e L, Л/0 * Л/,

Так как из физики следует, что на 1+ поле скоростей должно иметь особенность типа стока, а на I" - оно является гладким, то можно получить, что решение уравнения (9) имеет вид п tq-t если L является разомкнутой кривой и 77, (/)е Н* на L, и

2 к 2 если L является замкнутой кривой, где?72(0 е Н(а) на [о,/].

Если L является единичной окружностью, то решением уравнения (6), где в этом случае K(t0,t) = 0, должна быть функция, t„-t

2л- 2 т.е. в смысле главного значения по Коши должно быть верно равенство по

4л- ^ 2 2л 2, 4л Если теперь взять нетрадиционное условие непротекания , то есть на , то те же рассуждения приводят к уравнению

К(Мо)»ма =-й0{М0)пМа-£-щ{мфМч)-£з(мй-мч\ М0 е Ь. (12)

Оно отличается от уравнения (9) тем, что оно написано для любой точки М0 е X, в частности, и для точки Мч, но теперь в правой части появилась 8 - функция (обобщённая функция). Поэтому теперь интегральный оператор в равенстве (12) надо понимать в смысле теории псевдодиференциальных операторов [55,56]. В частности, если Ь единичная окружность, то уравнение (12) приобретает вид {П0{ма) = 0) з)

Поэтому из теории псевдодифферециальных операторов без исследования поведения касательных скоростей на V сразу следует, что решением уравнения (13) будет функция

14)

2л 2

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

Классически в работах С.М. Белоцерковского [12] метод дискретных вихрей формулируется для циркуляционного обтекания в виде т р

-1— = /(*<>)> / = 1, .л , где Г,- интенсивность / дискретного

1=1 ) ~ вихря. Затем в работах И.К.Лифанова вместо Г, используется представление в виде /(*,)•/г. Поскольку в работах Белоцерковского Г(. это есть дискретный вихрь равный завихренности на отрезке 1х0>|.,,х0(], то возникает целесообразность исследования квадратурной формулы для сингулярного интеграла следующего вида:

-1 X0j Х «=1 Х1 *„,,

В § 1 рассматривается сингулярный интеграл:

1М-1&Г <15) ' 'о по окружности Ь единичного радиуса с центром в начале координат; р(г) — функция класса Н(а), удовлетворяет условию Гёльдера степени а на Ь,

0<а<1.

На Ь выбирается два множества узлов (см. рис. 3) Е={^ к=1, .,п) и к=1,.,п}, таких, что к=1,.,п, разбивают окружность (рис. 3) на п равных частей, а точка Аж является серединой дуги где полагаем п+1~*ь Выбранные таким образом множества Е и Ео называются каноническим разбиением окружности Ь. (рис. 3 )

Теперь для сингулярного интеграла на окружности берется квадратурная формула:

Л 1 т~ J=l■■■>«■ (16)

-\h~hj ,„,

В этом случае справедлива теорема:

Теорема 2.1. Пусть (р($ удовлетворяет условию Н(а) на Ь. Тогда выполняется неравенство

1п П | I 1 j = l,.,n,

В §2 рассматривается сингулярный интеграл с ядром Гильберта

2т 0-0. = \ctg—±-<px(e)de, (17) О где q>x (в) - функция с периодом 2ж.

Квадратурная формула, аналогичная квадратурной формуле (16), для сингулярного интеграла с ядром Гильберта имеет вид: п Q — Q вы*\

Sn{^j) = t.ctSJ4rL \<P^)dO j = 1,(18) i=l ^ a 1 "Ok

Доказывается теорема:

Теорема 2.2. Пусть функция ^>,(0) еЯ на отрезке [0,2л] и имеет период 2к. Пусть точки Е={вь к=1,.,п} и Е0={ Ооь к=1,.,п} выбраны на [0,2л] следующим образом: вк+1 -вк = ~ = h, к = 1,., п -1, + 2п = вплЛ, в0к =вк+н/2,к = 1,.,п, т.е. точки tk=e'e" и tok=eie>ot, к=1,.,п, образуют каноническое разбиение L. Тогда справедливо следующее неравенство: |I{0Oj ) - (0Oj )| < о{па ln п), J = 1,., п.

В §3 рассматривается сингулярный интеграл с ядром Коши:

19)

L 0 где L = [a,b\ на действительной оси, а функция <p(t)e H* на L, т. е. р(*) = ~/-Z^} v > y/(t)eH(a) ца [a,b]. v ~а) \Ь~ч

Пусть точки а = /0,г„+1 =6 разбивают отрезок [а,б] на и + 1 равных частей длины Л = , а точка г07 является серединой отрезка 1^,/,+,], j = 0,1,.,«.

Точки множеств E={th к=1,.,п} и Ео={Щ, j=0,.,n} образуют каноническое разбиение отрезка [a,b\ с шагом h (см. рис 4).

00 toi to2 ton ,. to=a i - « " » ■ » - » " ■ i b=tn+1 tj 12 tn

Рис.4. Каноническое разбиение отрезка [a,Z>]. Квадратурная сумма для интеграла l(t0) в точке t0J следующая: лУ-ЁтЛ"!^- (20)

1 к

Доказывается теорема [32]: Теорема 2.4. Пусть <p(t)eH' на [а,ъ\ и множества Е и Е0 образуют каноническое разбиение этого отрезка . Тогда для любой точки t0J е Е0 справедливо неравенство: jf^-ZrVV» « — «п . ¿.I ' t ~ 'n . . г?(/0Д у = 0,1,

0) 4-1 ^0j где величина ) удовлетворяет неравенствам :

1) )< оДа* ) 0 Л, < 1, для всех точек /0, е [а + б, Ъ - б], где 8 >- 0;

2) ^ )о ^ А2 < 1, для всех точек е [а,Ь]. о

Далее рассматривается характеристическое уравнение

0, «.М). 'л »

21)

Известно, что решение индекса к =0,1,-1 этого уравнения задается формулой где у1 =1, у0 =у, =0,

Я" t-a 1

ЛЬ-^-аУ

Теорема 2.5. Пусть р(/)еЯ' на [а,б] и множества Е и Е0 образуют каноническое разбиение этого отрезка. Тогда между решениями систем линейных алгебраических уравнений

Я 1 'о* . .

1 'о/ h .

7=1,. п, 1 f°i / \ St-1- J =/У

-» foy h л 'о*

X \<pn{t)dt = C 1 '°Г ( \

-1 r0j у =1,.я-1, j = 0,1,.«, и решениями индекса к = 0,1,-1 уравнения (21) выполняется неравенство j 'о, т К ('У' в котором величина 6>„ ) удовлетворяет неравенствам:

1)0я(г4)£оДад,)> Для всех точек е [а + д,Ь-б], где 6>- 0;

2)^6к V1 ^ ), ^ х 0, для всех точек (к е [а,б]. 1

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

Вначале рассматривается сингулярное интегральное уравнение тг ' V — V л х0-х где /(х0)может быть гельдеровской функцией степени а т.е. /(х0)е Н(а) на [-1,1] или быть дельта - функцией.

Случай, когда /(х0)еЯ(а) был подробно рассмотрен [6], поэтому в данной работе рассматривается случай, когда в правой части стоит функция следующего вида

Я*о) = /*.<(*о) = Я• ¿(*о -Я), (-1,1), (23) где 8(х0 -д) = Нш 8Н (.х0 -д), й-м)

1 А' *0 6 И к д—,д + — 2 2.

0, к к

-и).

24)

Далее рассматривается метод численного решения уравнения (22) для различных индексов к =0,1,-1- На отрезке [-1,1] задаются два множества точек Е = {хк,к = 0,1,., и} и Е0 = {х07= 0,1,., и} удовлетворяющих условиям: 2 к+1 -хк\ = к, к = 1,.,я-1, к =--, п + 2

0>+1 ~ Х0) | = К ] = ОД,., п -1, н х\ хоо ~ 2' ~~ ** + 2' п' х00+\ = Н-д^, 1-х0л = к-д1Ь, О^р^д^, дгм <, р2 <+*>, где р! и р2 - заданные числа и точка д при любом п лежит в множестве Е0, д = *оу (см. [6;§1.3], там же указано, как эти множества построить).

Вначале рассмотрим решения индекса к =1. При моделировании тонкого разомкнутого профиля вихревым слоем в случае бесциркуляционного обтекания с отсосом внешнего потока, используя нетрадиционное условие непротекания, приходим к сингулярному интегральному уравнению (22), которое рассматриваем для решений класса индекса к =/, т.е. обращающиеся на обоих концах отрезка [-1,1] в бесконечность, а правая часть представляет собой дельта - функцию (см.(23)). Справедлива следующая

Лемма 3.1. Пусть в правой части уравнения (22) стоит функция 1д.я(хо), х0,£€(-1,1). Тогда решение системы линейных алгебраических уравнений

2>„(дг*)-Л» о, (25) и значения функции л 1 У1"?2 Я у(х) = —(26) связаны соотношением

И**) -гл** >1 * <х*)- *=* (2?) в котором величина вп(х^) удовлетворяет неравенствам:

1) для всех точек хк е [-\ + ё^-б\^>\щ + 6,1-д\ где д>0 сколь угодно мало, в„(хк)<С8-^, А,>-0, (28)

2) для всех точек хк е [-1,1]

0п(х,).А<СЛ\ 0, (29) где С& С - некоторые константы, не зависящие от п.

Чтобы показать, что функция у(х) в (26) является решением уравнения (22), сингулярный интеграл в нем рассмотрим как оператор в пространствах Н* (типа Соболевских с весом).

Напомним вначале понятие пространства Ь2р действительных функций, квадрат модуля которых интегрируем на [-1,1] с весом р = р{х), 1 т.е. таких функций /(*), х е [-1,1] , что $р(х)- ¡/(х)|2 <£с-< +а>, где р(х) >- 0 почти

-1 всюду на [-1,1]. Естественным образом вводится скалярное произведение функций в этом пространстве

-I

После этого пространство становится гильбертовым, и поэтому можно взять полную систему ортонормальных функций {Р„(х), п=0,1,.}, являющихся базисом в £2,/э т.е. хР(х)Рп{х)Рт{х)с1х = д:, где б" =0,п¥= т, и 8™ =1,/и, я = 0,1,., - символ Кронекера. Любая функция /(*) е Ь2р представляется рядом Фурье по системе этих функций, сходящимся по норме Ь2(Р: to = t Л»)РП to, /to = J'p(x)/(x)pn {x)dx n=0 -1 lift =lp(x)f2(x)dx = fj/2(n)^+co. 2" /1=0

Теперь можно дать следующее определение пространства Н* (типа соболевского) с весом р, связанного с пространством Ьг.р.

Определение [16; §4]. Пусть AeR - произвольное действительное число. Весовым пространством Нлр = Н1р ({Рп (х)}) на [-1,1] называется множество таких функций (или обобщенных функций дляЯ-чО) и(х), что функция оо п=0 й(п) = ]р{х)ы{х)рп {x)dx = Ш Рп to),

-I принадлежит пространству L,2tP, где п=тах{1,п}. Нормой функции и(х) в

2 V—«о 2Я 2 пространстве Я* назовем число ||м||я = я |м(я)| .Скалярным произведением функций и(х) и v(x) из пространства Нхр назовем число u,v) = 2]"on u(n)v(n) .Относительно этого скалярного произведения пространство Н* становится гильбертовым.

Рассматривая интеграл в (22) как оператор в пространствах Я* и Н*г, где рх = , * и р2= Vl-x2, показывается, что функция (26) является VI-х2 единственным решением (22).

Отметим, что справедливо следующее замечание: Замечание 3.1. Если рассмотреть гладкий слабоизогнутый тонкий профиль, то нетрадиционное выполнение условия непротекания приводит к тому, что правая часть в уравнении (22) будет представлять сумму гёльдеровской и дельта - функции: f{x0)+fSq{x0), где /(х0)е12рг. В этом случае все решения у(х) индекса к= 1 в паре пространств (Я* ;ЯД)Дк-1/2, даются формулой (ср. [36; (2.3.30)]) vM 1 уУь^Г-ЖК с 1 УГу q г\х) =--1-- I-+-7== +--,----.

7Г • "V1 — X -1 *-*о Wl-X2 n V1-X

X € (—1,1) , где функция х(х) удовлетворяет равенству j{x)dx = C.

При моделировании тонкого, разомкнутого профиля вихревым слоем в случае безударного обтекания с отсосом внешнего потока, используя нетрадиционное условие непротекания, приходим вновь к уравнению (22), которое рассматривается для решений индекса к =-1, при наличии в правой части функции fs,q(xo) в паре пространств Я*2 иЯ^ и показывается, что при этом должно выполняться условие разрешимости для индекса к =-1 . у(х);ш<ь=о (зо)

-1 л/ 1-х2 где /(х)- гельдеровская функция. Тогда из (30) следует, что Далее рассматривается уравнение

1 f ZW^L = /(Xe) + / (Xe)f х0 е (-1,1), (31)

71 м Х0 -X

Для численного решения уравнения (31) на основе идеи метода дискретных вихрей воспользуемся множествами Е = {хк, к = 1,. .,п} и Ео = {xoj.J - 0,1,., п}, построенными выше. Справедлива следующая

Лемма 3.3. Пусть функция f(х) принадлежит классу Н(а) на [-1,1] (см. [6; §1.1]), а функция fs,q(x) имеет вид (23). Тогда решение системы линейных алгебраических уравнений

-пп +1¿z2^ = /(*0j) + q¿h(*0,-я), j = 0,1,.,", (32)

Л Л XÜJ —Xk где yon - регуляризирующая переменная [6; §5.1], и значения функции

Кх).JEZ| я***. хе(1Л) (зз) связаны соотношением (27), в котором величина Qn(x¡) удовлетворяет неравенствам (28) и (29).

Наконец, моделируя тонкий, разомкнутый профиль вихревым слоем в случае циркуляционного обтекания с отсосом внешнего потока, и используя нетрадиционное условие непротекания, рассмотрим уравнение (22) для решений индекса к =0 в паре пространств ЯД и ЯД, где

1-х 1 + х

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

Е = {хк, к = ],. .,«/, Е0 = {х0и] = 0,1,., п}, построенные выше, аналогично леммам 3.1 и 3.3 с использованием теоремы 5.1.1 из [6; §5.1] доказывается

Лемма 3.5. Пусть в правой части уравнения (22) стоит функция/¿,ч(хо), х0 е (-1,1). Тогда решение системы линейных алгебраических уравнений у = !,.,„, (34)

Я *=1 j ~~ Хк и значение функции

Я,)-1^1114-2-, .(.Ц), (35) я V1 -нл:х-ц связаны соотношением (27), в котором величина вп(Хк) удовлетворяет неравенствам (28) и (29).

Если интеграл в уравнении (22) рассматривать как оператор в паре пространств то доказывается, что уравнение

1= (-1ДХ (36) я *~х х0 - X имеет единственное решение у(х), даваемое формулой (35).

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

I * + | д,о)х)кх)л = /(хо) + / (Хо)) х0 е(-1,1), (37) я м х0-х где функции /(х) и К(хо, х) принадлежат классу Н(а) на соответствующих множествах. Предполагается, что уравнение (37) имеет единственное решение в паре пространств , НД ). Пусть теперь оператор (37) есть оператор из Н* в НхРг для решений индекса к =1.

Решение у(х) ищется в виде

38)

Я л/\-х2 х-я Я"л/1 — X

Тогда для функции у/(х) получаем интегральное уравнение Фредгольма второго рода y/(x)-^K{x,T)y/(j)dT = Л(х,с) + /Ц(х), xQ е(-1,1), ^ х-х0

-, (39) л х-х0 « Vl-r2 (г —4Г> которое, в силу сделанного предположения, имеет единственное решение.

Аналогичные рассуждения сведения сингулярного интегрального уравнения к уравнению Фредгольма можно провести и для решений индексов к=-1, к=0. Только теперь при к=-1 условие разрешимости несколько усложняется для поиска числа Q.

Для численного решения уравнения (37) можно опять применить идеи метода дискретных вихрей и рассмотреть для нахождения решений индексов к = 1,-1,0 системы линейных алгебраических уравнений вида (25), (32) и (34) соответственно. Справедлива следующая

Теорема 3.1. Пусть в уравнении (37) функция /(*„) и ядро £(х0,х)принадлежат классу Н(а) на [-1,1] и [-1,1]* [-1,1] соответственно, а функция fs,q(x0) имеет вид (23) и пусть уравнение (23) имеет единственное решение в паре пространств Я*, для соответствующихТогда решения систем линейных алгебраических уравнений £K(XqjъХк).ГяМ.н =KXq)) + ftA(Хо.), j = i,.,(40)

Я 1 Х0 j ~Хк *=1

1 т^ГпМ-h + £K(x0J,xk)-yn(xk)-h =/(x0j) + fsq(x0J), j = 1,.,n, 71 *»1 X0j ~Xk 1

5>.(**)-й»С, (41)

A-l и соответствующие решения индекса к =0,1,-1 уравнения (37) связаны соотношением (27), в котором величина вп(хк) удовлетворяет неравенствам (28) и (29).

Если в задаче обтекания профиля моделировать профиль слоем диполей, то задача нахождения интенсивности слоя диполей сводится к следующему гиперсингулярному интегральному уравнению íi-^^Щт = f(xo), *0€(-1,1), (43) которое эквивалентно уравнению

1|1т,./(Хо), ,о£(и)> (44) я* 1 х0 -х при условии • (45)

Если в уравнении (43) справа стоит функция /д,ч(хо) (см.(37)), тогда уравнение (44) имеет своим решением функцию = ', х0 е (-1,1), (46) я V 1-х2 х-д а, следовательно, в силу равенства (45) уравнение (43) имеет своим решением функцию = -1п-,--, *0е(-1,1). (47)

Для численного нахождения решения (47) уравнения (43) в случае, когда в правой части стоит функция /д,ц(хо), опять воспользуемся идеей метода дискретных вихрей с использованием множества Е = {хк, к = 1,. .,п} и Е0 — {хор 3 — 0,1,., п}, построенных в § 1,гл. 3. Справедлива

Теорема 3.2. Пусть в правой части уравнения (43) стоит функция /д.ч(хо)> е (-1,1). Тогда решение системы линейных алгебраических уравнений

1 " я *=1

1 1 оу Хк+\ *0У Хк (*„;-<?), 3 = 1,2,., (48) и значения решения %(х), даваемые формулой (47) уравнения (43) связаны соотношением к(*о*)-8п(*ок)] * С* \ Хг > О, (49) где С не зависит от п.

В главе 4 рассматриваются схемы метода дискретных вихрей для сингулярных интегральных уравнений р ядром Гильберта первого рода в классе обобщенных функций и их обоснование.

В §1 рассматривается сингулярное интегральное уравнение с ядром Гильберта, которое мы получаем, моделируя телесный профиль вихревым слоем с отсосом внешнего потока, используя нетрадиционное условие непротекания

50) где в правой части есть дельта - функция (#0)=^(3> ~я\ 9а, Ч е [0;2;г]. произвольное действительное число, а дельта - функция б{вй-д) с носителем в точке # определяется условиями:

-«>-{*. V* • <51>

00, в0=д

В главе 1 показано, что уравнение (50), рассмотренное в обобщенных пространствах Нх, имеет решение только при выполнении условия:

2ло

При выполнении условия (53) решение уравнения (50) даётся формулой 1 2* где — \у{в\1в = С, в е [0;2я].

53)

54)

55)

Используя идеи метода дискретных вихрей, построим метод численного решения уравнения (50). Пусть отрезок [0;2л-] разбит двумя множествами точек £ = {<9*,£ = 1,2,.,«} и £0 = {0 ,у' = 1,2,.,«) таким образом,

2ж чтобы 0ы-0к=И = п к = 1,2,.,«; 0Л+1 =0,; 0оу =6»7+- и ?е£0, при этом д = 0ОЛ .

Теорема 4.1. Пусть в правой части уравнения (50) стоит сумма функций /(б>0)е н{а) и /¿ч(<90), [0,2л-]. Тогда между решением системы линейных алгебраических уравнений ./ = 1,2

2л- п

56) где

О, }*}ч

2 71

57) и решением уравнения (50), задаваемым формулой (54) при условии (55), выполняется соотношение гп(0к)-г(9к)<с^, 0<а<1. п

В §2 исследуется полное сингулярное интегральное уравнение, которое получается при моделировании обтекания произвольного телесного профиля вихревым слоем с отсосом внешнего потока с использованием нетрадиционного условия непротекания-. Задача нахождения интенсивности вихревого слоя сводится к решению сингулярного интегрального уравнения 2]к{в, ,в)у(е)с1в = Ж)+ Л.Д0о), в0,де [0,2л-]. (58)

Я о о

Его решение ищется в виде де[ 0,2л]

2л 2

Тогда для функции у/(в) уравнение (58) перейдет в уравнение

9. €[0,24 в-д

59)

60)

При этом условие разрешимости уравнения (58) в классе функций вида (59) получит вид

2<г,

61)

Для уравнения (58) предлагается следующий численный метод решения кЮЬ + ХФ»,.К)А = Ж,)+ К- -?), 1 я- *=1

Доказывается, что уп {вк)сходится в точках вк к точному решению.

Далее в §3 рассматривается" уравнение с гиперсингулярным периодическим ядром. Если в задаче обтекания профиля моделировать профиль слоем диполей, то задача нахождения слоя диполей сводится к следующему гиперсингулярному интегральному уравнению г\~~~в~в Л = №+&во ^ [0,2л]

ОбШ2 —

Уравнение (62) эквивалентно (в силу определения гиперсингулярного интеграла) уравнению

1 2} Iе*

2л $ 2

00 в 8'{вуе = г(ео)+Л,Д^о), в0е[0,2л] тдев0 е [о, 2л-]. Решение уравнения (62 )получаем в виде

62)

63) гИ = -|1п л I яп((<9-0„)/2)

БШ

0,12) л т{{в-д)!2) эш

Ь/2) С, д * 0.

64)

Для численного решения уравнения (62) вновь рассматриваем два множества точек Е и Е0, полагая д = е Е0. Используя идеи метода дискретных вихрей, построим метод численного решения уравнения (62). Заменим (62) следующей системой линейных алгебраических уравнений я *=1 2 1 \ к=\

Систему (65) можно записать в виде системы (56) относительно функции

2 л к=1 2

66)

4=1

Для однозначного нахождения ¿п (в0к) служат равенства к=\

Разрешая систему (66) и используя квадратурные формулы для интеграла с ядром Гильберта, имеем:

Лпп] 0-<а< 1.

Выводы:

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

1. Для сингулярного интеграла рассмотрены квадратурные формулы более отвечающие физической сущности метода дискретных вихрей и исследованы вопросы их сходимости.

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

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

Заключение.

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

1.Для сингулярного интеграла рассмотрены квадратурные формулы более отвечающие физической сущности метода дискретных вихрей и исследованы вопросы их сходимости.

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

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

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

2. Белоцерковский С.М., Скрипач Б.К. Табачников.В.Г. Крыло в нестационарном потоке газа,- М., Наука,-1987.-768 с.

3. Белоцерковский С.М. Ништ М.И. Отрывное и безотрывное обтекание тонких крыльев идеальной жидкостью.- М., Наука 1978.- 352 с.

4. Белоцерковский С.М. Лифанов И.К., Солдатов М.М. Метод дискретных особенностей в плоских задачах теории упругости / / ПММ.- 1988.- т. 47, вып. 5-с. 781-789.

5. Белоцерковский С.М., Лифанов И.К., Солдатов М.М. К численному решению задач теории упругости, ДАН СССР,- 1984, т. 277, №2-с.323-327.

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

7. Белоцерковский С.М., Лифанов И.К., Солдатов М.М. Метод дискретных особенностей в плоских задачах теории упругости с негладкой границей / / ПММ. -1987.-T. 51, вып. 2.- с.283-291.

8. Белоцерковский С. М., Лифанов И. К. О регуляризации численного решения сингулярных интегральных уравнений и нелинейных задач аэродинамики // Теория и методы решения неккоректно поставленных задач и их приложения. Новосибирск, 1978.- с. 173-179.

9. Белоцерковский С.М., Ништ М.И., ПономаревА.Т., Рысев О.В. Исследование парашютов и дельтопланов на ЭВМ.-М.: Машиностроение, 1987,-240 с.

10. Белоцерковский С.М., Лифанов И.К., Михайлов A.A. Расчет бесциркулярного обтекания произвольных тел / / Ученые записки ЦАГИ,-1987.-Т XVIII, №5-с.1-10

11. Белоцерковский С.М. Лифанов И.К. Некоторые сингулярные интегральные уравнения аэродинамики / / Дифференциальные уравнения.-1981 .-t.XVII, №9-с. 1539-1547

12. Березин И.С., Жидков Н. П., Методы вычислений, т. 1, 11.-М. Наука 1962.

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

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

15. Вайникко Г.М., Лебедева Н.В., Лифанов И.К. Численное решение сингулярого и гиперсингулярного интегрального уравнения первого рода на отрезке и дельта функция, Математический сборник, 2002.- т.193, №10, с. 1397-1410.

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

17. Габдулхаев Б.Г., Душков П.Н. О прямых методах решения сингулярных интегральных уравнений первого рода // Изв. вузов. Математика. -1973. №7.-с. 12-24

18. Габдулхаев Б.Г. Кубатурные формулы для многомерных сингулярных интегралов // Изв. вузов. Математика. 1975. - №4.-c.3-13

19. Габдулхаев Б.Г. Приближенное решение многомерных сингулярных интегральных уравнений, I // Изв. вузов. Математика.- 1975. -№7.-с.30-41

20. Габдулхаев Б.Г. Приближенное решение многомерных сингулярных интегральных уравнений, II // Изв. вузов. Математика. 1976. - №1.-с.30-41

21. Габдулхаев Б.Г., Онегов Л.А. Кубатурные формулы для сингулярных инп гралов // Изв. вузов. Математика.- 1976.-№7.-с. 100-105

22. Габдулхаев Б.Г. Конечномерные аппроксимации сингулярных интегралов и прямые методы решения особых интегральных и интегродифференциальных уравнений // Итоги науки и техники. Математ. анализ, т. 18- М.: ВИНИТИ АН CCCP.-1980.-c.251-307

23. Габдулхаев Б.Г. Оптимальные аппроксимации решений линейных задач Казань: Казан, ун-т. -1980.-232 с.

24. Габедава Г.В. Приближенное решение уравнения теории крыла методом Бубнова-Галеркина // Тр.мат.ин-та АН ГССР. -1974.-t.44.-с.52-56

25. Гагуа И.Б. О некотором применении кратных интегралов типа Коши // В кн.: Исследования по современным проблемам ТФКП. М.: Физматгиз. 1960.-с.345-352

26. Гандель Ю.В. К вопросу о решении одного класса парных уравнений математической физики // Харьков, университет. -1980, Деп. ВИНИТИ 13.04.81. -№ 1659-81.- 9 с.

27. Гандель Ю.В. О парных рядах Фурье некоторых смешанных краевых задач математической физики // ТФФА и их приложения, Харьков, Вища школа. 1982.-№38.-с. 15-18

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

29. Голубев В.В. Лекции по теории крыла, М, ГИТТЛ, 1949,- 480с.

30. Евсеева А.Ю. Об одной модификации метода дискретных вихрей численного решения сингулярных интегральных уравнений на отрезке. // Тр. IX МДОЗМФ 2000, - Орел, с.209-219

31. Гурса Э. Курс математического анализа, т.З- ч.2,-М.-Л., ГТТИ, 1934, -320с.

32. Кочин Н. Е., Кибель И. А., Розе Н.В., Теоретическая гидромеханника, ч. 1 -2.- М.- Физматиз, 1963.

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

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

35. Лифанов И.К.О сингулярных интегральных уравнениях с одномерными и кратными интегралами типа Коши // ДАН СССР.-1978.-т.239,№ 2.-е.265-268

36. Лифанов И.К. Формула Пуанкаре-Бертрана и сингулярные интегральные уравнения с кратными интегралами типа Коши // ДАН СССР. -1978. -т.243. № 1.-С.22-25

37. Лифанов И.К. О формулах обращения многомерных сингулярных интегралов// ДАН СССР-1979. -т.249, № 6.-с.1306-1309

38. Лифанов И.К. Топология кривых и численное решение СИУ-ий первого рода // IV Тираспольский симпозиум по общей топологии и ее приложениям.- 1979.-е. 82-85

39. Лифанов И.К. О методе дискретных вихрей // ПММ.-1979.-т.43, № 1.-с. 184-188

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

41. Лифанов И.К. О некорректности и регуляризации численного решения сингулярных интегральных уравнений первого рода // ДАН СССР. -1980.-t.255, № 5.-е. 1046-1050

42. Лифанов И.К. Квадратурные формулы и формула Пуанкаре-Бертрана для сингулярных интегралов // Сиб.мат.журн. -1980. t.XXI, № 6. -с.46-60

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

44. Лифанов И.К. О приближенном вычислении многомерных сингулярных интегралов // Семинар ин-та прикл. мат-ки им. И.Н.Векуа, Доклады.-1981 .-№ 15 .-с. 13-16.

45. Лифанов И. К., Тыртышников Е.Е.Теплицевы матрицы и сингулярные интегральные уравнения // сб. Вычислительные процессы и системы (под ред. Марчука Т.И.), М., Наука.-1990-вып.7 с.94-273.

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

47. Полтавский Л.Н., Математическое обоснование некоторых численных схем в аэродинамике. Докторская дисс. Москва, -1993.

48. Поляхов H.H. Об интегральном уравнении теории несущих поверхностей // Вести ЛГУ, Математика. Механика. Астрономия.- 1973,-№7.- с. 115-121.

49. Седов JI. И. Механика сплошной среды, ч. 1, М., Наука, -1973, 584с.

50. Belotserkovsky S.M. and Lifanov I.K., Metod of Diskrete Vortices. CRC Press, 1993.- 452p.

51. Gulyaev V. V., Lifanov I. K., Misko V. A. Mathematical model of the flow around airfll with energy high-lift devices // Russ. J. Numer. Anal. Math. Modelling. 1996. V. 11. №2. p. 155-1б"6.

52. Koenig D.G., Falarski M.D. Aerodinamic characteristics of a large scale model with a swept and augument jet flap// Nasa TM X62029.-1971

53. Vainikko G. Periodic integral and pseudodifferential equations.// Helsinki Univ. Technology, Institute of mathematics, 1997., Research report С13

54. Vainikko G. and Saranen J. Fast collocation solvers for integral equations on open arcs // J. Integral equations. Appl. 1999. — V.l l.No.l, p. 57-102.