Solving inhomogenous wave equation

  • Thread starter Thread starter wahoo2000
  • Start date Start date
  • Tags Tags
    Wave Wave equation
Click For Summary
SUMMARY

The discussion focuses on solving the inhomogeneous wave equation given by u_{tt} - u_{xx} = x with specific boundary conditions. The solution approach involves separating variables and determining the function S(x) that satisfies the equation. The participants identify the need to consider all cases for the eigenvalue λ, specifically λ < 0, λ = 0, and λ > 0, while correcting the initial assumptions regarding boundary conditions and the form of S(x). The final solution is contingent on correctly applying these principles to derive w(x,t) and ultimately u(x,t).

PREREQUISITES
  • Understanding of partial differential equations (PDEs)
  • Familiarity with boundary value problems
  • Knowledge of separation of variables technique
  • Concept of eigenvalues in differential equations
NEXT STEPS
  • Study the method of separation of variables in depth
  • Learn about eigenvalue problems in the context of PDEs
  • Explore the implications of different boundary conditions on solutions
  • Investigate the characteristics of inhomogeneous wave equations
USEFUL FOR

Mathematicians, physicists, and engineering students focusing on wave phenomena, particularly those dealing with partial differential equations and boundary value problems.

  • #31
wahoo2000 said:
Anyway.. if I keep the A
(ii) \Rightarrow w_{x}(1,t)=Akcos(kx)(Csin(kt)+Dcos(kt))=0 then cos(kx) should be zero because k and (Csin(kt)+Dcos(kt) are not zero, (and neither is A..), because if A=0 we have only trivial solution X(x)T(t)=0, and same if (Csin(kt)+Dcos(kt)=0, then per definition k is not zero since lambda<0.

Right, so what does cos(kx)=0 tell you about k?
 
Physics news on Phys.org
  • #32
that k=(n*Pi)/2...
so \lambda=-\frac{n^{2}\pi^{2}}{4}
 
Last edited:
  • #33
wahoo2000 said:
that k=(n*Pi)/2...
so \lambda=-\frac{n^{2}\pi^{2}}{4}

Close, but only true if n is odd. The best way to represent that is :

(ii) \Rightarrow k= \frac{(2n-1) \pi}{2} \quad \ni n=1,2,3 \ldots

\Rightarrow w(x,t)= \sum_{n=1}^{\infty} A_nsin \left( \frac{(2n-1) \pi}{2} x \right)(C_nsin \left(\frac{(2n-1) \pi}{2} t \right)+D_n cos \left( \frac{(2n-1) \pi}{2} t \right) )

And (iv) still implies that C=0,

\Rightarrow w(x,t)= \sum_{n=1}^{\infty} A_n D_n sin \left( \frac{(2n-1) \pi}{2} x \right) cos \left( \frac{(2n-1) \pi}{2} t \right)

So, what do you get from (iii) now?
 
  • #34
Well... that means that

<br /> \Rightarrow w(x,0)= \sum_{n=1}^{\infty} A_n D_n sin \left( \frac{(2n-1) \pi}{2} x \right) cos \left( 0\right) \Rightarrow \sum_{n=1}^{\infty} A_n D_n sin \left( \frac{(2n-1) \pi}{2} x \right)=\frac{x^{3}}{6}-\frac{x}{2}<br />
 
  • #35
wahoo2000 said:
Well... that means that

<br /> \Rightarrow w(x,0)= \sum_{n=1}^{\infty} A_n D_n sin \left( \frac{(2n-1) \pi}{2} x \right) cos \left( 0\right) \Rightarrow \sum_{n=1}^{\infty} A_n D_n sin \left( \frac{(2n-1) \pi}{2} x \right)=\frac{x^{3}}{6}-\frac{x}{2}<br />

Yes, and from this you can use Fourier's method to determine A_nD_n. Give it a shot; what do you get?
 
  • #36
Hmm, does that mean we are looking for Fourier sinus series of f(x)=x^3/6-x/2?
Is it A_{n}D_{n}=2\int^{1}_{0}(\frac{x^{3}}{6}-\frac{x}{2})(\frac{(2n-1)\pi}{2}x)dx
 
  • #37
wahoo2000 said:
Hmm, does that mean we are looking for Fourier sinus series of f(x)=x^3/6-x/2?
Is it A_{n}D_{n}=2\int^{1}_{0}(\frac{x^{3}}{6}-\frac{x}{2})(\frac{(2n-1)\pi}{2}x)dx

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

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

<br /> A_{n}D_{n}=2\int^{1}_{0}(\frac{x^{3}}{6}-\frac{x}{2})sin(\frac{(2n-1)\pi}{2}x)dx<br />
, 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:

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}}
 
  • #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:

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}}

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

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 2[(\frac{x^{3}}{6}-\frac{x^{2}}{2})(cos(\frac{(2n-1)\pi}{2}x)]^{1}_{0}<br /> <br /> +\frac{4}{(2n+1)\pi}\int^{1}_{0}(\frac{3x^{2}}{6}-\frac{2x}{2})(cos(\frac{(2n-1)\pi}{2}x)dx

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...

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

Do you mean:
<br /> A_{n}D_{n}=\frac{32}{(2n-1)^{4}\pi^{4}}*-(-1)^{n}<br /> ?
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.
<br /> 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 =
=-\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=
=\frac{32}{(2n+1)^{3}\pi{3}}[(sin(\frac{(2n-1)\pi}{2}x)]=\frac{32}{(2n-1)^{3}\pi^{3}}*-(-1)^{n}
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
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)<br /> and
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
 
  • #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.
<br /> 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 =
=-\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=
=\frac{32}{(2n+1)^{3}\pi{3}}[(sin(\frac{(2n-1)\pi}{2}x)]=\frac{32}{(2n-1)^{3}\pi^{3}}*-(-1)^{n}
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
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)<br /> and
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

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:

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

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

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

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:
<br /> 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<br />
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:
<br /> 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<br />
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
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<br /> \Rightarrow \sum_{m=0}^{\infty} \frac{1}{(2m+1)^{4}}=(1-\frac{4}{3})*(-\frac{\pi^{2}}{32})=\frac{\pi^{4}}{96}<br /> <br />
 
Last edited:
  • #52
Yep, and

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

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

\Rightarrow \sum_{n=1}^{\infty} \frac{1}{(2n-1)^{4}} = \sum_{m+1=1}^{\infty} \frac{1}{(2(m+1)-1)^{4}}<br /> = \sum_{m=0}^{\infty} \frac{1}{(2m+1)^{4}}=\frac{\pi^{4}}{96}
 
  • #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

  • · Replies 11 ·
Replies
11
Views
3K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 21 ·
Replies
21
Views
2K
Replies
4
Views
2K
  • · Replies 105 ·
4
Replies
105
Views
7K
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K