From Wikiversity
Jump to: navigation, search

My second wave of comments is color coded in green. Egm6322.s09 22:27, 2 April 2009 (UTC)

See my comments below.

After you made a correction for a section with a comment box, you want to put a comment in that same comment box on what you did. Egm6322.s09 13:49, 15 March 2009 (UTC)

The Classification of Second Order Partial Differential Equations[edit]


In the general form of PDE:

it can be classified by the value of .

PDE is hyperbolic.

PDE is parabolic.

PDE is elliptic.

For example:

In the equation:


This PDE is elliptic

Another way to find the type of equation (1):

preferred form

The equation (1) is an elliptic PDE.

Relationship Between Classifications and Transformations[edit]

Observations from the verification of PDE classifications:

1. Diffusion operator remains elliptic in polar coordinates.

Question:How about a different transformation of coordinate? Would classification remain the same?

2. Does classification make sense if it changes under transformation of coordinate?

The answer is no, because physics (e.g. distribution of temperature as a result of solution of heat equation) must remain the same regardless of how heat equation was solved (under different coordinate system).

Therefore, classification better remains the same under different coordinate system for it to make sense.

Egm6322.s09.three.liu 16:35, 24 April 2009 (UTC)

The Laplace Equation[edit] 17:58, 24 April 2009 (UTC)

The Laplace equation is the heat conduction equation with constant thermal conductivity and no heat generation.

The Laplace equation, symbolically:

Taking the Laplace equation in polar coordinates gives the following.

Features of an Axisymmetric Problem[edit] 17:58, 24 April 2009 (UTC)

Axisymmetric domain that includes the origin; 17:58, 24 April 2009 (UTC)
Axisymmetric domain that does not include the origin; 17:58, 24 April 2009 (UTC)

For an axis-symmetric problem, the theta terms drop out.

The equation then reduces to,

which is an ordinary differential equation (ODE).

In order to solve this ODE one should note that it may be rearranged.

Separating variables and integrating gives the solution:

where Ao and Bo are constants.

It should be noted that if the domain includes the origin (where r=0) then Ao must be zero for a finite solution.

Separation of Variables[edit] 17:58, 24 April 2009 (UTC)

Non-axisymmetric domain that includes the origin; 17:58, 24 April 2009 (UTC)
Non-axisymmetric domain that does not include the origin; 17:58, 24 April 2009 (UTC)

If the problem is not axis-symmetric, the Laplace equation may be solved using separation of variables.

Multiplying the Laplace equation by r2 gives:

... (a)

Observing that the equation has a portion that depends on r only and a part that depends on theta only, one may assume a solution of the form:

Thus the solution is a product of 2 functions: one which depends only on r [] and the other which only depends on theta [].

Plugging in this solution into (a) produces:

Dividing by and rearranging gives:

Move the figure to the left (instead of displaying it on the right) so to avoid blocking the hide/show link to open up the collapsible box. It is not possible to open up this collapsible box. Egm6322.s09 22:38, 2 April 2009 (UTC)

Figures moved to left. 22:38, 2 April 2009 (UTC)

And results in the following solution for .

If n=0:

In general, we want the solution to be periodic, such that:

The general form of the Laplace equation in polar coordinates takes the form that follows.

Another Axisymmetric Problem[edit]

What is "LP"? It is not clear how you arrived at the result, which is also wrong. Need more explicit explanation of the derivation, i.e., provide intermediate steps. I also mentioned to use our notation, not the notation by "LP". Egm6322.s09 15:08, 15 March 2009 (UTC)

My comment above had not been addressed; I did go over the above comment in class. Please take action. Egm6322.s09 22:38, 2 April 2009 (UTC)

Comment was addressed, but content was deleted. It has been re-posted, now with the correct result. 21:02, 6 April 2009 (UTC)

page 14.2


Pde hw4 1.svg


where T0 is a constant

Principle of Superposition[edit]

Problem P:


General Solution: Equation 2 p. 20-4.


Prob P1:

such that:

Prob P2:

such that:

For the Problem

Final Solution:

In general, for arbitrary function but periodic


for all and any K = constant

Pde hw4 2.svg

this is not periodic (not acceptable)

where: rn=an

Fourier Coefficients[edit]

