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

УДК 532.59

ПРИМЕНЕНИЕ ПОДХОДА ЛАГРАНЖА К РЕШЕНИЮ ОДНОМЕРНОЙ ЗАДАЧИ РАСПРОСТРАНЕНИЯ ВОЛН В ГАЗЕ В ЛИНЕЙНОЙ ПОСТАНОВКЕ

© 2009 В. В. Никонов, В. Г. Шахов Самарский государственный аэрокосмический университет

Рассматривается моделирование распространения одномерных волн в газе с помощью различных численных схем. Наряду с известными схемами предлагается схема решения задачи акустики, использующая подход Лагранжа к описанию движения среды. Результаты численного решения для трёх типов начальных условий сравниваются с аналитическим решением. Делаются выводы о применимости рассмотренных схем. Показано, что для предлагаемой схемы можно варьировать шагом по времени, не снижая точности решения.

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

1. Математическая формулировка задачи

Уравнения, описывающие одномерную задачу линейной акустики, получаются из уравнений Навье-Стокса при отбрасывании конвекционных и вязких членов. Закон изменения давления считается адиабатическим, в результате в размерных переменных будем иметь [1]

* и и _

+ Р --------г = 0.

и Х

л * *2 л *

ии с ир

----------1------------— = о

~ * 1 * ^ *

ОХ р их

(1)

Здесь и* - скорость, р* - плотность, х* -координата, Х* - время, с*2 = кр* / р* - квадрат скорости звука, р* - давление, к - показатель адиабаты.

Введя следующие безразмерные переменные: и = и*/и0, р = р*/р0, х = х*/Ь , Х = Х*Ь/и0, с = с* /и0, где Ь - характерный

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

5р ди _

— + — = 0, 5х и х

и и

ир

и Х

и х

(2)

где с2 = кР^(р0и2).

Систему (2), сделав следующую замену переменных

Wl = и + с0р,

(3)

можно привести к системе уравнении, аналогичных уравнениям конвекции (переноса)

[1]:

и Х

и w2

и Х

--с0

и х

= 0.

(4)

Система (4) имеет [1] следующее аналитическое решение:

и(Х, х) = 2 [и0 (х - с0Х) + и0(х + с0Х)] +

+ [р0(х - с0Х)-р0(х + с0Х) ],

p(t, х) = -1- [и0 (х - c0t) - u0(x + c0t)] + zc0

1 Г -| (5)

+ 2[p0(x - c0t) + p0(x + c0t)] ,

где нулевой верхний индекс указывает на значения переменных из начальных условий (НУ) задачи.

2. Численные схемы применяемых методов

2.1. Применение центральной конечно-разностной схемы для решения уравнений линейной акустики в переменных скорость-плотность

Будем использовать «шахматную» сетку, когда значения плотности находятся в центрах ячеек, а значения скорости - на границах между ячейками. Значения величин плотности и скорости в следующий момент времени при применении центральной конечно-разностной (ЦКР) схемы к системе (2) определяются следующим образом:

pj = pj-1 —(uj-1 — uj-1 )

Hi Hi h 05 1—0.5/j

tc2

ul+0.5 = U 1—0.5 -^(pl+l -pi-1), (6)

h

где D t - шаг по времени, h - шаг расчётной однородной сетки по пространственной координате.

Граничные условия (ГУ) на твердой стенке: если граница 1+0.5 ячейки совпадает со стенкой, то

uj-0.5 = 0. (7)

2.2. Применение метода Годунова для решения уравнений линейной акустики в переменных скорость-плотность

В методе Годунова [2] также используется «шахматная» сетка, только в центрах ячеек определяются «малые» переменные, а на границах ячеек - «большие» переменные. Значения величин плотности и скорости в следующий момент времени при применении метода Годунова находятся как

Если центр ячейки і+1 находится «внутри» тела, то ГУ на стенке определяются следующим образом:

<ї р:- =рі-1. (10)

При этом ГУ (10) подставляются в (9).

2.3. Применение схемы вперед по потоку для решения уравнений линейной акустики в форме конвективных уравнений

К системе уравнений (4) удобно применить конечно-разностную схему «вперёд по потоку» (ВП), в результате будем иметь

< = <-1 -^-К-' -

w!:=w,f +^с^2';;- W2І-'). с)

ГУ на стенке накладываются согласно (10), после чего они подставляются в (3).

2.4. Применение подхода Лагранжа для решения уравнений линейной акустики в форме конвективных уравнений

Систему уравнений (4) можно также решить с применением подхода Лагранжа к рассмотрению движения среды. Здесь предлагается использовать метод, применяемый для моделирования конвекции в методе «вихрь в ячейке» (ВЯ) [3]. При этом моделируется движение двух волн. Одна распрост-

13S

раняется вправо и переносит величину '1 . При этом в следующий момент времени

*:

'и

W'i = whi-1, хі = хі-1 + с0 дХ •

(12)

*: : 1 х: — х: 1 с Х

W9i = W9i , хі = хі с0 дХ

(13)

Wlk = w*iЛ(хк-1 - хі)

(14)

Л 0(2) =

Ч -0.5 < г/Ь < 0.5

