[cho99]

Hello dear users,

I'm very new to OpenFOAM, so I apologize for asking trivial questions. I already searched the forum for my problem but I didn't find an appropriate solution.

I try to simulate *just* a 2D axisymmetric water jet for different (low) Re numbers, and the results should look like results.png (this image was made during the experiment), but for me, both pressure and velocity fields looks a bit weird even though the solution converges (p.png and U.png at t = 3 s)

What am I doing wrong?

1) Did I choose the solver which doesn't fit in the case? (I tried pisoFoam and pimpleFoam)
2) The wrong turbulence model? (I tried RAS kEpsilon, and LES SpalartAllmarasDDES, both solutions didn't look good)

3) Or are there errors in the boundary conditions? Which are (for the LES case):

inlet (only 1 face, one can see it in the attachec U.png file):
p = zeroGradient
U = fixedValue / unifrom (0 -1 0) (this is for Re = 1000)
k = fixedValue / unifrom 0.5 (calculated with k = 0.5 u²)
nut = calculated / uniform 0
nuTilda = fixedValue / unifrom 0

outlet (All other faces except of 2 faces belongs to frontAndBack group):
p = fixedValue / uniform 0;
U = type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 -1 0);
k = type inletOutlet;
inletValue uniform 0.5,
value uniform 0.5;
nut = calculated / uniform 0
nuTilda = type inletOutlet;
inletValue uniform 0;
value uniform 0;

frontAndBack:
p = empty
U = empty
k = empty
nut = empty
nuTilda = empty

internalField:
p = uniform 0
U = unifrom (0 0 0)
k = uniform 0.5
nut = uniform 0
nuTilda = uniform 0


*I used no walls*


I would really appreciate any help.

Greetings,
Cho

[cho99]

Problem solved by visualizing the vorticity. div(U)