- #1
- 1,440
- 7
Hi guys. I am solving the axisymmetric free jet of an incompressible fluid. But I have troubles at r=0. Continuty equation can be written in cylindrical coordinates as:
1/r*d(rv)/dr + du/dz=0
v=radial velocity (v=0 at r=0)
u=axial velocity.
hz=delta(z)
hr=delta(r)
What happens at r=0?. I have to obtain a finite difference scheme in order to integrate the problem, so I have created a grid that has "i" index for "z" coordinate and "j" index for "r" coordinate. j=1 corresponds to the symmetry axe.
for j>1 I have no problem because r>0. But in r=0 there is a singularity point. How can I deal with this?. As you can see I know v(i,1)=0
1/r*d(rv)/dr + du/dz=0
v=radial velocity (v=0 at r=0)
u=axial velocity.
hz=delta(z)
hr=delta(r)
What happens at r=0?. I have to obtain a finite difference scheme in order to integrate the problem, so I have created a grid that has "i" index for "z" coordinate and "j" index for "r" coordinate. j=1 corresponds to the symmetry axe.
for j>1 I have no problem because r>0. But in r=0 there is a singularity point. How can I deal with this?. As you can see I know v(i,1)=0