VIV Square Cylinder
 

Hey, so Im doing a simulation of viv for a square cylinder and Im trying to replicate the results of this article: http://www.sciencedirect.com/science...89974613000832

but so far I havent gotten close. Everything not moving (mesh, residuals etc.) has been validated so the parts Im not certain of is the dynamics, like the UDF and meshing settings.

My UDF:

#include "udf.h"


DEFINE_SDOF_PROPERTIES(sdof_props,prop,dt,time,dti me)
{
real cg;
real k = 3.20886;
cg = DT_CG(dt)[1];
prop[SDOF_MASS] = 5;
prop[SDOF_LOAD_F_Y] = -k*cg;
prop[SDOF_ZERO_TRANS_X] = TRUE;
prop[SDOF_ZERO_ROT_X] = TRUE;
prop[SDOF_ZERO_ROT_Y] = TRUE;
}

where k is calculated as 4*pi^2*fn*m where fn=3.1875/Re where Re=25 and m=5.

The dynamic mesh settings Im using are:

Smoothing, Layering, Six DOF and implicit update. With use of the SIMPLE scheme, Least Squares Cell Based, Standard and QUICK.

Anyone got any ideas? Cheers