[MagnusJohannes]

Hi,

I've been searching through the threads and online. I am looking to write a UDF for blood flow through the aorta. The inlet velocity profile needs to be parabolic and transient. I have a polynomial for the pulse, which I have used in a UDF, but the profile needs to be parabolic as well.

I have tried modifying existing codes, but as a beginner my results haven't been good. Any ideas?

Thanks,

John

[MagnusJohannes]

I've pasted my code below.



/************************************************** *********************/
/* vinlet_udf.c
*/
/* UDFs for specifying time dependant velocity profile boundary condition
*/
/************************************************** *********************/
#include "udf.h"

DEFINE_PROFILE(unsteady_velocity, thread, position)

{

real y, x[ND_ND], rad;

real z;

face_t f;

begin_f_loop(f, thread)

{
real t = RP_Get_Real("flow-time");

F_CENTROID(x, f, thread);

y = x[1];

z = x[2];

rad = pow(y*y + z*z, 0.5);


F_PROFILE(f, thread, position) = 1538234.528 * t*t*t*t*t*t*t*t*t - 6233585.283 * t*t*t*t*t*t*t*t + 10657407.82 * t*t*t*t*t*t*t - 10008440.91 * t*t*t*t*t*t + 5636807.043 * t*t*t*t*t - 1954154.406 * t*t*t*t + 412312.4917 * t*t*t - 50109.49185 * t*t + 2903.215824 * t + 7.688636046; /* This is my pulsatile polynomial*/

}

end_f_loop(f, thread)