Partial Differential Equations result -- How to simplify trig series?

In summary: I hope, instead of starting at 0 and having to use 4k+1, 4k+3, 2k+1, 2k+3 and so on and so forth. It's just for the sake of time since the exam is tomorrow.Also I noticed that the last series could be further reduced tou(x, t) = \sum_{k=0}^{\infty} \frac{4}{(2k+1)^2}(-1)^{k}\exp(-(2k+1)^{2}t)\sin((2k+1)x)Which is the same asu(x, t) = \sum_{k=0}^{\infty
  • #1
AnotherParadox
35
3
Homework Statement
Solve the Boundary value problem
[itex]u_{t}=u_{xx}[/itex]



[itex]u(0, t) = u(\pi ,t)=0[/itex]



[itex]u(x, 0) = f(x)[/itex]

[tex]f(x)=\left\{\begin{matrix}

x; 0 < x < \frac{\pi}{2}\\

\pi-x; \frac{\pi}{2} < x < \pi

\end{matrix}\right.[/tex]
Relevant Equations
[tex]u(x, t) = \sum_{n=1}^{\infty} A_{n} \exp(- \lambda t)\sin(\sqrt{\lambda}\cdot x))[/tex]
[tex]A_{n}=\frac{2}{L}\int_{0}^{L}f(x)sin(\sqrt{\lambda}\cdot x)dx[/tex]
Solve the boundary value problem
Given
[itex]u_{t}=u_{xx}[/itex]

[itex]u(0, t) = u(\pi ,t)=0[/itex]

[itex]u(x, 0) = f(x)[/itex]
[tex]f(x)=\left\{\begin{matrix}
x; 0 < x < \frac{\pi}{2}\\
\pi-x; \frac{\pi}{2} < x < \pi
\end{matrix}\right.[/tex]

L is π - 0=π

λ = α2 since 0 and -α lead to trivial solutions
Let
[itex]u = XT[/itex]

[itex]X{T}'={X}''T[/itex]

[itex]\frac{{X}''}{X}=\frac{{T}'}{T}=-\lambda [/itex]

[itex]{X}''+\lambda =0[/itex]

[itex]X(0)=X(\pi)=0[/itex]

[itex]{T}'+\lambda T=0[/itex]

[itex]X(x)=C_{1}\cos (\sqrt{\lambda }\cdot x)+C_{2}\sin(\sqrt{\lambda }\cdot x)[/itex]

[itex]\sqrt{\lambda } \cdot \pi = n\pi[/itex]

[itex]\lambda =n ^{2}; n=1,2,3,..[/itex]

[itex]{T}'+n^{2}T=0[/itex]

(I'm not entirely sure how t got into this next step, I know T is actually T(t) but still not sure what it's doing in the next step with the exponential)

[itex]T(t)=C_{3}\exp (-n ^{2}t)[/itex]

[itex]u_{n}(x,t)=C_{2}C_{3} \exp(-n^{2}t)\sin(nx)=A_{n} \exp(-n^{2}t)\sin(nx)); n = 1, 2, 3, ..[/itex]
[tex]u(x, t) = \sum_{n=1}^{\infty} A_{n} \exp(-n^{2}t)\sin(nx))[/tex]
[tex]u(x,0)=f(x)[/tex]
[tex]A_{n}=\frac{2}{\pi}\left [\int_{0}^{\pi/2}x\sin(nx)dx + \int_{\pi/2}^{\pi}(\pi-x)\sin(nx)dx \right ][/tex]
[tex]=\frac{2}{\pi}\left [\left [ \frac{x(-\cos(nx))}{n} \right ]_{0}^{\pi/2} + \int_{0}^{\pi/2} \frac{\cos(nx)}{n}dx + \left [\frac{(\pi-x)(-\cos(nx))}{n} \right ] _{\pi/2}^{\pi} - \int_{\pi}^{\pi/2}\frac{\cos(nx)}{n}dx \right ][/tex]
[tex]=\frac{-1}{n}\cos (\frac{n \pi}{2}) + \frac{2}{\pi n^{2}}\sin(\frac{n \pi}{2}) + \frac{1}{n}\cos (\frac{n \pi}{2}) + \frac{2}{\pi n^{2}}\sin(\frac{n \pi}{2})[/tex]
[tex]=\frac{4}{\pi n^{2}}\sin(\frac{n \pi}{2})[/tex]
[tex]u(x, t) = \sum_{n=1}^{\infty} \frac{4}{\pi n^{2}}\sin(\frac{n \pi}{2}) \exp(-n^{2}t)\sin(nx))[/tex]

