Register to reply

Solving inhomogenous wave equation

by wahoo2000
Tags: equation, inhomogenous, solving, wave
Share this thread:
gabbagabbahey
#37
Sep21-08, 01:01 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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.
wahoo2000
#38
Sep21-08, 01:21 PM
P: 37
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 :)
gabbagabbahey
#39
Sep21-08, 01:24 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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.
wahoo2000
#40
Sep21-08, 02:08 PM
P: 37
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]
gabbagabbahey
#41
Sep21-08, 02:40 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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.
wahoo2000
#42
Sep21-08, 03:43 PM
P: 37
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.. :( :(
gabbagabbahey
#43
Sep21-08, 03:49 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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.
wahoo2000
#44
Sep21-08, 04:20 PM
P: 37
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]
wahoo2000
#45
Sep21-08, 04:49 PM
P: 37
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]
gabbagabbahey
#46
Sep21-08, 05:12 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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?
wahoo2000
#47
Sep21-08, 05:40 PM
P: 37
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]
gabbagabbahey
#48
Sep21-08, 06:00 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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?
wahoo2000
#49
Sep21-08, 06:37 PM
P: 37
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
[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 .....
gabbagabbahey
#50
Sep21-08, 06:51 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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
[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?
wahoo2000
#51
Sep21-08, 07:07 PM
P: 37
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]
gabbagabbahey
#52
Sep21-08, 07:11 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
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]
wahoo2000
#53
Sep21-08, 07:14 PM
P: 37
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.
gabbagabbahey
#54
Sep21-08, 07:25 PM
HW Helper
gabbagabbahey's Avatar
P: 5,003
Quote Quote by wahoo2000 View Post
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).


Register to reply

Related Discussions
Inhomogenous electrodynamics wave equation Classical Physics 3
An inhomogenous equation-help ! Calculus & Beyond Homework 17
Solving Wave Equation / Imaginary Numbers Calculus & Beyond Homework 1