[letzel]

Dear Foamers,

given a measured time series of patch-normal velocity Un, is it possible with groovyBC and timelines or lookuptables to set time-dependent inletOutlet boundary conditions on that patch?
I want to set up the following scenario:

• in case of inflow (say Un >= 0), set inflow boundary condition (value of Un uniform over patch),
• in case of outflow (say Un < 0), set zeroGradient pressure boundary condition.

Regards,
Marcus

[gschaider]

Originally Posted by letzel

Dear Foamers,

given a measured time series of patch-normal velocity Un, is it possible with groovyBC and timelines or lookuptables to set time-dependent inletOutlet boundary conditions on that patch?
I want to set up the following scenario:

• in case of inflow (say Un >= 0), set inflow boundary condition (value of Un uniform over patch),
• in case of outflow (say Un < 0), set zeroGradient pressure boundary condition.

Regards,
Marcus

Yep. Assuming you have a timeline Un defined

valueExpression "-normal()*Un";
fractionExpression "Un<0 ? 0 : 1";

[letzel]

Originally Posted by gschaider

Yep. Assuming you have a timeline Un defined

valueExpression "-normal()*Un";
fractionExpression "Un<0 ? 0 : 1";

Excellent. That means the U boundary condition on that patch is done.

Follow-up question: How do I have to set the corresponding pressure boundary condition on that patch? I would expect that a similar case distinction depending on the sign of Un needs to be implemented for p.

[gschaider]

Originally Posted by letzel

Excellent. That means the U boundary condition on that patch is done.

Follow-up question: How do I have to set the corresponding pressure boundary condition on that patch? I would expect that a similar case distinction depending on the sign of Un needs to be implemented for p.

That depends on your physics. But any switching here can be done in a similar fashion (of course you can use Un inside a pressure condition as well)

[letzel]

Sorry, the description of my scenario given above was not precise. Let me try to clarify. Basically I am planning to mimick "inletOutlet", the patch field type derived from mixed (table 5.4 in OF 2.2.2 user guide) using the patch-normal velocity component Un:

• in case of inflow (Un >= 0),
  • for pressure use zeroGradient, and
  • for U use valueExpression "-normal()*Un"
• in case of outflow (Un < 0),
  • for pressure use fixedValue, and
    
  • for U use zeroGradient.

My question is how to implement this pressure boundary condition with the help of fractionExpression "Un<0 ? 0 : 1" or
fractionExpression "Un<0 ? 1 : 0".

Whether or not this pressure boundary condition will work well in this case, this will be the next question...

[gschaider]

Originally Posted by letzel

Sorry, the description of my scenario given above was not precise. Let me try to clarify. Basically I am planning to mimick "inletOutlet", the patch field type derived from mixed (table 5.4 in OF 2.2.2 user guide) using the patch-normal velocity component Un:

• in case of inflow (Un >= 0),
  • for pressure use zeroGradient, and
  • for U use valueExpression "-normal()*Un"
• in case of outflow (Un < 0),
  • for pressure use fixedValue, and
    
  • for U use zeroGradient.

My question is how to implement this pressure boundary condition with the help of fractionExpression "Un<0 ? 0 : 1" or
fractionExpression "Un<0 ? 1 : 0".

Whether or not this pressure boundary condition will work well in this case, this will be the next question...

I'd say the second one. I'm not quite sure what your problem is: the meaning of the ?:-operator or fractionExpression. Because if you understand these two it should be quite clear

[letzel]

Thank you, Bernhard, the second one worked out fine.

As a follow-up question, I should now compare this with a conditional totalPressure boundary condition in case of outflow as follows:

• in case of inflow (Un >= 0),
  • for pressure use zeroGradient, and
  • for U use valueExpression "-normal()*Un"
• in case of outflow (Un < 0),
  • for pressure use totalPressure, and
    
  • for U use zeroGradient.

Does the groovyBC syntax allow for this? If so, how?

[gschaider]

Originally Posted by letzel

Thank you, Bernhard, the second one worked out fine.

As a follow-up question, I should now compare this with a conditional totalPressure boundary condition in case of outflow as follows:

• in case of inflow (Un >= 0),
  • for pressure use zeroGradient, and
  • for U use valueExpression "-normal()*Un"
• in case of outflow (Un < 0),
  • for pressure use totalPressure, and
    
  • for U use zeroGradient.

Does the groovyBC syntax allow for this? If so, how?

U we've already got covered I think.

p shouldn't be a big problem either. The catch is that you've got to specify the total pressure as an expression. You can't use the readymade BC