View Full Version : Solving PDEs using Fouries Series ???
Makveger
Mar3-11, 04:46 PM
Hello
I am trying to solve 2D Laplace's equation (\nabla2u) 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,....) ??
JJacquelin
Mar4-11, 08:18 AM
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)= ?
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.
JJacquelin
Mar6-11, 10:34 AM
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.
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).
JJacquelin
Mar6-11, 11:05 AM
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 completly fulfill the boundary conditions.
So, the answer of Makveger is a key point.
Yes, JJacquelin. You are right, I totally agree with you!
lapace or separation of variables works here.. either way is good.
vBulletin® v3.8.7, Copyright ©2000-2012, vBulletin Solutions, Inc.