A Implementing a weird-looking boundary condition (PDE/FDM)

204
10
So I have this problem, taken from Kraus's heat transfer book.

63ye6d.png


So deriving the computational molecule, the conditions for (3.251a), (3.251b) is a bit of a no brainer. The issue I am having is about the boundaries for (3.251c) and (3.251d). This is actually the first time I have seen this kind of boundary condition.

How do I deal with this? My hypothesis is to integrate the equations (lol) but the constants of integration stay around... I have no idea what to do at all. Bi and γ are constants. Thanks!
 
204
10
Oh, to add.

There's an analytical expression for the solution; and I am able to graph the resulting multivariate graph. I intend to do the numerical analysis as another way of solving it. The problem is, even if I implement the 'integrating the boundary condition' part or if I use the ghost point strategy, both of it does not work - the graph from the numerical solution is way too far from the analytical expression's graph.
 
19,464
3,894
You're solving this by finite difference, correct?
 
204
10
Yes sir. Finite difference.
 
19,464
3,894
$$\theta (1+\Delta R, Z)-\theta (1-\Delta R, Z)=-2(\Delta R) Bi\ \theta (1,Z)$$
so
$$\frac{\partial ^2 \theta}{\partial R^2}=2\frac{(\theta (1-\Delta R, Z)-\theta (1,Z))-(\Delta R) Bi\ \theta (1,Z)}{(\Delta R)^2}$$
or
$$\frac{\partial ^2 \theta}{\partial R^2}=2\frac{(\theta (1-\Delta R, Z)-\theta (1,Z))}{(\Delta R)^2}-2\frac{Bi\ \theta(1,Z)}{\Delta R}$$
or
$$\frac{\partial ^2 \theta}{\partial R^2}=2\frac{(\theta (1-\Delta R, Z)-[1+(\Delta R) Bi]\theta (1,Z))}{(\Delta R)^2}$$
 
204
10
Thanks for replying sir. Actually I did it already, but I seem to be getting erroneous results. Could you have my computational molecules for each boundary checked first?

I intend to implement a solution similar to Gauss-Seidel iteration; thus I kept on factoring out the 'center' molecule.
 
Last edited:
19,464
3,894
Sorry. I'll help with the finite differencing, but, as far as the debugging is concerned, you're on your own.
 
204
10
Oh, it's alright sir.

Apparently my computational molecules for the boundaries are correct; but the numerical solution is still given different values. I'm starting to think that the problem comes from the corners - problem is what should I do with the corners (0,1) and (1,1). What boundary should I invoke? Say for the upper left corner (0,1); is it 3.251b or 3.251d?
 
19,464
3,894
Oh, it's alright sir.

Apparently my computational molecules for the boundaries are correct; but the numerical solution is still given different values. I'm starting to think that the problem comes from the corners - problem is what should I do with the corners (0,1) and (1,1). What boundary should I invoke? Say for the upper left corner (0,1); is it 3.251b or 3.251d?
Both
 

Related Threads for: Implementing a weird-looking boundary condition (PDE/FDM)

  • Last Post
Replies
2
Views
3K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
2
Views
615
Replies
3
Views
545
Replies
2
Views
2K
Replies
6
Views
2K
Replies
1
Views
634
Replies
1
Views
2K

Hot Threads

Top