User:Egm6322.s09.mafia/HW4
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)
Contents
The Classification of Second Order Partial Differential Equations[edit]
Classifications[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:
...(1)
This PDE is elliptic
HW: What would be the classification of the diffusion operator in polar coordinates? Would it make sense if the classification changed along with transformation of coordinates? i.e If we transformed the diffusion operator from cartesian to polar coordinates, would its classification change? 

Egm6322.s09.Three.nav 13:36, 24 April 2009 (UTC) The diffusion operator D(.) is given by div(grad(.)). Recall, In cartesian coordinates, div(grad(u))= In polar coordinates, div(grad(u))= A PDE of the form (1) can be characterized by the nature of the determinant Elliptic if acb^{2}>0 Parabolic if acb^{2}=0 Hyperbolic if acb^{2}<0 Comparing the cartesian form of D(u)with (1), we see that a= 1, b=0, c=1 D(u) in cartesian coordinates is elliptic. In polar coordinates, the matrix X= is transformed and written as the matrix given by
To answer the first part of the HW, it would NOT and does not make sense if the classification changed along with the transformation. This is because no matter the type of coordinate system used, the physics of the problem remains the same. Hence if a different coordinate system was used, then a change in the nature of the transformed PDE would imply a change in the physics of the problem which is nonsensical. Hence transformation should not change classification Lets try to prove this. If (1) was transformed into another set of coordinates, the classification would change only if the determinant det() was different from det(X).
= (acb^{2})(C^{2}+S^{2}) But which is the same as det(X). Hence it is proved that Even if the PDE of form (1) in cartesian coordinates is transformed to another coordinate system (say polar), the classification remains constant. Consider the polar form of the diffusion operator D(u)=
Determinant= acb^{2}= >0 r. Hence the polar form of the diffusion operator is also 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]
Egm6322.s09.Three.ge 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]
Egm6322.s09.Three.ge 17:58, 24 April 2009 (UTC)
For an axissymmetric 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 A_{o} and B_{o} are constants.
It should be noted that if the domain includes the origin (where r=0) then A_{o} must be zero for a finite solution.
Separation of Variables[edit]
Egm6322.s09.Three.ge 17:58, 24 April 2009 (UTC)
If the problem is not axissymmetric, the Laplace equation may be solved using separation of variables.
Multiplying the Laplace equation by r^{2} 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:
HW: Why is the separation constant n^{2} positive? 

Egm6322.s09.Three.nav 13:38, 24 April 2009 (UTC) The sign of the separation constant n^{2} potentially determines the roles of the 'r' and ' coordinates. In the method used to solve Laplace's Equation here, a +ve sign implies periodicity in the coordinate. One of the applications of solving the Laplace's Equation is in electrostatics. Here the physics of the problem generally dictates periodicity in the direction. Hence the sign of n^{2} is positive. Also here the equation for 'r' = is a nonlinear equation and one that cannot, by its very nature ,dictate oscillation or periodicity. In such a case, we usually prefer having the oscillation dependence in the azimuthal () direction. 
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. Egm6322.s09.bit.la 22:38, 2 April 2009 (UTC)
And results in the following solution for .
If n=0:
HW: The solution to the Laplace Equation is called harmonic. Why? 

Egm6322.s09.Three.nav 13:38, 24 April 2009 (UTC) Consider , one of the equations obtained through the separation of variables. This resembles the equation of harmonic motion (the equation characterizing the motion of a simple harmonic oscillator). This is why the solution G() is called harmonic. 
In general, we want the solution to be periodic, such that:
Why do we want the solution to be periodic? 

Our solution resembles that of a harmonic oscillator which means it will complete one complete cycle of motion going from to . So when it completes this cycle, its solution at cannot be different from the solution at which is why we need the solution to be periodic. 
The general form of the Laplace equation in polar coordinates takes the form that follows.
Another Axisymmetric Problem[edit]
EGM6322.S09.TIAN 17:20, 24 April 2009 (UTC)
Substituting into the equation yields:

Multiplying out the equations and rearranging terms gives:
=

