Thermal properties temperature dependent

[RemDelaporteMathurin]

Coucou

[shimwell]

We would need to change the following functions in 3-MOTtiedup.py

calculate_specific_heat()
calculate_density_heat()
calculate_thermal_conductivity_heat()

[shimwell]

Perhaps as a first attempt we could use some simple SciPy interpolation as demonstrated here

https://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.interp1d.html#scipy.interpolate.interp1d

[RemDelaporteMathurin]

Could this be done in an external library ?

[RemDelaporteMathurin]

Hey @shimwell look what I found!

ft11_magnetostatics.py.zip

class Permeability(Expression):
def init(self, markers, **kwargs):
self.markers = markers def eval_cell(self, values, x, cell):
if self.markers[cell.index] == 0:
values[0] = 4pi1e-7 # vacuum
elif self.markers[cell.index] == 1:
values[0] = 1e-5 # iron (should really be 6.3e-3)
else:
values[0] = 1.26e-6 # copper

mu = Permeability(markers, degree=1)