Ok so far so good. An was found and put back into the general form of the solution.

Here's where it gets tricky
Apparently there is further simplification of this, now these next steps were provided to me by my instructor (as were the previous too) during lecture and I cannot make sense of them. I see what is supposed to happen, cos and sin follow a pattern on intervals of π or nπ where they evaluate to predictable values like 1,0,-1 etc. and we can ignore some values because they lead to all 0s or the trivial solution.

Apparently he continues by the following steps

[tex]
u(x, t) = \sum_{k=0}^{\infty} \frac{4}{\pi (4k+1)^{2}}\exp(-(4k+1)^{2}t)\sin((4k+1)x)) -\sum_{k=0}^{\infty} \frac{4}{\pi (4k+3)^{2}}\exp(-(4k+3)^{2}t)\sin((4k+3)x))[/tex]

Separated into two series and replaced n with 4k+1 and 4k+3, eliminated the first trig function with some version of previously stated pattern recognition and started the series at 0 instead of 1

Now I can't even verify if this is equivalent to the other series let alone figure out the pattern or logic to this to figure it out for other scenarios.

And someone stated that it can be reduced further into

[tex]
u(x, t) = \sum_{k=0}^{\infty} (-1)^{k}\exp(-(2k+1)^{2}t)\sin((2k+1)x))
[/tex]
If someone can please break it down and explain how this is done I will be very grateful as the exam is soon and I cannot seem to figure out how to follow the process to come up with these series reductions

Thank you!

Also it would be nice to know how t ends up in this step[itex]{T}'+n^{2}T=0[/itex]
[itex]T(t)=C_{3}\exp (-n ^{2}t)[/itex]

but not nearly as important.

Thanks again
 
Last edited:
Physics news on Phys.org
  • #2
What values can ##\sin(n\frac\pi 2)## take?
For each such value, what values of n produce that result?
 
Last edited:
  • Like
Likes AnotherParadox
  • #3
AnotherParadox said:
Also it would be nice to know how t ends up in this stepT′+n2T=0T′+n2T=0{T}'+n^{2}T=0
T(t)=C3exp(−n2t)T(t)=C3exp⁡(−n2t)T(t)=C_{3}\exp (-n ^{2}t)

but not nearly as important.
This is just solving a first order ODE in ##t##. Clearly the result must depend on ##t## unless ##T’ = 0## or you have a trivial solution.

Having taught and authored a book on this subject, I would say that if you have an exam soon then understanding this would be far more important than being able to simplify the resulting expressions.

Edit: Rereading this I see an ambiguity in my phrasing so just for clarity: I did not teach the book, I taught students ...
 
Last edited:
  • Like
Likes AnotherParadox
  • #4
haruspex said:
What values can ##\sin(n\frac\pi 2)## take?
For each such value, what values of n produce that result?

Thanks for the response

Here's my thought process

It can't take 0 since sin(0*pi/2) and sin(0*x) is 0 and that leads to the trivial solution.
It can take 1 since sin(1*pi/2) does not equal 0 it equals 1
It can't take 2 since sin(2*pi/2) = 0 and the trivial solution
It can take 3 since sin(3pi/2) = -1 and not zero
It can't take 4 since sin(4pi/4) = 0 and the trivial solution

So I know it can't take any even number since they all lead to 0 and thus the trivial solution

Odd numbers i.e. n = 1, 3, 5, 7.. lead to the pattern 1, -1, 1, -1

