Метод крупных частиц для исследования течений газовзвесей

УДК 519.63 + 532.529.5

DOI: 10.14529/^сг150210

МЕТОД КРУПНЫХ ЧАСТИЦ

ДЛЯ ИССЛЕДОВАНИЯ ТЕЧЕНИЙ ГАЗОВЗВЕСЕЙ

Ю.М. Ковалёв1, Е.А. Ковалёва2

1 Южно-Уральский государственный университет, г. Челябинск;

2 Челябинский государственный университет, г. Челябинск

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

Ключевые слова: численный метод, математическая модель, газовзвесь, законы сохранения, ударные волны, число Куранта.

Введение

Несмотря на наличие большого числа вычислительных пакетов и увеличение быстродействия вычислительной техники, разработка эффективных численных методов для решения задач в рамках новых математических моделей механики сплошных сред в настоящее время является актуальной задачей. Появление новых математических моделей, с одной стороны, связано с отсутствием в природе чистых веществ, что требует активного развития математических моделей многокомпонентных сред, достоверно описывающих физические процессы, применяемые в различных отраслях науки и техники. С другой стороны, развитие вычислительной техники позволяет получать численные решения для новых [1] все более сложных математических моделей многокомпонентных сред. Более того, есть такие проблемы, когда математическое моделирование является единственным средством предварительного изучения явлений (например, [2]). Адекватность математических моделей многокомпонентных сред физическим процессам предъявляет достаточно жесткие требования к математическим моделям: с одной стороны, уравнения сохранения должны быть инвариантны относительно преобразования Галилея [3], с другой стороны, должны выполняться законы сохранения для смеси [4]. В работах [1, 5] было показано каким образом можно выполнить оба эти условия.

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

1. Математическая модель газовзвеси

Рассмотрим одномерный плоский случай математической модели течения газа с твердыми частицами (аэровзвесь), которая описывается системой уравнений сохранения [5]. Данная система уравнений двухфазной аэровзвеси [5] без химических превращений имеет следующий вид:

£Р1,£Р1^1 = П др2 + др2у2 = 0 дп + дпу2 = _

дх ' дх ' дх 0

Р = Р1 (Р1, = Р2{р°2, Т2), = (р[, Т1), е2 = е2(р2, Т2) ,

Р\ = р!«1, Рг = Р°2а2 , а1 + а2 = 1, = ег + = 1, 2) ,

/ = па2р1 Са

Система уравнений

Здесь индексы 1, 2 относятся соответственно к газу и частицам; р°, а^ (£ = 1, 2) - истинные плотности и объемные содержания фаз; р¿, Еь - парциальная плотность, скорость, темпера-

тура, внутренняя и полная энергия 7-й фазы; р - давление, п - число частиц в единице объема смеси; су1и с2 - теплоемкости фаз: С0 - постоянная для нормирования внутренней энергии газовой фазы: - теплопроводность газовой фазы; R1 - универсальная газовая постоянная; Саи Ыи - коэффициент трения и число Нуссельта, определяемые числами Рейнольдса (Ке) и Прандтля (Рг) относительного движения фаз соответственно: к - показатель адиабаты Пуассона; ^ - диаметр частиц.

Уравнения (1.1) - уравнения неразрывности газа и частиц и уравнение сохранения числа частиц в единице объема смеси; (1.2) - уравнения импульса газа и частиц; (1.3) и (1.4) - уравнения сохранения внутренней энергии газа и частиц соответственно. (1.6) - уравнения, определяющие члены теплового (ф) и силового (/) взаимодействия между фазами: (1.7) - уравнения состояния фаз. В данной работе не рассматриваются более сложные уравнения состояния [8].

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

Умножая уравнение сохранения импульса газовой фазы на р1, а уравнение сохранения импульса конденсированной фазы на р2, получим уравнения сохранения кинетической энергии газа и частиц соответственно:

ХдолУл . до-,г2

V

которые после простых преобразований принимают следующий вид:

дрл^ дР1у^ др

+ -аТ^ = -

др2Ц~ др2у2^ д-р

+ = Тх + (19) Преобразуем левые части уравнений сохранения внутренней энергии газа (1.3) и частиц (1.4) к дивергентному виду. С учетом равенств (1.1) они могут быть представлены в виде:

др1е1 др1е1у1 ра± йгр1 „¿у... „ ^ „„ п 1т

1Т + = — + - ^ -

др2е2 др2е2у2 ра2 й2р°2 , .

—+^т = т^)—+щ.

