Уравнение состояния Мурнагана представляет собой зависимость между объемом тела и давлением, которому оно подвергается. Это одно из многих уравнений состояния, которые использовались в науках о Земле и физике ударных волн для моделирования поведения материи в условиях высокого давления. Своим названием он обязан Фрэнсису Д. Мурнагану [1] , который предложил его в 1944 году, чтобы отразить поведение материала в максимально широком диапазоне давлений и отразить экспериментально установленный факт: чем сильнее сжимается твердое тело, тем труднее его сжимать. дальше.
Уравнение Мурнагана при определенных предположениях выводится из уравнений механики сплошной среды . Он включает в себя два регулируемых параметра: модуль несжимаемости K 0 и его первую производную по давлению K 0 , оба измеренные при давлении окружающей среды. В общем случае эти коэффициенты определяются путем регрессии экспериментально полученных значений объема V в зависимости от давления P. Эти экспериментальные данные могут быть получены методом рентгеновской дифракции или ударными испытаниями. Регрессию также можно выполнить на значениях энергии как функции объема, полученных в результате расчетов ab-initio и молекулярной динамики .
Уравнение состояния Мурнагана обычно выражается следующим образом: Если уменьшение объема при сжатии невелико, т.е. для V / V0 более примерно 90%, уравнение Мурнагана может моделировать экспериментальные данные с удовлетворительной точностью. Более того, в отличие от многих предложенных уравнений состояния, оно дает явное выражение объема как функции давления V ( P ). Но диапазон ее применимости ограничен, а физическая интерпретация неадекватна. Однако это уравнение состояния продолжает широко использоваться в моделях твердых взрывчатых веществ. Из более сложных уравнений состояния наиболее используемым в физике Земли является уравнение состояния Берча-Мурнагана . В ударной физике металлов и сплавов другим широко используемым уравнением состояния является уравнение состояния Ми – Грюнайзена .
Изучение внутреннего строения Земли посредством познания механических свойств составляющих внутренних слоев планеты предполагает экстремальные условия; давление можно исчислить сотнями гигапаскалей , а температуру — тысячами градусов. Изучение свойств материи в этих условиях можно проводить экспериментально с помощью таких устройств, как ячейка с алмазной наковальней для статического давления, или подвергая материал ударным волнам . Это также привело к теоретической работе по определению уравнения состояния, то есть отношений между различными параметрами, определяющими в данном случае состояние материи: объемом (или плотностью), температурой и давлением.
Есть два подхода:
Различные авторы предложили десятки уравнений. [2] Это эмпирические зависимости, качество и актуальность которых зависят от их использования и могут оцениваться по разным критериям: количеству задействованных независимых параметров, физическому смыслу, который можно приписать этим параметрам, качеству экспериментальные данные и непротиворечивость теоретических предположений, лежащих в основе их способности экстраполировать поведение твердых тел при высоком сжатии. [3]
Обычно при постоянной температуре модуль объемного сжатия определяется следующим образом: Самый простой способ получить уравнение состояния, связывающее P и V, — это предположить, что K постоянен, то есть не зависит от давления и деформации твердого тела, тогда мы просто найдите закон Гука. В этом случае объем уменьшается экспоненциально с давлением. Это неудовлетворительный результат, поскольку экспериментально установлено, что по мере сжатия твердого тела его становится труднее сжимать. Чтобы пойти дальше, необходимо принять во внимание изменение упругих свойств твердого тела при сжатии.
Предположение Мурнагана состоит в том, чтобы предположить, что модуль объемного сжатия является линейной функцией давления: [1] Уравнение Мурнагана является результатом интегрирования дифференциального уравнения: Мы также можем выразить объем в зависимости от давления:
Однако Пуарье критикует это упрощенное изложение как недостаточно строгое. [4] Эту же зависимость можно показать и иначе, исходя из того, что несжимаемость произведения модуля и коэффициента теплового расширения не зависит от давления для данного материала. [5] Это уравнение состояния также является общим случаем более старого соотношения Политропа [6], которое также имеет соотношение постоянной степени.
В некоторых обстоятельствах, особенно в связи с расчетами ab initio, предпочтительным будет выражение энергии как функции объема, [7] которое можно получить путем интегрирования приведенного выше уравнения в соответствии с соотношением P = − dE / dV . Его можно записать в K 0 , отличный от 3,
Несмотря на свою простоту, уравнение Мурнагана способно воспроизвести экспериментальные данные для диапазона давлений, который может быть весьма большим, порядка К 0 /2. [8] Это также остается удовлетворительным, поскольку соотношение V / V 0 остается выше примерно 90%. [9] В этом диапазоне уравнение Мурнагана имеет преимущество по сравнению с другими уравнениями состояния, если кто-то хочет выразить объем как функцию давления. [10]
Тем не менее, другие уравнения могут дать лучшие результаты, и несколько теоретических и экспериментальных исследований показывают, что уравнение Мурнагана неудовлетворительно для многих задач. Таким образом, в той степени, в которой отношение V / V 0 становится очень низким, теория предсказывает, что K 'достигнет 5/3, что является пределом Томаса-Ферми. [10] [11] Однако в уравнении Мурнагана K ′ является постоянным и устанавливается равным начальному значению. В частности, значение K 0 = 5/3 в некоторых ситуациях становится несовместимым с теорией. Фактически, при экстраполяции поведение, предсказанное уравнением Мурнагана, довольно быстро становится маловероятным. [10]
Независимо от этого теоретического рассуждения, опыт ясно показывает, что К ' уменьшается с давлением, или, другими словами, что вторая производная модуля несжимаемости К ″ строго отрицательна. Теория второго порядка, основанная на том же принципе (см. следующий раздел), может объяснить это наблюдение, но этот подход все еще неудовлетворителен. Действительно, это приводит к отрицательному объемному модулю в пределе, когда давление стремится к бесконечности. Фактически, это неизбежное противоречие, какое бы полиномиальное разложение ни было выбрано, потому что всегда будет доминирующий член, который расходится до бесконечности. [3]
Эти важные ограничения привели к отказу от уравнения Мурнагана, которое В. Хользапфель называет «полезной математической формой, не имеющей никакого физического обоснования». [12] На практике анализ данных сжатия осуществляется с использованием более сложных уравнений состояния. Наиболее часто используемым в научном сообществе является уравнение Берча-Мурнагана второго или третьего порядка по качеству собираемых данных. [13]
Наконец, очень общим ограничением этого типа уравнений состояния является их неспособность учитывать фазовые переходы, вызванные давлением и температурой плавления, а также множественные переходы твердое тело-твердое, которые могут вызвать резкие изменения плотности и модуля объемного сжатия. в зависимости от давления. [3]
На практике уравнение Мурнагана используется для выполнения регрессии набора данных, где получают значения коэффициентов K 0 и K 0 . Полученные коэффициенты и зная соотношение объема к условиям окружающей среды, позволяют нам в принципе рассчитать объем, плотность и модуль объемного сжатия при любом давлении.
Набор данных в основном представляет собой серию измерений объема для различных значений приложенного давления, полученных в основном с помощью дифракции рентгеновских лучей. Также можно работать с теоретическими данными, вычисляя энергию для разных значений объема методами ab initio, а затем регрессируя эти результаты. Это дает теоретическое значение модуля упругости, которое можно сравнить с экспериментальными результатами.
В следующей таблице перечислены некоторые результаты для различных материалов с единственной целью проиллюстрировать некоторые численные анализы, выполненные с использованием уравнения Мурнагана, без ущерба для качества полученных моделей. Учитывая критику физического смысла уравнения Мурнагана, высказанную в предыдущем разделе, к этим результатам следует относиться с осторожностью.
Чтобы улучшить модели или избежать критики, изложенной выше, было предложено несколько обобщений уравнения Мурнагана. Обычно они заключаются в отказе от упрощающего предположения и добавлении еще одного регулируемого параметра. Это может улучшить качество утонченности, но также привести к усложнению выражений. Также поднимается вопрос о физическом смысле этих дополнительных параметров.
Возможная стратегия состоит в том, чтобы включить дополнительный термин P 2 в предыдущую разработку, [17] [18] требуя, чтобы . Решение этого дифференциального уравнения дает уравнение Мурнагана второго порядка: где . Находится естественным путем в уравнении первого порядка, приняв . Развитие до порядка выше 2 в принципе возможно [19] , но ценой добавления регулируемого параметра для каждого члена.
Можно привести и другие обобщения: