How do you solve the heat equation?
2 Answers
(to remember: q equals m-cat)
Explanation:
specific heat is the amount of heat required to raise the temperature of one gram of substance by one degree Celsius. It is different for every element/compound and is usually given.
Just in case this was a partial differential equations question instead... the heat equation as I know it is:
#kgrad^2u - (delu)/(delt) = 0# #" "" "" "# #0 < x < L# ,#t > 0#
#u(0,t) = u(L,t) = 0# #" "" "# #t > 0#
#u(x,0) = f(x)# #" "" "" "" "# #0 < x < L# where:
#grad^2# is the second order derivative in all dimensions considered. Its form depends on what coordinates you are in.#k# is a constant describing the extent of heat conduction.#u# is the function describing the current temperature.#(delu)/(delt)# is the change in temperature over time.#grad^2u# is the change in temperature throughout space.
In one dimension cartesian, it would be:
#bb(k(del^2u)/(delx^2) - (delu)/(delt) = 0)#
To solve this equation, we would apply separation of variables. Since there is a change in temperature,
#k(del^2(phi(x)G(t)))/(delx^2) = (del(phi(x)G(t)))/(delt)#
#kG(t)(d^2phi(x))/(dx^2) = phi(x)(dG(t))/(dt)#
#1/(phi(x))(d^2phi(x))/(dx^2) = 1/(kG(t))(dG(t))/(dt) = -lambda#
Since the two sides of the equation each depend on only one variabe now, they are each equal to the same constant, known as the "separation constant". We call it
Now we have two ordinary differential equations to solve.
#1/(phi(x))(d^2phi(x))/(dx^2) = -lambda#
#1/(kG(t))(dG(t))/(dt) = -lambda#
The first:
#(d^2phi(x))/(dx^2) + lambdaphi = 0#
We assume the solution
#color(green)(phi(x) = Acos(sqrt(lambda)x) + Bsin(sqrt(lambda)x))#
The second equation gives:
#(dG(t))/(dt) = -lambdakG#
#int 1/GdG = -int lambdakdt#
#lnG(t) = -lambdakt + C#
#color(green)(G(t) = Ae^(-lambdakt))# where
#A = e^C# .
Now that we have two ordinary differential equations solutions, we apply boundary conditions to the space equation to get:
#u(0,t) = u(L,t) = 0#
#=> phi(0) = Acos(sqrt(lambda)0) + Bsin(sqrt(lambda)0) = 0#
Since
#phi(x) = Bsin(sqrt(lambda)x)#
The other boundary condition gives:
#phi(L) = Bsin(sqrt(lambda)L) = 0#
To get this,
#color(blue)(u_n(x,t) = sum_(n=1)^(oo) B_n sin((npix)/L)e^(-((npix)/L)^2kt))# where
#B_n = 2/Lint_(0)^(L) phi(x)sin((npix)/L)# .