Level 5 — MasteryPower, Thermal & Reliability

Power, Thermal & Reliability

90 minutes60 marksprintable — key stays hidden on paper

Level 5 — Mastery (cross-domain: physics + mathematics + coding) Time limit: 90 minutes Total marks: 60

Notation: dynamic power Pdyn=αCV2fP_{dyn}=\alpha C V^2 f; static (leakage) power Pstat=VIleakP_{stat}=V\,I_{leak} with Ileak=I0eVth/(nVT)I_{leak}=I_0\,e^{-V_{th}/(n\,V_T)}, thermal voltage VT=kT/qV_T=kT/q. Steady-state die temperature T=Tamb+PθJAT=T_{amb}+P\,\theta_{JA}.


Question 1 — DVFS, energy-optimal operating point (24 marks)

A core executes a fixed workload of WW clock cycles. Its power model is P(V,f)=αCV2fdynamic+VIleakstatic,P(V,f)=\underbrace{\alpha C V^2 f}_{\text{dynamic}}+\underbrace{V I_{leak}}_{\text{static}}, and, in the region of interest, the maximum stable frequency scales with voltage as f=β(VVth).f = \beta\,(V - V_{th}).

(a) The task must finish in a deadline of DD seconds. Show that operating at the lowest frequency that still meets the deadline (i.e. f=W/Df = W/D) minimises the dynamic energy per task, and derive an expression for that dynamic energy EdynE_{dyn} in terms of α,C,W,β,Vth\alpha, C, W, \beta, V_{th} and DD. (7)

(b) Static energy per task is Estat=Pstat(W/f)E_{stat}=P_{stat}\cdot(W/f). Explain, using your EdynE_{dyn} and EstatE_{stat} expressions, why lowering ff below some point increases total energy per task, and state the qualitative condition (in words + inequality on the derivative) that defines the energy-optimal frequency. (6)

(c) Take numbers: αC=2×109 F\alpha C = 2\times10^{-9}\ \mathrm{F}, β=5×109 Hz/V\beta = 5\times10^{9}\ \mathrm{Hz/V}, Vth=0.3 VV_{th}=0.3\ \mathrm{V}, Ileak=0.8 AI_{leak}=0.8\ \mathrm{A} (assume approximately voltage-independent for this part), W=3×109W = 3\times10^{9} cycles. Compute the total energy per task at f1=1 GHzf_1=1\ \mathrm{GHz} and at f2=2 GHzf_2=2\ \mathrm{GHz}. Which is more energy efficient, and by what percentage? (8)

(d) State one physical mechanism that eventually breaks the "lower ff always saves dynamic energy" assumption at very low voltage. (3)


Question 2 — Thermal budget, throttling & dark silicon (22 marks)

A chip has TDP =95 W=95\ \mathrm{W}, package thermal resistance θJA=0.35 C/W\theta_{JA}=0.35\ \mathrm{^\circ C/W}, ambient Tamb=35 CT_{amb}=35\ \mathrm{^\circ C}, and a throttle trip point Tmax=100 CT_{max}=100\ \mathrm{^\circ C}.

(a) Compute the steady-state junction temperature at TDP. Is the cooling solution adequate? Compute the maximum sustainable power PmaxP_{max} before the trip point is reached. (6)

(b) The chip has 16 identical cores. Each active core dissipates 9 W9\ \mathrm{W}. Uncore/leakage floor is 12 W12\ \mathrm{W} regardless. Compute the maximum number of cores that can run simultaneously without throttling. Explain how this quantitatively illustrates the dark silicon problem. (6)

(c) A transient burst raises power to 130 W130\ \mathrm{W}. The die has heat capacity Cth=0.9 J/CC_{th}=0.9\ \mathrm{J/^\circ C}. Model temperature by CthdTdt=PTTambθJAC_{th}\dfrac{dT}{dt}=P-\dfrac{T-T_{amb}}{\theta_{JA}}. Starting from the steady state of part (a), derive T(t)T(t) and compute the time until the throttle trips at 100 C100\ \mathrm{^\circ C}. (10)


