How to Solve the 2D Laplace Equation with Fourier Series in a Grounded Square?

[humo90]

I don't know how to build up Fourier series for this problem. I tried several times but no results. Can anybody help me.

Consider the 2-d Laplace equation for the potential V(x,y) in the region inside
an empty square whose sides lie on the lines x = 0, x = a , y = 0, and y = a, respectively.
The sides at x = a , y = 0, and y = a are grounded, and
V(0,y) = C_{1}sin\frac{(\pi)y}{a} + C_{2}sin\frac{2(\pi)y}{a}
where C1 and C2 are constants

(a) Find V(x,y) everywhere inside the square.
Hint: use the trivial solution for X(x) with X(x)=0 at x=a

(b) Calculate the potential for C1=100V, C2=50V in the middle of the box. Do you find
any dependence on the box size?

 

[runnergirl]

Well I took a stab at it, and I think I have an idea where the solution should go. If you the standard Laplace's Equation: \nabla^2V=0. Where in this case we are dealing with cartesian coordinates.

I'm hoping you know the general solutions are going to be sines and cosines for the y-direction and sinh and cosh for the x-direction (if you don't know why I can explain).

So our V(x,y) = \sum_{i=1}^\infty A_n*sin\left(\frac{n\pi y}{a}\right)\left[cosh\left(\frac{n\pi x}{a}\right)+sinh\left(\frac{n\pi x}{a}\right)\right].

Now if you apply the boundary condition properly and use the orthogonality property, you should get the series in only two terms, when n = 1 and n = 2.

I hope I'm approaching this correctly, but let me know if it does make sense.

 

[vela]

Show us your work. Show/explain how you solved for X(x) and Y(y) before applying the boundary conditions and what applying the three V=0 boundary conditions got you. It doesn't make sense to try apply the last condition if you don't have the solution correct up to that point.

 

[runnergirl]

Sorry about that, I'm new to posting so I'm unsure how much I should post of a solution.

Ok so separation of variables (in cartesian coordinates):
\frac{\partial^2 V}{\partial x^2} + \frac{\partial^2 V}{\partial y^2} = 0
where V(x,y) = X(x)Y(y).

Substituting:
\frac{1}{X}\frac{d^2X}{dx^2}+\frac{1}{Y}\frac{d^2Y}{dy^2} = 0

If you set it equal to a constant say \alpha^2 then you're two differential equations will be:
\frac{d^2X}{dx^2} - \alpha X = 0
\frac{d^2Y}{dy^2} + \alpha Y = 0
And solutions:
X(x) = Acosh(\alpha x)+Bsinh(\alpha x)
Y(y) = A_1cos(\alpha y)+B_1sin(\alpha y)

Since Y(0)=Y(a)=0\implies A_1 = 0 and \alpha = \frac{n\pi}{a}

So since we're left with the sine term and with the boundary conditions given for X(x) doesn't allow us to get rid of any of those terms the full solution would be:
V(x,y) = \sum_{i=1}^\infty A_n*sin\left(\frac{n\pi y}{a}\right)\left[cosh\left(\frac{n\pi x}{a}\right)+sinh\left(\frac{n\pi x}{a}\right)\right]

So then when you apply the last one and use orthogonality that should give you only two terms n=1 and n=2. If you need me to explain how I deduced how I only got two terms I can.

 

[vela]

My post was meant for humo90.

You can't quite set A=B=1 yet.

One trick you can use in this type of problem when you have a boundary where the function vanishes is to write the solution in the form
X(x) = A\cosh\alpha(x-a)+B\sinh\alpha(x-a)You don't really need to do this. It just makes the solution look nicer.