From Wikiversity
Jump to: navigation, search

R 9.1 - Coefficients for heat conduction on a sphere[edit]







1. Find
2. Compute the first 5 non-zero coefficients
3. Plot Exact Function, Fourier-Legendre series with 3 Terms and 5 Terms


Problem 1.



Converting to











Problem 2.
Substitute j into Eq 9.1.11 to Eq 9.1.15

Problem 3.
Use Equation(9.1.10),Expand the fourier-legendre function


Egm6322.s12.sungsik 01:54, 7 February 2012 (UTC) Egm6322.s12.sungsik 02:05, 10 February 2012 (UTC)

R 9.2 - Coefficients for general solution of heat conduction on a sphere[edit]


Boundary condition on the surface of the sphere:


Fourier-Legendre series for the boundary condition in 9.2.1:


General solution for heat conduction on a sphere:



Show that



Solved on my own

The purpose of this problem is to show that the coefficients of the general formulation of the heat conduction problem on a sphere (obtained through solving the Laplace equation via separation of variables) are equal to the coefficients obtained from expanding the Fourier-Legendre series subject to the boundary condition of Eq. 9.2.1. As this condition is valid when r=1, Eq. 9.2.3 becomes:


Applying the boundary conditions and comparing to the Fourier-Legendre expansion of 9.2.2:


Recall from class notes that the only nonzero indices of are j=0,2,4. This is due to the orthogonality and oddness/eveness of the Legendre polynomials. This is covered in the lecture notes but will be reviewed here for clarity.

We know that the Legendre polynomials, by definition, are a subset of all polynomials:


Any set of polynomials of degree n are spanned by the same number of orthogonal polynomials (for our purposes, Legendre polynomials will suffice). This also indicates that a polynomial of degree m > n is orthogonal to the subspace as no linear combination of degree n polynomials can produce a new polynomial of a higher order. This is written as:


The boundary condition of the problem, given in Eq. 9.2.1, is a 4th order polynomial; thus, when computing the coefficients of the Fourier-Legendre expansion in Eq. 9.2.2 according to the formula


the numerator on right hand side will be zero for j>4 due to the orthogonality of the polynomials (the inner product of orthogonal functions is zero). Thus,


It should also be noted that the boundary condition is an even function whereas Legendre polynomials of order 2k+1 for k = 0,1,2,... are odd functions; thus, since the inner product of odd and even functions are zero,


From Eq. 9.2.10 and 9.2.11, the only surviving coefficients are . Thus, after evaluating these coefficients using Eq. 9.2.9, Eq. 9.2.6 becomes:


As the expansion of the left hand side of Eq. 9.2.12 is in terms of the same Legendre polynomials as the right hand side, it is clear that the coefficients must be identical. For , the values of are zero and for the coefficients must be the same as given in the class notes. Thus,


R*9.3 Identities of Legendre Equation[edit]


We are given the 1st and 2nd homogeneous solutions of the Legendre Equation from p.41-4, 41-5, 7-1, and 37-4.


Derive the following identities for the 2nd homogeneous solutions of the Legendre Equation.

Also find the expression for the in terms of the inverse hyperbolic tangent of x

From lecture notes Mtg 43b


Solved on our own

The Legendre differential equation was defined in lecture page 7-1 for EGM 6321.

King (2003) very clearly derives the solutions for these equations in pages 5-6 and 31-33. First assert the method of Frobenius. Choose a power series solution of the form below.

Substitute the term above into Legendre's equation to get a series representation.

Rearrange the equation to simplify the grouping of terms below.

Further simplification gives a form that allows derivation of homogeneous solutions.

The series definitions gives an indicial equation of the form c(c-1) = 0. Apply Frobenius General Rule 1 on p. 14 of King (2003) - "If the indicial quation has two distinct roots, whose difference is an integer, and one of the coefficients of becomes indeterminate on putting , both solutions can be generated by putting in the recurrence relation." Therefore, the solution of Legendre's equation is given in expanded form below.

The simplified form of the soltuion is given below. This equation gives rise to the Legendre polynomials.

In the solution, the Legendre polynomials are given by the following summation.

The first three terms are listed as .

These terms provide the basis for finding the Legendre functions: . The key is to utilize the reduction of order formula. For a homogeneous differential equation, we have the definition below.

If a solution exists, then a 2nd solution is given by this integral.

If we apply this equation to Legendre polynomials, then the following form is given for the differential equation.

The astute mathematician must note that the integral for the inner integral must be derived first in the reduction of order formula.

So the Legendre function can be derived by using substitution variables of t. In general, the equation below gives the nth Legendre function for the nth Legendre polynomial.

For n=0, the following derivation occurs for .

Similarly, for n=1, can be derived.

Lastly, for n=2, the integration gives .

Note that in each Legendre function there exists a term that can be simplified into the inverse hyperbolic tangent. Apply the definition from Wolfram for inverse hyperbolic tangent at

