Comment Jane/Lonny Eachus goes Sky Dragon Slayer (Score 1) 708
Now that we've agreed on the inner shell temperature of ~149.9F, let's take the last step. Calculate the enclosed source temperature.
Draw a boundary just inside the inner surface of the enclosing shell. Because nothing in the boundary is changing with time, power in = power out. The same constant electrical power flows in as before the shell was added. Net radiative power flows out from the source to the enclosing shell's inner surface.
As before, that net radiative power is described by Wikipedia’s equation which accounts for areas and view factors.
#Completely surrounded by shell with finite conductivity.
var('sigma T_c T_h A_c A_h F_hc power epsilon_h epsilon_c')
eq1 = power == sigma*(T_h^4 - T_c^4)/((1-epsilon_h)/(epsilon_h*A_h) + 1/(A_h*F_hc) + (1-epsilon_c)/(epsilon_c*A_c))
soln4 = solve(eq1.subs(T_c=338.629929346551,power=15028.4258648090,sigma=5.670373e-8,epsilon_h=0.11,epsilon_c=0.11, F_hc=1, A_h=510.064471909788, A_c=511.185932522526),T_h)
soln4[0].rhs().n()
... Please explain what calculations you are using where, because I find it hard to tell the Sage-formatted calculations apart. [Jane Q. Public, 2014-09-08]
The first line "var('sigma..." declares my variables.
The line "eq1 = power == sigma..." is my "power in = power out" equation using Wikipedia's equation for net radiative power.
The next line plugs in all the relevant variables and solves it for the enclosed source temperature T_h.
The last line displays the answer.
So I've described my method for calculating the enclosed source temperature from start to finish. Before I post that final answer, can we agree with my method? If not, could you please describe your method?