blockDebyeEoS

Energy density:

e = 3 n k T D₃(x)

where n is the number density of atoms, D₃(x) is the third Debye function, x=T_D/T, and T_D is the Debye temperature:

T_D = [h c_s / (2 k)] (6 n / π)^1/3 

c_s being the sound speed in the material (value obtainable via the parameter, *.c_s).

With proper values of T_D (in K), the atomic mass (in Da), and the mass density (as parameters), this block will provide the energy density (in J m^-3) as a function of the input temperature (in K): e(T). The specific heat capacity, c_v(T) (in J kg^-1 K^-1), is also calculated from the derivative wrt T:

c_v = [e/T - 9 n k (B(x)-D₃(x))] / ρ = n k [12 D₃(x) - 9 B(x)] / ρ