[Tobi]

Hi all,

I have a question to the totalPressure. First the simple test case I made. It is a 2d rectangle with inlet left and outlet right; top bottom wall and of course empty patches for front and back.

The velocity at inlet and outlet is adjusted by the pressure.
The pressure at the outlet is totalPressure with p0 = 0 and the inlet is -5 or +5 as fixed value.

Behavior

Okay if I set a fixed pressure drop (pressure gradient) of inlet and outlet (both fixedValue), the fluid will accelerate till infinity because we will keep the pressure gradient, that will accelerate the fluid more and more. That is clear.

If we fix the inlet pressure and use totalPressure at the outlet, then (as far as I understand) the pressure at the faces should be adjusted due to the fluxes and the totalPressure we set.

    1. incompressible subsonic:
        \f[
            p_p = p_0 - 0.5 |U|^2
        \f]

So far so good.

Working boundary condition set-up

inlet p fixedValue < 0
outlet p totalPressure p0 = 0

In this set-up the flow is entering at the outlet and leaving at the inlet and hence the pressure at the faces are adjusted till a balance between pressure and fluxes fulfill the totalPressure condition.

NOT working boundary condition set-up

inlet p fixedValue > 0
outlet p totalPressure p0 = 0

In this set-up the flow is entering at the inlet and leaving at the outlet but the pressure at the outlet is kept at 0 Pa (like the total pressure).


Reason for the "NOT working boundary condition set-up"

The reason why the second case is not working as I expect is due to the fact that this line in the source file:

operator==(p0p - 0.5*(1.0 - pos(phip))*magSqr(Up));

will be at least

operator==p0p;

because pos(phip) is 1 if the flux is directed outside and 0 if the flux is directed inside.

My question

Is there a reason for the direction (only working if the flux is going inside that patch?). I made a test where I switched from:

operator==(p0p - 0.5*(1.0 - pos(phip))*magSqr(Up));

to

operator==(p0p - 0.5*(-1.0*sign(phip))*magSqr(Up));

that is actual working for the inflow and outflow.

For the inflow I get the same result as the normal one (that is clear)
For the outflow I get another result that looks okay but the velocities are total different.


At the moment I am only searching for the reason why outflow is not valid (or forbidden)?

Thanks in advance,
Tobi


Answer:
Its clear. TotalPressure acts as fixedValue for outflow and for inflow we adjust the pressure due to fluxes.