# Continuity equation question

1. Jul 14, 2004

### Clausius2

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

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 simmetry 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

2. Jul 23, 2004

### mika

1/r*d(vr)/dr+du/dz=0 =>

1/r*(v+r*dv/dr)+du/dz = v/r + dv/dr +du/dz

at r=0, v=0 ->

dv/dr + du/dz = 0

does this help?

Last edited: Jul 23, 2004
3. Jul 23, 2004