User:Egm6322.s14.team2/report3

=Report 3 PEA 2 Spring 2014 Team 2=

Problems are given in lecture notes

Problem Statement
sec.67a, Pb-67.1 (equilibrium of line defects, Legendre diff. eq., Legendre polynomials)

Provide a detailed justification for the following relations:

and thus,

$$ P(u) $$ satisfies the differential equation,

Find
Provide a detailed justification for the following: (a):

(b):

(c):

(d): $$ P(u) $$ satisfies the differential equation,

(e):

Solution

 * {| style="width:100%" border="0"

On our honor, we did this assignment on our own, without looking at the solutions in previous semesters or other online solutions.
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * }

Given at equilibrium,

($$) is the definition of S(u), a summation of N terms.

Substituting ($$)for the N-1 terms (excluding j=i) in ($$) and evaluating the behavior of the function as u approaches ui results in ($$).

P(u) is defined in equation ($$) and the derivative is displayed in ($$).

Substituting ($$) and ($$) into ($$), ($$) was established.

Using the chain rule the following equation can be defined using ($$).

($$) is the derivative of S as $$ u \to u_i $$.

Substituting ($$) and ($$) into ($$), multiplying the last term by (u-ui)/(u-ui), and evaluating the behavior of the system as $$ u \to u_i $$ yields P''(u)/P'(u) approaches 0 as $$ u \to u_i $$.

Evaluating the product in ($$), a power series with constant coefficients that are products of the ui positions can be expressed as follows:

And the first and second derivatives can also be expressed as a power series as follows:

Given that P''(u) is proportional to $$\frac{P(u)}{u^{2}-1}$$,

Substituting the power series expressions of the solutions into the above equation:

Considering just the final term in the summation:

The above equation simplifies to c=N(N-1). Substituting the coefficient, c, back into ($$)

Now, using the transformation for P as a function of Q, and the first and second derivatives:

Substituting the transformation into the differential equation in ($$)

Problem Statement
sec.67b, Pb-67.2 (equilibrium of line defects, Hermite diff. eq., Hermite polynomials)

Equation ($$) is given:

Find a suitable change of variable to show that ($$) ((6) p.68-5) can be written as the Hermite differential equation of the form:

For the 2 cases of $$ N = 7,8 $$ particles, find the non-dimensional equilibrium positions and plot them.

Find
(a): A suitable change of variable to show that ($$) can be written as the Hermite differential equation of the form of ($$)

Solution

 * {| style="width:100%" border="0"

On our honor, we did this assignment on our own, without looking at the solutions in previous semesters or other online solutions.
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * }

Part (a): Find Suitable Change of Variable
For this problem we seek to determine a suitable change of variables that will give us ($$) from ($$). The variable of interest is changing from u to x in this problem. We need to find some sort of relation between the two variables that will gives us what we want. We can observe the differences between the two forms of the differential equations and conclude that a factor of 2 on the zeroth and first derivative terms is an observable difference. We can also observe that the transformed differential equation is of the same order as the original, therefore the change of variables does not include a derivative term. It must therefore only contain the independent variables u and x.

Using the transformation $${\rho}(u)={\lambda}u$$ and the first and second derivatives below in the function P:

Make the substitution $$ u{\lambda}=\rho (u) $$ and then divide the equation by $${\lambda}^{2}$$ to get.

Examining equation ($$) in comparison to the the Hermite differential equation $$, $${\lambda}=\frac{1}{\sqrt{2}}$$ and $${\rho}(u)=\frac{u}{\sqrt{2}}$$. Using y instead of P(u) and x instead of $$ \rho (u)$$, the equation takes the form of the Hermite differential equation above.

Plot Solution for 2 Cases
Here we are asked to plot the solution to the Hermite differential equation for the cases when $$ N = 7,8 $$.

Here the differential equation takes the following form when $$ N=7 $$.

The equation for the coefficients of the Hermite polynomial can be obtained by using a power series solution to the even or odd Hermite differential equation:

Therefore, the 7th and 8th order Hermite polynomials, respectively, are:



The roots from Wolfram Alpha using the functions HermiteH[7,x] and HermiteH[8,x] are:

This equation is a second order linear differential equation that is homogenous with variable coefficients. A solution is sought in the form of a finite power series. It is finite in this case because the original polynomial that we are solving for is a finite polynomial i.e. it only contains terms up to $$ u^{N} $$

Wolfram Alpha gives the solution

Problem Statement
sec.70a, Pb-70.1 (Kepler’s equation, acceleration vector in planar motion using Euler formula)

Show,

Find
(a): ($$)

Solution

 * {| style="width:100%" border="0"

On our honor, we did this assignment on our own, without looking at the solutions in previous semesters or other online solutions.
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * }

Using polar coordinates, the unit vectors that form a basis in the orbital plane and their first order time derivatives are:

The position vector, r, in this basis is:

The time derivative of the position vector is the velocity:

The second derivative with respect to time is the acceleration:

After combining like terms and factoring a $$\frac{1}{r}$$ out of the $$e_{\theta}$$ term, the equation above simplifies to:

Recognizing the derivative of $$r^{2}\dot{\theta}$$ and performing the integration, Kepler's equation takes the form:

Problem Statement
A general First Order Linear Differential Equation with Varying Coefficients (L1-ODE-VC) has the following form:

Find
(a)

(b)

Solution

 * {| style="width:100%" border="0"

On our honor, we did this assignment on our own, without looking at the solutions in previous semesters or other online solutions.
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * }

