 104
 14
 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 guasseidel methods (or maybe a different method, not the point):
(1/s)⋅(∂/∂s)[(s/ρ)(∂ψ/∂s)] + (1/s^{2})⋅(∂/∂Φ)[(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/s_{j})⋅(1/Δs^{2})⋅[(s_{j+1}/ρ_{j+1,l}){ψ_{j+1,l}  ψ_{j,l}}  (s_{j1}/ρ_{j1,l}){ψ_{j,l}  ψ_{j1,l}}] + (1/s_{j}^{2})⋅(1/ΔΦ^{2})⋅[(1/ρ_{j,l+1})⋅{ψ_{j,l+1}  ψ_{j,l}}  (1/ρ_{j,l1})⋅{ψ_{j,l}  ψ_{j,l1}}]
Cell Centered Discretization (see cell centered grid attached):
(1/s_{j})⋅(1/Δs^{2})⋅[(s_{j+3/2}/ρ_{j+1,l}){ψ_{j+1,l}  ψ_{j,l}}  (s_{j1/2}/ρ_{j1,l}){ψ_{j,l}  ψ_{j1,l}}] + (1/s_{j+1/2}^{2})⋅(1/ΔΦ^{2})⋅[(1/ρ_{j+1,l})⋅{ψ_{j,l+1}  ψ_{j,l}}  (1/ρ_{j1,l})⋅{ψ_{j,l}  ψ_{j,l1}}]
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?
(1/s)⋅(∂/∂s)[(s/ρ)(∂ψ/∂s)] + (1/s^{2})⋅(∂/∂Φ)[(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/s_{j})⋅(1/Δs^{2})⋅[(s_{j+1}/ρ_{j+1,l}){ψ_{j+1,l}  ψ_{j,l}}  (s_{j1}/ρ_{j1,l}){ψ_{j,l}  ψ_{j1,l}}] + (1/s_{j}^{2})⋅(1/ΔΦ^{2})⋅[(1/ρ_{j,l+1})⋅{ψ_{j,l+1}  ψ_{j,l}}  (1/ρ_{j,l1})⋅{ψ_{j,l}  ψ_{j,l1}}]
Cell Centered Discretization (see cell centered grid attached):
(1/s_{j})⋅(1/Δs^{2})⋅[(s_{j+3/2}/ρ_{j+1,l}){ψ_{j+1,l}  ψ_{j,l}}  (s_{j1/2}/ρ_{j1,l}){ψ_{j,l}  ψ_{j1,l}}] + (1/s_{j+1/2}^{2})⋅(1/ΔΦ^{2})⋅[(1/ρ_{j+1,l})⋅{ψ_{j,l+1}  ψ_{j,l}}  (1/ρ_{j1,l})⋅{ψ_{j,l}  ψ_{j,l1}}]
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

8 KB Views: 7

8.8 KB Views: 8