How do I solve the Laplace equation using Robbin's Boundary Conditions?

  • #1
shreddinglicks
216
6
Homework Statement
Equation

(d^2 u)/(dx^2 )+(d^2 u)/(dy^2 )=0

boundary conditions

u(0,y)=0
u(a,y)=50
u(x,0)=0
u^' (x,b)=-h[u(x,b)-70]
Relevant Equations
(d^2 u)/(dx^2 )+(d^2 u)/(dy^2 )=0
I've tried a few things. I did one method to try to accomplish the removal of the -70 in the derivative boundary condition. It came out as below. When plotting it however it gave a solution that didn't make sense.

1730136669757.png


1730136684280.png
 
Last edited:
Physics news on Phys.org
  • #2
What is the boundary condition at [itex]x = a[/itex]? What do you mean by the notation [itex]u'[/itex], since [itex]u[/itex] is a function of 2 varaibles? Do you mean [itex]\frac{\partial u}{\partial x}[/itex]? Is [itex]h[/itex] a constant?

If that is correct, then I would set [tex]
u(x,y) = \frac{70h}{1 + ah}x + v(x,y)[/tex] where [itex]v[/itex] satisfies [tex]\begin{split}
v(0,y) = 0, \\
\frac{\partial v}{\partial x}(a,y) + hv(a,y) &= 0, \\
v(x,0) &= -\frac{70h}{1 + ah}x, \\
v(x,b) &= 50 - \frac{70h}{1 + ah}x. \end{split}[/tex]
 
  • #3
pasmith said:
What is the boundary condition at [itex]x = a[/itex]? What do you mean by the notation [itex]u'[/itex], since [itex]u[/itex] is a function of 2 varaibles? Do you mean [itex]\frac{\partial u}{\partial x}[/itex]? Is [itex]h[/itex] a constant?

If that is correct, then I would set [tex]
u(x,y) = \frac{70h}{1 + ah}x + v(x,y)[/tex] where [itex]v[/itex] satisfies [tex]\begin{split}
v(0,y) = 0, \\
\frac{\partial v}{\partial x}(a,y) + hv(a,y) &= 0, \\
v(x,0) &= -\frac{70h}{1 + ah}x, \\
v(x,b) &= 50 - \frac{70h}{1 + ah}x. \end{split}[/tex]
Yes, exactly as you said. du/dx and h is a constant.

I realized I made an error.

The conditions are actually

1730136193713.png

I'll fix it also on the original post.

It's what you have but y instead of x.

1730138509262.png

1730137038134.png


The problem is, this gives an answer where for all y when u(0,y)=0 when x should be zero, x does not equal zero.

The funny thing is if I manipulate the equation by multiplying the first term by x/a. I can plot a solution that makes sense.

1730137022111.png
 

Attachments

  • 1730136289863.png
    1730136289863.png
    2.4 KB · Views: 1
  • 1730136397561.png
    1730136397561.png
    4.8 KB · Views: 1
  • 1730136876131.png
    1730136876131.png
    2 KB · Views: 3
Last edited:
Back
Top