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

ISSN 1992-6502 (Print)_

2015. т. 19, № 2 (68). с. 179-182

Ъюьшм QjrAQllQj

ISSN 7775-7789 (Online) http://journal.ugatu.ac.ru

удк 519.6

Модификация методов численного интегрирования для решения осесимметричных задач

1 2 3

В. П. Житников , О. Р. Зиннатуллина , Н. И. Житникова

1 zhitnik@ugatu.ac.ru, 2 olga_zr@mail.ru, 3 zhitnikova_ni@mail.ru ФГБОУ ВПО «Уфимский государственный авиационный технический университет» (УГАТУ) Поступила в редакцию 25 февраля 2015 г.

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

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

ВВЕДЕНИЕ

Вычисление продольной и радиальной составляющих напряженности осесимметричного поля сводится к численному интегрированию аналитической функции комплексного переменного / {2) по формулам Г. Н. Положего [1] в узловых точках коллокации 2(ск), к = 1, ..., N

Im°f f [z(а)]_—_

0 л1(z(a) - z(ак ))(z() - z(ak)) '

(1)

1 Im Y f [[№()- Re z(ok ))) (2)

Im z(ak ) 0 V()- z(ak ))()- z(ak ^ '

0

Подынтегральные функции (1), (2) содержат особенности, которые требуют применения специальных численных методов. Обычно численное интегрирование производится с помощью квадратурной формулы наивысшей алгебраической степени точности (по формуле Гаусса [2]).

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

1. ОЦЕНКА ПОГРЕШНОСТИ ВЫЧИСЛЕНИЯ ИНТЕГРАЛОВ ПОЛОЖЕГО ПРИ ПРЯМОМ ПРИМЕНЕНИИ КВАДРАТУРНЫХ ФОРМУЛ ГАУССА

При вычислении интегралов Г. Н. Положего (1), (2) отрезок интегрирования с е [0, стах] разбивается на п отрезков узловыми точками ст. (рис. 1, а).

Рис. 1. Разбиение отрезка интегрирования

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

а

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

Рассмотрим интеграл, подобный (1),

(с- *о)(а + ¿о)

(3)

при ¿0 = 2к = с2. Такой случай вызывается необходимостью вычисления значений интегралов в каждой узловой точке, в том числе и при с=с2 (рис. 1, б).

Производя замену переменных

, х + 1 , к ,

t = к-, т = — ах

2 2

(при этом отрезок интегрирования перейдет в [-1, 1]), получим следующий интеграл

к

I

-1

ах

\х+! - 2к т кх+1 + 2к

2

2

1 ах %

= — 11-= — 1—

-^(3 - х)(х + 5) 6

(4)

Таким образом, значение представленного интеграла не зависит от величины шага к. Согласно [2], не зависит от к и погрешность вычисления интеграла. Сходимость к точному решению при к ^ 0 отсутствует.

Ухудшение сходимости может также вызываться тем, что увеличение точности в алгоритме решения общей задачи достигается путем уменьшения шага разбиения. При этом особенность подынтегральной функции при z=z0 учитывается только при интегрировании на ближайшем к особой точке отрезке. При уменьшении шага уменьшается и величина этого отрезка (рис. 1, а). На большей части отрезков особенность не учитывается.

2. ЧИСЛЕННЫЕ МЕТОДЫ ВЫЧИСЛЕНИЯ ИНТЕГРАЛОВ ПОЛОЖЕГО

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

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

По аналогии с методом Гаусса потребуем выполнения следующих условий

1

А1 + А2 + А3 = | х

0

ах

0 л](ь - х)(х + а)

А^1 + А2 х2 + Азхз =

ах

= А), (5)

= А

0 л\(Ъ - х)(х + а)

А^2 + А2 + А3 х2 = | х2

ах

где

0 л1Ф - х)(х + а)

.а - Ь + 2 .а - Ь В0 = агс81п--агс81п-,

а + Ь а + Ь

В1 =,[аЪ - 1)(а + 1) - ^В0.

= В2

В 2 =1 3

(3 а-Ь - 2- 3 ~1аЬ +

+

(а + Ь)2 (а - Ь)2 8 4

В0.

Решение этой системы 3 линейных относительно Aj уравнений запишется в виде

А3 =

Ао =

В2 - (х1 + х2 )В1 + х1х2 В0

(х3 - х2 )(х3 - х1) В1 - х1В0 -(х3 - х1 ))3

