Example: Heater on two layers

We will compute the temperature rise 2ω of a heater made of gold on a thin film deposited atop a crystalline silicon substrate system. Here, the substrate acts as a semi-infinite medium, and we SI units are used.

The complete code can be found here.

Load modules

using ThinFilmsTools
using Plots
pyplot()

Input data

We define the parameters needed to construct the structures of the layers.

# Half-width of heater line [m]
b = (12.5/2)*1e-6
# Length of heater line [m]
l = 1.0e-3
# Range of frequencies [Hz]
f = exp10.(LinRange(0, 9, 1500))
# Power [W]
p = 0.030^2*22.11
# Heater thermal resistance
ρh = [0.0 0.0]
# Interface thermal resistances
thresistances = [0.0 0.0]

Wrap information into structures

We construct the system of layers where the heater is at the top and the substrate at the bottom.

layers = [
   LayerTOM(310.0, 1.0, 0.2e-6, 2.441e6), # heater
   LayerTOM(0.1, 1.0, 1.0e-6, 2320*700*0.1), # specimen_1
   LayerTOM(160.0, 1.0, 525.0e-6, 2320*700.0), # substrate
]
hgeometry = HeaterGeometry(b, l, ρh)
source = Source(p, f)

Solve

sol = three_omega(layers, hgeometry, source, thresistances)

Plot results

plot(TOMPlot(), sol)
gui()