Équation d'état de Mie-Grüneisen

L' équation d'état de Mie-Grüneisen est une équation d'état qui relie la pression et le volume d'un solide à une température donnée[1],[2]. Elle est utilisée pour déterminer la pression dans un solide parcouru par une onde de choc. La relation de Mie-Grüneisen est une forme spéciale du modèle de Grüneisen qui décrit l'effet de la modification du volume d'un réseau cristallin sur ses propriétés vibrationnelles. Plusieurs variantes de l'équation d'état de Mie-Grüneisen sont utilisées; elles reposent toutes sur la définition du paramètre de Grüneisen. Celui-ci peut-être définit à partir de grandeurs thermo-mécaniques macroscopique du matériau ou à partir des considérations reposant sur la physique statistique.

Une définition possible du paramètre de Grüneisen est la grandeur sans dimension suivante

Γ = α K T ρ c V {\displaystyle \Gamma ={\frac {\alpha K_{T}}{\rho c_{V}}}}

α {\displaystyle \alpha } est le coefficient de dilatation thermique, K T {\displaystyle K_{T}} est le module d'élasticité isostatique, ρ {\displaystyle \rho } la masse volumique du solide et c v {\displaystyle c_{v}} la capacité thermique massique à volume constant. Γ {\displaystyle \Gamma } est en général constant pour un matériau donné, et de l'ordre de l'unité. Cette définition est associée aux travaux de Grüneisen (1926) qui s'appuie sur le modèle d'Einstein (1907) des vibrations atomiques dans un solide et du modèle de Debye (1912) pour décrire le comportement de la capacité thermique des solides en fonction de la température. En utilisant les définitions de α = 1 V ( V T ) P {\displaystyle \alpha ={\frac {1}{V}}\left({\frac {\partial V}{\partial T}}\right)_{P}} et K T = V ( P V ) T {\displaystyle K_{T}=-V\left({\frac {\partial P}{\partial V}}\right)_{T}} ainsi que l'identité mathématique ( V T ) P ( P V ) T ( T P ) V = 1 {\displaystyle \left({\frac {\partial V}{\partial T}}\right)_{P}\left({\frac {\partial P}{\partial V}}\right)_{T}\left({\frac {\partial T}{\partial P}}\right)_{V}=-1} , on obtient que α K T = ( P T ) V {\displaystyle \alpha K_{T}=\left({\frac {\partial P}{\partial T}}\right)_{V}} . En intégrant par rapport à la température cette dernière relation à volume constant on obtient une relation entre la pression et l'énergie interne:

P P r e f = α K T d T = Γ ρ c v d T = Γ ( U U r e f ) / V {\displaystyle P-P_{ref}=\int \alpha K_{T}dT=\Gamma \int \rho c_{v}dT=\Gamma (U-U_{ref})/V}

P r e f {\displaystyle P_{ref}} et U r e f {\displaystyle U_{ref}} sont respectivement la pression et l'énergie interne d'un état de référence (qui peut être l'état dans lequel la température est de 0K) et V {\displaystyle V} le volume. Cette relation constitue la forme générale d'une équation d'état de Mie-Grüneisen.

On peut ré-arranger la définition du paramètre de Grüneisen, en faisant apparaître le volume: Γ V = α K T m c V = α K T C V {\displaystyle {\frac {\Gamma }{V}}={\frac {\alpha K_{T}}{mc_{V}}}={\frac {\alpha K_{T}}{C_{V}}}} et en première approximation, on considère que Γ V = Γ 0 V 0 {\displaystyle {\frac {\Gamma }{V}}={\frac {\Gamma _{0}}{V_{0}}}} est une constante, où Γ 0 {\displaystyle \Gamma _{0}} et V 0 {\displaystyle V_{0}} sont des constantes évaluées à pression nulle. Pour tenir compte des variations possible du paramètre de Grüneisen, on peut développer des modèles avec par exemples une correction d'ordre 1: Γ = Γ 0 V V 0 + a ( 1 V V 0 ) {\displaystyle \Gamma =\Gamma _{0}{\frac {V}{V_{0}}}+a(1-{\frac {V}{V_{0}}})}

De manière encore plus générale, lorsque l'on considère un solide soumis à de forte variation de pression (par exemple du à l'impact d'un projectile à très haute vitesse, ou obtenu à l'aide d'une impulsion laser), on peut utiliser la théorie des ondes de choc pour établir une équation d'état de Mie-Gruneisen généralisée. Plus spécifiquement, l'équation de Mie-Grüneisen la plus utilisée dans les logiciels de simulation numérique de thermo-mécanique est donnée par[3]

P ( μ , E ) = ρ 0 C 0 2 [ 1 + ( 1 Γ 0 2 ) μ a 2 μ 2 ] [ 1 ( S 1 1 ) μ S 2 μ 2 1 + μ S 3 μ 3 ( 1 + μ ) 2 ] 2 + ( Γ 0 + a μ ) E {\displaystyle P(\mu ,E)={\frac {\rho _{0}C_{0}^{2}\left[1+(1-{\frac {\Gamma _{0}}{2}})\mu -{\frac {a}{2}}\mu ^{2}\right]}{\left[1-(S_{1}-1)\mu -S_{2}{\frac {\mu ^{2}}{1+\mu }}-S_{3}{\frac {\mu ^{3}}{(1+\mu )^{2}}}\right]^{2}}}+(\Gamma _{0}+a\mu )E}

μ = ρ ρ 0 1 {\displaystyle \mu ={\frac {\rho }{\rho _{0}}}-1} , C 0 = K T ρ 0 {\displaystyle C_{0}={\sqrt {\frac {K_{T}}{\rho _{0}}}}} et a {\displaystyle a} le coefficient de la correction d'ordre 1 à la variation en volume du paramètre de Grüneisen et E {\displaystyle E} l'énergie interne volumique. Cette expression est uniquement valable pour un solide comprimé, c'est a dire pour μ > 0 {\displaystyle \mu >0} . Pour un matériau détendu, on doit utiliser la relation

P ( μ , E ) = ρ 0 C 0 2 μ + ( Γ 0 + a μ ) E {\displaystyle P(\mu ,E)=\rho _{0}C_{0}^{2}\mu +(\Gamma _{0}+a\mu )E}

Voir également

Notes

  • (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Mie–Grüneisen equation of state » (voir la liste des auteurs).

Références

  1. Roberts, J. K., & Miller, A. R. (1954). Heat and thermodynamics (Vol. 4). Interscience Publishers.
  2. Burshtein, A. I. (2008). Introduction to thermodynamics and kinetic theory of matter. Wiley-VCH.
  3. Lukyanov, A.A., « An equation of state for anisotropic solids under shock loading », Eur. Phys. J. B, vol. 64,‎ , p. 159-164 (lire en ligne Accès payant)
  • icône décorative Portail de la physique