Egm6322.s09.Three.nav 13:39, 24 April 2009 (UTC)

How do we derive the fourier coeffcients C0, An and Bn?

An orthogonal basis: {1, cosm, sinm} is used.

As defined above, two functions f(), g() are orthogonal if

Eg. Consider f()= cos () and g()= cos()

For more details, here is a link that explains the math in greater detail[1]. Similarly using '1' as the basis function and simplifying, we get

Using this knowledge, we proceed to determine the fourier coefficients C0, An and Bn, using the orthogonal basis {1, cosm, sinm}

Consider the equation -----(1)

Say we have the Boundary condition T(r=a, )= ()

Then, ----(2)

At n=0,

Multiplying through by {1} and integrating over [0,2]

Solving for An in Equation (2), multiply through by cos(m) and integrate over [0,2].

Evaluating each term in the resultant equation,

Term on the LHS= \int_{0}^{2\pi}T^{*}\left(\theta \right)cosm\theta d\theta

First term on RHS=

= 0, by definition of orthogonality

Second term on RHS=

(by definition of orthogonality)

Hence it is seen that though it was assumed m , simplification of the second term in (2) determines that m=n.

Third term on RHS=

= 0, by definition of orthogonality

(2) multiplied through by cos() or cos() and integrated over [0,2] reduces to

We solve similarly for Bn in (2). Multipling through with sin() and integrating over [0,2], it is seen that the first and second terms = 0, by definition. Third term reduces to . Then the modified Equation(2) becomes

Find the Fourier coefficients C0,An,Bn

Due to orthogonality of Fourier basis function

But you were asked to derive the above Fourier coefficients (and therefore the orthogonality property of the Fourier basis functions). Egm6322.s09 20:11, 15 March 2009 (UTC)

Edits to this section were made my Navya on March 27th. Either due to the work of a vandal, or a failure of wikiversity, they were removed.

They are now being placed in by Andrew Lapetina --Egm6322.s09.lapetina 20:56, 27 March 2009 (UTC)

A Third Example of the Laplace Equation[edit]

Egm6322.s09.bit.sahin 16:31, 24 April 2009 (UTC)

A domain which is a quadrant of annulus is subjected a boundary conditions that

temperature at is

temperature at is

Domain and boundary conditions

where and are given constants. Also, the boundaries at and are insulated which means no heat flow at these boundaries. According to the Fourier's law:

here denotes heat flux tensor. Relevant to the insulated conditions,

on .

How so? In general, express in polar coordinates then deduce . Such approach is important when the insulated boundaries do not coincide with the Egm6322.s09 20:11, 15 March 2009 (UTC)

Necessary changes were made Egm6322.s09.bit.sahin 16:16, 10 April 2009 (UTC)

The general solution of the Laplace Eq. is

Eliminating Terms Based on BC[edit]

We can eliminate the following terms that do not satisfy the boundary conditions:

1) , indipendent of

2) , , cannot satisfy the

To show the second one, let's differentiate the general solution with respect to

Since and at , and must be zero to satify the condition that .

Thus the solution has the following form

Using the boundary conditions, we have

Since the only term in the boundary condition is the term with , boundary conditions can only be satisfied for n=4, all other and must be zero. Then we have,

Eventually the solution for the temperature distribution is

The figure below shows the plot of the solution for a given data:

, , ,

MATLAB Plot of Function

MATLAB Plot of Function

--EGM6322.S09.TIAN 17:33, 24 April 2009 (UTC)

You were asked to plot in both polar coordinates and in cartesian coordinates; the figure shown is a plot in cartesian coordinates; I moved this figure to the left for a better presentation. Also provide the matlab codes used to create these plots. Egm6322.s09 15:08, 15 March 2009 (UTC)

My comment above had not been addressed in this updated version. Please take action. Egm6322.s09 22:38, 2 April 2009 (UTC)

Now there it is. EGM6322.S09.TIAN 21:39, 9 April 2009 (UTC)

The Power Law[edit]

--Egm6322.s09.lapetina 02:03, 17 April 2009 (UTC)

The Power Law is a very common relationship in the universe. It is defined as:


The Power Law is observed in classical physics, biology, economics, and many other natural and social sciences. The exponent dominates the nature of the equation. In electrostatics and gravitation, , while in Stefan-Boltzmann equations, . More can be found on the Power Law here.

