[MFWilliams]

Hello,


Is it possible to implement a time varying velocity boundary conditon using GroovyBC with a polynomial or lookup table?


I have been having great trouble trying to implement a time varying velocity BC. The aim is to simulate inlet velocity from three identical heart beats in a row and I have started trying GroovyBC to limited success.


I have managed to create a boundary condition for one heart cycle using GroovyBC and a polynomial but am not sure how to repeat it for three cycles.



I have now been looking at this post to try implement a lookup velocity table, but this one does not seem to be time varying.

Originally Posted by dferrando

Thank you very much. I do that but I have another Error that I don't understand:

--> FOAM FATAL IO ERROR: 
attempt to read beyond EOF

file: /home/diego/OpenFOAM/diego-4.1/run/windAroundBuildings/0/U.boundaryField.inlet.lookuptables.valueColumns at line 74.

    From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
    in file db/IOstreams/Tstreams/ITstream.C at line 82.

FOAM exiting

My lookuptables and expresion:

    inlet
    {
        type            groovyBC;
        lookuptables
        (
            {
            name profileX;
            fileName "$FOAM_CASE/0/U1.csv";
            outOfBounds clamp;
            readerType csv;
            hasHeaderLine 2;
            timeColumn 0;
            valueColumns 1;
            separator " "; 
            }        
        );
        lookuptables
        (
            {
            name profileY;
            fileName "$FOAM_CASE/0/U1.csv";
            outOfBounds clamp;
            readerType csv;
            hasHeaderLine 2;
            timeColumn 0;
            valueColumns 2;
            separator " "; 
            }        
        );
        lookuptables
        (
            {
            name profileZ;
            fileName "$FOAM_CASE/0/U1.csv";
            outOfBounds clamp;
            readerType csv;
            hasHeaderLine 2;
            timeColumn 0;
            valueColumns 3;
            separator " "; 
            }   
       );

    valueExpression "vector(profileX(pos().z), profileY(pos().z), profileZ(pos().z))";
    }

Thanks Bernhard!

Diego.

Please can someone tell me if it is possible to implement a time varying velocity lookup table GroovyBC?



Thank you for your help in advance
Merlin