[marhamat]

Hi guys
I'm studing the effect of adding body force(F) to oodles solver.
Where F is exponential function:
F= A*exp(B*(x^2+y^2+z^2)
that i force it to my case geometry center for only 1 iteration . A & B are constants.
My study case is a channel with uniform inlet velocity(=10m/s).and the sides are symmetryplane.
......
so i changed oodles solver:

if (momentumPredictor)
{
solve(UEqn == -fvc::grad(p)+Fcoefficient*F);

}

Where Fcoefficient changes between 0&1 for iteration control.so it is 1 for only 1 iteration

You can see some figures that shows the vorticty variation with time changing in flow :
:
:
:
:
:
:
How i can explan the variation of vorticity distribution by time changing?
Can you explan me more detailes in this field ?

Have a nice weekend
Marhamat