I Help discretizing this PDE (stream function)

Summary
I have attempted two different methods of discretization for a PDE on a 2D annulus grid in cylindrical coordinates. I would appreciate anyone pointing out problems or suggesting better methods. This is part of research I am conducting as an undergrad.
I have a PDE that I want to solve for a stream function ψ(j,l) by discretizing it on a 2D annulus grid in cylindrical coordinates, then solving with guas-seidel methods (or maybe a different method, not the point):

(1/s)⋅(∂/∂s)[(s/ρ)(∂ψ/∂s)] + (1/s2)⋅(∂/∂Φ)[(1/ρ)(∂ψ/∂Φ)]

Where s and Φ are cylindrical coordinates, ρ is mass density (known through other means)

I've attempted to discretize this equation in a cell centered manner (which is what my professor suggested) but I also discretized it with ψ defined at the vertexes because it is easier for me to think about (I am fairly new to this type of numerical equation solving).

Vertex Centered Discretization (please see vertex centered grid attached):

(1/sj)⋅(1/Δs2)⋅[(sj+1j+1,l){ψj+1,l - ψj,l} - (sj-1j-1,l){ψj,l - ψj-1,l}] + (1/sj2)⋅(1/ΔΦ2)⋅[(1/ρj,l+1)⋅{ψj,l+1 - ψj,l} - (1/ρj,l-1)⋅{ψj,l - ψj,l-1}]

Cell Centered Discretization (see cell centered grid attached):

(1/sj)⋅(1/Δs2)⋅[(sj+3/2j+1,l){ψj+1,l - ψj,l} - (sj-1/2j-1,l){ψj,l - ψj-1,l}] + (1/sj+1/22)⋅(1/ΔΦ2)⋅[(1/ρj+1,l)⋅{ψj,l+1 - ψj,l} - (1/ρj-1,l)⋅{ψj,l - ψj,l-1}]

In each case I am attempting to approximate the outside differential terms by finding difference in the slope of ψ at either end of Δs for each iteration.

Does this make sense to anyone familiar with this type of numerical problem solving? Are there any obvious flaws, instabilities or oversights I've missed?
 

Attachments

Oh, I should probably mention that I am only giving the terms which I need help discretizing. The rest of the equation is straightforward to discretize and is mostly constant terms.
 
19,178
3,787
Shouldn't the mass density be constant? If not, please describe the exact problem being solved.
 
No, the density depends on s and Φ, the entire equation is:

(1/s)⋅(∂/∂s)[(s/ρ)(∂ψ/∂s)] + (1/s2)⋅(∂/∂Φ)[(1/ρ)(∂ψ/∂Φ)] - 2Ω + ρ(c0 + c1ψ) = 0

Where Ω is frame rotation rate and c0, c1 are arbitrary constants (from the first and second term of a taylor expansion).
 
19,178
3,787
$$\frac{1}{s}\frac{\partial}{\partial s}\left(\frac{s}{\rho}\frac{\partial \psi}{\partial s}\right)=\frac{1}{s_i}\frac{\left[\frac{s_{i+1/2}}{\rho_{i+1/2}}(\psi_{i+1,j}-\psi_{i,j})-\frac{s_{i-1/2}}{\rho_{i-1/2}}(\psi_{i,j}-\psi_{i-1,j})\right]}{(\Delta s)^2}$$
 

Want to reply to this thread?

"Help discretizing this PDE (stream function)" You must log in or register to reply here.

Related Threads for: Help discretizing this PDE (stream function)

Replies
3
Views
514
Replies
12
Views
3K
Replies
2
Views
756
  • Posted
Replies
7
Views
1K
  • Posted
Replies
3
Views
2K

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving

Hot Threads

Top