Plane channel flow
 

Hello,

I'm a very begginer in OpenFOAM. I'm just trying to simulate a laminar plane channel flow using icoFoam.

I got the conditions from the analytical solution and fvSchemes and fvSolution from Cavity tutorial.
My problem is that Courant number keeps growing until explode.

The setup is attached.

Thanks in advance.

Hi,

Here are two test-cases that may help you:
http://openfoamwiki.net/index.php/Ma...leSimple2DFlow
http://openfoamwiki.net/index.php/Bl...Flow_Benchmark

Regards,
Jose Santos

Obrigado, Jose Santos!

I've got it to work. I had a problem in the BC and my time step was too high.

Hi Rodrigo,

I try to solve the same case and I think I also have some problems with my boundary conditions, especially with the one for the outlet. Could you please post your boundary conditions? That would be great! Thanks a lot!

Hello

I also started with this simple case to collect experiments in OF.

Copied the cavity tutorial and generated new mesh with new BC and initial conditions to solve with icoFoam a 2D channel case.

Now I am wondering - in folder .../case/0/ the dimension definition for pressure p?

-> dimensions [ 0 2 -2 0 0 0 0] -> this is no pressure!?

Anybody an idea about this background?

Thanks & regards

Hi,

p is the kinematic pressure. If you multiply it by the density, you'll get Pascal units (kg/ms^2).

Regards,
Jose Santos

Hello.

You have my setup attached.
Hope it helps.

@santos

O.K. p is the kinematic pressure, I understand this. But what should I do if I don't want to set-up a dynamic pressure but a static pressure? If I just change the dimension, I get an error compiling it, which says, that the dimensions don't match. So, what do I have to do if I want to initialize my flowfield with a static pressure?

Hello Sven,

p is the static pressure divided by density (this is what Santos means with kinematic pressure).
You should not change units to set a static pressure, just divide the value by the constant density and set the resulting value.