0, г/Ь <-0.5 Уг/Ь > 0.5. (15)

W'

W,

= w;:5 0.5(хі-1 + хі) < хк-1 < 0.5(хі + хі+1),

= w2:, 0.5(х:_'+х:)<хк-1 <0.5(х:+х:+1).

(16)

Рассмотрим ГУ на стенке для правой волны. Если ячейка при своём движении (12) оказывается «внутри» тела, то она отражается от его поверхности, и её координата определяется как

Другая волна распространяется влево и переносит величину '2, причём

х: = 2хЬ - хі - с0 д^

(17)

Здесь w1*: и w2: обозначают величины W' и w2 на искажённой после перемещения сетке. Для того, чтобы получить величины W'

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

где хь - координата стенки тела. При этом

правая волна превращается в левую волну и в точке с координатой (17) оказывается уже величина

^: = -иҐ - с0рҐ

(18)

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

Равенство (18) получается при подстановке ГУ (10) в (3). После этого полученная в соответствии с (18) величина перераспределяется в ячейки расчётной сетки согласно (14) или (16). ГУ на стенке для левой волны получается аналогично.

3. Тестирование методов

Описанные выше четыре численных схемы тестировались на трёх задачах, отличающихся НУ и ГУ

В первой тестовой задаче в точке НУ соответствуют однородному потоку

и0 = 1, р0 = 1,

(19)

Также возможен другой подход к перераспределению значений '1 и '2 в ячейки сетки. При этом временно строится искажённая после перемещения (12) или (13) частиц сетка и определяются границы её ячеек. Так для 1-ой искажённой ячейки в момент времени ] границы равны

{0.5(х1-1 + х1),0.5(х1+ х1+1)} .

Тогда значения '1 и '2 на однородной основной сетке в момент времени ] определятся как

а с момента времени 1 > 0 в точке х = 0 ставится непроницаемая перегородка.

НУ второй тестовой задачи задавались следующим образом:

0 Г\

и = 0,

р0 =■

1, х < 0.5

1.0002, х > 0.5.

(20)

а твёрдые границы отсутствовали.

В третьей тестовой задаче твёрдые границы отсутствовали, а НУ задавались как

и0 =

1, 0.4 < х < 0.6 0, х < 0.4 У х > 0.6.

р0 = 1.

(21)

В тестовых расчетах использовались конечные области. В первой задаче область

принималась равной х е , а во второй и третьей задачах области принимались равными х е [ 0.1].

Численное решение сравнивалось с аналитическим (5). Аналитическое решение первой задачи (19) имеет следующий вид:

и(1:,х) =

1,

х > М |х| <

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

Результаты, показанные на рисунках 1-5, получены после 20 шагов по времени

(1 = 5.76195268-10-4) для с0 = 347.105,

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

Погрешность решения определялась следующим образом:

Р(х,х) =

1, х > е01

1—

0 < х < е^

1 +-----, — еД < х < 0.

(22)

і ь

^ = ка ,

е0

(23)

называемому критерием Куранта-Фридриха-Леви [4], с величиной коэффициента ка = 1. Заметим, что для предложенной схемы, ис-

8„ = тах

и, — и

и

•100%.

(24)

В результате расчетов выяснилось, что наилучшие результаты для первых трёх рассмотренных методов получаются, если шаг по времени удовлетворяет выражению

где иех - значение скорости аналитического решения.

После 20 шагов по времени ошибка численного решения (24) составляет менее 0.002 % для результатов, показанных на рисунках 1, 3, 5. Графики численного решения, приведенные на рисунках 2 и 4, носят «пилообразный» характер, и ошибка (24) составляет 100 %.

При изменении шага сетки погрешность решения не меняется, что объясняется линейностью задач и выбором шага по времени в виде (23).

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

для ка, равных 1, 2, 4 и 20. При этом по-

0

0

1.005

1.004

1.003

1.002

1.001

1

0.999

0.998

0.997

0.996

0.995

-0.6 -0.4 -0.2 0 0.2 0.4 0.6

-0.6 -0.4 -0.2 0 0.2 0.4 0.6

Рис. 1. Распределение плотности и скорости в задаче (19) для всех рассмотренных методов:

- численное решение,____________________________- аналитическое решение

Р 1.0003 1.00025 1.0002 1.00015 1.0001 1.00005 1

0.99995

0.9999

0 0.2 0.4 0.6 0.8 1х

Рис. 2. Распределение плотности и скорости в задаче (20) для метода, использующего ЦКР схему в прямой постановке:

- численное решение, - аналитическое решение

Р 1.0003 1.00025 1.0002 1.00015 1.0001 1.00005 1

0.99995

0.9999

и 0.01

0

-0.01

-0.02

-0.03

-0.04

0.2 0.4 0.6 0.8 1

0.2 0.4 0.6 0.8 1

,х ч, , х

Рис. 3. Распределение плотности и скорости в задаче (20) для метода Годунова; метода, использующего ВП схему, и метода, использующего подход Лагранжа:

- численное решение, - аналитическое решение

1.005

1.004

1.003

1.002

1.001

1

0.999

0.998

