Code for simple heat transfer modeling can be found here: https://git.ligo.org/voyager/mariner40//blob/master/CryoEngineering/qil_simple_heat_transfer.ipynb
My original 1D cooldown script modeled conductive cooling along the copper braid as: Pcool = k * A/L * (T  T_set); where Pcool is the cooling power, k is the thermal conductivity of copper, A is the crosssectional area of the braid, and L is the length of the braid. I changed the code to instead use the tabulated power vs. temperature points for the CH104 coldhead, taken from Paco's script qil_heat_estimate.ipynb. The first figure compares the interpolated curves from the tabulated values (at 50Hz and 60Hz operation), to the original conductive transfer model (50K setpoint). The original conductive powertemp relationship is linear, which overestimates the cooling power at high temperatures. Switching to the tabulated points results in more realistic model. Moving forward, I intend to use the 50Hz interpolated curve.
The script considers radiative heating to the coldplate from the the chamber bottom (rough aluminum) and the outer shield (coated in aluminum foil). It assumes over a long period of time that the inner shield and coldplate temperatures are equal.
The second figure shows the results of this model alongside the actual coolddown data extracted from the CTC100. It is clear that the model is not accounting for additional radiative heat sources that would explain the slower cooldown and higher final temperature. Adding in model complexity is my current focus.
