Numerical PDE's II - Circular Domain

Nusc
Messages
752
Reaction score
2

Homework Statement


Approximate the solution of:
<br /> <br /> \frac{\partial^{2}u }{{\partial r^{2} }} + \frac{1}{r}\frac{\partial u}{\partial r} + \frac{1}{r^{2}}\frac{\partial^{2}u}{\partial \theta^{2}} = 0, 1&lt;r&lt;3, 0&lt;\theta&lt;\pi,<br /> <br /> u(1,\theta) = 1-\cos(\theta),u(2,\theta)=u(r,0)=0 u(r,\pi)=4-2r<br /> <br />

Homework Equations


<br /> <br /> \delta r = \frac{r_{outer} - r_{inner}}{N}<br /> \delta \theta = \frac{\theta_{end} - \theta_{start}}{M}<br /> <br />
n,m are positve integers

<br /> \delta r_i = r_{inner} + i\delta r<br />
<br /> \delta theta_j = \theta_{start} + j\delta \theta<br />

The Attempt at a Solution


<br /> \frac{u_{i+1,j} - 2U_{i,j} + U_{i-1,j}}{(\delta r)^2} + \frac{1}{r_i}\frac{u_{i+1,j}-u_{i-1,j}}{2\delta r} +\frac{u_{i,j+1} - 2U_{i,j} + U_{i,j-1}}{(\delta \theta)^2} = f_{i,j}<br />

Multiply this equation by (\delta r)^2<br />
<br /> u_{i+1,j} - 2u_{i,j} + u_{i-1,j} + \frac{1}{r_i}\frac{u_{i+1,j}-u_{i-1,j}}{2}\delta r + \frac{1}{r_i}\frac{\delta r}{\delta \theta}(u_{i+1,j}-2u{i,j}+u{i,j-1} = f_{i,j}(\delta r)^2<br /> <br />

What do I do next>
 
Last edited:
Physics news on Phys.org
And the question is?
 
Nevermind, I got it.

But does anyone know if such a problem with those given BC's can be applied elsewhere?
When would you model something with a circular annulus?
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top