Here is what I did.

1. Open FEATool Multiphysics and select 2D Navier Stokes

2. Go to Geometry and create 0.1 x 0.5 (I tried to use SI units. So these are 0.1 m diamter and 0.5 m length.

3. Create 0.001 Meshes both for Subdomain and Boundary Grids

4. Go to Equation and set density to 1,000 and mu to 9.47E-05 Pa s. These are Si units.

5. Go to Boundary and set inlet velocity, u0, to 2 for the left side (Domain 4), set Outflow/pressure and p0=0 for the right side (Domain 2), and set Wall/no-slip for the top and the bottom (Regions 1 and 3).

6. Go to Solve and run by hitting “=”.

Result 1.png is what I got. The velocity is obviously too high and it is a wrong result. Result 2.png is the result I got when I ran analysis with density of 1, not 1,000. Result 2 is the result what I am looking for. The difference between result 1 and result 2 is created by the difference in densities. When I use 1,000 for density to meet the SI unit system, the answer is obvisouly wrong. However, when I use 1 for density the result seems right to me. So I suspect there is something wrong the units for density I used.

Result1.pngResult2.png