close

Вход

Забыли?

вход по аккаунту

?

Математическое обоснование численного интегрирования несобственных интегралов встречающихся в аэродинамике.

код для вставкиСкачать
НАУЧНЫЙ ВЕСТНИК МГТУ ГА
серия Аэромеханика и прочность
2007
№ 111
УДК 533.6:681.3
МАТЕМАТИЧЕСКОЕ ОБОСНОВАНИЕ ЧИСЛЕННОГО ИНТЕГРИРОВАНИЯ
НЕСОБСТВЕННЫХ ИНТЕГРАЛОВ, ВСТРЕЧАЮЩИХСЯ В АЭРОДИНАМИКЕ
(краткое сообщение)
В.Ю. СМИРНОВ
Статья представлена доктором технических наук, профессором Барзиловичем Е.Ю.
Рассматривается численное интегрирование встречающихся в задачах аэродинамики несобственных интегралов в
смысле конечного значения по Адамару. Приводится выражение для интегралов по Адамару, не содержащее производных
под знаком интеграла, которое удобно для численного интегрирования. Приводится также процедура численного интегрирования несобственных интегралов.
Математические модели аэродинамической интерференции воздушного судна и грузов на внешней
подвеске содержат уравнения с несобственными интегралами. В частности, несобственные интегралы
возникают при использовании формулы Био-Савара-Лапласа. Как правило, они являются несобственными интегралами второго рода, которые не существуют ни по Риману, ни в смысле главного значения по
Коши, так как содержат сильную особенность. Такие интегралы впервые были описаны Адамаром [1].
Например, интеграл, возникающий при вычислении скорости в контрольной точке, расположенной внутри вихревого потока, имеет вид [3]:
0, 5 b j

K 
a (s )
f (s )

I 1 = ∫ 2 1+
ds = ∫
ds,
2


2
2
− 0,5b j s
− 0 ,5 b j s
a
(
s
)
+
b
(
s
)


2
где K = cos(ψ j − ψ i ) / sin χ1 j , a (s) = ξ i − ξ j − s cos χ1 j , b(s) = s sin χ1 j .
0,5b j
(1)
Интеграл (1) является частным случаем несобственного интеграла вида
b
*∫
a
v(x)
xn
dx ,
(2)
где n – некоторое натуральное число; v(x) – функция, которая имеет в некоторой окрестности точки 0
производную порядка n и интегрируема по Риману на отрезке [a, b].
На основе доказательства теорем существования [5, 6] можно определить интеграл в смысле конечного значения по Адамару вида (2) следующим образом:
b
*∫
a
v(x)
xn
n − 2 1 v ( i ) (0) 1− ( −1) n − i −1
v(x)
dx = lim ( ∫ n dx − ∑
).
ε→0 L x
ε n −i −1
i = 0 i! n − i −1
(3)
С помощью интегрирования по частям интеграл в смысле конечного значения по Адамару типа
b
*∫
a
v(x )
(x − x 0 ) n
dx
(4)
можно представить в следующем виде [5]:
b
*∫
a
v(x )
(x − x 0 )
n −2
1
v (i) (x )
n − i −1
i = 0 ( n −1)( n − 2)...( n − i −1) ( x − x 0 )
dx = − ∑
n
1  ( n −1)
ν
( x ) ln x − x 0
( n −1)! 
b
a
b
− ∫ ln x − x 0 ν
a
(n)

( x )dx .