The inverse of the exponent function is the logarithm.

Application of the Power Law[edit]

The thermal conductivity of solids is summarized in the following graph [5].

Plotting here is on a log-log scale. As a result, the slopes appear linear, rather than exponential.

Log-log vs. Semi-log Plotting --Egm6322.s09.lapetina 17:53, 24 April 2009 (UTC)

The exponential varying thermal conductivity of solids is very important for solving Fourier's Law:

where is the heat flux.

From this equation, we can find the units of in the following fashion:



If we consider the thermal conductivity as where is temperature, we can find the heat flux at any given temperature using the power law by the following equation:


This equation can be solved over a given domain if and are known.

Using data from the aforementioned graph, we see that for diamond,  :


so .

The Wave Equation and String Vibration[edit]

--Egm6322.s09.lapetina 02:04, 17 April 2009 (UTC)

Free Body Diagram of String Vibration

The Wave Equation can be studied by examining the physics of string vibration in one spatial dimension. An excellent book on this topic is The Theory of Sound by Lord Rayleigh.

I never mentioned the book by Rossing et al. I mentioned the classic The Theory of Sound by Lord Rayleigh; see the Lecture plan. Egm6322.s09 15:08, 15 March 2009 (UTC)

Correction made. --Egm6322.s09.lapetina 20:57, 27 March 2009 (UTC)

In the accompanying free-body diagram for Case 1, forces in the direction are :

is very small here, therefore .

Inserting this into the equation results in:


We can neglect second order terms here, since they are very small and nearly cancel, leaving:

, suggesting is constant throughout the string.

In the direction:

where .

Where is small, we can simplify the first two terms as:

multiplied by the constant .

Therefore, for an infinitely small length of string ,

Sign problem in ; inconsistent with Taylor series expansion. The key for the derivation to be independent of any figure of free-body diagram is to treat the tension function and the slope angle as algebraic quantities so that and . Note the plus signs in both cosine terms in , and the plus signs in both sine terms in . This systematic approach is unlike what you used to see in undergraduate courses such as statics, dynamics, etc. See further comments in class. Egm6322.s09 15:08, 15 March 2009 (UTC)

Explanation of sign switch demonstrated in figure is articulated in the text of the original problem.

Alternative solution provided. Image addition will come in near future.

--Egm6322.s09.lapetina 14:33, 16 March 2009 (UTC)


Here, , is mass/unit area.

Wave Equation In 1-D Space (actually a 2-D (x,t) Problem)[edit]

--EGM6322.S09.TIAN 17:28, 24 April 2009 (UTC)


In this case,

Unsteady Heat Equation[edit]

--EGM6322.S09.TIAN 17:24, 24 April 2009 (UTC)

1-D space[edit]

Here, is heat conductivity, is heat source, is heat capacity.

We assume is constant.

2-D space[edit]

We assume is constant here.

General to 3 independent variables (x,y,z)


In 2-D case,

Solution of unsteady heat equation without heat source, which means in polar coordinate.[edit]

Separation of Variables.


Plug into :

The solution for problem 5.12 in Selvadurai (2000) is missing. Egm6322.s09 15:15, 15 March 2009 (UTC)

The statement of the boundary value problem was incomplete (1st question). Do complete this problem in an updated version of report R4. Egm6322.s09 22:38, 2 April 2009 (UTC)

The problem has been completed. 19:09, 10 April 2009 (UTC)


  1. [1]
  5. Simon Phillpot, February 2009 Lecture, Materials Science and Engineering, University of Florida

Signatures[edit] 03:49, 6 March 2009 (UTC) 19:34, 5 March 2009 (UTC)

--Egm6322.s09.three.liu 19:55, 5 March 2009 (UTC)

--EGM6322.S09.TIAN 21:05, 5 March 2009 (UTC)

--Egm6322.s09.lapetina 22:43, 5 March 2009 (UTC)

Egm6322.s09.Three.nav 16:05, 6 March 2009 (UTC)Egm6322.s09.Three.nav

Egm6322.s09.bit.sahin 18:04, 6 March 2009 (UTC) 18:08, 6 March 2009 (UTC)

Egm6322.s09.bit.gk 19:20, 6 March 2009 (UTC)