Уравнение состояния Ми-Грюнайзена

страница в разработке

Уравнение состояния Ми-Грюнайзена

При больших давлениях и температурах принято представлять давление [math]p[/math] в конденсированном веществе в виде суммы "холодной" и "тепловой" компонент:

[math]p = p_0 + p_T, p_T = p - p_0[/math]

Холодная компонента, часто называемая "холодной кривой" (cold curve), обусловлена деформированием кристаллической решетки, а вторая - тепловыми колебаниями атомов. Иными словами, холодное давление зависит только от объема, а тепловое - от объема и тепловой энергии [math] E_T [/math]:

[math]p = p_0(V) + p_T(V,E_T)[/math]

Тепловая энергия - часть внутренней энергии твердого тела, обусловленная тепловым движением атомов. В первом приближении тепловая энергия равна [math] c_V T [/math]. На практике часто предполагается линейная связь теплового давления и тепловой энергии:

[math] p = p_0(V) + \frac{\varGamma(V)}{V} E_T[/math]

Данное уравнение называют уравнением состояния Ми-Грюнайзена, а функцию [math]\varGamma(V)[/math] - коэффициентом Грюнайзена.

Уравнение состояния для кристаллов простой структуры

[math] p_0 = \frac{1}{2V_0d\theta^d}\sum_{k=1}^n N_k\varPhi_k A_k^2,~~~~\varGamma = -\frac{\sum_{k=1}^n N_k((d+2)\varPhi'_k A_k^2 + 2\varPhi''_k A_k^4 )}{d\sum_{k=1}^n N_k (d\varPhi_k +2\varPhi'_k A_k^2)} [/math]

где [math]k[/math] - номер координационной сферы, [math]n[/math] - их число, [math]N_k[/math] - число атомов на [math]k[/math]-ой координационной сфере, [math] A_k = \rho_k R \theta[/math] - радиус координационной сферы, [math] \rho_k=A_k/A_1 [/math] - безразмерные константы решетки, [math]R[/math] - радиус первой координационной сферы в отсчетном положении, [math]\varPhi^{(n)}_k = \varPhi^{(n)}(A_k^2)[/math].

Холодная кривая для потенциалов Леннарда-Джонса, Ми, Морзе

Потенциал Леннарда-Джонса [math] \PI(r) =D[(\frac{a}{r}\right)^{12}-2(\frac{a}{r})^{6}] [/math] [math] \hence p_0 = \frac{6MD}{dV_0\theta^{d}}(\theta^{-12}-\theta^{-6}) [/math]

Потенциал Ми

$$%\be{}
  \PI(r) =\frac{D}{n-m}\left[m\left(\frac{\DS a}{\DS
  r}\right)^{n}-n\left(\frac{\DS a}{\DS r}\right)^{m} \right]
  \hence
  p_0 =\frac{m n
  MD}{2d(n-m)V_0\theta^{d}}\,\(\theta^{-n}-\theta^{-m}\)%.
$$%\ee

%\item

Потенциал Морзе

$$%\be{}
  \PI(r) = D\[e^{2\alpha(a-r)}-2e^{\alpha(a-r)}\]
  \hence
  p_0 = \frac{\alpha a MD}{d V_0\theta^{d-1}}\,\[e^{2\alpha a(1-\theta)}-e^{\alpha a(1-\theta)}\]%.
$$%\ee

%\end{itemize}

Здесь $D$ --- энергия связи, $a$ --- длина связи,

$\alpha$ --- параметр, характеризующий ширину потенциальной ямы; $m, n$ --- параметры потенциала Ми.

Статьи

• Кривцов А. М., Кузькин В. А. Получение уравнения состояния идеальных кристаллов простой структуры // Механика твёрдого тела. — 2011. — № 3.

Ссылки

• Статья про уравнение Ми-Грюнайзена в Википедии