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

УДК 536.71

МОДЕЛИРОВАНИЕ УДАРНОГО СЖАТИЯ И ТЕПЛОВОГО РАСШИРЕНИЯ ПЯТИ МЕТАЛЛОВ1

В.Ф. Куропатенко2, С.Ю. Филатов3

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

Ключевые слова: уравнение состояния вещества; симплекс-метод; тепловое расширение; ударное сжатие.

Введение

В настоящее время в литературе описано большое количество уравнений состояния веществ (УРС) от очень простых до очень сложных [1-3]. При изучении поведения сплошных сред, которые подвергаются воздействию динамических нагрузок, применяется математическое моделирование, выполняющееся на электронных вычислительных машинах (ЭВМ). Число компонентов в сложных моделях многокомпонентных сред с химическими реакциями и фазовыми переходами может достигать нескольких десятков, а то и сотен. При расчете смеси в таких моделях затраты машинного времени могут растягиваться на значительное время даже с использованием современных ЭВМ. Предлагаемый малопараметрический УРС значительно сократит время расчета на ЭВМ сложных физических процессов.

Уравнения на поверхности сильного разрыва

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

р(В - и ) = Ро (В - и0), (1)

р.( в - и )и - Р = Ро •( В - ио )ио - Ро, (2)

Ро = 0. (3)

В уравнениях (1)-(3) величины без индекса характеризуют состояние за разрывом, В - скорость ударной волны. Величины с индексом, такие как ро - плотность, ио - массовая скорость, Ро - давление, Ео - удельная внутренняя энергия, описывают состояние вещества перед разрывом. Рассмотрим ударную волну в покоящемся веществе. Согласно [1] уравнение состояния рассматривалось в предположении, что величины Ро и Ео пренебрежимо малы по сравнению с Р и Е . Уравнения (1)-(3) принимают вид при Ро = о, Ео = о и ио = о :

р.(В - и )-Ро В = о, РоВи = Р, (4)

E =1DU-2

. Р )

(5)

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

2 Куропатенко Валентин Федорович - доктор физико-математических наук, профессор, главный научный сотрудник, Российский Федеральный Ядерный Центр - Всероссийский Научно-Исследовательский Институт Технической Физики им. академика Е.И. Заба-бахина, г. Снежинск.

E-mail: v.f.kuropatenko@rambler.ru

3 Филатов Сергей Юрьевич - научный сотрудник, Российский Федеральный Ядерный Центр - Всероссийский Научно-Исследовательский Институт Технической Физики им. академика Е.И. Забабахина, г. Снежинск.

E-mail: phil@chel.ru

Куропатенко В.Ф., Моделирование ударного сжатия

Филатов С.Ю. и теплового расширения пяти металлов

Система, состоящая из трех уравнений, (4), (5) содержит пять величин Р, р, Е , и, В. Если любые две из них получены экспериментальным путем, то остальные величины можно найти из выражений (4) и (5). Точка на ударной адиабате определяется полностью. Зависимость между экспериментально измеренными В и и известна уже более пятидесяти лет и описывается линейным соотношением

В = С0 + Ь ■и. (6)

На основе большого количества обработанных экспериментальных данных в справочнике [2] содержится информация по В (и) соотношениям.

Давление Р определяется калорическим уравнением состояния (УРС) вида Р = Р(р,Е) в том случае, если р и Е - независимые термодинамические переменные. Сложные и трудоемкие современные УРС рассмотрены в [3], [4], однако, для экспресс-расчетов достаточно использовать простые УРС.

В [1] рассмотрено уравнение состояния вида

Р = Рх (р) + РТ (р,5), Е = Ех (р) + Ет (р,5), (7)

где

Рг =

р°с°2

(х- -1), Рт = р°С°2 • f (5)• х^, (8)

"X 1 n -1

n

n i ,„ 1 \ Л с

= f isV-y. (9)

Y-1

, = _cl (хш+-1), , = am xi-

nn

\

C° - скорость звука перед ударной волной при Р = Р°, р = р° , S - энтропия, n = const, Y= const, х = р°/р .

Зависимости Рх (EX, х) и Рт (ET, х) следуют из (8) и (9)

Рх = (i-1)Mx. + р°с°2 • Гi- 1, Рт J-l-nML, (1°)

х V х J х

УРС получим путем подстановки (1°) в (7), которое принимает следующий вид

Р = (y- 1)р°Ех"1 +ф(х), (11)

n-Y Y-1 Y

ф( х )=" ' х~П +-

(п -1) п (п -1) X п Согласно [1] уравнение (11) называется «уравнением с согласованными у и п », если у= п. Таким образом, уравнение (11) примет вид