Из уравнений неразрывности газовой и конденсированной фаз (1.1) легко получить следующие равенства:

о да1 , да^Ут. „ ¿2Р° _ ° (да2 , да2у2\

а2~= р2\дГ + ~дГГ

Подставляя данные выражения в уравнения (1.10) и (1.11) соответственно, получим:

др1е1 др1е1у1 гда1 да1у1^ „¿у... ,, ^ „„ п

1Т" + = + ~) + - ^ - пЧ>

+ = - +

дЬ дх ^ \ дЬ дх ) ^ у '

Ковалёв Ю.М., Ковалёва Е.А.

Метод крупных частиц для исследования течений газовзвесей

В случае несжимаемости конденсированной фазы уравнения сохранения внутренней энергии газовой (1.3) и конденсированной (1.4) фаз, легко преобразуются к виду:

^ + +Ц?-) + nf^v1-v2)-щ,

+ =

Для получения уравнение сохранения полной энергии смеси просуммируем левые и правые части уравнений (1.8), (1.9), (1.14), (1.15). В результате получим уравнение сохранения полной энергии смеси в виде

д

Система уравнений (1.1), (1.2), (1.5)-(1.7), (1.14)—(1.16) представляет собой замкнутую систему уравнений для описания течений газовзвесей, инвариантную относительно преобразования Галилея.

= 0,

дуЛ = др

де-, Л1Т = Р( дХ

де7 = nq,

2. Модификация метода крупных частиц для расчета течений газовзвеси

В соответствии с идеологией метода крупных частиц [6] систему законов сохранения газовзвеси (1.1), (1.2), (1.5)—(1.7), (1.14)—(1.16) на эйлеровом этапе можно представить следующим образом:

%=0.

+ a-^-) + nf

д

Учитывая несжимаемость конденсированной фазы (р2 = const), запишем уравнения (2.1), (2.3), (2.5) в более удобном для представления на эйлеровом этапе виде:

dPi п ° да2 п дп _ о даг _ «1-^=0, р2 — =0, -=0, Pl — = °,

Pi ^ = -р

дЕ1 + п ££2 , п d

Подставляя уравнение состояния газовой фазы (1.7) в уравнение (2.7) получим следующее базовое соотношение для определения давления на эйлеровом этапе

+ - *) - "")

Используя явные разностные представления для равенства (2.9), легко получить выражения для определения предварительных значений давления на новом т + 1 временном слое на границах i - 1/2 и i + 1/2 для ячеек i - 1, i и i + 1:

ххт+1

.

Pi+1/2 — ;