Question 3 — Electromigration & voltage droop (14 marks)

(a) Black's equation gives median-time-to-failure MTTF=AJneEa/(kT)\mathrm{MTTF}=A\,J^{-n}\,e^{E_a/(kT)}, with n=2n=2, Ea=0.7 eVE_a=0.7\ \mathrm{eV}, k=8.617×105 eV/Kk=8.617\times10^{-5}\ \mathrm{eV/K}. A wire has MTTF=10\mathrm{MTTF}=10 years at T1=350 KT_1=350\ \mathrm{K} and fixed current density. Find the MTTF at T2=380 KT_2=380\ \mathrm{K}. State the reliability implication for DVFS/turbo. (8)

(b) During a fast current step ΔI=20 A\Delta I=20\ \mathrm{A} over Δt=2 ns\Delta t = 2\ \mathrm{ns} the power-delivery network has effective inductance L=15 pHL=15\ \mathrm{pH}. Compute the Ldi/dtL\,di/dt voltage droop. If the supply is 1.0 V1.0\ \mathrm{V} with a ±5%\pm5\% tolerance, does this violate the rail? State how a decoupling capacitor of C=100 nFC=100\ \mathrm{nF} mitigates it (mechanism, one sentence). (6)

Answer keyMark scheme & solutions

Question 1

(a) (7 marks) Time to finish WW cycles at frequency ff is t=W/ft=W/f. Dynamic energy Edyn=Pdynt=αCV2fWf=αCV2WE_{dyn}=P_{dyn}\cdot t=\alpha C V^2 f\cdot\frac{W}{f}=\alpha C V^2 W. (2) — energy depends only on VV, not directly on ff for a fixed cycle count. Since f=β(VVth)V=Vth+f/βf=\beta(V-V_{th})\Rightarrow V=V_{th}+f/\beta, lower ff ⇒ lower VV ⇒ lower EdynE_{dyn} (monotone in VV). (2) The lowest ff meeting deadline DD is f=W/Df=W/D (finishing exactly on time). (1) Substituting: Edyn=αCW(Vth+WβD)2.E_{dyn}=\alpha C W\left(V_{th}+\frac{W}{\beta D}\right)^{2}. (2)

(b) (6 marks) Estat=PstatW/f=VIleakW/fE_{stat}=P_{stat}\cdot W/f = V I_{leak}W/f. As f0f\to0, t=W/ft=W/f\to\infty, so the core leaks for longer and EstatE_{stat}\to\infty. (3) Total E=Edyn(f)+Estat(f)E=E_{dyn}(f)+E_{stat}(f): EdynE_{dyn} decreases with ff but EstatE_{stat} increases as 1/f1/f-ish (leakage integrated over longer time). The optimum is where the marginal saving in dynamic energy equals the marginal rise in static energy: dEtotdf=0,i.e. dEdyndf=dEstatdf.\frac{dE_{tot}}{df}=0,\quad \text{i.e. } \frac{dE_{dyn}}{df}=-\frac{dE_{stat}}{df}. (3)

(c) (8 marks) Voltage at each frequency: V=Vth+f/βV=V_{th}+f/\beta.

  • f1=1×109f_1=1\times10^9: V1=0.3+(1×109)/(5×109)=0.3+0.2=0.5 VV_1=0.3+ (1\times10^9)/(5\times10^9)=0.3+0.2=0.5\ \mathrm{V}.
  • f2=2×109f_2=2\times10^9: V2=0.3+(2×109)/(5×109)=0.3+0.4=0.7 VV_2=0.3+ (2\times10^9)/(5\times10^9)=0.3+0.4=0.7\ \mathrm{V}.

Edyn=αCV2W=2×109V23×109=6V2E_{dyn}=\alpha C V^2 W = 2\times10^{-9}\cdot V^2\cdot 3\times10^9 = 6\,V^2 J. Estat=VIleakW/f=V0.83×109/fE_{stat}=V I_{leak}\,W/f = V\cdot0.8\cdot 3\times10^9/f.

