Example: Heater on a substrate (single layer)

We will compute the temperature rise 2ω of a simple heater made of gold on a crystalline 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] = current^2 * resistance
p = 0.030^2*22.11
# Heater thermal resistance
ρh = [0.0 0.0]
# Interface thermal resistances
thresistances = [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(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; int_limit=5.0e6)

Plot results

plot(TOMPlot(), sol)
gui()