+ «2™+1/2(^ + 1 - - ^(<1/2^1/2^ +

Дл 1,1+1/2

+ 2/^1/2

"1,1+1/2

Здесь Д£ - шаг по времени, Дх - шаг по пространству. Полученные значения давления используются для определения промежуточных величин скоростей на эйлеровом этапе:

ЖУ^-пГТГД*,

рГ^ + ПГ/ТД^

vl,i — vm - Pi + 1/2

v2,i — vm -V2,i Pi + 1/2

Промежуточные значения скоростей конденсированной и газовой фаз на границах ячеек определяются как средние арифметические от их значений в двух соседних ячейках:

«ffi/2 =

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

¿и1 = + -^пГчГЫ

и полной энергии смеси

_m ?т+1 | _т Рт+1 _ _т т?т , _т г-т

Pl,iil,i + Pl,i 2,i = Pl,ihl,i + P2,t^2,£ -

^„n fim+1 ~m+l n ftm+l ^m+l Л At

lal,i + l/2yl,i + l/2Pi + l/2 al,i-l/2Vl,i-l/2Pi-l/2j ^

—

На этапе Лагранжа и заключительном этапе метода крупных частиц для каждой фазы были использованы формулы, приведенные в монографии О.М. Белоцерковского и Ю.М. Давыдова [6].

Заключение

1. Тестирование предложенной модификации метода крупных частиц проводилось на решении задач о распространении ударных волн в «замороженной» газовзвеси [9] и в облаке газовзвеси [10].

2. Было показано, что применение на этапе Эйлера уравнений (2.10)-(2.15) более эффективно, чем применение метода крупных частиц [6] и модифицированного метода крупных частиц [11] прирешении задач о распространении ударных волн в «замороженной» газовзвеси [9] и в облаке газовзвеси [10] даже при отсутствии искусственной вязкости и большем числе Куранта.

Авторы выражают свою благодарность профессору В. Ф. Куропатенко за полезные обсуждения

и интерес к работе.

Работа выполнена при поддержке РФФИ грант №13 - 01 - 00072.

Литература

1. Куропатенко, В.Ф. Новые модели механики сплошных сред / В.Ф. Куропатенко // Инженерно-физический журнал. - 2011. - Т. 84, № 1. - С. 74-92.

2. Гришин, А.М. Об усилении ударных волн при их взаимодействии с фронтом лесного пожара /А.М. Гришин, Ю.М. Ковалев //Доклады Академии наук. - 1990. - Т. 312, № 1. - С. 50-54.

3. Ковалев, Ю.М. Математическая модель газовзвеси с химическими превращениями в приближении парных взаимодействий / Ю.М. Ковалев, Е.Е. Пигасов // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». - 2014. - Т. 7, № 3. - С. 40-49.

4. Ковалев, Ю.М. Математический анализ уравнений сохранения двухфазных смесей /Ю.М. Ковалев, Е.А. Ковалева // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». - 2014. - Т. 7, № 2. - С. 29-37.

5. Ковалев, Ю.М. Анализ возможности применения некоторых численных методов для решения задач механики многокомпонентных сред / Ю.М. Ковалев, Е.А. Ковалева // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2014. - Т. 14, № 1. - С. 57-62.

6. Белоцерковский, О.М. Метод крупных частиц в газовой динамике / О.М. Белоцерковский, Ю.М. Давыдов. - М.: Наука, 1982. - 392 с.

7. Гришин, Ю.А. Новые схемы метода крупных частиц и использование их для оптимизации газовоздушных трактов двигателей /Ю.А. Гришин //Математическое моделирование. - 2002. -Т. 14, № 8. - С. 51-55.

8. Ковалев, Ю.М. Уравнения состояния и температуры ударного сжатия кристаллических ВВ / Ю.М. Ковалев // Физика горения и взрыва. - 1984.- Т. 20, № 2. - С. 102-107.

9. Кругликов, Б.С. Ослабление воздушных ударных волн экранирующими решётками / Б.С. Кругликов, А.Г. Кутушев // ФГВ. - 1988. - № 1. - С. 115-117.

Ковалёв Ю.М., Ковалёва Е.А.

Метод крупных частиц для исследования течений газовзвесей

10. Кругликов, Б.С. Ослабление воздушных ударных волн слоями запыленного газа и решетками /Б.С. Кругликов, А.Г. Кутушев //ПМТФ. - 1988. - № 1. - С. 51-57.

11. Ивандаев, А.И. Численное исследование нестационарных волновых течений газовзвесей с выделением границ двухфазных областей и контактных разрывов в несущем газе / А.И. Ивандаев, А.Г. Кутушев // Численные методы в механике сплошных сред. - 1983. - Т. 14, № 6. - С. 47-60.

Ковалев Юрий Михайлович, д-р физ.-мат. наук, профессор, заведующий кафедрой вычислительной механики сплошных сред, Южно-Уральский государственный университет, г. Челябинск; yum_kov@mail.ru.

Ковалева Елена Адамовна, канд. физ.-мат. наук, доцент кафедры математических методов в экономике, Челябинский государственный университет, г. Челябинск; ea_kov@mail.ru.

Поступила в редакцию 18 марта 2015 г.

DOI: 10.14529/ctcr150210

METHOD OF LARGE PARTICLES FOR RESEARCH OF CURRENTS OF GAS-SUSPENSIONS

Yu.M. Kovalev, South Ural State University, Chelyabinsk, Russian Federation, yum_kov@mail.ru, E.A. Kovaleva, Chelyabinsk State University, Chelyabinsk, Russian Federation, ea_kov@mail.ru

Modification of a method of large particles in the annex to researches of currents of gas-suspensions was developed. It is shown that this modification of method of large particles allows to carry out calculations of behavior of shock waves in gas-suspensions without introduction in an explicit form of artificial viscosity. It allowed to eliminate distortions of a physical picture of a current of the gas-suspension connected with existence of the ostsillyation taking place at distribution of shock waves in non-uniform environments. It was shown that the obvious option of modification of a method of large particles when using of the offered modification for carrying out calculations of distribution of shock waves in gas-suspensions with large numbers of Courant can be used. It allows to reduce considerably time of calculation of a task and to avoid carrying out the difficult iterative procedures inherent in implicit differential schemes. It is shown that the offered modification is effective and allows to carry out calculations of shock waves in gas-suspensions with large numbers of Courant.

Keywords: numerical method, mathematical model, gas-suspensions, conservation laws, shock waves, number of Courant.

References

1. Kuropatenko V.F. [New Models of Mechanics of Continuous Environments]. Engineering Physical Journal, 2011, vol. 84, no. 1, pp. 74-92.

2. Grishin A.M., Kovalev Yu.M. [About Strengthening of Shock Waves at their Interaction with the Front of Forest Fire]. Reports of Academy of Sciences, 1990, vol. 312, no. 1, pp. 50-54.

3. Kovalev Yu.M., Pigasov E.E. [Mathematical Model of a Gas-Suspension with Chemical Transformations in Approach of Pair Interactions]. Bulletin of the South Ural State University. Ser. Mathe-matic Modelling and Programming, 2014, vol. 7, no. 3, pp. 40-49.

4. Kovalev Yu.M., Kovaleva E.A. [Mathematical Analysis of the Equations of Preservation Two-Phase Mixes]. Bulletin of the South Ural State University. Ser. Mathematic Modelling and Programming, 2014, vol. 7, no. 2, pp. 29-37.

5. Kovalev Yu.M., Kovaleva E.A. [Analysis of Possibility of Application of Some Numerical Methods for the Solution of Problems of Mechanics of Multicomponent Environments]. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2014, vol. 14, no. 1, pp. 57-62.

6. Belotserkovskiy O.M., Davydov Yu.M. Metod krupnykh chastits v gazovoy dinamike [Method of Large Particles in the Gas Dynamics]. Moscow, Nauka Publ., 1982. 392 p.

7. Grishin Yu.A. [New Schemes of a Method of Large Particles and Use them for Optmization of Air-Gas Paths of Engines]. MathematicModelling, 2002, vol. 14, no. 8, pp. 51-55.

8. Kovalev Ju.M. [Equations of a State and Temperature of Shock Compression Crystal Explosives]. Physics of Burning and Explosion, 1984, vol. 20, no. 2, pp. 102-107.

9. Kruglikov B.S., Kutushev A.G. [Weakening of Air Shock Waves by the Shielding Lattices]. Physics of Burning and Explosion, 1988, no. 1, pp. 115-117.

10. Kruglikov B.S., Kutushev A.G. [Weakening of Air Shock Waves by Layers of Dusty Gas and Lattices]. Applied Mechanics and Technical Physics. 1988, no. 1, pp 51-57.

11. Ivandaev A.I., Kutushev A.G. [Numerical Research of Non-Stationary Wave Currents of Gas-Suspensions with Allocation of Borders of Two-phase Areas and Contact Gaps in the Bearing Gas]. Numerical Methods in Mechanics of Continuous Environments, 1983, vol. 14, no. 6, pp. 47-60.

Received 18 March 2015

БИБЛИОГРАФИЧЕСКОЕ ОПИСАНИЕ СТАТЬИ

Ковалёв, Ю.М. Метод крупных частиц для исследования течений газовзвесей / Ю.М. Ковалёв, Е.А. Ковалёва // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». -2015. - Т. 15, № 2. - С. 91-96. Ш1: 10.14529/Ссг150210

REFERENCE TO ARTICLE

Kovalev Yu.M., Kovaleva E.A. Method of Large Particles for Research of Currents of Gas-Suspensions. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2015, vol. 15, no. 2, pp. 91-96. (in Russ.) DOI: 10.14529/ctcr150210