Р = (п - 1)роЕ ■ х"1 +роСо2 (х"1 -1), (12)

где р0, С0 и п - постоянные величины.

Далее рассмотрим, при каких условиях уравнения (4), (5), (12) согласуются с линейной зависимостью В (и). Зависимость В (и) получим, путем исключения Р, Е и х в выражениях (4),

(5) и (12).

в = n±i U + А

4 V

C°2 +| ^U | . (13)

n +1 -

2

Сравнивая выражение (13) с выражением (6), можно увидеть, что эти зависимости различа-

п + 1

ются. Наименьшие отличия находятся в области 4 и С0 .

Так как предположение о постоянстве п в УРС (12) приводит к нелинейной зависимости В (и) в выражении (13), которая, в свою очередь, не согласуется с линейной зависимостью

В (и) (6), сделаем наоборот. Линейную зависимость В (и) (6) примем за основу, Р и Е выра-

зим через С0 и Ь из выражений (4) - (6). Полученные выражения подставим в (11) и найдем уравнение в котором п вдоль ударной адиабаты не может быть постоянной величиной

п = 2Ь-(2 - Ь-(1 - х))-1, (14)

Уравнение состояния

Поведение вещества вместо УРС (12) будем описывать уравнением, в котором справедливо предположение, что п не зависит от х

Р = (п(х)- 1)-А>кх~1Е + РокС0к ф(X). (15)

В уравнении (15) p0K, C0K - плотность и скорость звука в точке P = 0, T = 0, х = 1.

В [5-7] использовалось УРС типа (15), где величина п зависит от х. В уравнении (14) величина п линейно зависит от х в том диапазоне D и U, в котором справедлива зависимость D(и) (6). Из выражений (4), (5) и (15) следует зависимость п (х) (14), именуемая далее экспериментальной.

Функция п(х) должна иметь максимум в окрестности х ~ 1 согласно [2] и [3]. Таким образом, если х = 0 и п = п0 , то справедливо утверждение, что nm > п0 . Будем использовать простую функцию вида

п (х ) = п0 + (^ - п0 )-;-ту, (16)

2

2.(2 2 Г ах +(х -^)

16 х2т

где а-

\2 '

4-(^ - п0 )

Зная, что при х = 0 и х = ^ п (х) = п0 , а также при х = хm п (х) = nm , из (16) найдем производную п (х) в виде

dn ( ^ - п0 )•2ах2 (^ - х2)

I ^ - п0 )• 2ах (^ - х

Лх ( 2 / 2 2 \2 I т +(х -^) I

(17)

Производная п (х) равна нулю при х = 0, х = ^ и х = хm. То есть, с большой долей достоверности вещество можно считать газом с у = п0 при бесконечно большом сжатии х ~ 0 и бесконечно большом разрежении. Значение п0 , как и значения nm, хm, подберем из условия наилучшего описания экспериментальных данных, т.к. область применимости малопараметрического УРС ограничена конечным сжатием.

Функцию ф(х) из (15) возьмем в виде, максимально близком к выражению (12)

1 - х

Ф( х ) = хп0+Т . (18)

Разделим давление и энергию на холодные и тепловые составляющие для определения температуры и теплоемкости:

P = PX (х) + PT (х^), E = EX (х) + ET (х,T). (19)

В соответствии с [7] зависимость ET (х^) возьмем в виде

AT 2

^ =Щ+¥ ■ <20)

где A - индивидуальная характеристика вещества, для простых веществ близкая к 3Л/Ц; Я -универсальная газовая постоянная; ¡1 - молекулярная масса. Дифференцируя выражение (20), получим теплоемкость при постоянном х

Куропатенко В.Ф., Филатов С.Ю.

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

Cv =

AT •(2в(х) + T)

И х) + т Г

В соответствии с [7] выражение для теплового давления будет иметь вид

Рт =

A dQ(х) T2

в(х) dх в(х) + T

Зная, что справедливо уравнение

Р°к

