Approximation of boundary value problem using finite differences

  • Thread starter ataylor
  • Start date
  • #1
1
0

Homework Statement


A hot fluid is flowing through a thick-walled cylindrical metal tube at a constant temperature of 450C. The cylinder wall has an inner radius of 1 cm and an outer radius of 2 cm and the surrounding temperature is 20C. The temperature distribution u(r) in the metal is defined by the differential equation:

[tex]r\frac{d^2 u}{dr^2} + \frac{du}{dr} = 0 [/tex]

With the boundary conditions:

[tex]u\left(1\right) = 450 [/tex]

[tex]\frac{du}{dr}\left(2\right) = -K\left(u-20\right)[/tex]

where K=1.

Approximate the boundary value problem by creating a matrix problem using the finite difference method.

Homework Equations





The Attempt at a Solution


First, the interval [tex]1\leq r \leq 2[/tex] is divided by the number of iterations, N, to give the step length, h.

The derivatives of all the values inbetween the boundaries can be approximated by using finite differences as follows:

[tex]r\frac{u_{N-1}-2u_{N}+u_{N+1}}{h^2} + \frac{u_{N+1}-u_{N-1}}{2h} = 0[/tex]

Which after rearrangement gives:

[tex]r\frac{u_{N-1}-2u_{N}+u_{N+1}}{h^2} = -\frac{u_{N+1}-u_{N-1}}{2h}[/tex]

This can be represented as a linear matrix problem Ax = B by using the left-hand coefficients to create a tridiagonal N-by-N matrix (A), the coefficients being:

[tex]1\frac{r_{N}}{h^2}, -2\frac{r_{N}}{h^2}, 1\frac{r_{N}}{h^2}[/tex]


and by using the right-hand values to create a n-by-1 column matrix (B). This problem should be solveable using gaussian elimination to give the matrix X, which contains the values of u at all N iteration points.

The first and last value of the B matrix also have to be corrected due to the boundary values, which themselves are outside the range of the matrix.

When solving similar problems in the past, the boundary conditions have always been given as two function values at different points i.e u(1) and u(2) and usually there has been some way of calculating values for the B matrix. The problem I'm having here is finding a way to use the second boundary value, since it is given as -(u-20), with u presumably being u(2) which is unknown and finding a way to generate the matrix B, i.e the values of du/dr.
 

Answers and Replies

Related Threads on Approximation of boundary value problem using finite differences

Replies
8
Views
663
Replies
4
Views
17K
  • Last Post
Replies
4
Views
2K
  • Last Post
Replies
5
Views
3K
  • Last Post
Replies
5
Views
1K
  • Last Post
Replies
10
Views
2K
  • Last Post
Replies
3
Views
1K
Replies
4
Views
820
  • Last Post
Replies
6
Views
5K
  • Last Post
Replies
2
Views
672
Top