How can vacuum heating be calculated for substrates in a chamber?

5 replies · 3K views · last reply
ektov_konstantin
Messages
5
Reaction score
0
TL;DR
I have a physically correct approach but I get incorrect result and need some help!
Hello, everyone!

I am calculating the time of substrates' heating in the vacuum chamber.

It is a simple mathematical model to calculate roughly how much time does it take to heat polycor substrates in the vacuum chamber.

At least all the tube heater`s power heats the substrates (this is because the geometry of the vacuum chamber and beacuse all inner surface is polished).

So our goal is to get the time when conductive heating will be equal to the radiative heat losses.

Look at the scheme.

Heat balance diagram: P*dt arrives from the tube heater, radiative losses epsilon sigma0 S T^4 dt leave the surface, and cmdT is stored in the Polycor substrates.


So the differential equation is .
Differential equation: cmdT minus epsilon times sigma_0 times S times T to the fourth power equals P times dt.


And dividing the variables we get
Equation t(T_0) = integral from 293 K to T_0 of cm divided by (P + epsilon sigma_0 S T^4) dT.
, where T_0 is the Temperature when conductive heating equals radiative losses.

If P=800 W, I have this dependence...

Line graph of t(temperature, P_TEH) in seconds against temperature in degrees Kelvin from about 15,000 to 100,000; the curve rises steeply from 0.65310 and levels off near 0.65319, with a small vertical spike near 52,000.


and as you can see, it is not physically right beacuse of two reasons. First - why there is a horisontal asymptote (it must be vertical!) ? Second - the heating is very-very fast...

For instance, I have calculated the heating time if there are no radiative losses.

t = c*m*(T_0-293)/W

and for W=800 W, T_0=623 K (technological process` temperature) we get t=0,205 seconds - and it is rubbish too I think...

Additionaly I add the constants and variables that I have used in my MathCAD app.

Screenshot of MathCAD worksheet listing constants: c_Al2O3 := 750 J/kg*degC, rho_Al2O3 := 3.76e3 kg/m^3, epsilon := 0.1, sigma_0 := 5.67e-8 W/(m^2*K^4), S_Al2O3 := 60*48e-6 = 2.88e-3 m^2, substrates := 23 items, delta := 2e-3 mm, and P_TEH := 800 W.


MathCAD expression: m_Al2O3 := rho_Al2O3 · S_Al2O3 · delta · substrates = 0.498 kg
 
  1. Is ALL of the thermal energy of the heater focused on the work piece?
    • How?
  2. What about conductive losses from the heater thru its mounting?
  3. Is the work piece suspended and not in contact with a support that gives a path for conductive losses?
    • How?

The graph shows a temperature of 1×105K (100 000 Kelvin); is that realistic? (Have you obtained a Star for a heat source?)

Graphs generally have the independent variable (Time) on the X-axis, and the dependent variable (Temperature) on the Y-axis; that makes them easier for most of us to interpret.

Anxiously awaiting clarifications of the above!
Tom
 
Tom.G said:
  1. Is ALL of the thermal energy of the heater focused on the work piece?
    • How?
  2. What about conductive losses from the heater thru its mounting?
  3. Is the work piece suspended and not in contact with a support that gives a path for conductive losses?
    • How?

The graph shows a temperature of 1×105K (100 000 Kelvin); is that realistic? (Have you obtained a Star for a heat source?)

Graphs generally have the independent variable (Time) on the X-axis, and the dependent variable (Temperature) on the Y-axis; that makes them easier for most of us to interpret.

Anxiously awaiting clarifications of the above!
Tom
Hello, Tom!

Thank you for your reply.

Imagine, that this is abstract mathematical model.
I extremely agree with you that...
1. There are conductive losses in the heater mounting and to the polycor substrates basement (this is round plate)
2. Not all the power from the heater heats the substrates but if I have to make this assumption to ease my calculation. It is affordable because all inner surface of the chamber is polished.

Later I will make the heat transfer calculation in ANSYS application. It will consider all radiation heating and conductive losses.

BUT my mathematical calculation gives unphysical and illogical results.

Just imagine that this is an abstract body... And losses are radiative. Temperature must come to stationary value. But it doesn't. What I get is that time (!) becomes static. Time stops. The question is why?!

My differential equation is solved correct - I have rechecked it.
 
Lnewqban said:
Your equation seems to be correct.
The calculations may have unit's problems if the result is not logical.
Your diagram should look something like this:

View attachment 297387
Yes. That`s right. Thank you for your answer. But I`ve checked units. They are all in SI... May I show you my MathCAD file?
 
ektov_konstantin said:
Yes. That`s right. Thank you for your answer. But I`ve checked units. They are all in SI... May I show you my MathCAD file?
You are welcome.
Please, do.