Рт =( n ( х )--et ,

Ет

получим характеристическую функцию в( х), которая связана с n (х) уравнением

d \пв(х) n (х)-1

dх

Зависимость в(х) следует из выражений (24) и (16)

в(х) = 0° • х1-п° •

х2 + B • х1

2 . 1 2

х + — • хт

B

(21)

(22)

(23)

(24)

(25)

где B =

2 + ( nm - n° ) 2-( nm - n° ) '

Дополнительно, сделаем такое условие, чтобы изобара Р = 1° 4 ГПа проходила через точку, характеризующую нормальное состояние ( Р = 1°-4 ГПа, T = 293 °К, р = р° , СР = СР° ), а также через точку плавления при Р = 1°-4 ГПа ( T = Тпл , СР = СРпл , р = рпл). При описании зависимости СР (T) при Р = const следует выразить СР через T и х. С этой целью воспользуемся уравнением из [8]

T •

СР = Су —-

дР

дт

х

р°к •

дР дх

(26)

Зная, что

дР

х

V дт Jх

= ° продифференцируем выражение (23). Производная

дР

дт

запишется

х

в виде

дРт

~дТг

(n (х)- 1)^р°к

• С

(27)

х

Так как выражение (19) справедливо, то производная I — I будет иметь вид

дР

дР

дх

дх

dFx +1 дР

т dх V дх j т

(28)

Уравнение состояния вещества определяется совокупностью уравнений (15)-(17) и (25) и содержит 7 параметров: р0К, С0К, А, п0, пт, хт и в0 . Численные значения этих параметров определяются таким образом, чтобы наилучшим образом описать поведение ударной адиабаты состояния вещества в точке, характеризующей нормальное состояние Р0 = 10-4 ГПа, Т0 = 293 °К,

р0 =р(Р0,Т0 ) , СР0 = СР (р0,т0) и в точке плавления Рт = 10-4 ГП^ Т1 = Тт , р =рт (Тпл, Р0 ) , СР1 = СР (Тпл > Р0 ) .

х

2

х

Результаты расчетов

Расчеты выполнены для нескольких простых веществ - металлов. При помощи оригинальной версии симплекс-метода выполнен подбор основных параметров уравнения состояния. В качестве целевой функции выбрана сумма квадратичных разностей между расчетными и экспериментальными значениями теплового расширения вещества и теплоемкости при постоянном давлении. Для рассмотренных материалов получено удовлетворительное совпадение с экспериментальными данными. На рис. 1 показано сравнение безразмерного давления для ударных адиабат рассчитываемых материалов. Величина безразмерного давления определяется по формуле (29):

P

п=Р0СТ <29)

В качестве сравниваемых величин давлений используются экспериментальные, взятые из [2], и теоретические, определяемые из выражения (15).

Зависимость теплоемкости от температуры приведена на рис. 2, экспериментальные значения теплоемкости взяты из справочника [9]. В таблицу сведены результаты расчетов подбора основных параметров УРС.

П(х)

X

Рис. 1. Сравнение экспериментальных и теоретических зависимостей П(х) на ударной адиабате

Результаты подбора основных параметров уравнения состояния

Величины Вещества

А1 (алюминий) Аи (золото) (магний) РЬ (свинец) Си (медь)

«0 1,62 1,515 1,291 1,719 1,49

пт 1,933 2,564 1,944 2,027 2,997

Хт 1,908 2,485 2,613 1,66 2,484

Рок >г/см3 2,908 19,668 1,833 11,501 9,014

Сок, км/с 5,642 3,402 4,724 2,102 3,904

АЛО6, кДж/г 920,39 140,65 1106,14 127,6 394,94

0О,К 15,012 8,004 43,141 39,198 1,466

Куропатенко В.Ф., Филатов С.Ю.

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

Ср,джш*1д uoü

12D0

Í000

300

600

400

200

0

О 200 «О ООО BD0 1ОО0 1400 1600

Температура. К

Рис. 2. Сравнение экспериментальных и теоретических зависимостей Ср(Т) Литература

1. Забабахин, Е.И. Некоторые вопросы газодинамики взрыва / Е.И. Забабахин. - Снежинск: РФЯЦ-ВНИИФ, 1977. - 203 с.

2. Экспериментальные данные по ударно-волновому сжатию и адиабатическому расширению веществ / Р.Ф. Трунин, Л.Ф. Гударенко, М.В. Жерноклетов, Г.В. Симаков. - Саров: РФЯЦ-ВНИИЭФ, 2001.-439 с.

3. Жарков, В.И. Уравнение состояния твердых тел при высоких давлениях и температурах / В.И. Жарков, В .А. Калинин. - М.: Наука, 1968. - 310 с.

4. Куропатенко, В.Ф. Уравнение состояния компонентов низкотемпературной плазмы / В.Ф. Куропатенко // Энциклопедия низкотемпературной плазмы: сб. науч. тр. - М.: Янус-К, 2008. - Серия Б, Том VII. - Ч. 2. - С, 436-450.

5. Куропатенко, В.Ф. Уравнение состояния продуктов детонации конденсированных ВВ / В.Ф. Куропатенко // Численные методы механики сплошных сред. - Новосибирск, ВЦ СО АН СССР, ИТПМ СО АН СССР, 1977. - Т. 8, №6. - С, 68-71.

6. Куропатенко, В.Ф. Уравнение состояния продуктов детонации плотных ВВ / В.Ф. Куропатенко // Физика горения и взрыва. - 1989. - № 6. - С. 112-117.

