[cardio]

I'm trying to stimulate flow in the Left ventricular region of the heart during diastole for which I have written the following UDF :


#include <stdio.h>
#include "udf.h"
#include "dynamesh_tools.h"

#define pi 3.141592654
static real r_lv=0.023, flowrate_off=0.8;
static real force, dv ,flowrate ,t_use;

DEFINE_CG_MOTION(Velociraptor,dt,vel,omega,time,dt ime)
{
Thread *t=DT_THREAD ((Dynamic_Thread *)dt);
face_t f;
real *AA;
/*reset velocities*/
NV_S(vel, =, 0.0);
NV_S(omega, =, 0.0);
if (!Data_Valid_P())
return;
/*t=DT_THREAD (dt);*/
force = 0.0;
begin_f_loop (f, t)
{AA = F_AREA_CACHE (f, t);
force += F_P (f,t) * AA[0];}
end_f_loop(f,t);
t_use = 0.639638*CURRENT_TIME/0.6 - 0.039638; /* revise later */
flowrate = -99.865 *pow(t_use, 4);
flowrate = flowrate + 119.124 *pow(t_use, 3);
flowrate = flowrate - 43.389 *pow(t_use, 2);
flowrate = flowrate + 4.3084 *t_use;
flowrate = flowrate + 0.1951;
flowrate =(flowrate + 0.052) * 0.3/0.3748;

flowrate = - flowrate /1000.0 * flowrate_off; /* /1000.: unit factor */
vel[1] = flowrate / pi*r_lv*r_lv ;
}

The problem I'm facing is that the geometry is not getting updated.It updates for the first time step,but then that remains constant:

Updating mesh to time 1.00000e-03 (step = 00001)...
Dynamic Mesh Statistics:
Minimum Volume =2.91284e-011
Maximum Volume =2.54127e-009
Maximum Cell Skew =8.51365e-001 (cell zone 4)
Maximum Face Skew =4.59499e-001 (face zone 19)
done.
Updating mesh to time 2.00000e-03 (step = 00002)...
Dynamic Mesh Statistics:
Minimum Volume =2.91284e-011
Maximum Volume =7.82814e-009
Maximum Cell Skew =8.51365e-001 (cell zone 4)
Maximum Face Skew =4.59499e-001 (face zone 19)
done.
Updating mesh to time 3.00000e-03 (step = 00003)...
Dynamic Mesh Statistics:
Minimum Volume =2.91284e-011
Maximum Volume =7.82815e-009
Maximum Cell Skew =8.51365e-001 (cell zone 4)
Maximum Face Skew =4.59499e-001 (face zone 19)
done.

Can someone tell me what the problem might be ?

[lgm]

Hi, i'm having the same problem
Do you resolve this? thanks