Session 5 Heat Conduction in Cylindrical and Spherical Coordinates I
|
|
- Jocelin Ellis
- 6 years ago
- Views:
Transcription
1 Session 5 Heat Conduction in Cylindrical and Spherical Coordinates I 1 Introduction The method of separation of variables is also useful in the determination of solutions to heat conduction problems in cylindrical and spherical coordinates. A few selected examples will be used for illustration. Recall that the simplest form of the heat equation in cylindrical coordinates (r, φ, z) is ρc p T t = 1 r r (kr T r )+ 1 r 2 And in spherical coordinates (r, φ, θ) ρc p T t = 1 r 2 r (kr2 T r )+ 1 r 2 sin θ φ (k T φ )+ z (k T z )+g T (k sin θ θ θ )+ 1 r 2 sin 2 θ φ (k T φ )+g 2 Fundamental Solutions to Steady State, Unidimensional Problems As in the case of Cartesian coordinates, analytical solutions are readily obtained for unidirectional problems in cylindrical and spherical coordinates. Solutions to steady unidimensional problems can be readily obtained by elementary methods as shown below. 2.1 Cylindrical Coordinates Consider the infinite hollow cylinder with inner and outer radii r 1 and r 2, respectively. At steady state the heat equation in cylindrical coordinates with azimuthal symmetry becomes d dt (r dr dr )=0 1
2 the general solution of which is T (r) =Aln r + B where the constants A and B must be determined from the specific boundary conditions involved. This represents the steady state loss of heat through a cylindrical wall. Two particular sets of boundary conditions are now investigated in detail. Case 1.- Assume the heat flux at r = r 1, q 1 is given (Neumann boundary condition), while the temperature at r = r 2 is specified as T 2 (Dirichlet boundary condition). Therefore, at r = r 1 one has q(r 1 )=q 1 = k dt dr r=r 1 Substitution into the general solution given above yields A = q 1r 1 k Now, from the boundary condition at r = r 2 one has The required solution is then B = T 2 + q 1r 1 k ln r 2 T (r) = q 1r 1 k ln(r 2 r )+T 2 Finally, the unknown temperature at r = r 1 is given by T 1 = q 1r 1 k ln(r 2 r 1 )+T 2 Note that since ln(r 2 /r 1 ) > 0, if q 1 > 0 (i.e. heat flows into the cylindrical shell through the inner radius) then T 1 >T 2. Case 2.- This case is identical to the first except that boundary conditions are reversed relative to r 1 and r 2, i.e. one has a Dirichlet boundary condition at r = r 1 (Known temperature T 1 ) and a Neumann boundary condition at r = r 2 (Known heat flux q 2 ). Proceeding as before one gets the solution and the unknown temperature at r = r 2 is T (r) = q 2r 2 k ln(r 1 r )+T 1 T 2 = q 2r 2 k ln(r 1 r 2 )+T 1 Since ln(r 1 /r 2 ) < 0, if q 2 < 0 is (i.e. heat flows into the cylindrical shell through the outer radius) then T 2 >T 1. 2
3 2.2 Spherical Coordinates Consider a hollow sphere with inner and outer radii r 1 and r 2, respectively. At steady state the heat equation in spherical coordinates with azimuthal and poloidal symmetry becomes the general solution of which is d dt (r2 dr dr )=0 T (r) = A r +B where the constants A and B must be determined from the specific boundary conditions involved. This represents the steady state loss of heat through a spherical shell. As in the previous section, two particular sets of boundary conditions are now investigated in detail. Case 1.- Assume the heat flux is given at r = r 1 as q 1 (Neumann boundary condition) while the temperature is specified as T 2 at r = r 2 (Dirichlet boundary condition). Therefore, at r = r 1 one has q(r 1 )=q 1 = k dt dr r=r 1 yielding A = q 1r1 2 k The second boundary condition then leads to so that the desired solution is B = T 2 q 1r 2 1 k and the unknown temperature at r = r 1 is 1 r 2 T (r) = q 1r 2 1 k (1 r 1 r 2 )+T 2 T 1 = q 1r 2 1 k ( 1 r 1 1 r 2 )+T 2 Note that the since (1/r 1 1/r 2 ) > 0, if q 1 > 0 (i.e. heat flows into the spherical shell through the inner radius) then T 1 >T 2,. Case 2.- Now assume instead that T (r 1 )=T 1 is given (Dirichlet) and q(r 2 )=q 2 is also known (Neumann). With these A = q 2r 2 2 k 3
4 and B = T 1 q 2r 2 2 k 1 r 1 so that the desired solution is and the unknown temperature at r = r 2 is T (r) = q 2r 2 2 k (1 r 1 r 1 )+T 1 T 2 = q 2r 2 2 k ( 1 r 2 1 r 1 )+T 1 Note also that since (1/r 2 1/r 1 ) < 0, if q 2 > 0 (i.e. heat flows into the spherical shell through the outer radius) then T 2 >T 1. 3 Steady State Multidimensional Problems Many multidimensional problems can be solved by simple extension of the separation of variables method. As an example consider the problem of steady state heat conduction in a short cylinder (radius r 0,heightL). The goal is to determine the steady state temperature field T (r, z) inside the cylinder subject to specific conditions on its boundaries. The governing equation in this case is 2 T r + 1 T 2 r r + 2 T z =0 2 Let the boundary conditions at r = r 0 and z = L be and at z =0, T(R, z) =T(r, L) =0 T(r, 0) = T 0 An additional but implicit requirement is that the solution remain bounded at r = 0. Now assuming the required solution is of the form T (r, z) =R(r)Z(z) and proceeding as in the Cartesian case, yields the following set of ODEs d dr (r dr dr )+λ2 rr =0 4
5 and 2 Z z 2 λ2 Z =0 where λ 2 is the separation constant representing the eigenvalues of the Sturm-Liouville problem for R, i.e. λ n are the roots of J 0 (λ n r 0 )=0 for n =1,2,... Solving the above, incorporating the boundedness consition as well as the first two boundary conditions and performing a linear combination of the eigensolutions gives T (r, z) = R n (r)z n (z) = a n J 0 (λ n r)sinh(λ n [L z]) Finally, substituting the boundary condition at z = 0 requires that T 0 = [a n sinh(λ n L)]J 0 (λ n r) which is the Fourier-Bessel series representation of T 0 and leads directly to the Fourier coefficients a n = 1 2 sinh(λ n L) r0j (λ n r 0 ) r0 0 T 0 J 0 (λ n r)rdr 4 Analytical Solutions to Transient Problems Analytical solutions to various transient problems in cylindrical and spherical coordinates can be obtained using the method of separation of variables. Several examples are presented below. 4.1 The Quenching Problem for a Cylinder with Fixed Temperature at its Boundary Consider the quenching problem where a long cylinder (radius r = b) initially at T = f(r) whose surface temperature is made equal to zero for t>0. The heat equation in this case is 1 r r (r T r )= 1 α T t 5
6 subject to T =0 at r = b and T r =0 at r =0. According to the separation of variables method we assume a solution of the form T (r, t) = R(r)Γ(t). Substituting this into the heat equation yields 1 d dr (r rr dr dr )= 1 dγ αγdt This equationonly makes sense when both terms are equal to a negative constant λ 2.The original problem has now been transformed into two boundary value problems, namely with general solution dγ dt + αλ2 Γ=0 Γ(t) =Cexp( αλ 2 t) where C is a constant and r 2 d2 R dr + r dr 2 dr + r2 λ 2 R =0 subject to R(b) = 0 and necessarily bounded at r =0. The last equationis a special case of Bessel s eqnarray*, the only physically meaningful solution of which has the form R(r) =A J 0 (λr) where A is another constant and J 0 (z) is the Bessel function of first kind of order zero of the argument given by J 0 (z) =1 z2 (1!) z4 2 (2!) 2 2 z6 4 (3!) Since R(b) = 0 this requires that J 0 (λb) = 0 which defines the eigenvalues and eigenfunctions for this problem. The eigenvalues are thus given as the roots of J 0 (λ n b)=0 6
7 Specifically, the first four roots are: λ 1 b =2.405, λ 2 b =5.520, λ 3 b =8.654, and λ 4 b = A particular solution is then T n (r, t) =R n (r)γ(t) =A n J 0 (λ n r) exp( αλ 2 nt) where A n = A nc and n =1,2,3,... Superposition of particular solutions yields the more general solution T (r, t) = T n (r, t) = R n (r)γ(t) = A n J 0 (λ n r) exp( αλ 2 nt) All is left is to determine the A n s. For this we use the given initial condition, i.e. T (r, 0) = f(r) = A n J 0 (λ n r) This is the Fourier-Bessel series representation of f(r) and one can use the orthogonality property of the eigenfunctions to write b b b rj 0 (λ m r)f(r)dr = A n rj 0 (λ m r)j 0 (λ n r)dr = A n rj0 2 (λ mr)dr = = b2 A m 2 [J 0 2 (λ mb)+j1 2 (λ mb)] = b2 A m 2 J 1 2 (λ mb) where J 1 (z) = dj 0 (z)/dz is the Bessel function of first kind of order one of the argument. Therefore, A n = so that the required solution is T (r, t) = 2 b 2 2 b rj b 2 J1 2 0 (λ n r)f(r)dr (λ n b) 0 J 0 (λ n r) b J1 2 (λ n b) exp( αλ2 nt) r J 0 (λ n r )f(r )dr 0 An important special case is obtained when f(r) =T i = constant. The required solution then becomes T (r, t) = 2T i b Exercise: Derive the above result. J 0 (λ n r) λ n J 1 (λ n b) exp( αλ2 nt) 7
8 4.2 The Quenching Problem for a Cylinder with Convective Heat Losses at its Boundary Many problems involving more complex boundary conditions can also be solved using the separation of variables method. As an example consider the quenching problem where a long cylinder (radius r = b) initially at T = f(r) is exposed to a cooling medium at zero temperature which extracts heat uniformly from its surface. The heat equation in this case is subject to at r = b, and subject to 1 r r (r T r )= 1 α k T r = ht T t T r =0 at r =0. Separation of variables T (r, t) =R(r)Γ(t) gives in this case T (r, t) = 2 b 2 m=1 β 2 exp( αβ m t) mj 0 (β m r) (βm 2 +( h k )2 )J0(β 2 m b) b 0 r J 0 (β m r )f(r )dr A special case of interest is when the initial temperature is constant = T i and the surrounding environment is at a non-zero temperature = T. In this case the above equationreduces to T (r, t) =T + 2 b (T i T ) m=1 1 J 1 (β m b)j 0 (β m r) β m J0 2 (β m b)+j1(β 2 m b) e β2 m αt where the eigenvalues β m are obtained from the roots of the following transcendental eqnarray* with Bi = hb/k. Exercise: Derive the above result. βbj 1 (βb) BiJ 0 (βb) =0 8
9 4.3 The Quenching Problem for a Sphere with Fixed Temperature at its Boundary Consider the quenching problem where a sphere (radius r = b) initially at T = f(r) whose surface temperature is made equal to zero for t>0. The heat equation in this case is subject to 1 r 2 r (r2 T r )= 1 T α t at r = b and T =0 T r =0 at r =0. According to the separation of variables method we assume a solution of the form T (r, t) = R(r)Γ(t). Substituting this into the heat equation yields 1 d dr r 2 (r2 R dr dr )= 1 dγ αγdt Again, this equationonly makes sense when both terms are equal to a negative constant λ 2. The original problem has now been transformed into two boundary values problems, namely with general solution dγ dt + λ2 αγ=0 where C is a constant and Γ(t) =Cexp( αλ 2 t) 1 d dr r 2 (r2 dr dr )+λ2 R=0 subject to R(b) = 0 and necessarily bounded at r = 0. The general solution of the last problem is R(r) =A sin(λr) r + B cos(λr) r = A sin(λr) r 9
10 where A and B are constants and B = 0 since the temperature must be bounded at r =0. Moreover, the boundary condition at r = b yields the eigenvalues and the eigenfunctions λ n = nπ b R n (r) = A n r sin(λ nr) = 1 r sin(nπr b ) for n =1,2,3,... A particular solution is then T n (r, t) =R n (r)γ(t) = A n r sin(λ nr) exp( αλ 2 nt) where A n = A nc and n =1,2,3,... Superposition of particular solutions yields the more general solution A n T (r, t) = T n (r, t) = R n (r)γ(t) = r sin(λ nr) exp( αλ 2 nt) To determine the A n s we use the given initial condition, i.e. T (r, 0) = f(r) = A n r sin(λ nr) this is again a Fourier series representation with the Fourier coefficients given by A n = 2 b b 0 f(r )sin( nπr b )r dr An important special case results when the initial temperature f(r) =T i = constant. In this case the Fourier coefficients become A n = T ib nπ cos(nπ) = T ib nπ ( 1)n Exercise: Derive the above result. 4.4 The Quenching Problem for a Sphere with Convective Heat Losses at its Boundary Consider a sphere with initial temperature T (r, 0) = F (r) and dissipating heat by convection into a medium at zero temperature at its surface r = b. The heat conduction equation in 1D spherical coordinates is 2 T r + 2 T 2 r r = 1 T α t 10
11 In terms of the new variable U(r, t) = rt(r, t) the mathematical formulation of the problem is subject to 2 U r 2 = 1 U α t U(0,t)=0, r =0 and U r +(h k 1 )U =0, b r = b U(r, 0) = rf(r), t =0 This is just like heat 1D conduction from a slab so the solution is T (r, t) = 2 β e αβ2 mt m 2 +(h/k 1/b)2 rm=1 b(βm 2 +(h/k 1/b) 2 )+(h/k 1/b) sin(β mr) b r F (r )sin(β m r )dr and the eigenvalues β m are the positive roots of r =0 β m b cot(β m b)+b(h/k 1/b) =0 Consider now as another example a hollow sphere a r b intially at F (r) and dissipating heat by convection at both its surfaces via heat transfer coefficients h 1 and h 2. Introducing the transformation x = r a, the problem becomes identical to 1D heat conduction from a slab and the solution is where T (r, t) = 1 r m=1 e αβ2 m t 1 b N(β m ) R(β m,r) r F(r )R(β m,r )dr r =a R(β m,r)=β m cos(β m [r a]) + (h 1 /k +1/a)sin(β m [r a]) and the eigenvalues β m are the positive roots of tan(β m [b a]) = β m([h 1 /k +1/a]+[h 2 /k 1/b]) β 2 m [h 1/k +1/a][h 2 /k 1/b] 11
12 5 Non-homogeneous Problems Nonhomogeneous problems for the cylinder or sphere can be solved using the same approach used to solve similar problems in Cartesian coordinates. Consider a long cylinder initially at T = F (r) inside which heat is generated at a constant rate g 0 and whose boundary r = b is subjected to T = 0. The mathematical statement of the problem consists of the heat eqnarray* 2 T r + 1 T 2 r r + 1 k g 0 = 1 T α t with T =0atr=band t>0andt =F(r)att= 0. The problem can be split into a steady state problem giving T s (r) =(g 0 /4k)(b 2 r 2 ) and a homogeneous transient problem giving T h (r, t). The solution to the original problem becomes Exercise: Solve the above problem. T (r, t) =T s (r)+t h (r, t). 6 Transient Temperature Nomographs: Heisler Charts The solutions obtained for 1D nonhomogeoeus problems with Neumann boundary conditions in cylindrical and spherical coordinate systems using the method of separation of variables have been collected and assembled in the form of transient temperature nomographs by Heisler. As in the Cartesian case, the charts are a useful baseline against which to validate one s own analytical or numerical computations. The Heisler charts summarize the solutions to the following three important problems. The first problem is the 1D transient homogeneous heat conduction in a solid cylinder of radius b from an initial temperature T i and with one boundary insulated and the other subjected to a convective heat flux condition into a surrounding environment at T. Introduction of the following nondimensional parameters simplifies the mathematical formulation of the problem. First is the dimensionless distance R = r b next, the dimensionless time then the dimensionless temperature τ = αt b 2 θ(x, τ) = T(r, t) T T i T 12
13 and finally, the Biot number Bi = hb k With the new variables, the mathematical formulation of the heat conduction problem becomes subject to at R =0, at R =1,and 1 θ (R R R R )= θ τ θ R =0 θ + Biθ =0 R θ=1 in 0 R 1 for τ =0. As the second problem consider the cooling of a sphere (0 r b) initially at a uniform temperature T i and subjected to a uniform convective heat flux at its surface into a medium at T with heat transfer coefficient h. In terms of the dimensionless quantities Bi = hb/k, τ = αt/b 2, θ =(T(r, t) T )/(T i T )andr=r/b, the mathematical statement of the problem is 1 θ R 2 (R2 R R )= θ τ in 0 <R<1, subject to θ R =0 at R =0,and θ + Biθ =0 R at R =1,and θ=1 in 0 R 1, for τ =0. As before, the solutions to the above problems are well known and are readily available in graphical form (Heisler charts). 13
14 7 Solution of Transient Multidimensional Problems by Product Superposition As in the Cartesian coordinates case, the solutions obtained for unidimensional systems can often be combined by product superposition in order to obtain solutions to multidimensional problems. Specifically, the solution to the problem of unsteady state conduction in a finite cylinder (radius r 0, height 2L) which is initially at temperature T i and is subsequently exposed to a quenching environment at temperature T, can be readily written down in dimensionless form as follows [ T (r, z, t) T T i T ] finitec ylinder =[ T(r, z, t) T T i T ] 2Ls lab [ T(r, z, t) T ] infinitec ylinder T i T 14
3.3 Unsteady State Heat Conduction
3.3 Unsteady State Heat Conduction For many applications, it is necessary to consider the variation of temperature with time. In this case, the energy equation for classical heat conduction, eq. (3.8),
More informationConduction Heat Transfer. Fourier Law of Heat Conduction. x=l Q x+ Dx. insulated x+ Dx. x x. x=0 Q x A
Conduction Heat Transfer Reading Problems 10-1 10-6 10-20, 10-48, 10-59, 10-70, 10-75, 10-92 10-117, 10-123, 10-151, 10-156, 10-162 11-1 11-2 11-14, 11-20, 11-36, 11-41, 11-46, 11-53, 11-104 Fourier Law
More informationChapter 4: Transient Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University
Chapter 4: Transient Heat Conduction Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Objectives When you finish studying this chapter, you should be able to: Assess when the spatial
More informationConduction Heat Transfer. Fourier Law of Heat Conduction. Thermal Resistance Networks. Resistances in Series. x=l Q x+ Dx. insulated x+ Dx.
Conduction Heat Transfer Reading Problems 17-1 17-6 17-35, 17-57, 17-68, 17-81, 17-88, 17-110 18-1 18-2 18-14, 18-20, 18-34, 18-52, 18-80, 18-104 Fourier Law of Heat Conduction insulated x+ Dx x=l Q x+
More informationChapter 5 Time-Dependent Conduction
Chapter 5 Time-Dependent Conduction 5.1 The Lumped Capacitance Method This method assumes spatially uniform solid temperature at any instant during the transient process. It is valid if the temperature
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationReview: Conduction. Breaking News
CH EN 3453 Heat Transfer Review: Conduction Breaking News No more homework (yay!) Final project reports due today by 8:00 PM Email PDF version to report@chen3453.com Review grading rubric on Project page
More informationHeat Transfer in a Slab
Heat Transfer in a Slab Consider a large planar solid whose thickness (y-direction) is L. What is the temperature history of the slab if it is suddenly brought into contact with a fluid at temperature
More informationTRANSIENT HEAT CONDUCTION
TRANSIENT HEAT CONDUCTION Many heat conduction problems encountered in engineering applications involve time as in independent variable. This is transient or Unsteady State Heat Conduction. The goal of
More informationMECH 375, Heat Transfer Handout #5: Unsteady Conduction
1 MECH 375, Heat Transfer Handout #5: Unsteady Conduction Amir Maleki, Fall 2018 2 T H I S PA P E R P R O P O S E D A C A N C E R T R E AT M E N T T H AT U S E S N A N O PA R T I - C L E S W I T H T U
More informationIntroduction to Heat and Mass Transfer. Week 8
Introduction to Heat and Mass Transfer Week 8 Next Topic Transient Conduction» Analytical Method Plane Wall Radial Systems Semi-infinite Solid Multidimensional Effects Analytical Method Lumped system analysis
More informationUnsteady State Heat Conduction in a Bounded Solid How Does a Solid Sphere Cool?
Unstead State Heat Conduction in a Bounded Solid How Does a Solid Sphere Cool? We examined the cooling a sphere of radius R. Initiall the sphere is at a uniform temperature T 0. It is cooled b convection
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationPHY752, Fall 2016, Assigned Problems
PHY752, Fall 26, Assigned Problems For clarification or to point out a typo (or worse! please send email to curtright@miami.edu [] Find the URL for the course webpage and email it to curtright@miami.edu
More informationChapter 4 TRANSIENT HEAT CONDUCTION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 4 TRANSIENT HEAT CONDUCTION LUMPED SYSTEM ANALYSIS Interior temperature of
More informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationChapter 10: Steady Heat Conduction
Chapter 0: Steady Heat Conduction In thermodynamics, we considered the amount of heat transfer as a system undergoes a process from one equilibrium state to another hermodynamics gives no indication of
More informationNumerical Heat and Mass Transfer
Master degree in Mechanical Engineering Numerical Heat and Mass Transfer 02-Transient Conduction Fausto Arpino f.arpino@unicas.it Outline Introduction Conduction ü Heat conduction equation ü Boundary conditions
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More information1 Introduction. Green s function notes 2018
Green s function notes 8 Introduction Back in the "formal" notes, we derived the potential in terms of the Green s function. Dirichlet problem: Equation (7) in "formal" notes is Φ () Z ( ) ( ) 3 Z Φ (
More informationIntroduction to Heat and Mass Transfer. Week 8
Introduction to Heat and Mass Transfer Week 8 Next Topic Transient Conduction» Analytical Method Plane Wall Radial Systems Semi-infinite Solid Multidimensional Effects Analytical Method Lumped system analysis
More informationChapter 22 Gauss s Law. Copyright 2009 Pearson Education, Inc.
Chapter 22 Gauss s Law 22-1 Electric Flux Electric flux: Electric flux through an area is proportional to the total number of field lines crossing the area. 22-1 Electric Flux Example 22-1: Electric flux.
More informationHEAT CONDUCTION USING GREEN S FUNCTIONS
HEAT CONDUCTION USING GREEN S FUNCTIONS Preface to the first edition Preface to the second edition Author Biographies Nomenclature TABLE OF CONTENTS FOR SECOND EDITION December 2009 Page viii x xii xiii
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationq x = k T 1 T 2 Q = k T 1 T / 12
Conductive oss through a Window Pane q T T 1 Examine the simple one-dimensional conduction problem as heat flow through a windowpane. The window glass thickness,, is 1/8 in. If this is the only window
More informationUniversity of Rome Tor Vergata
University of Rome Tor Vergata Faculty of Engineering Department of Industrial Engineering THERMODYNAMIC AND HEAT TRANSFER HEAT TRANSFER dr. G. Bovesecchi gianluigi.bovesecchi@gmail.com 06-7259-727 (7249)
More informationIntroduction to Heat and Mass Transfer. Week 7
Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms
More informationBOUNDARY PROBLEMS IN HIGHER DIMENSIONS. kt = X X = λ, and the series solutions have the form (for λ n 0):
BOUNDARY PROBLEMS IN HIGHER DIMENSIONS Time-space separation To solve the wave and diffusion equations u tt = c u or u t = k u in a bounded domain D with one of the three classical BCs (Dirichlet, Neumann,
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More informationChapter 2: Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering, Hashemite University
Chapter : Heat Conduction Equation Dr Ali Jawarneh Department of Mechanical Engineering, Hashemite University Objectives When you finish studying this chapter, you should be able to: Understand multidimensionality
More informationQuestions Chapter 23 Gauss' Law
Questions Chapter 23 Gauss' Law 23-1 What is Physics? 23-2 Flux 23-3 Flux of an Electric Field 23-4 Gauss' Law 23-5 Gauss' Law and Coulomb's Law 23-6 A Charged Isolated Conductor 23-7 Applying Gauss' Law:
More informationChapter 3: Transient Heat Conduction
3-1 Lumped System Analysis 3- Nondimensional Heat Conduction Equation 3-3 Transient Heat Conduction in Semi-Infinite Solids 3-4 Periodic Heating Y.C. Shih Spring 009 3-1 Lumped System Analysis (1) In heat
More informationConsider a volume Ω enclosing a mass M and bounded by a surface δω. d dt. q n ds. The Work done by the body on the surroundings is
The Energy Balance Consider a volume enclosing a mass M and bounded by a surface δ. δ At a point x, the density is ρ, the local velocity is v, and the local Energy density is U. U v The rate of change
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep Copyright
More informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationSteady and unsteady diffusion
Chapter 5 Steady and unsteady diffusion In this chapter, we solve the diffusion and forced convection equations, in which it is necessary to evaluate the temperature or concentration fields when the velocity
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationE. not enough information given to decide
Q22.1 A spherical Gaussian surface (#1) encloses and is centered on a point charge +q. A second spherical Gaussian surface (#2) of the same size also encloses the charge but is not centered on it. Compared
More informationFIND: (a) Sketch temperature distribution, T(x,t), (b) Sketch the heat flux at the outer surface, q L,t as a function of time.
PROBLEM 5.1 NOWN: Electrical heater attached to backside of plate while front surface is exposed to convection process (T,h); initially plate is at a uniform temperature of the ambient air and suddenly
More informationQUESTION ANSWER. . e. Fourier number:
QUESTION 1. (0 pts) The Lumped Capacitance Method (a) List and describe the implications of the two major assumptions of the lumped capacitance method. (6 pts) (b) Define the Biot number by equations and
More informationwhere u 0 is given. In all our applications it is possible to subtract a known function v(x, t) from u(x, t) so that w(x, t) = u(x, t) v(x, t)
CHAPTER PDE Partial Differential Equations in Two Independent Variables D. An Overview Drawing on the Sturm-Liouville eigenvalue theory and the approximation of functions we are now ready to develop the
More informationThe temperature of a body, in general, varies with time as well
cen58933_ch04.qd 9/10/2002 9:12 AM Page 209 TRANSIENT HEAT CONDUCTION CHAPTER 4 The temperature of a body, in general, varies with time as well as position. In rectangular coordinates, this variation is
More information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More informationPrinciples of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points)
Principles of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points) 1. Are the following statements true or false? (20 points) a. Thermal conductivity of a substance is a measure
More informationConsider a volume Ω enclosing a mass M and bounded by a surface δω. d dt. q n ds. The Work done by the body on the surroundings is.
The Energy Balance Consider a volume Ω enclosing a mass M and bounded by a surface δω. δω At a point x, the density is ρ, the local velocity is v, and the local Energy density is U. U v The rate of change
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep
More informationChapter 2: Steady Heat Conduction
2-1 General Relation for Fourier s Law of Heat Conduction 2-2 Heat Conduction Equation 2-3 Boundary Conditions and Initial Conditions 2-4 Variable Thermal Conductivity 2-5 Steady Heat Conduction in Plane
More informationChapter 22 Gauss s Law. Copyright 2009 Pearson Education, Inc.
Chapter 22 Gauss s Law Electric Flux Gauss s Law Units of Chapter 22 Applications of Gauss s Law Experimental Basis of Gauss s and Coulomb s Laws 22-1 Electric Flux Electric flux: Electric flux through
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationChapter 21: Gauss law Tuesday September 13 th. Gauss law and conductors Electrostatic potential energy (more likely on Thu.)
Chapter 21: Gauss law Tuesday September 13 th LABS START THIS WEEK Quick review of Gauss law The flux of a vector field The shell theorem Gauss law for other symmetries A uniformly charged sheet A uniformly
More informationTime-Dependent Conduction :
Time-Dependent Conduction : The Lumped Capacitance Method Chapter Five Sections 5.1 thru 5.3 Transient Conduction A heat transfer process for which the temperature varies with time, as well as location
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
More informationMa 530. Partial Differential Equations - Separation of Variables in Multi-Dimensions
Ma 530 Partial Differential Equations - Separation of ariables in Multi-Dimensions Temperature in an Infinite Cylinder Consider an infinitely long, solid, circular cylinder of radius c with its axis coinciding
More informationHomework 7 Solutions
Homework 7 Solutions # (Section.4: The following functions are defined on an interval of length. Sketch the even and odd etensions of each function over the interval [, ]. (a f( =, f ( Even etension of
More informationLecture 8 Multiple Choice Questions :
Lecture 8 Multiple Choice Questions : 1. A point charge -3Q lies at the centre of a conducting shell of radius 2R. The net charge on the outer surface of the shell is -3Q Zero +1.5 Q d. +3Q 2. Two identical
More informationProcess of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion
Brian England Objective Process of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion Cylindrical Coordinates Initial and Boundary Conditions Methodologies and
More informationSolutions to PS 2 Physics 201
Solutions to PS Physics 1 1. ke dq E = i (1) r = i = i k eλ = i k eλ = i k eλ k e λ xdx () (x x) (x x )dx (x x ) + x dx () (x x ) x ln + x x + x x (4) x + x ln + x (5) x + x To find the field for x, we
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More informationUnsteady State Heat Conduction in a Bounded Solid
Unsteady State Heat Conduction in a Bounded Solid Consider a sphere of radius R. Initially the sphere is at a uniform temperature T. It is cooled by convection to an air stream at temperature T a. What
More informationOne dimensional steady state diffusion, with and without source. Effective transfer coefficients
One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport
More informationPhysics 114 Exam 1 Fall 2016
Physics 114 Exam 1 Fall 2016 Name: For grading purposes (do not write here): Question 1. 1. 2. 2. 3. 3. Problem Answer each of the following questions and each of the problems. Points for each question
More information2 Which of the following represents the electric field due to an infinite charged sheet with a uniform charge distribution σ.
Slide 1 / 21 1 closed surface, in the shape of a cylinder of radius R and Length L, is placed in a region with a constant electric field of magnitude. The total electric flux through the cylindrical surface
More informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More informationPhysics 6303 Lecture 8 September 25, 2017
Physics 6303 Lecture 8 September 25, 2017 LAST TIME: Finished tensors, vectors, and matrices At the beginning of the course, I wrote several partial differential equations (PDEs) that are used in many
More informationDUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs
DUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs MATH 4354 Fall 2005 December 5, 2005 1 Duhamel s Principle
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationPDE and Boundary-Value Problems Winter Term 2014/2015
PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 6 Saarland University 17. November 2014 c Daria Apushkinskaya (UdS) PDE and BVP lecture 6 17. November 2014 1 / 40 Purpose of Lesson To show
More informationUNIVERSITY OF SOUTH CAROLINA
UNIVERSITY OF SOUTH CAROLINA ECHE 460 CHEMICAL ENGINEERING LABORATORY I Heat Transfer Analysis in Solids Prepared by: M. L. Price, and Professors M. A. Matthews and T. Papathansiou Department of Chemical
More informationPARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS
PARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS NAKHLE H. ASMAR University of Missouri PRENTICE HALL, Upper Saddle River, New Jersey 07458 Contents Preface vii A Preview of Applications and
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More information4. Analysis of heat conduction
4. Analysis of heat conduction John Richard Thome 11 mars 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Conduction 11 mars 2008 1 / 47 4.1 The well-posed problem Before we go further with
More information1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used?
1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used?. During unsteady state heat transfer, can the temperature
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationSAMPLE FINAL EXAM SOLUTIONS
LAST (family) NAME: FIRST (given) NAME: ID # : MATHEMATICS 3FF3 McMaster University Final Examination Day Class Duration of Examination: 3 hours Dr. J.-P. Gabardo THIS EXAMINATION PAPER INCLUDES 22 PAGES
More informationPhysics 116C Solutions to the Practice Final Exam Fall The method of Frobenius fails because it works only for Fuchsian equations, of the type
Physics 6C Solutions to the Practice Final Exam Fall 2 Consider the differential equation x y +x(x+)y y = () (a) Explain why the method of Frobenius method fails for this problem. The method of Frobenius
More informationGauss s Law. 3.1 Quiz. Conference 3. Physics 102 Conference 3. Physics 102 General Physics II. Monday, February 10th, Problem 3.
Physics 102 Conference 3 Gauss s Law Conference 3 Physics 102 General Physics II Monday, February 10th, 2014 3.1 Quiz Problem 3.1 A spherical shell of radius R has charge Q spread uniformly over its surface.
More information# Points Score Total 100
Name: PennID: Math 241 Make-Up Final Exam January 19, 2016 Instructions: Turn off and put away your cell phone. Please write your Name and PennID on the top of this page. Please sign and date the pledge
More informationChapter 2: Heat Conduction Equation
-1 General Relation for Fourier s Law of Heat Conduction - Heat Conduction Equation -3 Boundary Conditions and Initial Conditions -1 General Relation for Fourier s Law of Heat Conduction (1) The rate of
More informationChapter 3: Steady Heat Conduction
3-1 Steady Heat Conduction in Plane Walls 3-2 Thermal Resistance 3-3 Steady Heat Conduction in Cylinders 3-4 Steady Heat Conduction in Spherical Shell 3-5 Steady Heat Conduction with Energy Generation
More informationChapter 3 Steady-State, ne- mens onal C on uction
Chapter 3 Steady-State, One-Dimensional i Conduction 3.1 The Plane Wall 3.1.1 Temperature Distribution For one-dimensional, steady-state conduction in a plane wall with no heat generation, the differential
More informationd Wave Equation. Rectangular membrane.
1 ecture1 1.1 2-d Wave Equation. Rectangular membrane. The first problem is for the wave equation on a rectangular domain. You can interpret this as a problem for determining the displacement of a flexible
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationGauss s Law & Potential
Gauss s Law & Potential Lecture 7: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Flux of an Electric Field : In this lecture we introduce Gauss s law which happens to
More information1 Conduction Heat Transfer
Eng690 - Formula Sheet 2 Conduction Heat Transfer. Cartesian Co-ordinates q x xa x A x dt dx R th A 2 T x 2 + 2 T y 2 + 2 T z 2 + q T T x) plane wall of thicness 2, x 0 at centerline, T s, at x, T s,2
More informationTransient Heat Conduction in a Circular Cylinder
Transient Heat Conduction in a Circular Cylinder The purely radial 2-D heat equation will be solved in cylindrical coordinates using variation of parameters. Assuming radial symmetry the solution is represented
More informationFOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS
fc FOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS Second Edition J. RAY HANNA Professor Emeritus University of Wyoming Laramie, Wyoming JOHN H. ROWLAND Department of Mathematics and Department
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6
Lectures on Nuclear Power Safety Lecture No 6 Title: Introduction to Thermal-Hydraulic Analysis of Nuclear Reactor Cores Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More informationGauss s Law. Chapter 22. Electric Flux Gauss s Law: Definition. Applications of Gauss s Law
Electric Flux Gauss s Law: Definition Chapter 22 Gauss s Law Applications of Gauss s Law Uniform Charged Sphere Infinite Line of Charge Infinite Sheet of Charge Two infinite sheets of charge Phys 2435:
More informationThe Laplacian in Polar Coordinates
The Laplacian in Polar Coordinates R. C. Trinity University Partial Differential Equations March 17, 15 To solve boundary value problems on circular regions, it is convenient to switch from rectangular
More information13 Spherical Coordinates
Utah State University DigitalCommons@USU Foundations of Wave Phenomena Library Digital Monographs 8-204 3 Spherical Coordinates Charles G. Torre Department of Physics, Utah State University, Charles.Torre@usu.edu
More informationLecture 11: Fourier Cosine Series
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce Not to be copied, used, or revised without eplicit written permission from the copyright owner ecture : Fourier Cosine Series
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationName: Math Homework Set # 5. March 12, 2010
Name: Math 4567. Homework Set # 5 March 12, 2010 Chapter 3 (page 79, problems 1,2), (page 82, problems 1,2), (page 86, problems 2,3), Chapter 4 (page 93, problems 2,3), (page 98, problems 1,2), (page 102,
More informationLast Time. Magnetic Torque Magnetic Dipole in a B-Field: Potential Energy
Last Time Magnetic Torque Magnetic Dipole in a B-Field: Potential Energy 1 Today Magnetic Force Motors and Generators Gauss' Law Force on a Magnetic Dipole df df Idl B IBdl df IBdl sin df df If θ =, i.e.
More informationIf there is convective heat transfer from outer surface to fluid maintained at T W.
Heat Transfer 1. What are the different modes of heat transfer? Explain with examples. 2. State Fourier s Law of heat conduction? Write some of their applications. 3. State the effect of variation of temperature
More information