Solving PDEs using Fouries Series ?

In summary, Solving PDEs using Fouries Series can be difficult because the solution must include both sinusoidal and Heaviside functions.
  • #1
Makveger
17
0
Solving PDEs using Fouries Series ?

Hello
I am trying to solve 2D Laplace's equation ([tex]\nabla[/tex]2u) using Fourier series using these boundary conditions for a square domain of length L:
u(x, 0) = 0
u(0, y) = 0
u(L, y) = 0
u(x,L) = Uo

After solving the 2 ODEs(separating variables method) the solution is in this form (using the boundary conditions except the last one):
u(x,y) = E*Sin(nπx/L)*Sinh(nπy/L)

And I'm stuck here,after using the 4th boundary conditions to convert the above function into a series
How can I convert this formula into a series??
How can the solution be the sum of all the values of u for all the values of n(Where n =1,2,3,...) ??
 
Physics news on Phys.org
  • #2


u(0, y) = 0 and u(x,L) = Uo
Then, can you clearly state u(0,L)= ?

u(L, y) = 0 and u(x,L) = Uo
Then, can you clearly state u(L,L)= ?
 
  • #3


I think the sinuous functions in x direction are the basis you choose. However, for the y direction, it must be both sinuous and co-sinuous functions. When you get the matrix representation in y direction, I think we need to impost the boundary condition on both sides. The simple way would be change the first row and the last row by the boundary condition in the so called physical space.
 
  • #4


Roughly I agree with the proposal of Kuan. But I am afraid that it will be not so simple.
In fact, around the points of (0 , L) and (L , L), a special beaviour occurs. That is why I asked Makvenger to clarify the bounding conditions on these particular points. Up to now, Makvenger gave no answer. So I think that it would be premature to go further in the mathematical development.
Without more relevant information, the solution have to include not only sinusoidal fuctions, but also Heaviside functions.
 
  • #5


JJacquelin, I see what you mean! For me, I simply ignored these 2 points at (0 , L) and (L , L). No matter what conditions they are, the Fourier series is going to diverge pointwisely at at (0 , L) and (L , L).
 
  • #6


OK. Kuan, that's right.
But, is it allowed to ignore what appened close to whese points ?
If this is of no importance regarding to the physical problem, OK, the solution expessed only in terms of sinusoidal functions will be sufficient.
If not, the solution has to be more complicated, in order to completely fulfill the boundary conditions.
So, the answer of Makveger is a key point.
 
  • #7


Yes, JJacquelin. You are right, I totally agree with you!
 
  • #8


lapace or separation of variables works here.. either way is good.
 
  • #9
Hey can somebody tell me the Fourier transform of e^-(pi*t^2)?
 
  • #11
Thanks for your help. But can you please help me out with the mathematics. I have tried it but got struck in the integral.
 
  • #12
Makveger said:
Hello
I am trying to solve 2D Laplace's equation ([tex]\nabla[/tex]2u) using Fourier series using these boundary conditions for a square domain of length L:
u(x, 0) = 0
u(0, y) = 0
u(L, y) = 0
u(x,L) = Uo

After solving the 2 ODEs(separating variables method) the solution is in this form (using the boundary conditions except the last one):
u(x,y) = E*Sin(nπx/L)*Sinh(nπy/L)

And I'm stuck here,after using the 4th boundary conditions to convert the above function into a series
How can I convert this formula into a series??
How can the solution be the sum of all the values of u for all the values of n(Where n =1,2,3,...) ??

You haven't written your ##u(x,y)## correctly. You have ##X_n(x)=\sin\frac{n\pi x}{L}## and ##Y_n(y) =\sinh\frac {n\pi y}{L}##. So your prospective solution is$$
u(x,y) = \sum_{n=1}^\infty c_nX_n(x)Y_n(y) = \sum_{n=1}^\infty c_n\sin\frac{n\pi x}{L}
\sinh\frac {n\pi y}{L}$$Now apply your fourth BC ##u(x,L)=U_0##:$$
U_0 =\sum_{n=1}^\infty c_n\sin\frac{n\pi x}{L}
\sinh\frac {n\pi L}{L}$$Now use your half range sine expansion to get the constants. Notice that what you usually call ##a_n## in a sine expansion is ##a_n = c_n\sinh(n\pi)## in this problem.
 
  • #13
JJacquelin said:

For anyone using wolframalpha or mathematica to compute Fourier transforms, be aware that the default convention is to use the unitary-transformation convention: factors of ##1/\sqrt{2\pi}## outside both the transform and the inverse. Typically in physics the factors are 1 in front of the transform and ##1/2\pi## in front of the inverse.

(Overall the difference doesn't matter as long as you use a consistent convention, or unless you want to compare to a particular solution which uses a particular convention).
 

1. What is a Fourier series and how is it used to solve PDEs?

A Fourier series is a mathematical tool used to represent a periodic function as a sum of sines and cosines. It can be used to solve partial differential equations (PDEs) by transforming a PDE into a simpler problem involving ordinary differential equations (ODEs) that can then be solved using Fourier series techniques.

2. What are the advantages of using Fourier series to solve PDEs?

One advantage of using Fourier series to solve PDEs is that it allows for the solution to be expressed in terms of a series of simple functions, rather than a complex equation. This can make the solution easier to work with and interpret. Additionally, Fourier series can be used to solve a wide range of PDEs, making it a versatile tool for scientists and engineers.

3. Can Fourier series be used to solve PDEs with non-periodic boundary conditions?

Yes, Fourier series can still be used to solve PDEs with non-periodic boundary conditions by using a technique called Fourier transforms. This involves extending the domain of the function to infinity and then taking the inverse Fourier transform to obtain the solution in the original domain.

4. Are there any limitations to using Fourier series to solve PDEs?

One limitation of Fourier series is that it can only be applied to linear PDEs. Additionally, some PDEs may have solutions that are not well-represented by Fourier series, particularly if the solution has sharp discontinuities or singularities.

5. How can I determine if a PDE can be solved using Fourier series?

A PDE can be solved using Fourier series if it is linear and has constant coefficients. Additionally, the boundary conditions must be periodic or can be made periodic through a change of variables. If these conditions are met, then Fourier series is a viable method for solving the PDE. Otherwise, alternative methods may need to be used.

Similar threads

Replies
4
Views
1K
  • Differential Equations
Replies
2
Views
2K
  • Differential Equations
Replies
1
Views
2K
  • Differential Equations
Replies
2
Views
2K
  • Differential Equations
Replies
3
Views
2K
  • Differential Equations
Replies
7
Views
332
  • Differential Equations
Replies
2
Views
3K
  • Differential Equations
Replies
22
Views
2K
  • Differential Equations
Replies
4
Views
2K
Replies
1
Views
1K
Back
Top