At f1f_1: Edyn=6(0.25)=1.5E_{dyn}=6(0.25)=1.5 J; Estat=0.50.83=1.2E_{stat}=0.5\cdot0.8\cdot3=1.2 J; E1=2.7E_1=2.7 J. (3) At f2f_2: Edyn=6(0.49)=2.94E_{dyn}=6(0.49)=2.94 J; Estat=0.70.81.5=0.84E_{stat}=0.7\cdot0.8\cdot1.5=0.84 J; E2=3.78E_2=3.78 J. (3)

f1f_1 is more efficient. Difference =(3.782.7)/3.78=28.6%=(3.78-2.7)/3.78=28.6\% lower energy at f1f_1. (2) (Accept "E2E_2 is (3.782.7)/2.7=40%(3.78-2.7)/2.7=40\% higher than E1E_1" if stated relative to E1E_1.)

(d) (3 marks) At very low voltage leakage dominates and total energy rises (race-to-idle failure); also VV cannot drop below VthV_{th} (or the near-threshold region) where the f=β(VVth)f=\beta(V-V_{th}) model breaks down and circuit reliability/variation and subthreshold conduction rise sharply. (Any one valid mechanism.)


Question 2

(a) (6 marks) T=Tamb+PθJA=35+95×0.35=35+33.25=68.25 CT=T_{amb}+P\theta_{JA}=35+95\times0.35=35+33.25=68.25\ \mathrm{^\circ C}. (2) Well below 100 C100\ \mathrm{^\circ C} ⇒ cooling adequate. (2) Pmax=(TmaxTamb)/θJA=(10035)/0.35=185.7 WP_{max}=(T_{max}-T_{amb})/\theta_{JA}=(100-35)/0.35=185.7\ \mathrm{W}. (2)

(b) (6 marks) Power with nn cores: P=12+9nPmax=185.7P=12+9n\le P_{max}=185.7. (2) 9n173.7n19.39n\le173.7\Rightarrow n\le19.3; but only 16 cores exist ⇒ all 16 can run: P=12+144=156 W<185.7P=12+144=156\ \mathrm{W}<185.7. (2) Here thermal ceiling permits all cores → not thermally limited. To illustrate dark silicon, if θJA\theta_{JA} or per-core power were higher (e.g. per-core 12 W: 12+12n185.7n14.512+12n\le185.7\Rightarrow n\le14.5, so only 14 of 16 usable), the remaining cores must stay dark. The point: with more cores than the power/thermal budget can sustain, a fraction must be powered off ("dark silicon"). (2)

(c) (10 marks) ODE: CthT˙=P(TTamb)/θJAC_{th}\dot T=P-(T-T_{amb})/\theta_{JA}. Let τ=CthθJA=0.9×0.35=0.315 s\tau=C_{th}\theta_{JA}=0.9\times0.35=0.315\ \mathrm{s}. (2) Steady-state target for burst: T=Tamb+PθJA=35+130×0.35=35+45.5=80.5 CT_\infty=T_{amb}+P\theta_{JA}=35+130\times0.35=35+45.5=80.5\ \mathrm{^\circ C}? — wait, 130×0.35=45.5130\times0.35=45.5, so T=80.5 CT_\infty=80.5\ \mathrm{^\circ C}.

Since T=80.5 C<100 CT_\infty=80.5\ \mathrm{^\circ C}<100\ \mathrm{^\circ C}, the throttle never trips — the burst settles below the trip point. (3) Solution form: T(t)=T+(T0T)et/τT(t)=T_\infty+(T_0-T_\infty)e^{-t/\tau} with T0=68.25 CT_0=68.25\ \mathrm{^\circ C}: T(t)=80.5+(68.2580.5)et/0.315=80.512.25et/0.315.T(t)=80.5+(68.25-80.5)e^{-t/0.315}=80.5-12.25\,e^{-t/0.315}. (3) Setting T=100T=100 has no solution (T<80.5T<80.5 always). Conclusion: trip time = ∞ (never trips); the transient is thermally safe. (2) (Full credit for correctly deriving T=80.5 CT_\infty=80.5\ \mathrm{^\circ C} and concluding no trip; recognising the trick is the mastery test.)