Therefore, the following definitions are derived for 9.3.

Note 4 tilde command by Manuel Steele below.

Egm6322.s12.team2.steele.m2 08:24, 10 February 2012 (UTC)

R 9.4 -Verify the second homogeneous solutions of Legendre equation[edit]


Given the following equations




Verify the expressions for as shown on the following equations





Solved on my own

The Legendre polynomial is as follows:


where .


Then from 9.4.1 we have

Although , the fact that 1-1=0 makes this term meaningless. We'd better define , if

R 9.5 - Legendre Function[edit]


Show that


Show that



is odd or even depending on the index n





When n=even
From Eq(9.5.1)

when n=Odd


Egm6322.s12.sungsik 02:22, 7 February 2012 (UTC)

R 9.6 - Equalities in spherical coordinates[edit]


Definition of : Boundary condition on the surface of the sphere:


Cosine law:



Show the following equalities:




Solved on my own

The first inequality in Eq. 9.6.3 can be shown trivially simply by substituting the definition of from Eq. 9.6.1 into the equation. This immediately yields the solution.


Showing the equality in Eq. 9.6.4 requires somewhat more effort. It is first necessary to define the coordinates of points P and Q in spherical coordinates, as shown in the class notes:


The distance d between two points in Cartesian space is given by the three-dimensional Pythagorean theorem:


Substituting the spherical coordinates in Eq. 9.6.6 into Eq. 9.6.7 yields an expression for :




Expanding and collecting like terms yields:


Again collecting like terms in the term on the far right yields:


Recognizing the term inside the parenthesis as the sum of cosines identity shown in Eq. 9.6.12:


Eq. 9.6.11 can now be written as:


Comparing this expression for with the one given by the law of cosines in Eq. 9.6.2 indicates that the term in the brackets must be equal to cos for the equality to hold true:


R9.7 Identities of Binomial Theorem[edit]


We are given the following Binomial Equation from 43-21.


Derive the following identities for the Binomial Theorem.

From lecture notes Mtg 43c


Solved on our own

The first step is to write the general equation of the Binomial Theorem.

Note that the Binomial Theorem can be extended to infinity when .

The biniomial coefficent can be expanded as follows.

Note that if we set y=1 then the equation can be simplified.

The Binomial Equation was adapted by Isaac Newton for real numbers beyond non-negative integers. The key is to factor (n-k)! from the binomial coefficient. Then if r=-s, we get the following form of the Binomial Equation.

Note that if s = 1/2, then the equation can be further simplified.

Realize that the binomial coefficient is a fraction, so the numerator and denominator can be multiplied by 2. This will give us the form for the solution.

Therefore, the given steps derive the desired form below where we let i=k.

Manuel Steele 4 tilde command signature

Egm6322.s12.team2.steele.m2 08:26, 10 February 2012 (UTC)

R* 9.8 -Generating Legendre polynomials[edit]


Part 1

Recurrence relation RR2:


with and

Part 2



and the binomial theorem:





in terms of condition 1, and 2 respectively. and compare the result


Solved on my own

Part 1

From 9.8.1, we have


with and then

Part 2


From the 9.8.3 and 9.8.4 we have

when computing for , we can see that the highest order is ,thus

Since the highest order we consider is 6,then


So, using two different methods, we obtained the same result.

Contributing Members[edit]

Team Contribution Table
Problem Number Assigned To Solved By Typed By Proofread By
9.1 sungsik sungsik sungsik Matt Shields
9.2 Matt Shields Matt Shields Matt Shields Lang Xia
9.3 Manuel Steele Manuel Steele Manuel Steele
9.4 Lang Xia Lang Xia Lang Xia Manuel Steele
9.5 sungsik sungsik sungsik Manuel Steele
9.6 Matt Shields Matt Shields Matt Shields Lang Xia
9.7 Manuel Steele Manuel Steele Manuel Steele
9.8 Lang Xia Lang Xia Lang Xia

Members Signatures[edit]

Egm6322.s12.team2.Xia, I solved R9.4 and R*9.8 without others assistance and proofread R9.2 and R9.6. 10:55, 10 February 2012 (UTC)
Egm6322.s12.sungsik 02:07, 10 February 2012 (UTC)

References for Report 9[edit]

King, A.C., J. Billingham and S.R. Otto. "Differential Equations: Linear, Nonlinear, Ordinary, Partial." New York, NY: Cambridge University Press, 2003.

Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 43a) Mtg 43a 24 Jan 2011.

Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 43a) Mtg 43b 24 Jan 2011.

Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 43c) Mtg 43c 31 Jan 2012.

Vu-Quoc, L. Class Lecture: Principles of Engineering Analysis. University of Florida, Gainesville, FL, (Meeting 43a) Mtg 43d 31 Jan 2011.