Assignment Solutions of Partial Difierential Equations

1 Assignment Solutions of Partial Differential EquationsWeijiu LiuDepartment of MathematicsUniversity of Central Arkansas201 Donaghey Avenue, Conway, AR 72035, USA1 Assignment Derive the heat equation for a rod assuming constant thermal properties with variablecross-sectional areaA(x) assuming no byAthe the cross-sectional quantities: Thermal energy densitye(x, t) = the amount of thermal energy per unit volume. Heat flux (x, t) = the amount of thermal energy flowing across boundaries per unitsurface area per unit time. Heat sourcesQ(x, t) = 0. Temperatureu(x, t). Specific heatc= the heat energy that must be supplied to a unit mass of a substanceto raise its temperature one unit.

2 Mass density (x) = mass per unit volume. Fourier s Law: the heat flux is proportional to the temperature gradient = K0 u.(1)Conservation of heat energy:Rate of change of heat energy in time = Heat energy flowing across boundaries per unittime + Heat energy generated insider per unit time heat energy =e(x, t)A(x) x. Heat energy flowing across boundaries per unit time = (x, t)A(x) (x+ x, t)A(x+ x).Then t[e(x, t)A(x) x] = (x, t)A(x) (x+ x, t)A(x+ x).Dividing it by xand letting xgo to zero giveA(x) e t= A(x) x (x) A x.(2)Heat energy per unit mass =c(x)u(x, t) A x. Soe(x, t)A(x) x=c(x)u(x, t) A(x) x,and thene(x, t) =c(x)u(x, t).

3 2It then follows from Fourier s law thatc A(x) u t=A(x) x(K0 u x)+K0 A x u x.(3)and then the heat equationA(x) u t=k(A(x) 2u x2+ A x u x),(4)wherek=K0c is the thermal Consider a thin one-dimensional rod without source of thermal energy whose lateralsurface is not insulated. Letw(x, t) dente the heat energy flowing out of the lateral sidesper unit surface area per unit time. Assume thatw(x, t) is proportional to the temperaturedifference between the rodu(x, t) and a known outside temperature (x, t). Derive theequation for the byAthe the cross-sectional area, andPthe lateral quantities: Thermal energy densitye(x, t) = the amount of thermal energy per unit volume.

4 Heat flux (x, t) = the amount of thermal energy flowing across boundaries per unitsurface area per unit time. Temperatureu(x, t). Specific heatc= the heat energy that must be supplied to a unit mass of a substanceto raise its temperature one unit. Mass density (x) = mass per unit volume. Fourier s Law: the heat flux is proportional to the temperature gradient = K0 u.(5)Conservation of heat energy:Rate of change of heat energy in time = Heat energy flowing across boundaries per unittime + Heat energy generated insider per unit time heat energy =e(x, t)A x. Heat energy flowing across boundaries per unit time = (x, t)A (x+ x, t)A.

5 Heat energy flowing out of the lateral sides per unit time =w(x, t)P x= [u(x, t) (x, t)]h(x)P x, whereh(x) is a t[e(x, t)A(x) x] = (x, t)A(x) (x+ x, t)A(x+ x) [u(x, t) (x, t)]h(x)P it byA xand letting xgo to zero give e t= x PA[u(x, t) (x, t)]h(x).(6)Heat energy per unit mass =c(x)u(x, t) A x. Soe(x, t)A x=c(x)u(x, t) A x,and thene(x, t) =c(x)u(x, t) .It then follows from Fourier s law thatc u t= x(K0 u x) PA[u(x, t) (x, t)]h(x).(7)42 Assignment Determine the equilibrium temperature distribution for a one-dimensional rod withconstant thermal properties with the following sources and boundary conditions:(a)Q= 0, u(0) = 0, u(L) =T.

6 (f)Q=K0x2, u(0) =T, u (L) = (a) Equilibrium satisfiesu (x) = 0,whose general solution isu=c1+ boundary conditionu(0) = 0 impliesc1= 0 andu(L) =Timpliesc2=T/Lso thatu=Tx/L.(f) In equilibrium,usatisfiesu (x) = Q/K0= x2,whose general solution (by integrating twice) isu= x4/12 +c1+ boundary conditionu(0) =Tyieldsc1=T,whileu (L) = 0 yieldsc2=L3/3. Thusu= x4/12 +L3x/3 + Suppose u t= 2u x2+x, u(x,0) =f(x), u x(0, t) = , u x(L, t) = 7.(a) Calculate the total thermal energy in the one-dimensional rod (as a function of time).(b) From part (a), determine a value of for which an equilibrium exists.

7 For this valueof , determine limt u(x, t).Solution. (a) Integrating the equation, we obtain:ddt L0u(x, t)dx= L0( 2u x2+x)dx= u x L0+12L2= 7 + intfrom 0 tot, we obtain the total thermal energy L0u(x, t)dx= L0f(x)dx+(7 +12L2)t.(8)(b) In order for an equilibrium to exist,(7 +12L2)tmust be 0. So = 7 + equilibrium satisfies (x) +x= general solution (after integrating twice) is = 16x3+c1+ boundary condition yieldsc2= 7 + = 16x3+c1+(7 +12L2) u(x, t) = (x),using (8), we obtain L0f(x)dx= L0u(x, t)dx= limt L0u(x, t)dx= L0 (x)dx= L0( 16x3+c1+(7 +12L2)x)dx= 124L4+c1L+12(7 +12L2) it givesc1= L0f(x)dx 72L2 524L4L,and then = 16x3+(7 +12L2)x+ L0f(x)

8 Dx 72L2 For conduction of thermal energy, the heat flux vector is = K0 inaddition the molecules move at an average velocityV, a process called convection, then = K0 u+c the corresponding equation for heat flow, including both conductionand convection of thermal energy (assuming constant thermal properties with no sources). quantities: Thermal energy densitye(x, t) = the amount of thermal energy per unit volume. Heat flux (x, t) = the amount of thermal energy flowing across boundaries per unitsurface area per unit Heat sourcesQ(x, t) = heat energy per unit volume generated per unit time.

9 Temperatureu(x, t). Specific heatc= the heat energy that must be supplied to a unit mass of a substanceto raise its temperature one unit. Mass density (x) = mass per unit of heat energy:Rate of change of heat energy in time = Heat energy flowing across boundaries per unittime + Heat energy generated insider per unit time heat energy = Re(x, t)dV. Heat energy flowing across boundaries per unit time = ndS. Heat energy generated insider per unit time = RQ(x, t)dV= t Re(x, t)dV= divergence theorem give t Re(x, t)dV= R then e t= .(9)Heat energy per unit volume =c(x)u(x, t).

10 Soe(x, t) =c(x)u(x, t) .It then follows thatc u t= (c uV) + (K0 u).(10)and then the convection-diffusion equation u t+ (uV) =k 2u,(11)wherek=K0c is called the thermal Assignment (d) Find the eigenvalues and the corresponding eigenfunctions of the eigenvalue prob-lem d2 dx2= , (0) = 0,d dx(L) = 0.(i) If >0, =c1cos( x) +c2sin( x). (0) = 0 impliesc1= 0, whiled dx(L) = 0impliesc2 cos( L) = 0. Thus L= 2+n (n= 1,2, ). Then the eigenvaluesare n=( 2+n )2/L2and the corresponding eigenfunctions are n= sin(( 2+n )xL)(n= 1,2, ).(ii) If = 0, =c1+c2x. (0) = 0 impliesc1= 0, whiled dx(L) = 0 impliesc2= 0.