Part a
The equation in 1 can be written as follows:

The general solution for the Euler Integrating Factor Method is:

And $$h(x)$$ is:

($$k_1$$ is arbitrary, it scales the integrating factor by a constant but has no effect on the final solution. Here we set it to 0 for ease of computation.)

Plugging $$a_0$$ into $$h(x)$$ and integrating gives:

Plugging $$h(x)$$ into the general solution we get:

Integrating this gives the solution:

Part b
The equation 2 can be rewritten as:

($$) is divided by $$(x^2+1)$$ in order to have the coefficient of the leading derivative be unity.

$$h(x)$$ is then found by the following equation:

($$) then becomes:

Which after evaluating the integral can be simplified to:

Problem Statement
Redo commented problems from Report 2.

Given Information
The following equation is given for ds2:

xi is written in spherical coordinates in ($$) through ($$).

Problem Statement
Show that the infinitesimal length ds2 in ($$) can be written in spherical coordinates, ($$), as given in Pb.R*7.3 on pp.39-[1,3] of sec.39 of the notes.

Solution

 * {| style="width:100%" border="0"

On our honor, we did this assignment on our own, without looking at the solutions in previous semesters or other online solutions.
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * }

The derivatives of xi are given in ($$) through ($$).

($$) was developed by substituting the squared derivatives in the summation of ($$) and combining like terms.

Using the Pythagorean Trigonometric Identity, ($$) reduces to ($$).

The magnitude of the tangent vector h are then identified as follows:

The definition of the Laplace operator, $${\triangle}u$$, in cartesian coordinates:

The Laplace operator can be expressed in spherical coordinates from equation (2) from the class notes 39-2,

After substituting the values of h from ($$) through ($$), the Laplacian takes the form:

Factoring the $$cos{\theta}$$ from the first term, the partial derivative with respect to r, and the $$sec{\theta}$$ from the third term, the partial derivative with respect to $${\varphi}$$, since both values are fixed in the partial derivatives, the Laplacian takes the form:

Given Information
The following equation is given for P2(x), the Legendre Polynomial of degree n=2:

Problem Statement
Using variation of parameters, show that $$ is a second solution to the Legendre differential equation in $$ as given in Pb.R*6.11 on pp.37-4 in sec.37 of the notes.

Solution

 * {| style="width:100%" border="0"

On our honor, we did this assignment on our own, without looking at the solutions in previous semesters or other online solutions.
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * }

The Legendre differential equation of degree n=2 is given in equation 5-3.

Using reduction of order, some function, v, times P2(x) and the first two derivatives of the product are substituted into the Legendre differential equation in equation 5-7. Furthermore, since P2(x) is a solution to the Legendre equation the coefficient of the first order term, v, reduces to zero. Since the lowest order term cancels out, the second order differential equation is reduced to a first order differential equation where w=v' and w'=v'' in Equation 5-8 which can be solved through separation of variables.

Through partial fraction decomposition, the integral is

Through integration

Using properties of the natural logarithm and then inverting it using the exponential function, equation 5-9 becomes

Substituting v' back into Equation 5-10 and again using partial fraction decomposition

Using integration by parts

where u and v are defined as follows:

Using the properties of the natural logarithm to combine the first two terms and performing the simple integration on the last term:

Substituting v into Equation 5-4 then identifying the new solution as Q2

R2.6: Testing Gauss Laguerre quadrature
The Bessel function integration test performed by Dr. Burkardt's code test the following integral:

$$ e^{-2} \int^{\infty}_{2} \frac{\sin{x-1}}{\sqrt{x (x-2)}} dx $$

and has the following raw output:

7                 0.162669                 1         0.19313       0.0304616                 2       0.0346675        0.128001                 4       0.0367188         0.12595                 8       0.0395037        0.123165                16       0.0970831       0.0655858                32        0.100708       0.0619605                64        0.107105       0.0555637

The top row cow is the problem number, and the exact solution. The first column of numbers is the number of integration points, the second column is the approximated answer, and the third column is the absolute error.

This code did not test a Bessel integral of the form given in problem R2.7.

Problem Statement
Test the following integrals using Gauss-Laguerre quadrature.

Where $$ \epsilon= 0 $$

Find
Plot the convergence of the integrals in:

(a) ($$)

(b) ($$)

as the number of integration points are increased.

Solution

 * {| style="width:100%" border="0"

On our honor, we did this assignment on our own, without looking at the solutions in previous semesters or other online solutions.
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * style="width:92%; padding:10px; border:2px solid #8888aa" |
 * }

Neither integral given converges, nor has a principle value for $$\epsilon=0$$. However, the exercise is carried out for a variety of decreasing values of $$\epsilon$$ ranging from $$\epsilon=1$$ to $$\epsilon=10^{-6}$$

The Gauss-Laguerre quadrature rule was given in Equation 7.7 Report2

Plots were generated using the following numpy script.

Source:

Part b


In both Figure 6.1 and Figure 6.2, the quadrature converge quickly to the correct value for larger values of $$\epsilon$$, but fail to converge to the correct answer for smaller $$\epsilon$$

=Contributing Team Members=


 * Cameron Stewart Solved problem 6, part b of problem 4, and the part of problem 5 pertaining to R2.7. Reviewed problem 3.
 * Elizabeth Bartlett Solved problems 1, 2, 3, and the portion of problem 5 related to R2.4 and R2.5. Reviewed problem 4.
 * Kevin Frost Solved problem 4 part a.
 * Christopher Neal Solved problem 2.

=References=