Question 3

(a) (8 marks) Same JJ, so MTTF ratio depends only on temperature: MTTF2MTTF1=exp ⁣[Eak(1T21T1)].\frac{\mathrm{MTTF}_2}{\mathrm{MTTF}_1}=\exp\!\left[\frac{E_a}{k}\left(\frac1{T_2}-\frac1{T_1}\right)\right]. (3) Eak=0.7/8.617×105=8124.6 K\frac{E_a}{k}=0.7/8.617\times10^{-5}=8124.6\ \mathrm{K}. 13801350=0.00263160.0028571=2.2556×104 K1\frac1{380}-\frac1{350}=0.0026316-0.0028571=-2.2556\times10^{-4}\ \mathrm{K^{-1}}. Exponent =8124.6×(2.2556×104)=1.8326=8124.6\times(-2.2556\times10^{-4})=-1.8326. Ratio =e1.8326=0.1599=e^{-1.8326}=0.1599. (3) MTTF2=10×0.1599=1.60\mathrm{MTTF}_2=10\times0.1599=1.60 years. (1) Implication: raising temperature (turbo/high V,f) sharply shortens interconnect lifetime — reliability budget limits sustained boost. (1)

(b) (6 marks) Droop =Ldi/dt=15×1012×(20/2×109)=15×1012×1010=0.15 V=150 mV=L\,di/dt=15\times10^{-12}\times(20/2\times10^{-9})=15\times10^{-12}\times10^{10}=0.15\ \mathrm{V}=150\ \mathrm{mV}. (3) Tolerance is ±5%\pm5\% of 1.0 V=±50 mV1.0\ \mathrm{V}=\pm50\ \mathrm{mV}. 150 mV50 mV150\ \mathrm{mV}\gg50\ \mathrm{mV}violates the rail. (2) Decap mechanism: the capacitor stores charge locally and supplies the fast transient current, holding VV up during the di/dtdi/dt event before the (inductive) main supply responds. (1)

[
  {"claim":"Q1c E1 total energy = 2.7 J", "code":"V1=Rational(1,2); Edyn1=6*V1**2; Estat1=V1*Rational(8,10)*3; result=(Edyn1+Estat1==Rational(27,10))"},
  {"claim":"Q1c E2 total energy = 3.78 J", "code":"V2=Rational(7,10); Edyn2=6*V2**2; Estat2=V2*Rational(8,10)*Rational(15,10); E2=Edyn2+Estat2; result=(E2==Rational(378,100))"},
  {"claim":"Q1c f1 is ~28.6% lower energy than f2", "code":"E1=Rational(27,10); E2=Rational(378,100); frac=(E2-E1)/E2; result=(abs(float(frac)-0.2857)<0.001)"},
  {"claim":"Q2a Pmax = 185.71 W", "code":"Pmax=(100-35)/Rational(35,100); result=(abs(float(Pmax)-185.714)<0.01)"},
  {"claim":"Q2c burst steady state 80.5C never reaches 100C", "code":"Tinf=35+130*Rational(35,100); result=(Tinf==Rational(805,10) and float(Tinf)<100)"},
  {"claim":"Q3a MTTF2 = 1.60 years", "code":"import sympy as sp; Ea=0.7; k=8.617e-5; ratio=sp.exp((Ea/k)*(1/380-1/350)); m2=10*ratio; result=(abs(float(m2)-1.60)<0.02)"},
  {"claim":"Q3b droop = 0.15 V", "code":"droop=15e-12*(20/2e-9); result=(abs(droop-0.15)<1e-9)"}
]