Which is a similar pattern as (-1)^n which is for n = 0, 1, 2, 3,: 1, -1, 1, -1

So this can be (-1)^(n) which is for n = 0, 1, 2, 3: To match that part of the expression

Evaluating the rest of the series 2k+1 is needed instead of n is needed so that it's 1,3,5,7 for the 0,1,3,4.. I see how this logic was applied to come up with the 4k+1 and 4k+3 so that checks out too now

So far so good however I think he may of made a mistake and I would like confirmation of this

So he changes

[tex]u(x, t) = \sum_{n=1}^{\infty} \frac{4}{\pi n^{2}}\sin(\frac{n \pi}{2}) \exp(-n^{2}t)\sin(nx)[/tex]

to
[tex]
u(x, t) = \sum_{k=0}^{\infty} \frac{4}{\pi (4k+1)^{2}}\exp(-(4k+1)^{2}t)\sin((4k+1)x) -\sum_{k=0}^{\infty} \frac{4}{\pi (4k+3)^{2}}\exp(-(4k+3)^{2}t)\sin((4k+3)x)[/tex]

I think this checks out and the series match however on the next reduction I think a few things are left out because he shows
[tex]
u(x, t) = \sum_{k=0}^{\infty} (-1)^{k}\exp(-(2k+1)^{2}t)\sin((2k+1)x)
[/tex]

Which looks even better but I believe it's missing the front part and it should be[tex]
u(x, t) = \sum_{k=0}^{\infty} \frac{4}{\pi (2k+1)^2} (-1)^{k}\exp(-(2k+1)^{2}t)\sin((2k+1)x)
[/tex]

Or even more

[tex]u(x, t) = \frac{4}{\pi} \sum_{k=0}^{\infty} \frac{1}{ (2k+1)^2} (-1)^{k}\exp(-(2k+1)^{2}t)\sin((2k+1)x)[/tex]

If anyone can confirm this and/or tell me about other short-cuts/quick verifying tricks/ and ways of thinking that would be much appreciated.

Orodruin said:
This is just solving a first order ODE in ##t##. Clearly the result must depend on ##t## unless ##T’ = 0## or you have a trivial solution.

Having taught and authored a book on this subject, I would say that if you have an exam soon then understanding this would be far more important than being able to simplify the resulting expressions.

Edit: Rereading this I see an ambiguity in my phrasing so just for clarity: I did not teach the book, I taught students ...

Np, I understood what you meant, thanks for the response.

Here's my work