b
a
+
(5)
Несобственные интегралы, применяемые в математических моделях аэродинамики и в других
задачах математической физики, вычисляются, как правило, с помощью численных методов. Однако
при таких вычислениях возникают существенные трудности.
81
Математическое обоснование численного интегрирования несобственных интегралов...
Например, формула для вычисления несобственного интеграла вида (5) имеет существенный недостаток, который состоит в том, что при вычислении на ЭВМ необходимо считать интегралы от
производных функции v(x) до n-го порядка включительно. Эта процедура является чрезвычайно трудоемкой.
Поэтому хотелось бы иметь выражение для интеграла в смысле конечного значения по Адамару
типа (4), которое не содержало бы интегралов от производных.
Для того чтобы избавиться от этого недостатка, преобразуем формулу (5) так, чтобы она содержала только значения производных в точке x0 и на концах отрезка.
Для этого вычислим следующий интеграл:
n −1 v ( i ) ( x
v(x ) − ∑
b
i!
i =0
I=∫
(x − x 0 ) n
a
=−
0 )( x − x 0 )
i
n −1 v ( i ) ( x )( x − x ) i
1 b
1
0
0
dx = −
(v(x) − ∑
)d
=
∫
n −1 a
i!
( x − x 0 ) n −1
i =0
n −1 v ( i ) ( x )( x − x ) i
1 
1
0
0
(
v
(
x
)
−
)
∑
n −1 
i!
( x − x 0 ) n −1
i =0
b
a
n −1 v ( i ) ( x
0 )( x − x 0 )
i =0
(i −1)!
v ′( x ) − ∑
b
−∫
i −1
( x − x 0 ) n −1
a
dx.
Проинтегрировав по частям n-1 раз, получим
n −2
n −1 v ( i ) ( x )( x − x ) i − j
1
1
1
0
0
I= −∑
...
( v ( j) ( x ) − ∑
)
+
 j=0 n −1 n − j−1
(i − j)!
( x − x 0 ) n − j−1
i= j
b
1
+
( v ( n −1) ( x ) − v ( n −1) ( x 0 )) 

