[elobb]

I am trying to model flow through a channel, Womersley number one, with a sinusoidal pressure gradient. I saw elsewhere on this forum that using a udf with a scheme file to implement it at each timestep is the way forward, though I am having trouble programming the scheme file.

This is the udf and scheme file I started with..

udf...

#include "udf.h"

DEFINE_EXECUTE_AT_END(execute)
{
real pressure, t;
t=CURRENT_TIME;
pressure=sin(t);; /* The pressure gradient can be a function of time. */
RP_Set_Real("flag", pressure); /* This will store the pressure gradient value in a variable called flag*/
}

scheme file...

(rp-var-define 'flag 0 'real #f)
(define print
(lambda()
(define nu (rpgetvar 'flag))
(rpsetvar 'periodic/pressure-derivative nu)
)
)

Any help would be much appreciated

Regards