[tex]{T}'+n^2T=0[/tex]
[tex]{T}'=-n^2T[/tex]
[tex]\frac{{T}'}{T}=-n^2[/tex]
[tex]\ln\left ( \frac{{T}'}{T} \right )=-n^2[/tex]
[tex]\int \left ( \frac{{T}'}{T} \right )=\int-n^2dt[/tex]

I'll stop here since I figured it out. the -n^2 has to be integrated with respect to t which multiplies t in. Though it makes me wonder why not with respect to n. I suppose for the boundary conditions it needs to be t? I kind of see what is being said with T'=/=0

I agree it's more important to understand this ODE but I can easily write in the correct step and receive full credit but he counts off significantly for not reducing the series if I understand his grading correctly.
 
Last edited:
  • #5
AnotherParadox said:
however on the next reduction I think a few things are left out because he shows
[tex]
u(x, t) = \sum_{k=0}^{\infty} (-1)^{k}\exp(-(2k+1)^{2}t)\sin((2k+1)x)
[/tex]

Which looks even better but I believe it's missing the front part and it should be[tex]
u(x, t) = \sum_{k=0}^{\infty} \frac{4}{\pi (2k+1)^2} (-1)^{k}\exp(-(2k+1)^{2}t)\sin((2k+1)x)
[/tex]
I agree, he seems to have dropped the factor ##\frac 4{\pi}##.
 
  • Like
Likes AnotherParadox
  • #6
haruspex said:
I agree, he seems to have dropped the factor ##\frac 4{\pi}##.
Also the term for 1/(2k+1)^2, right?
 
  • #7
AnotherParadox said:
Also the term for 1/(2k+1)^2, right?
Ah, yes... and that's worse.
 
  • Like
Likes AnotherParadox
  • #8
AnotherParadox said:
Though it makes me wonder why not with respect to n.
Because the derivative is with respect to t while n is a constant.

AnotherParadox said:
I agree it's more important to understand this ODE but I can easily write in the correct step and receive full credit but he counts off significantly for not reducing the series if I understand his grading correctly.
This is really bad pedagogy if it is correct.
 
  • Like
Likes AnotherParadox
  • #9
Orodruin said:
Because the derivative is with respect to t while n is a constant.

I like this reasoning

Here's some work I thought up with considering this more

[tex]{T}'+n^2T=0[/tex]
[tex]\frac{{dT}}{dt}+n^2T=0[/tex]
[tex]\frac{{dT}}{dt}=-n^2T[/tex]
[tex]\frac{\frac{{dT}}{dt}}{T}=-n^2[/tex]
[tex]\int \left ( \frac{dT}{T} \right )=\int-n^2dt[/tex]

from there it's ln -> exp

Let me know if this checks out
Orodruin said:
This is really bad pedagogy if it is correct.
I agree. I wasn't aware of how bad it was though. To be fair it is not a math major course, it's an engineering elective (still important imho). I could very well leave out a lot of these steps entirely and not be counted off however the reduction portion could be anywhere to 10%-30% of the question worth. He often teaches from memory (impressively) and makes mistakes that don't get corrected until I look over the class notes on my own which can be awfully difficult at times (wish he would check /confirm his answers from a record). Still a decent teacher, just a bit aged and dated in his ways. Also different standards for non math majors I suppose.
 
Last edited:
  • #10
AnotherParadox said:
Let me know if this checks out
Yes. The more standard approach is to assume ##T = C e^{rt}## and inserting into the differential equation to solve for ##r##, but what you did is certainly fine and does not depend on an ansatz.
 
  • Like
Likes AnotherParadox

1. What are Partial Differential Equations (PDEs)?

Partial Differential Equations (PDEs) are mathematical equations that involve multiple independent variables and their partial derivatives. They are used to describe physical phenomena such as heat transfer, fluid flow, and electromagnetic fields.

2. How are PDEs different from ordinary differential equations (ODEs)?

PDEs involve partial derivatives, while ODEs only involve ordinary derivatives. This means that PDEs are used to model systems with multiple variables, while ODEs are used to model systems with only one variable.

3. What is the result of solving a PDE?

The result of solving a PDE is a mathematical expression that describes the behavior of the system being modeled. This expression can be used to make predictions and analyze the system's behavior under different conditions.

4. How do I simplify a trigonometric series in a PDE?

To simplify a trigonometric series in a PDE, you can use techniques such as the Fourier series or the Laplace transform. These methods allow you to express the series in terms of simpler trigonometric functions or algebraic equations.

5. Can PDEs be solved analytically?

Some PDEs can be solved analytically, meaning that an exact solution can be found using mathematical techniques. However, many PDEs do not have analytical solutions and must be solved numerically using computational methods.

Similar threads

  • Calculus and Beyond Homework Help
Replies
3
Views
419
  • Calculus and Beyond Homework Help
Replies
1
Views
221
  • Calculus and Beyond Homework Help
Replies
1
Views
539
  • Calculus and Beyond Homework Help
Replies
11
Views
366
  • Calculus and Beyond Homework Help
Replies
2
Views
919
  • Calculus and Beyond Homework Help
Replies
1
Views
348
  • Calculus and Beyond Homework Help
Replies
1
Views
445
  • Calculus and Beyond Homework Help
Replies
1
Views
288
  • Calculus and Beyond Homework Help
Replies
16
Views
567
  • Calculus and Beyond Homework Help
Replies
6
Views
762
Back
Top