0.997

0.996

0.995

и

0.2 0.4 0.6 0.8 1

1

0.8

0.6

0.4

0.2

0

0.2 0.4 0.6 0.8 1

Рис. 4. Распределение плотности и скорости в задаче (21) для метода, использующего ЦКР схему в прямой постановке:

- численное решение, - аналитическое решение

1.005

1.004

1.003

1.002

1.001

1

0.999

0.998

0.997

0.996

0.995

и

0.2 0.4 0.6 0.8 1

Рис. 5. Распределение плотности и скорости в задаче (21) для метода Годунова; метода, использующего ВП схему, и метода, использующего подход Лагранжа:

- численное решение, - аналитическое решение

грешность численного решения незначительно уменьшалась при увеличении ка. Это объясняется меньшим количеством выполняемых арифметических операций (в том числе и перераспределения в ячейки сетки). Для рассмотренных задач применение подходов (14) и (16) распределения значений в ячейки сетки показало одинаковые результаты.

В заключение можно сделать следующие выводы:

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

2. Метод Годунова, метод, использующий конечно-разностную схему вперед по потоку, и предложенный метод, использующий подход Лагранжа, адекватно моделируют распространение волн в одномерных задачах линейной акустики.

3. Предложенный метод позволяет варьировать в широких пределах шагом по времени.

Библиографический список

1. Годунов, С.К. Уравнения математической физики [Текст] / С.К. Годунов. - М.: Наука, 1971. - 416 с.

2. Годунов, С.К. Численное решение многомерных задач газовой динамики [Текст] / С.К. Годунов, А.В. Забродин, М.Я. Иванов, АН. Крайко, ГП. Прокопов. -М.: Наука, 1976. - 400 с.

3. Никонов, В.В. Моделирование двумерного ламинарного пограничного слоя с помощью метода “вихрь в ячейке” [Текст] / В.В. Никонов, В.Г. Шахов // Вестник СГАУ

- Самара. - 2006, № 3 (11). - С. 25 - 30.

4. Ferziger, J. Computational methods for fluid dynamics [Text] / J. Ferziger, M. Peric. -3 rev. ed. - Springer-Verlag. - 2002. - 423 p.

References

1. Godunov, S.K. Mathematical physics equations [Text] / S.K. Godunov. - M.: Science.

- 1971. - 416 p.

2. Godunov, S.K. Numerical solution of multidimensional problems of gas dynamics [Text] / S.K. Godunov, A.V. Zabrodin, M.Ya. Ivanov, A.N. Krayko, G.P. Prokopov. - M.: Science. - 1976. - 400 p.

3. Nikonov, V.V. Two-dimensional laminar boundary layer simulation by Vortex-in-cell method [Text] / V.V. Nikonov, V.G. Shakhov // Vestnik (Bulletin) SSAU. - Samara. - 2006, N 3

(11). - pp. 25 - 30.

4. Ferziger, J. Computational methods for fluid dynamics [Text] / J. Ferziger, M. Peric. -3 rev. ed. - Springer-Verlag. - 2002. - 423 p.

APPLYING LAGRANGE’S APPROACH TO SOLVING A ONE-DIMENSIONAL

PROBLEM OF WAVE PROPAGATION IN GAS IN LINEAR FORMULATION

© 2009 V. V. Nikonov, V. G. Shakhov Samara State Aerospace University

Modelling of one-dimensional wave propagation in gas using various numerical schemes is discussed. Alongside the known schemes a scheme of solving an acoustics problem is proposed using Lagrange’s approach to the description of environment motion. The results of numerical solution for three types of initial conditions are compared to those of analytical solution. Conclusions are drawn as to the applicability of the schemes discussed. The possibility of varying by step over the time without reducing the accuracy of the solution is shown.

Perfect gas, one-dimensional wave, numerical method, acoustics problem, “upstream” scheme, Godunov’s method, Lagrange’s approach, analytical solution.

Информация об авторах

Никонов Валерий Владимирович, кандидат технических наук, инженер научно-технологического парка «Авиатехнокон» Самарского государственного аэрокосмического университета. Область научных интересов: вихревые методы, прямое численное моделирование несжимаемых и сжимаемых течений, пограничный слой. Е-mail: v_nikonov@mail.ru.

Шахов Валентин Гаврилович, кандидат технических наук, профессор, заведующий кафедрой аэрогидродинамики Самарского государственного аэрокосмического университета. Область научных интересов: теория пограничного слоя, турбулентность, численные методы, аэродинамика летательных аппаратов. Е-mail: shakhov@ssau.ru.

Nikonov, Valery Vladimirovitch, engineer of the Science-and-Technology Park “Aviatechnokon”, Samara State Aerospace University, candidate of technical sciences. Area of research: vortex methods, direct numerical modelling of compressible and incompressible flows, boundary layer. Е-mail: v_nikonov@mail.ru.

Shakhov, Valentin Gavrilovitch, head of aerohydrodynamics department of Samara State Aerospace University, professor, candidate of technical sciences. Area of research: boundary layer theory, turbulence, numerical methods, aircraft aerodynamics. Е-mail: shakhov@ssau.ru.