[vitorspadeto]

Can you tell me if my code is right?
My intention is to apply a pulsating pressure at the inlet with the FUNCION: P = A + A*sin(pi*time) , WHERE A=100 Pa (unit in Pa) .... NOTE that the units are [0 2 -2 0 0 0 0] on the description of p file.



This is my code:

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    
    outlet
    {
        type fixedValue;
        value uniform 0;
    }
    

inlet
    {
     type            groovyBC;
     variables "var=100/1050;";
     valueExpression "(var+(var*sin(pi*time())))"; //*normal();
    }
 

wall
    {
        type extrapolatedPressure;
        value uniform 0;
    }

1- I divided "var" by 1050kg/m³(that is density of fluid). Am I right? WHERE A=100 Pa (unit in Pa) .... NOTE that the units are [0 2 -2 0 0 0 0]


2- I don't use "(var+(var*sin(pi*time())))*normal()" because I believe that pressure is not a vector. Am I right?


Best Regards,
Vitor Spadeto Venturin

[gschaider]

Originally Posted by vitorspadeto

Can you tell me if my code is right?
My intention is to apply a pulsating pressure at the inlet with the FUNCION: P = A + A*sin(pi*time) , WHERE A=100 Pa (unit in Pa) .... NOTE that the units are [0 2 -2 0 0 0 0] on the description of p file.



This is my code:

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    
    outlet
    {
        type fixedValue;
        value uniform 0;
    }
    

inlet
    {
     type            groovyBC;
     variables "var=100/1050;";
     valueExpression "(var+(var*sin(pi*time())))"; //*normal();
    }
 

wall
    {
        type extrapolatedPressure;
        value uniform 0;
    }

1- I divided "var" by 1050kg/m³(that is density of fluid). Am I right? WHERE A=100 Pa (unit in Pa) .... NOTE that the units are [0 2 -2 0 0 0 0]


2- I don't use "(var+(var*sin(pi*time())))*normal()" because I believe that pressure is not a vector. Am I right?

Looks OK. Why are you asking? To check whether the BC does what you want without running the solver you can use replayTransientBC (which comes with swak)