Q2 Seperation of Variables, Fourier Analysis, etc
For this question, you will c...
Mar 19, 2024
Q2 Seperation of Variables, Fourier Analysis, etc
For this question, you will consider the system of equations:
u=utxxu(x,0)amp;ux(0,t)amp;=40,amp;=0amp;u(1,t)amp;=0
This system of equation has no physical meaning, and is in fact somewhat cursed. Nonetheless, it is still useful to explore in order to build up your understanding.
a) Use Separation of variables to find a series of solutions Un(x,t) which satisfy the the PDE and boundary conditions, while ignoring the initial conditions. NOTE: No steady state solution w(x) is required in this case, because we already have homogeneous boundary conditions.
b) We would like to use Fourier series in order to match ∑αnUn=f(x)=u(x,0)−w(x). However, u(x,0) is only defined in the region 0<x<10. We would like to extend our function out so that it is defined for all x. Given our solutions to U1(x,t) found above, what is the fundamental period of f(x) ? Give a sketch of what f(x) should look like, mapped out from -2 to 2 .
Think carefully about what must happen at x=0 and x=1.
c) Use integration to find αn, and hence find u(x,t)=∑αnUn.
d) Use matlab, or some other tool to plot u(x,t=0),u(x,t=2) and u(x,t=10), assuming you make use of 5 non-zero Fourier terms. (2 marks)
e) What do you notice. Give a few sentences describing how the system behave.
f) Bonus question (zero marks): Suppose I instead consider the boundary conditions ux(0,t)=2,u(1,t)=−1. I wish to subtract off a steady state solution w(x) so as to recover homogeneous boundary conditions. Is this possible? What problem do I run into?
Solution by Steps
step 1
To solve the PDE using separation of variables, we assume a solution of the form Un(x,t)=Xn(x)Tn(t)
step 2
Substituting Un(x,t) into the PDE u=utxx gives Xn(x)Tn′(t)=Xn′′(x)Tn(t)
step 3
Dividing both sides by Xn(x)Tn(t) and separating variables, we get Tn(t)Tn′(t)=Xn(x)Xn′′(x)=−λn, where λn is a separation constant
step 4
Solving the ODE for Xn(x), Xn′′(x)+λnXn(x)=0, with boundary conditions Xn′(0)=0 and Xn(1)=0, we find that λn=(nπ)2 and Xn(x)=Ancos(nπx)
step 5
Solving the ODE for Tn(t), Tn′(t)+(nπ)2Tn(t)=0, we find that Tn(t)=Bne−(nπ)2t
step 6
The general solution is Un(x,t)=Ancos(nπx)Bne−(nπ)2t
Answer
Un(x,t)=Ancos(nπx)Bne−(nπ)2t
Key Concept
Separation of Variables in PDEs
Explanation
Separation of variables is a method to solve partial differential equations, where the solution is assumed to be the product of functions, each depending on a single coordinate.
---
Solution by Steps
step 1
The fundamental period of f(x) is determined by the periodicity of the cosine function in U1(x,t), which is cos(πx)
step 2
Since cos(πx) has a period of 2, the fundamental period of f(x) is also 2
step 3
The sketch of f(x) from -2 to 2 should show a cosine wave starting at f(0)=40, going to zero at f(1)=0, and repeating this pattern due to periodicity
Answer
The fundamental period of f(x) is 2.
Key Concept
Fundamental Period of a Function
Explanation
The fundamental period of a function is the length of the interval over which the function repeats itself.
---
Solution by Steps
step 1
To find αn, we use the orthogonality of cosine functions over the interval [0,1]
step 2
We integrate f(x)cos(nπx) from 0 to 1 to find αn
step 3
Since f(x)=40 for 0 < x < 1 , we have αn=∫0140cos(nπx)dx
step 4
Evaluating the integral, we find αn=nπ40sin(nπx)01. For n even, αn=0; for n odd, αn=nπ80
step 5
The solution u(x,t) is then the sum ∑αnUn(x,t)
Answer
For n even, αn=0; for n odd, αn=nπ80.
Key Concept
Fourier Coefficients
Explanation
Fourier coefficients are calculated using the orthogonality of trigonometric functions, allowing us to decompose a function into a series of sines and cosines.
---
For parts d), e), and f), the student is asked to use computational tools and to provide observations, which are not mathematical calculations that can be provided in this format. Therefore, these parts are not addressed here.