[FMDenaro]

Originally Posted by IbizaFR

Thank you FMDenaro.

So for all cells I check the max value of divergence and if not satisfied according to the maximum permitted value then I continue to new iteration.
Is it right?

do not forget that you have a system of equations therefore you have to satisfy continuity and momentum equation simultaneously

[IbizaFR]

Thank you FMDenaro.
What is the checking condition for the momentum?

[FMDenaro]

The residuals of the two momentum equations

[IbizaFR]

Do you mean the absolute change in u,v from one itetation to the next one?

[FMDenaro]

Originally Posted by IbizaFR

Do you mean the absolute change in u,v from one itetation to the next one?

I mean the residuals of the equations. The difference between two iterations is not the residual.

[IbizaFR]

Thank you.
The check has to be done in the provisional velocities or to the corrected in SIMPLE?

[IbizaFR]

Something else.
Do I have to use periodic boundary conditions mandatorily?

[FMDenaro]

Originally Posted by IbizaFR

Something else.
Do I have to use periodic boundary conditions mandatorily?

It would be better

[IbizaFR]

Can I avoid them with other way?

[FMDenaro]

Originally Posted by IbizaFR

Can I avoid them with other way?

The theoretical solution requires two plates estending infinitely in both directions.

Using fictious inlet/outlet BC.s requires that you prescribe the exact solution at inflow and a developed condition at outlet. I never tried such BC.s but I suppose they will drive to the exact solution at the steady state.

[IbizaFR]

So at the inflow and outflow I can give linear distribution of u and p=0 at both ends. Is it right?

[IbizaFR]

Thank you FMDenaro,

How can I implement periodic boundary conditions in this problem? Can you suggest me a book to study?

[IbizaFR]

Thank you FMDenaro,

Because I have not dealt with periodic boundary conditions before can you suggest me a book to study about this matter?

[FMDenaro]

Originally Posted by IbizaFR

Thank you FMDenaro,

Because I have not dealt with periodic boundary conditions before can you suggest me a book to study about this matter?

You cannot find a complete book about the periodic boundary conditions... it is a numerical issue that is simply implemented in a code by using a link between boundary nodes.
I give you an example in 1D for a periodicity over a lenght L. You define the nodes 1,2...N in such a way that x(i)=(i-1)*h+Xs. This way x(N)=Xs+L.
Now, you have that nodes 1 and N are linked by periodicity. You compute the equation either in node 1 or N and prescribe the computed value at the node you do not compute.

[IbizaFR]

Thank you,

Let me ask you:
What is Xs?
Do I have to alter all the momentum and pressure equations in SIMPLE at two ends or only the u,v,p at these ends?

[FMDenaro]

Originally Posted by IbizaFR

Thank you,

Let me ask you:
What is Xs?
Do I have to alter all the momentum and pressure equations in SIMPLE at two ends or only the u,v,p at these ends?

Xs is only the general value of the position of the left point, it could be also 0. The periodicity involves all the variables.

[IbizaFR]

Originally Posted by FMDenaro

Xs is only the general value of the position of the left point, it could be also 0. The periodicity involves all the variables.

I cann't understand if I write x(N)=Xs+L with Xs=0 so X(N)=L, how the code will understand that the start and finish node are the same?

[FMDenaro]

Originally Posted by IbizaFR

I cann't understand if I write x(N)=Xs+L with Xs=0 so X(N)=L, how the code will understand that the start and finish node are the same?

Solve all the equations at i=1 and set the obtained values at i=N

[IbizaFR]

So all the loops will not calculate anything at the i=N and set the calculated values at i=1 to i=N. Is it right this?

Thank you

[FMDenaro]

Originally Posted by IbizaFR

So all the loops will not calculate anything at the i=N and set the calculated values at i=1 to i=N. Is it right this?

Thank you

yes, correct