Section 10.2 Heat Conduction in a Cylinder with an Electrical Heat Source Following the same steps as BS&L used. > restart; > A:=r-> 2*Pi*r*L; area of conduction > Q:=r-> A(r)*qr(r); Heat conducted across cylinder of radius r in the direction of increasing r: qr is the flux of energy. > Q(r); the rate of thermal energy into cylindrical shell at r: eq. 10.2-2 > Q(r + dr); rate of thermal energy out across cyllindrical surface at r+dr: eq. 10.2-3 > V:=r->A(r)*dr; Volume of shell > Epro:=(r)->V(r)*Se; Thermal production rate by electrical current: eq. 10.2-4 > Epro(r); > LHS:=limit((Q(r+dr)-Q(r))/dr,dr=0); Left hand side of eq. 10.2-5 > de:=simplify(r*(LHS=Epro(r)/dr)/A(r)); energy balance, out - in = energy production by electrical dissipation giving 10.2-6 > sol:=dsolve({de,qr(0)=finite},qr(r)); Solving 10.2-6 with BC 10.2-8 to get 10.2-9 > assign(sol); Giving qr(r) a value > qr:=unapply(qr(r),r); Making qr a function of r > deFL:=qr(r)=-k*D(T)(r); Fourier's Law > sol:=dsolve({deFL,T(R)=T0},T(r)); solving the de with the bc 10.2-12 to get 10.2-13 > assign(sol);T:=unapply(T(r),r); > dT:=simplify(T(r)-T0); > dToverCon:=simplify((dT/(Se*R^2/(4*k))));dTmax:=simplify(T(0)-T0);dTavg:=simplify(int((T(r)-T0)*r,r=0...R)/int(r,r=0...R)); QatR:=Q(R); Our answers agree with eqs. 10.2-13, 14, 15 and 16
>