[IbizaFR]

Hello to everyone,

I am writing SIMPLE code for couette Laminar flow with the existence of pressure.

I have read similar topics for this problem on this forum but it is not clear yet.

My question is about:

1) How can I apply the given pressure gradient in the SIMPLE algorithm? Maybe in the provisional velocities equations u*,v* can I replace the term p(i+1,j)-p(i,j) with the given pressure gradient?

2) Can you suggest me the correct boundary conditions at inlet, outlet, walls for the horizontal velocity, vertical velocity and pressure?

Thank you.

[AliE]

Originally Posted by IbizaFR

Hello to everyone,

I am writing SIMPLE code for couette Laminar flow with the existence of pressure.

I have read similar topics for this problem on this forum but it is not clear yet.

My question is about:

1) How can I apply the given pressure gradient in the SIMPLE algorithm? Maybe in the provisional velocities equations u*,v* can I replace the term p(i+1,j)-p(i,j) with the given pressure gradient?

2) Can you suggest me the correct boundary conditions at inlet, outlet, walls for the horizontal velocity, vertical velocity and pressure?

Thank you.

Hello,

You have to employ periodic bcs at inlet/outlet and add the pressure gradient as body force. In practice, you have to add to your momentum equation a term like dpdx*vol, where vol is the cell volume.

For the pressure equation, set periodic bcs for inlet/outlet, ehile simple zero gradient for solid walls.

[IbizaFR]

Thank you very much.

So for the x-momentum I add this term dp/dx*vol. In the y-momentum?

what are periodic bcs? Can you state about the velocities in boundaries?

[AliE]

Your source term is naturally a vector. Since you are simulating a couette flow only x-momentum is non zero. Thus let's say that you want to impose a gradient of 0.25, the add at each cell 0.25*volume. The sign will determine the flow direction.

For the moving wall you can set a standard (1 0 0) fixed velocity while (0 0 0) for the other. Periodic bc means that you inlet velocity will take the outlet velocity and
Pressure as boundary condition. Immagine to have a rigid translation of your channel in order to produce an infinite channel. Those conditions are used to simulate something repeting infinitly such as couette or channel flows. In practice you must consider the bundary faces as internal faces.

[IbizaFR]

I understand.
So in the x-momentum I add lets say 0.25*volume and in the y-momentum nothing.

For the velocities at the moving wall I have u=1 (dimensionless), v=0 and at the stationary wall u=v=0.

For the inlet what is the relationship that connects with outlet velocity and pressure?
Thank you.

[AliE]

Yes, everithing is understood as non dimensional. In order to create a couple of periodic boindaries you have to adjust a bit the data structure of you mesh. For examole yuo can create a couple of "ghost" columns cointaing the values of velocity and pressure at I=1 and I=NI. If your code is not parallel you can choose other simple options, but you have to make sure that the outlet cells know the inlet cells values and vice versa. Just to give you an example, if you want to calculate the face velocity at outlet by linear interpolation it would be equalt to uf=(u1+uN)/2 where u1 and uN are the velocity values at the cell centers for the first and the last sets of cells along y-direction in your channel. The same applies for other physical quantities. Your coeffs will become the same as the internal face case. As you can see you should "cheat" your code in order to mimic an infinite channel.

[IbizaFR]

Thank you AliE,

Before continuing in order to make it crystal clear is it right for the walls to take θp/θn=0 where n is the normal direction?

[AliE]

Sure. Pressure is periodic at inlet/oultlet and zero grad at solid walls. Also look at this paper which is very usefull to you:

Implementation of boundary conditions in the finite-volume pressure-based method—Part I: Segregated solvers. By Darwish et al.

[IbizaFR]

Thank you.

So for the inlet and outlet I have to relate the velocity and the pressure in these two positions. I am right?

[AliE]

Yes think the inlet as the neoghbour of the outlet.

[IbizaFR]

Thank you AliE,

Lets take the simple Couette flow where dp/dx=0.
So I have to put zero all the terms regarding the pressure gradient in the x and y momentum provisional velocities equations. Is it correct this?
In this case do I have in my code pressure correction equation since there is not any pressure?

[AliE]

Sure. You have to see the pressure eq as sonething that enforce and guarantee mass consevation (continuity eq). In you simple case however the pressure field will ve constan.

[IbizaFR]

Do I solve pressure correction equation in this case?

[AliE]

Yes solve it in any case.

[IbizaFR]

I solve it and the pressure correction calculates the velocity correction?

[AliE]

Yes this is thr standard alghoritm. If the pressure correction will be constant than its gradient will be zero e there will be no correction. But this is a parcicular case and in any case if you have problem with bcs your code will diverge highlithing the error.

[IbizaFR]

Thank you,

So I solve pressure correction equation then calculate velocity correction and then I enforce constant presssure everywhere e.g. p(1:nx,1:ny,n+1)=1 where nx, ny the x,y nodes and n the itelation level. Is it correct?

[AliE]

No. Solve momentum, solve pressure correction, correct velocity field and pressure field (p=p+p'). The pressure is a result of the computation.

[IbizaFR]

How will I enforce constant pressure field everywhere since I have only simple couette flow?

[AliE]

You are misunderstanding a bit what you are going to do.
You have not to impose a constant pressure. The pressure is a reuslt of your calculation and it is expected to be constant in your case (and only in your particular case). So simply follow the steps of alghoritm. If you have a pressure gradient, yhen add it as body force.