def cv(T):
'''
Thermal conductivity paralell to fibers
source: Runyan and Jones 2008 (arxiv:0806.1921)
'''
if (not isinstance(T, np.ndarray)):
T = np.array(T)
if (np.any(T > 80)):
print 'Warning: Temperature too high (valid range is 2 to 80 K)'
if (np.any(T < 2)):
print 'Warning: Temperature too low (valid range is 2 to 80 K)'
coeff = [
-1.37737, -3.40668, 20.5842, -53.1244, 73.2476, -57.6546, 26.1192,
-6.34790, 0.640331
]
if len(np.shape(T)) == 0:
# T is a scalar
if T > 80:
return bf.NIST_10(80., coeff)
else:
return bf.NIST_10(T, coeff)
inds = np.where(T > 80.)[0]
val = np.zeros(np.shape(T))
val[inds] = bf.NIST_10(80., coeff)
val[T < 80] = bf.NIST_10(T[T < 80], coeff)
return val
def cv_normal(T, verbose=True):
'''
Thermal conductivity in W / m K
'''
if (not isinstance(T, np.ndarray)):
T = np.array(T)
if np.any(T > 300) and verbose:
print 'Warning: Temperature too high (valid range is .3 to 300 K)'
if np.any(T < .3) and verbose:
print 'Warning: Temperature too low (valid range is .3 to 300 K)'
if np.any(T < 10) and np.any(T > 10) and vebose:
print 'Warning: Using two different fits! Check your results'
lowT = np.where(T < 10)
highT = np.where(T >= 10)
k = np.zeros(np.shape(T))
if len(lowT) > 0:
alpha = 12.8e-3
Beta = 2.41
gamma = -9.21
n = .222
k[lowT] = bf.kt(T[lowT], alpha, Beta, gamma, n)
if len(highT) > 0:
coeff = [
-4.1236, 13.788, -26.068, 26.272, -14.663, 4.4954, -.6905, .0397
]
k[highT] = bf.NIST_10(T[highT], coeff)
return k
def cv_warp(T, verbose=True):
if np.any(T > 300) and verbose:
print 'Warning: Temperature too high (valid range is 12 to 300 K)'
if np.any(T < 12) and verbose:
print 'Warning: Temperature too low (valid range is 12 to 300 K)'
coeff = [
-2.64827, 8.80228, -24.8998, 41.1625, -39.8754, 23.1778, -7.95635,
1.48806, -0.11701
]
return bf.NIST_10(T, coeff)
def cv(T):
if (not isinstance(T, np.ndarray)):
T = np.array(T)
if (np.any(T > 300)):
print 'Warning: Temperature too high (valid range is 1 to 300 K)'
if (np.any(T < 1)):
print 'Warning: Temperature too low (valid range is 1 to 300 K)'
coeff = [
22.401433, -141.1343, 394.95461, -601.1537, 547.83202, -305.9969,
102.38656, -18.81023, 1.4576882
]
return bf.NIST_10(T, coeff)
def cv(T):
'''
Thermal conductivity paralell to fibers
source: Runyan and Jones 2008 (arxiv:0806.1921)
'''
if (not isinstance(T, np.ndarray)):
T = np.array(T)
if (np.any(T > 300)):
print 'Warning: Temperature too high (valid range is 1 to 300 K)'
if (np.any(T < 1)):
print 'Warning: Temperature too low (valid range is 1 to 300 K)'
coeff = [0.07918, 1.0957, -0.07277, 0.08084, 0.02803, -0.09464,
0.04179, -0.00571]
return bf.NIST_10(T, coeff)
def cv(T):
'''
'''
if (not isinstance(T, np.ndarray)):
T = np.array(T)
if (np.any(T > 300)):
print 'Warning: Temperature too high (valid range is 4 to 300 K)'
if (np.any(T < 4)):
print 'Warning: Temperature too low (valid range is 4 to 300 K)'
coeff = [
2.7380, -30.677, 89.430, -136.99, 124.69, -69.556, 23.320, -4.3135,
0.33829
]
return bf.NIST_10(T, coeff)