7. Куропатенко, В.Ф. Моделирование отклика веществ на динамическое воздействие / В.Ф. Куропатенко // Химическая физика. - 2002. - Т. 21, № 10. - С. 46-54.

8. Куропатенко, В.Ф. Модели механики сплошных сред / В.Ф. Куропатенко. - Челябинск: Челябинский государственный университет. - 2007. - 302 с.

9. Физико-химические свойства элементов. Справочник / под ред. Самсонова Г.В. - Киев: Наукова думка, 1965. - 808 с.

Поступила в редакцию 26 июня 2014 г.

• At (расчет) 1 ■ Al ¡экспеэнмент: 2 - Au (расчет) * Al (эксперимент) ^^^ Си [расчет! ■ Си (эксперимент) - Me) (расчет) 4 * Mg (эксперимент) 5 — Pb (расчет) « Pb (экперммент) -

4 4 4

/ i Щ 1

■ m • ш 3

• • 2 -

и! ■ „"í---- 5

Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2014, vol. 6, no. 3, pp. 30-36

MODELLING OF SHOCK COMPRESSION AND THERMAL EXPANSION OF FIVE METALS

V.F. Kuropatenko^, S.Yu. Filatov2

The method of building equations of state in the area of low pressures is offered. The method is based on linear dependence of shock wave velocity on a substance velocity and thermal expansion data at constant pressure. Linear dependence of shock wave velocity on mass velocity leads to the dependence of Gruneisen parameters on the density. This data is used together with thermal expansion data at constant pressure. The parameters for the equation of state are computed for five metals with the help of the original version of Simplex method. A sum of quadratic differences between calculated and experimental values is selected as a target function. The theoretical values have agreed with the experimental one.

Keywords: equation of state; Simplex method; thermal expansion; shock compression.

References

1. Zababakhin E.I. Nekotorye voprosy gazodinamiki vzryva (Some issues of gas dynamics of explosion). Snezhinsk, RFYaTs-VNIIF, 1977. 203 p. (in Russ.).

2. Trunin R.F., Gudarenko L.F., Zhernokletov M.V., Simakov G.V. Eksperimental'nye dannye po udarno-volnovomu szhatiyu i adiabaticheskomu rasshireniyu veshchestv (Experimental data about shock wave compression and adiabatic expansion of substances). Sarov: RFYaTs-VNIIEF, 2001. 439 p. (in Russ.).

3. Zharkov V.I., Kalinin V.A. Uravnenie sostoyaniya tverdykh tel pri vysokikh davleniyakh i tem-peraturakh (Equation of state for solids at high pressure and temperatures). Moscow, Nauka, 1968. 310 p. (in Russ.).

4. Kuropatenko V.F. Uravnenie sostoyaniya komponentov nizkotemperaturnoy plazmy (Equation of state for components of low-temperature plasma). Entsiklopediya nizkotemperaturnoy plazmy. Chast 2. (Encyclopaedia of low-temperature plasma. Part 2). Moscow, Yanus-K Publ., 2008. Series of B. Vol. VII. pp. 436-450. (in Russ.).

5. Kuropatenko V.F. Chislennye metody mekhaniki sploshnykh sred. Novosibirsk, VTs SO AN SSSR Publ., ITPM SO AN SSSR Publ., 1977. Vol. 8, no. 6. pp. 68-71. (in Russ.).

6. Kuropatenko V.F. Fizika goreniya i vzryva. 1989. no. 6. pp. 112-117. (in Russ.).

7. Kuropatenko V.F. Khimicheskayafizika. 2002. Vol. 21, no. 10. pp. 46-54. (in Russ.).

8. Kuropatenko V.F. Modeli mekhaniki sploshnykh sred (Models of continuum mechanics). Chelyabinsk, Chelyabinskiy gosudarstvennyy universitet Publ., 2007. 302 p. (in Russ.).

9. Samsonov G.V. (ed.) Fiziko-khimicheskie svoystva elementov. Spravochnik (Physical and chemical properties of elements. Handbook). Kiev, Naukova dumka, 1965. 808 p. (in Russ.).

Received 26 June 2014

1 Kuropatenko Valentin Fedorovich is Dr. Sc. (Physics and Mathematics), Professor, Senior Researcher, Russian Federal Nuclear Center -Zababakhin All-Russian Scientific Research Institute of Technical Physics, Snezhinsk.

E-mail: vfkuropatenko@rambler.ru

2 Filatov Sergey Yurjevich is Resercher, Russian Federal Nuclear Center - Zababakhin All-Russian Scientific Research Institute of Technical Physics, Snezhinsk.

E-mail: phil@chel.ru