х2 - х1

А1 = В0 - А2 - А3 . Значения узлов квадратурной формулы

возьмем из [2] при п=3: ~1 = -

х3 =

~2 = 0

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

В. П. Житников, О. Р. Зиннатуллина и др. • Модификация методов

181

нашей квадратурной формулы:

х + 1

х =

х1 =

1 -

, х2 = — , х3 = — 2 2 3 2

1 +

Для оценки погрешности и ее составляющих использовался метод численной фильтрации, изложенный в [3, 4]. На рис. 2, 3 оценки погрешности представлены на графике в виде зависимости - ^ А (десятичного логарифма относительной погрешности) от числа узлов кол-локации п в логарифмической шкале.

Рис. 2. Оценка погрешности первого метода

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

Результаты исследования этого метода (рис. 2) показывают, что при постоянном шаге он обладает 3-м порядком точности и имеет только нечетные составляющие.

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

1 -с т-1 ~ 1 с т ~

2 С / (с)ас+ 2 I / (с)ас =

— с т с т-1

1 ст ~ 1 ст ~ ст ~

= 2 | /(-с)ас + 2 | /(с)ас= | /(с)ас. (6)

ст-1 Ст-1 Ст-1

Коэффициенты А]- удовлетворяют системе 3 линейных алгебраических уравнений

ах

А1 + А2 + А3 = |{ х + , =

1 2 3 01 2 ) 4(Ь - х)(х + а)

= В0, (7)

А х1 + I + Аз! х2 +

а - Ь

+ А3 (х3 +

+ ■

- Ь V 1Г а - Ь

АЛ х1 +

= I! х

ёх

2 ) Л (Ь - х)(х + а)

= В 2

а-Ь 2

+ АЛ х2 +

а - Ь 2

