[farich_d]

Hello, im doing some research about hydrokinetic turbine (transient and 2dimensional), where i attach moving flap in each hydrofoil and control its movement using UDF. The flap moving relatively to the rotation domain.

I always getting "negative cell volume" error while previewing mesh motion immediately (at early time step).

Then i tried to move only the flap without rotate the rotating domain,
and the preview mesh motion can running smoothly and shows how the flap move without getting "negative cell volume" error.

But when i rotate the rotating domain, while the flap is also moving, i keep getting "negative cell volume" error.

here is UDF

#include "udf.h"

DEFINE_CG_MOTION(bukatutup,dt,cg_vel,cg_omega,time ,dtime)

{
real t = CURRENT_TIME;
real pi;
pi = 3.141592654;
NV_S (cg_omega, =, 0.0);

if ((t > 0) && (t <= (0.5)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = 1.125*pi;
}

if ((t > (0.5)) && (t <= (1)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = -1.125*pi;
}

if ((t > (1)) && (t <= (1.5)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = 1.125*pi;
}

if ((t > (1.5)) && (t <= (2)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = -1.125*pi;
}

if ((t > (2)) && (t <= (2.5)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = 1.125*pi;
}

if ((t > (2.5)) && (t <= (3)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = -1.125*pi;
}

if ((t > (3)) && (t <= (3.5)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = 1.125*pi;
}

if ((t > (3.5)) && (t <= (4)))
{
cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = -1.125*pi;
}

}

DEFINE_CG_MOTION(domainrot,dt,cg_vel,cg_omega,time ,dtime)
{
real omega_z;

omega_z = 1.570796;

cg_vel[0] = 0.0;
cg_vel[1] = 0.0;
cg_vel[2] = 0.0;

cg_omega[0] = 0.0;
cg_omega[1] = 0.0;
cg_omega[2] = omega_z;
}