One can see clearly that the first five terms cancel. This leaves the following equation:
Which, after some manipulation, results in:
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 reposted, now with the correct result. Egm6322.s09.Three.ge 21:02, 6 April 2009 (UTC)
page 14.2
Application:
Example:
where T_{0} is a constant
Principle of Superposition[edit]
Egm6322.s09.xyz 16:29, 4 April 2009 (UTC)
The governing PDE is given as . This PDE is linear and therefore the solution can be expressed using the Principle of Superposition:
i.e the solution = solution for = constant solution for
The proof of the linearity of the and operators was presented by Team Mafia in R2. For completeness, the relevant portions of the proof are presented again below:
note: is linear
and is linear
note: is linear because it is another differential operator. Let and
The proof of linearity of each operator within the PDE yields the conclusion that the entire PDE is also linear. The Principle of Superposition is applicable to linear PDEs. Applying the Principle of Superposition, the original PDE can be split into two separate parts such that the temperature is
The solutions for and constitutes two separate problems that satisfy the following:
such that
and
such that
Problem P:
PDE:
General Solution: Equation 2 p. 204.
Superposition:
Prob P1:
such that:
Prob P2:
such that:
Egm6322.s09.xyz 16:33, 4 April 2009 (UTC)
represents the solution to the first portion of the temperature profile. It simply states that the temperature is constant. The general form of the solution to the heat equation is given as (see #Separation of Variables):
In order for this general form to converge to , all the coefficients need to be equal to zero. i.e.
For the Problem
Egm6322.s09.xyz 17:23, 24 April 2009 (UTC)
For problem P2:
because this is an axisymmetric problem with the origin ( r = 0 ) within the domain. (see Axisymmetric Problems)
. By inspection of the general form (see above) for the solution to these types of problems, the solution needs to converge to be a function of for . For these boundary conditions and , therefore the coefficients of these terms must be equal to zero.
The expression for is then
Egm6322.s09.xyz 16:33, 4 April 2009 (UTC)
At the boundary (r = a), the expression for is rewritten as:
such that
The resulting temperature profile is a function of only. Upon inspection, the term should be forced to go to zero. Otherwise the equality would not be satisfied. Therefore, the coefficient for all
Substitution of r = a into the equal for , yields:
For the case where:
For the case where:
The equality DOES NOT hold for all values where . Taking the case where as an illustrative example, the resulting expression would be:
This illustrates that for all
Final Solution:
In general, for arbitrary function but periodic
i.e.
for all and any K = constant
this is not periodic (not acceptable)
The general solution of the Laplace equation:
Where the partial differential equation governing the problem is:
Applying the boundary conditions explained in the previous homework
The resulting solution is:
where: r^{n}=a^{n}
Fourier Coefficients[edit]
Egm6322.s09.Three.nav 13:39, 24 April 2009 (UTC)
How do we derive the fourier coeffcients C_{0}, A_{n} and B_{n}?
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 C_{0}, A_{n} and B_{n}, 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 A_{n} 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 B_{n} 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 C_{0},A_{n},B_{n}
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)
Jean Baptiste Joseph Fourier was a French mathematician and physicist best known for first devising Fourier Series and their applications to heat flow problems. He was born in Auxerre and orphaned at the age of 9. He became a chair at the École Polytechnique at the height of his career. Fourier died in Paris. More information is available here
To show that is non linear
Let,
be an operator,such that is linear with respect to if,
Therefore, in the present problem ,assuming 2D, we have,
\therefore,we can see that
is non linear
Egm6322.s09.bit.gk 20:41, 24 April 2009 (UTC)
We can see from the free body diagram ,equating the horizontal forces ,we have
and equating the vertical forces, we have,
where is the transverse load acting on the membrane
When is small ,we can assume,
the above equation becomes,
Let this equation be
but, is the slope the membrane with respect to the and axes
where the displacement is given as
and
and
Substituting the slopes in ,we have ,
Let this be equation
We have Taylor series expansion as ,
substituting taylor series expansion in equation ,we have
Egm6322.s09.bit.gk 20:39, 24 April 2009 (UTC)
Egm6322.s09.xyz 16:34, 4 April 2009 (UTC)
The following is the definition of Orthogonal Functions as presented on the Wolfram MathWorld website^{[2]}:
"Two functions and are orthogonal over the inverval with weighting function if
If, in addition,
and
the functions and are said to be orthonormal"
Additional information can be found at Wiki Orthogonality^{[3]} and Dictionary.com^{[4]}
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
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 .
The boundaries are kept insulated which means that there is no heat flow at and . So,
and
Since at insulated surfaces we have
So, we obtain that
at and
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:
, , ,
EGM6322.S09.TIAN 17:33, 24 April 2009 (UTC)
function [x] = plot(A,b,x0)
% figure out T
[r,theta] = meshgrid(1:0.05:2,0:pi/40:pi/2);
a=1;
b=2;
Ta=5;
Tb=20;
for i=1:21;
for j=1:21;
T(i,j)=( (a^4 * b^4 *Tb / (b^8  a^8)) * (r(i,j) ^4 / a^4  a^4 / r(i,j) ^4)  (a^4 * b^4 *Ta / (b^8  a^8))* (r(i,j) ^4 / b^4  b^4 / r(i,j) ^4))*cos(4*theta(i,j));
end;
end;
% plot Polar
figure(1);
surf(r,theta,T); colorbar ;
title('Polar Coordinate');
Xlabel('R');
Ylabel('Theta');
Zlabel('T');
% plot Cartesian
for i=1:21;
for j=1:21;
x(i,j)=r(i,j)*cos(theta(i,j));
y(i,j)=r(i,j)*sin(theta(i,j));
end;
end;
figure(2)
surf(x,y,T)
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 StefanBoltzmann 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 loglog scale. As a result, the slopes appear linear, rather than exponential.
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:
while
Therefore:
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, :
while
so .
Egm6322.s09.lapetina 02:04, 17 April 2009 (UTC)
For , we want to find . We can estimate for diamond using the slope of graphite parallel to layers:
.
Extrapolating backwards shows its value is .
Therefore, for all ,
The Wave Equation and String Vibration[edit]
Egm6322.s09.lapetina 02:04, 17 April 2009 (UTC)
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 freebody 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 ,
Egm6322.s09.lapetina 02:04, 17 April 2009 (UTC)
Case 4: Vertical Reflection
This is the most trivial of the other three cases. In the direction, this changes only the sign on . Because this value is squared, we are left with:
Case 2: Modified Geometry
For cases two and three, we have a slightly more complicated situation, as both ends of the string point down and up, respectively. This negates our ability to simplify the equation in the same way for the string as shown. Looking first at the accompanying figure, we can our case on the left side.
However, recall the string is infinitesimally small. As seen in the accompanying image, we can simply examine a fraction of the infinitesimally small string , and encounter the same geometry. In the image the sum of forces in the direction on the left side of the image can be expressed as:
while the right side (showing a fraction of the infinitesimally small string) can be expressed as:
This can be simplified using the assumption:
,
leaving:
.
Neglecting second order terms here, since they are very small and nearly cancel, leaves:
, suggesting is constant throughout the string, which by definition continues to .
This means that even for geometries where the ends of the strings point in the same direction, .
In the direction, we start with:
.
However, for a fraction of the string , (i.e., the modified version shown in the right side of the image), this equation changes to:
.
Now, where is small, we can again simplify the first two terms as:
multiplied by the constant .
This leaves us with the same equation as for Cases One and Four:
Case Three: Vertical Reflection of Case Two
Our signs are merely switched again, as they were between Cases One and Four.
Essentially, no matter what the geometry of the string (so long as it is simple and does not cross, and is always differentiable), the equation of motion is:
.
This is because the original derivation is for a nearly linear string, and the can always be viewed as linear, as:
.
An alternative means of solving this problem is to view and as algebraic quantities rather than physical entities. The accompanying image shows the only free body diagram needed.
Here, we see that:
(1) (2)
For equation 1, we can assume is small, making , leading to:
.
Therefore:
This changes equation 2 to:
where, because is small:
higher order terms.
This leaves us with:
Sign problem in ; inconsistent with Taylor series expansion. The key for the derivation to be independent of any figure of freebody 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)
Egm6322.s09.lapetina 02:05, 17 April 2009 (UTC)
In two dimensions, the vertical position of a membrane of density with constant thickness (and corresponding constant strength and flexibility) and loading can be described as:
, where and are Cartesian positions, and is time. With no loading, the membrane rests in the plane
The stress in the membrane, is constant throughout. For all time and for all and , the membrane surface is differentiable by and , and small angle approximations can be made.
As shown in the image, we will examine a small square portion of width and height .
For this small portion, There are five forces acting on it, the tensions, , and the load . The tensions all act normal to the edge of the small portion, and at an angle to the plane, , respectively.
The tensions are applied at coordinates, respectively.
Newton's Second Law in the vertical direction is then:
Using the small angle approximation,
which is approximately equal to the derivative of the membrane in the direction of the tension force at each location.
Therefore:
.
The magnitudes of are equivalent to the constant stress multiplied by the area over which it is applied, such that:
Let us define:
and .
Substituting these values into Newton's Law and dividing by and leaves us with :
.
Taking the limits as and approach zero, we are left with a Poisson Equation:
where is the mass per unit area.
This is equivalent to:
.
For steady state, and a massless membrane, the equation becomes:
Applications[edit]
Here, , is mass/unit area.
Wave Equation In 1D Space (actually a 2D (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)
1D space[edit]
Here, is heat conductivity, is heat source, is heat capacity.
We assume is constant.
2D space[edit]
We assume is constant here.
General to 3 independent variables (x,y,z)
Here,
In 2D case,
Solution of unsteady heat equation without heat source, which means in polar coordinate.[edit]
Separation of Variables.
Here,
Plug into :
The solution for problem 5.12 in Selvadurai (2000) is missing. Egm6322.s09 15:15, 15 March 2009 (UTC)
An unloaded weightless membrane roof over an annular enclosure. The outer circular boundary is r=b, and the inner circular boundary, r=a, is subject to the following displacement:
where and are constants.
i) Formulate the Boundary value problem:
The boundary conditions are:
ii) Develop an expression for the membrane, given that:
where are constants.
Solution: Using boundary condition 1, the equation becomes:
Clearly,
and,
Thus the equation simplifies to:
Using the 3rd boundary condition:
One can clearly see that since the coefficient of the sine term is 1, all n terms that are not equal to one must be zero.
And,
It can also be seen that since there is no cosine term in the boundary equation,
Thus the equation reduces to:
Where,
Reapplying boundary condition three gives:
From this equation it is clearly evident that,
Egm6322.s09.bit.gk 20:59, 24 April 2009 (UTC) These equations and the last two boundary conditions:
Give 4 equation to solve for 4 unknown coefficients.
The resulting system of equations may be expressed as a matrices:
Using the matrix to solve for the constants one finds that:
iii) Calculate the resultant force and moment necessary to maintain the inner rigid disk shaped region in the displaced position.
To find the net force necessary to maintain the inner disk in its position, one must assume the membrane has an effective spring constant of k. The force is then k multiplied by the displacement w. To get the net force we integrate around the edge of the disk.
To find the moment necessary to keep the disk in its tilted position requires a bit more consideration. The moment results from an unbalanced force applied at some distance y. It should be noted that since the term is the same all along the boundary, it does not contribute to the moment necessary to maintain the disk's position.
Where the value of F is already given as:
To get the total moment the sum of the discrete should be taken. By increasing the number of and simultaneously decreasing their magnitude, the summation becomes on integral.
Thus,
This integral must be taken over the entire circle. Thus the moment required is:
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. Egm6322.s09.Three.ge 19:09, 10 April 2009 (UTC)
References[edit]
 ↑ [1]
 ↑ http://mathworld.wolfram.com/OrthogonalFunctions.html
 ↑ http://en.wikipedia.org/wiki/Orthogonal/
 ↑ http://dictionary.reference.com/browse/orthogonality
 ↑ Simon Phillpot, February 2009 Lecture, Materials Science and Engineering, University of Florida
Signatures[edit]
Egm6322.s09.xyz 03:49, 6 March 2009 (UTC)
Egm6322.s09.Three.ge 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) Egm6322.s09.bit.la 18:08, 6 March 2009 (UTC)
Egm6322.s09.bit.gk 19:20, 6 March 2009 (UTC)