+ А3 (х3 +

+ а - Ь = |/ + а - Ь V ах = в + = ¿Г + д/(Ь - х)(х + а) = '

4 :

Здесь

.а - Ь + 2 .а - Ь В0 = агс81п-;--агс81п

а+Ь

а - Ь ¡-г а - Ь + 2 ( )(-^ (а + Ь)

В2 = ^аЬ--4-# - 1)(а + 1) + B0,

а+Ь

В 4 =

(а - Ь)3 + 3 (а + Ь)2 (а - Ь)

аЬ 32

(а - Ь + 2)3 + -2 (а + Ь)2 (а - Ь + 2)

-,.У(Ь - 1)(а +1) + 3 (±*£ В0

32 128 0'

Решение этой системы уравнений запишется

в виде

С1 =|х1 + Ь] , С2 = {х2 + ^у^) =

С3 =| х3 +

а - Ь

2

= В2 - С1В0 -(С3 - С1 )В2

А2 = С2 - С1 ,

А = В4 -(С1 + С2 )В2 + С1С2 В0 3 ((3 - С2 )(( - С1) ' А1 = В0 - А2 - А3.

Значения узлов квадратурной формулы возьмем, как и выше, из [2] при п=3. Результаты исследования этого метода (рис. 3) показывают,

что он обладает 3-^-м порядком точности, порядок остальных составляющих увеличивается последовательно на 1.

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

2

2

2

1

1

2

а

2

2

4

4

2

-IgA

о _,_,_,_,_

1 2 3 4 J lg п 6

Рис. 3. Оценка погрешности второго метода

ЗАКЛЮЧЕНИЕ

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

Эти методы были использованы для решения осесимметричных задач электрохимического формообразования. Применение этих методов позволило исследовать тонкие эффекты формообразования обрабатываемой поверхности [5].

СПИСОК ЛИТЕРАТУРЫ

1. Положий Г. Н. Обобщение теории аналитических функций комплексного переменного. Киев: Киевск. ун-т, 1965. 442 с. [[ G. N. Polozhiy, Generalization of analytical functions of complex variable theory, (in Russian). Kiev: Kiev University, 1965. ]]

2. Крылов В. И. Приближенное вычисление интегралов. М.: Наука. 1967. 255 с. [[ V. I. Krylov, Approximate integrals calculation, (in Russian). Moscow: Nauka, 1979. ]]

3. Житников В. П., Шерыхалина Н. М. Моделирование течений весомой жидкости с применением методов многокомпонентного анализа. Уфа: Гилем, 2009. 336 с. [[ V. P. Zhitnikov, N. M. Sherykhalina, Modeling of weighted fluid flows using the methods of multi-component analysis, (in Russian). Ufa: Gilem, 2009. ]]

4. Житников В. П., Шерыхалина Н. М. Применение многократной фильтрации при численном решении задач методами теории функций комплексного переменного // Вычислительные технологии. 2013. Т. 18, № 1. С. 15-24. [[ V. P. Zhitnikov, N. M. Sherykhalina, "Using repeated filtration for problems numerical solution by complex variable theory methods," (in Russian), in Computational technologies, vol. 18, no. 1, pp. 15-24, 2013. ]]

5. Житников В. П., Зиннатуллина О. Р., Пореч-ный С. С., Шерыхалина Н. М. Особенности установления предельных решений нестационарных осесимметричных задач Хеле-Шоу // ПМТФ. 2009. Т. 50, № 4. С. 87-99. [[ V. P. Zhitnikov, O. R. Zinnatullina, S. S. Porechny, N. M. She-

rykhalina, "Determining the limiting solutions of nonstationary axisymmetric Hele-Shaw problems," in Journal of Applied Mechanics and Technical Physics, vol. 50, no. 4, pp. 617-627, 2009. ]]

ОБ АВТОРАХ

ЖИТНИКОВ Владимир Павлович, проф. каф. высокопро-изв. выч. технол. и систем. Дипл. инж.-физ. (МФТИ, 1973). Д-р физ.-мат. наук по мех. жидкости, газа и плазмы (Ка-занск. ун-т, 1993). Засл. деят. науки РБ. Иссл. в обл. волн. течений жидкости, э/хим. формообразования, числ.-аналит. методов.

ЗИННАТУЛЛИНА Ольга Рифовна, доц. каф. выч. мат. и кибернетики. Дипл. инж.-с/аналитик (УГАТУ, 2002). Канд. физ.-мат. наук по мех. жидкости, газа и плазмы (БашГУ, 2006). Иссл. в обл. моделирования э/хим. формообразования, числ.-аналит. методов, оценки погрешн. и достоверн. числ. результатов.

ЖИТНИКОВА Наталья Ивановна, доц. каф. выч. мат. и кибернетики. Дипл. инж.-мех. (УАИ, 1972). Канд. техн. наук по дин. и прочн. машин, приборов и аппаратуры (УГАТУ, 1994). Иссл. в обл. мех. оболочек, числ. методов.

METADATA

Title: Numerical methods of integration modification for ax-

isymmetric problems-solution. Authors: V. P. Zhitnikov, O. R. Zinnatullina, N. I. Zhitnikova Affiliation:

Ufa State Aviation Technical University (UGATU), Russia. Email: 1 zhitnik@ugatu.ac.ru. Language: Russian.

Source: Vestnik UGATU, vol. 19, no. 2 (68), pp. 179-182, 2015.

ISSN 2225-2789 (Online), ISSN 1992-6502 (Print). Abstract: The adaptive numerical methods are suggested for effective solution of the problems of axisymmetric fields modeling using the integral transformations of analytical functions of complex variable. The methods are based on quadrature formulas of Gaussian type with varying node points location, which is determined by integrated function singularities. With the help of numerical filtration of results obtained with different nods number the efficiency of methods suggested is shown. Key words: axisymmetric field; conformal mapping; integral transformation; adaptive quadrature formula; numerical filtration.

About authors:

ZHITNIKOV, Vladimir Pavlovich, Prof., Dept. of Common Sciens. Dipl. Engineer-physicist (Moscow Physical-Technical Inst., 1973). Cand. of Phys.-Math. Sci. (MIPT, 1984), Dr. of Phys.-Math. Sci. (KSU, 1993). ZINNATULLINA, Olga Rifovna, Docent, Dept. of Computer Sciens. Dipl. Engineer-system master (UGATU, 2002). Cand. of Phys.-Math. Sci. (BGU, 2006). ZHITNIKOVA, Natalya Ivanovna, Docent, Dept. of Computer Sciens. Dipl. Engineer-mech. (UGATU, 1972). Cand. of Techn. Sci. (UGATU, 1994).