-
Notifications
You must be signed in to change notification settings - Fork 41
/
SLB_2005.py
123 lines (103 loc) · 3.04 KB
/
SLB_2005.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
# This file is part of BurnMan - a thermoelastic and thermodynamic toolkit for the Earth and Planetary Sciences
# Copyright (C) 2012 - 2017 by the BurnMan team, released under the GNU
# GPL v2 or later.
"""
SLB_2005
^^^^^^^^
Minerals from Stixrude & Lithgow-Bertelloni 2005 and references therein
"""
from __future__ import absolute_import
from ..classes import mineral_helpers as helpers
from ..classes.mineral import Mineral
class stishovite(Mineral):
def __init__(self):
self.params = {
"formula": {"Si": 1.0, "O": 2.0},
"equation_of_state": "slb3",
"V_0": 14.02e-6,
"K_0": 314.0e9,
"Kprime_0": 4.4,
"G_0": 220.0e9,
"Gprime_0": 1.6,
"molar_mass": 0.0601,
"n": 3,
"Debye_0": 1044.0,
"grueneisen_0": 1.34,
"q_0": 2.4,
"eta_s_0": 5.0,
}
Mineral.__init__(self)
class periclase(Mineral):
def __init__(self):
self.params = {
"formula": {"Mg": 1.0, "O": 1.0},
"equation_of_state": "slb3",
"V_0": 11.24e-6,
"K_0": 161.0e9,
"Kprime_0": 3.8,
"G_0": 131.0e9,
"Gprime_0": 2.1,
"molar_mass": 0.0403,
"n": 2,
"Debye_0": 773.0,
"grueneisen_0": 1.5,
"q_0": 1.5,
"eta_s_0": 2.8,
}
Mineral.__init__(self)
class wuestite(Mineral):
def __init__(self):
self.params = {
"equation_of_state": "slb3",
"formula": {"Fe": 1.0, "O": 1.0},
"V_0": 12.06e-6,
"K_0": 152.0e9,
"Kprime_0": 4.9,
"G_0": 47.0e9,
"Gprime_0": 0.7,
"molar_mass": 0.0718,
"n": 2,
"Debye_0": 455.0,
"grueneisen_0": 1.28,
"q_0": 1.5,
"eta_s_0": 0.8,
}
Mineral.__init__(self)
class mg_perovskite(Mineral):
def __init__(self):
self.params = {
"formula": {"Mg": 1.0, "Si": 1.0, "O": 3.0},
"equation_of_state": "slb3",
"V_0": 24.45e-6,
"K_0": 251.0e9,
"Kprime_0": 4.1,
"G_0": 175.0e9,
"Gprime_0": 1.7,
"molar_mass": 0.1000,
"n": 5,
"Debye_0": 1070.0,
"grueneisen_0": 1.48,
"q_0": 1.4,
"eta_s_0": 2.6,
}
Mineral.__init__(self)
class fe_perovskite(Mineral):
def __init__(self):
self.params = {
"formula": {"Fe": 1.0, "Si": 1.0, "O": 3.0},
"equation_of_state": "slb3",
"V_0": 25.48e-6,
"K_0": 281.0e9,
"Kprime_0": 4.1,
"G_0": 138.0e9,
"Gprime_0": 1.7,
"molar_mass": 0.1319,
"n": 5,
"Debye_0": 841.0,
"grueneisen_0": 1.48,
"q_0": 1.4,
"eta_s_0": 2.1,
}
Mineral.__init__(self)
mg_bridgmanite = mg_perovskite
fe_bridgmanite = fe_perovskite