Hi everyone, I'm currently working on coding a finite element heat problem in Matlab that is as follows: I have a cross section of a beam made of an isotropic material (concrete) that is placed in a furnace, and the outside temperature is raised at a constant rate. I am interested in graphically displaying the temperature distribution at specific times to show how fast the interior of the cross section heats up. I have been able to make a Matlab model that works for boundary conditions in which the outside nodes of the cross section are held at a specific temperature, but I am not familiar enough with the heat problem to know what "force" vector I should characterize the outside temperature change as. The book that I am using (The Finite Element Method for Engineers by Huebner, talks about specified surface heating as a boundary condition, and characterizes this quantity as [Rq]. However, I am unsure if this is the proper quantity for my problem, as I am haven't taken thermodynamics and have only dealt with steady state heat problems in diff eq up to this point. Should I be dealing with heat flux? I've scoured Google,but have been unable to find any other paper that talks about specified surface heating. Conceptually, I am just unsure of what raising the ambient temperature around the cross section is classified as. Thanks!