( n −1)!
b
− ∫ v ( n ) ( x ) ln x − x 0 dx.
a
a
Окончательно имеем следующее выражение для интеграла в смысле конечного значения по
Адамару типа (4):
n −1
b
*∫
a
n −2
−∑
j= 0
v(x ) − ∑
b
v(x)
(x − x 0 )
n
i=0
dx = ∫
v ( i ) ( x 0 )( x − x 0 ) i
i!
(x − x 0 ) n
a
dx −
n −1 v ( i ) ( x )( x − x ) i − j
1
1
0
0
∑
( n −1)...( n − j −1) i = j
(i − j)!
( x − x 0 ) n − j−1
1
−
v ( n −1) ( x 0 ) ln x − x 0
( n −1)!
b
−
(6)
a
b
.
a
В качестве иллюстрации рассмотрим интеграл, встречающийся в задачах численного моделирования аэродинамической интерференции воздушного судна и грузов на внешних подвесках [3, 4], вида
b
v ′( x )
∫ x −x
a
dx .
(7)
0
Имеем
b
I=∫
v(x ) − v(x 0 )
(x − x 0 )
a
2
b
Отсюда следует, что
v(x) − v(x 0 )
1
dx = − ∫ ( v ( x ) − v ( x 0 ))d
=
x −x 0
x −x 0
a
v ′( x )
∫ x−x
a
b
0
b
dx = ∫
a
v(x ) − v(x 0 )
(x − x 0 ) 2
dx +
a
v ′( x )
dx .
a x−x0
b
+∫
v ( b) − v ( x 0 ) v ( a ) − v ( x 0 )
−
.
b−x 0
a−x0
c
v ′( x )
v ( x ) − v (0)
2 v (0) v (c) − v ( −c)
dx = ∫
dx −
Например, ∫
+
. Тогда
2
c
c
x
−c x
−c
c
b
82
В.Ю. Смирнов
c
∗∫
−c
v(x )
x
2
dx =
c
∫
v ( x ) − v ( 0)
x
−c
2
dx −
2 v (0)
.
c
Полученное выражение не содержит производных функции v(x), но интеграл в правой части существует только в смысле главного значения, поэтому для вычислений лучше взять интеграл
b
v ′( x )
∫ x − x dx = ∫ v ′( x )d ln x − x 0 = v ′( x )ln x − x 0
0
a
a
b
b
a
b
− ∫ v ′′( x )ln x − x 0 dx .
(8)
a
Рассмотрим
b
∫
v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 )
a
=−
=
b
dx = − ∫ ( v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 ))d
(x − x 0 ) 2
v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 )
x−x0
v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 )
x−x0
a
b
a
b
b
+∫
v ′( x ) − v ′( x 0 )
x−x0
a
dx =
+ ( v ′( x ) − v ′( x 0 )) ln x − x 0
a
1
=
x−x0
b
b
− ∫ v ′′( x ) ln x − x 0 dx.
a
a
Отсюда
∫
v ′′( x ) ln x − x 0 dx = − ∫
b
v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 )
a
a
(x − x 0 ) 2
b
−
v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 )
b
x−x0
a
dx −
+ ( v ′( x ) − v ′( x 0 )) ln x − x 0
b
.
a
Подставляя полученное выражение в соотношение (8), найдем представление для интеграла
вида (7)
b
v ′( x )
b
dx
=
∫ x−x
∫ v ′( x )d ln x − x 0 = v ′( x ) ln x − x 0 a +
0
a
a
b
v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 )
+∫
dx +
(x − x 0 ) 2
a
v ( x ) − v ( x 0 ) − v ′( x 0 )( x − x 0 ) b
+
− ( v ′( x ) − v ′( x 0 )) ln x − x 0
x−x0
a
b
(9)
b
.
a
Выражение (9) для интеграла (7), также как и выражение (6) для интеграла (4) не содержит под
знаком интеграла производных функции v(x) кроме значений производных в некоторых известных и
фиксированных точках, и интеграл существует в обычном смысле.
Таким образом, для несобственного интеграла получено выражение, не содержащее производных под знаком интеграла, которое удобно для численного интегрирования.
Представление несобственных интегралов в виде (6) и (9) позволяет реализовать процедуру их вычисления на основе численных методов. Эта процедура заключается в следующем:
– одним из численных методов вычисления определенных интегралов (прямоугольников, трапеций,
Симпсона, Котеса и т.д.) вычисляется значение собственного интеграла, входящего в правую часть формул (6) или (9). Причем значение интеграла вычисляется вне некоторой δ-окрестности особой точки, т.е.
вне отрезка [x0-δ; x0+δ], чтобы исключить возникновение программного деления на ноль;
– полное значение несобственного интеграла вычисляется по формулам (6) или (9).
Значение δ определяется в процессе вычисления интеграла из условия, что знаменатель подынтегрального выражения меньше некоторой заданной величины. Проведенные методические исследования [3, 4] показали, что целесообразно выбирать эту величину равной 10-6.
Полученные выше результаты и процедура численного интегрирования несобственных интегралов носят универсальный характер, так как могут найти применение не только в аэродинамике, но и в
других областях науки, где встречаются сингулярные интегралы.
83
Математическое обоснование численного интегрирования несобственных интегралов...
ЛИТЕРАТУРА
1. Адамар Ж. Задачи Коши для линейных уравнений с частными производными гиперболического типа. – М.: Наука,
1978.
2. Лифанов И.К. Метод сингулярных интегральных уравнений и численный эксперимент. – М.: ТОО "Янус", 1995.
3. Смирнов В.Ю. Расчет нестационарных аэродинамических характеристик грузов на дозвуковых скоростях полета //
Установки и системы управления авиационным вооружением. – М.: ВВИА им. Н.Е. Жуковского, 1994.
4. Смирнов В.Ю. Расчет линейных стационарных и нестационарных аэродинамических характеристик АСП на дозвуковых и сверхзвуковых скоростях полета // Установки и системы управления авиационным вооружением. – М.: ВВИА им.
Н.Е. Жуковского, 1994.
5. Смирнов В.Ю. Математическое обоснование вычисления несобственных интегралов, встречающихся в аэродинамике, для целых степеней особенности // Научный Вестник МГТУ ГА. Сер. Аэромеханика и прочность. – 2006. № 97 (1).
С. 76 – 79.
6. Смирнов В.Ю. Математическое обоснование вычисления несобственных интегралов, встречающихся в аэродинамике, для дробных степеней особенности // Научный Вестник МГТУ ГА. Сер. Аэромеханика и прочность. – 2006. № 97 (1).
С. 80 – 82.
MATHEMATICAL SUBSTANTIATION OF NUMERICAL INTEGRATION OF IMPROPER
INTEGRALS THAT ARE FOUND IN AERODYNAMICS
Smirnov V.J.
The numerical integration of improper integrals in terms of final value by Hadamard that are found in aerodynamics is considered. Expression for integrals on the Hadamard that doesn’t contain integrand derivatives that is convenient for a numerical integration is obtained.
Procedure of a numerical integration of improper integrals is obtained also.
Сведения об авторе
Смирнов Владимир Юрьевич, 1963 г.р., окончил МГУ (1985) и ВВИА им. проф. Н.Е. Жуковского (1987),
кандидат технических наук, доцент, заведующий кафедрой МАИ, автор более 50 научных работ, область научных интересов – аэродинамическая интерференция самолета и грузов, эксплуатация сложных технических систем, компьютерные технологии.
1/--страниц
Пожаловаться на содержимое документа