Solving inhomogenous wave equation

In summary: The boundary conditions for w(x,t) are:w(0,t)=1, w_{x}(1,t)=0w(x,0)=u(x,0)-S(x)=1-\frac{x^{3}}{6}+\frac{x}{2}, w_{t}(x,0)=0So, for w(0,t)=1 we get:w(0,t)=\frac{A}{B}sin(0)+cos(0)=\frac{A}{B}=1And for w_{x}(1,t)=0 we get:w_{x}(x,t)=\frac{A}{B}kcos(kx)-ksin(kx) ==> w_{x}(1
  • #36
Hmm, does that mean we are looking for Fourier sinus series of f(x)=x^3/6-x/2?
Is it [tex]A_{n}D_{n}=2\int^{1}_{0}(\frac{x^{3}}{6}-\frac{x}{2})(\frac{(2n-1)\pi}{2}x)dx[/tex]
 
Physics news on Phys.org
  • #37
wahoo2000 said:
Hmm, does that mean we are looking for Fourier sinus series of f(x)=x^3/6-x/2?
Is it [tex]A_{n}D_{n}=2\int^{1}_{0}(\frac{x^{3}}{6}-\frac{x}{2})(\frac{(2n-1)\pi}{2}x)dx[/tex]

Yes, you can evaluate the integral by using integration by parts 3 times.
 
  • #38
I think I made a typing error there, of course it should be

[tex]
A_{n}D_{n}=2\int^{1}_{0}(\frac{x^{3}}{6}-\frac{x}{2})sin(\frac{(2n-1)\pi}{2}x)dx
[/tex]
, right? I'll attack it right now :)
 
  • #39
wahoo2000 said:
I think I made a typing error there, of course it should be

[tex]
A_{n}D_{n}=2\int^{1}_{0}(\frac{x^{3}}{6}-\frac{x}{2})sin(\frac{(2n-1)\pi}{2}x)dx
[/tex]
, right? I'll attack it right now :)

Yes, I just assumed you meant to put the sin in there.
 
  • #40
Ok..puh, now I have tried.. I am not too confident it is OK though :( so much to write if I write everything, but I can show how I did and thought if this is completely wrong..
This is what I got:

[tex]A_{n}D_{n}=-\frac{2}{(2n+1)^{2}\pi^{2}}(-(-1)^{n})-\frac{16}{(2n+1)^{2}\pi^{2}}-\frac{64}{(2n+1)^{4}\pi^{4}}[/tex]
 
  • #41
wahoo2000 said:
Ok..puh, now I have tried.. I am not too confident it is OK though :( so much to write if I write everything, but I can show how I did and thought if this is completely wrong..
This is what I got:

[tex]A_{n}D_{n}=-\frac{2}{(2n+1)^{2}\pi^{2}}(-(-1)^{n})-\frac{16}{(2n+1)^{2}\pi^{2}}-\frac{64}{(2n+1)^{4}\pi^{4}}[/tex]

This is incorrect. Show me what you did when you applied integration by parts the first time.
 
  • #42
first time was [tex]2[(\frac{x^{3}}{6}-\frac{x^{2}}{2})(cos(\frac{(2n-1)\pi}{2}x)]^{1}_{0}

+\frac{4}{(2n+1)\pi}\int^{1}_{0}(\frac{3x^{2}}{6}-\frac{2x}{2})(cos(\frac{(2n-1)\pi}{2}x)dx[/tex]

But I saw now when writing this that I used x^2/2 instead of x/2 in my calculations.. :( :(
 
  • #43
wahoo2000 said:
first time was [tex]2[(\frac{x^{3}}{6}-\frac{x^{2}}{2})(cos(\frac{(2n-1)\pi}{2}x)]^{1}_{0}

+\frac{4}{(2n+1)\pi}\int^{1}_{0}(\frac{3x^{2}}{6}-\frac{2x}{2})(cos(\frac{(2n-1)\pi}{2}x)dx[/tex]

But I saw now when writing this that I used x^2/2 instead of x/2 in my calculations.. :( :(

Try it with x/2 and show me what you get.
 
  • #44
Here is the first step.. if you are online, just tell me if that step i correct before I do it all again...

[tex]A_{n}D_{n}=2[(\frac{x^{3}}{6}-\frac{x}{2})(cos(\frac{(2n-1)\pi}{2}x)]^{1}_{0}

+\frac{4}{(2n+1)\pi}\int^{1}_{0}(\frac{3x^{2}}{6}-\frac{1}{2})(cos(\frac{(2n-1)\pi}{2}x)dx
[/tex]
where the first term, [tex]2[(\frac{x^{3}}{6}-\frac{x}{2})(cos(\frac{(2n-1)\pi}{2}x)]^{1}_{0}=0[/tex]
 
  • #45
I worked it through anyway... all terms except the last became zero.
I now got [tex]A_{n}D_{n}=\frac{32}{(2n-1)^{3}\pi^{3}}*-(-1)^{n}[/tex]
 
  • #46
wahoo2000 said:
I worked it through anyway... all terms except the last became zero.
I now got [tex]A_{n}D_{n}=\frac{32}{(2n-1)^{3}\pi^{3}}*-(-1)^{n}[/tex]

Do you mean:
[tex]
A_{n}D_{n}=\frac{32}{(2n-1)^{4}\pi^{4}}*-(-1)^{n}
[/tex] ?
If so, then your right. What are w(x,t) and u(x,t) then?
 
  • #47
Hmm.. no, actually i got ^3... I will show you how I did and then you can see where it went wrong. I don't write out the terms that went zero.
[tex]
A_{n}D_{n}=+\frac{4}{(2n+1)\pi}\int^{1}_{0}(\frac{3x^{2}}{6}-\frac{1}{2})(cos(\frac{(2n-1)\pi}{2}x)dx = [/tex]
[tex]=-\frac{16}{(2n+1)^{2}\pi{2}}\int^{1}_{0}(x*sin(\frac{(2n-1)\pi}{2}x)dx = \frac{32}{(2n+1)^{3}\pi{3}}\int^{1}_{0}(cos(\frac{(2n-1)\pi}{2}x)dx=[/tex]
[tex]=\frac{32}{(2n+1)^{3}\pi{3}}[(sin(\frac{(2n-1)\pi}{2}x)]=\frac{32}{(2n-1)^{3}\pi^{3}}*-(-1)^{n}[/tex]
edit: I think I lost one inner derivative there at the end...but it doesn't make sense to me anyway..in that case it would be 64/((2n+1)^4*pi^4)
:/


Hope we can find the mistake. I go on with "your" AnDn meanwhile.

Then it means that
[tex] w(x,t)=-\frac{32}{\pi^{4}} \sum_{n=1}^{\infty} \frac{(-1)^{n}}{(2n+1)^{4}}sin(\frac{(2n-1)\pi}{2}x)cos(\frac{(2n-1)\pi}{2}t)
[/tex] and
[tex]u(x,t)=w(x,t)+S(x)= w(x,t)=-\frac{32}{\pi^{4}} \sum_{n=1}^{\infty} \frac{(-1)^{n}}{(2n+1)^{4}}sin(\frac{(2n-1)\pi}{2}x)cos(\frac{(2n-1)\pi}{2}t)-\frac{x^{3}}{6}+\frac{x}{2}+1[/tex]
 
  • #48
wahoo2000 said:
Hmm.. no, actually i got ^3... I will show you how I did and then you can see where it went wrong. I don't write out the terms that went zero.
[tex]
A_{n}D_{n}=+\frac{4}{(2n+1)\pi}\int^{1}_{0}(\frac{3x^{2}}{6}-\frac{1}{2})(cos(\frac{(2n-1)\pi}{2}x)dx = [/tex]
[tex]=-\frac{16}{(2n+1)^{2}\pi{2}}\int^{1}_{0}(x*sin(\frac{(2n-1)\pi}{2}x)dx = \frac{32}{(2n+1)^{3}\pi{3}}\int^{1}_{0}(cos(\frac{(2n-1)\pi}{2}x)dx=[/tex]
[tex]=\frac{32}{(2n+1)^{3}\pi{3}}[(sin(\frac{(2n-1)\pi}{2}x)]=\frac{32}{(2n-1)^{3}\pi^{3}}*-(-1)^{n}[/tex]
edit: I think I lost one inner derivative there at the end...but it doesn't make sense to me anyway..in that case it would be 64/((2n+1)^4*pi^4)
:/


Hope we can find the mistake. I go on with "your" AnDn meanwhile.

Then it means that
[tex] w(x,t)=-\frac{32}{\pi^{4}} \sum_{n=1}^{\infty} \frac{(-1)^{n}}{(2n+1)^{4}}sin(\frac{(2n-1)\pi}{2}x)cos(\frac{(2n-1)\pi}{2}t)
[/tex] and
[tex]u(x,t)=w(x,t)+S(x)= w(x,t)=-\frac{32}{\pi^{4}} \sum_{n=1}^{\infty} \frac{(-1)^{n}}{(2n+1)^{4}}sin(\frac{(2n-1)\pi}{2}x)cos(\frac{(2n-1)\pi}{2}t)-\frac{x^{3}}{6}+\frac{x}{2}+1[/tex]

In your second line you multiplied by 4 instead of 2 and in your 3rd line you missed a factor of 2/(2n-1)pi also, they should all be (2n-1)'s not(2n+1)'s. This gives:

[tex]u(x,y)= -\frac{32}{\pi^{4}} \sum_{n=1}^{\infty} \frac{(-1)^{n}}{(2n-1)^{4}}sin(\frac{(2n-1)\pi}{2}x)cos(\frac{(2n-1)\pi}{2}t)-\frac{x^{3}}{6}+\frac{x}{2}+1[/tex]

or if you make the substitution n=m+1 you get:

[tex]u(x,y)= \frac{32}{\pi^{4}} \sum_{m=0}^{\infty} \frac{(-1)^{m}}{(2m+1)^{4}}sin(\frac{(2m+1)\pi}{2}x)cos(\frac{(2m+1)\pi}{2}t)-\frac{x^{3}}{6}+\frac{x}{2}+1[/tex]

Now, what is u(1,0) according to the boundary conditions? What is it according to the above formula?
 
  • #49
Ah, I see.. The missed inner derivative I found myself, but I squared once instead of multiplying by 2.. Time to go to bed soon :rolleyes:
[tex]
u(1,0)= -\frac{32}{\pi^{4}} \sum_{n=1}^{\infty} \frac{(-1)^{n}}{(2n-1)^{4}}*-(-1)^{n}*1-\frac{1}{6}+\frac{1}{2}+1
[/tex]
but I don't have u(1,0) in my BC?

EDIT: and (-1)^n*-(-1)^n = 1 ...
 
Last edited:
  • #50
wahoo2000 said:
Ah, I see.. The missed inner derivative I found myself, but I squared once instead of multiplying by 2.. Time to go to bed soon :rolleyes:
[tex]
u(1,0)= -\frac{32}{\pi^{4}} \sum_{n=1}^{\infty} \frac{(-1)^{n}}{(2n-1)^{4}}*-(-1)^{n}*1-\frac{1}{6}+\frac{1}{2}+1
[/tex]
but I don't have u(1,0) in my BC?

EDIT: and (-1)^n*-(-1)^n = 1 ...

You have u(x,0)=1 don't you? Why wouldn't u(1,0)=1 too?:wink:
 
  • #51
Aha.. but then
[tex] u(1,0)= -\frac{32}{\pi^{4}} \sum_{m=0}^{\infty} \frac{(-1)^{m}}{(2m+1)^{4}}*-(-1)^{m}*1-\frac{1}{6}+\frac{1}{2}+1=1
\Rightarrow \sum_{m=0}^{\infty} \frac{1}{(2m+1)^{4}}=(1-\frac{4}{3})*(-\frac{\pi^{2}}{32})=\frac{\pi^{4}}{96}

[/tex]
 
Last edited:
  • #52
Yep, and

[tex]\sum_{n=1}^{\infty} \frac{1}{(2n-1)^{4}} =\frac{\pi^{4}}{96}[/tex]

can be used to answer your question by substituting n=m+1 :

[tex]\Rightarrow \sum_{n=1}^{\infty} \frac{1}{(2n-1)^{4}} = \sum_{m+1=1}^{\infty} \frac{1}{(2(m+1)-1)^{4}}
= \sum_{m=0}^{\infty} \frac{1}{(2m+1)^{4}}=\frac{\pi^{4}}{96}[/tex]
 
  • #53
Thank you SO much for your patience. Tomorrow I am going through this problem from the start again to make sure I get each step.
 
  • #54
wahoo2000 said:
Thank you SO much for your patience. Tomorrow I am going through this problem from the start again to make sure I get each step.

Your welcome. The critical point of the problem was the part about the boundary conditions for w(x,t). Choosing the right boundary conditions (by choosing boundary conditions for S) makes the solution much easier. In general, you should try to set the conditions on S such that w(x,t) vanishes at one endpoint (x=0 in this case).
 

Similar threads

  • Calculus and Beyond Homework Help
Replies
11
Views
689
  • Calculus and Beyond Homework Help
Replies
2
Views
756
  • Calculus and Beyond Homework Help
Replies
21
Views
758
  • Calculus and Beyond Homework Help
Replies
1
Views
383
  • Calculus and Beyond Homework Help
Replies
7
Views
135
  • Calculus and Beyond Homework Help
Replies
14
Views
239
  • Calculus and Beyond Homework Help
Replies
7
Views
512
  • Calculus and Beyond Homework Help
Replies
5
Views
560
  • Calculus and Beyond Homework Help
Replies
4
Views
636
  • Calculus and Beyond Homework Help
Replies
1
Views
797
Back
Top