[milad1990]

Hi all

I am studying wake downstream of a cylinder using LES. I am calculating LIFT and DRAG coefficient using both "forcecoeff" function and "force" function but I can not get similar results. This is my "force" and "force coefficient" functions :


forces_object
{
type forces;
functionObjectLibs ("libforces.so");

enabled true;

//writeControl outputTime;
writeControl adjustableRunTime;
writeInterval 0.001;

patches ("cylinder");

pName p;
Uname U;

////Density only for incompressible flows
rho rhoInf;
rhoInf 1000;

////Centre of rotation
CofR (0 0 0);
}
///////////////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////////////
forceCoeffs1
{
// Mandatory entries
type forceCoeffs;
libs ("libforces.so");
patches ("cylinder");


// Optional entries

// Field names
p p;
U U;
rho rhoInf;
rhoInf 1000;

// Reference pressure [Pa]
pRef 0;

// Include porosity effects?
porosity no;

// Store and write volume field representations of forces and moments
writeFields no;

// Centre of rotation for moment calculations
CofR (0 0 0);

// Lift direction
liftDir (0 1 0);

// Drag direction
dragDir (1 0 0);

// Pitch axis
pitchAxis (0 0 1);

// Freestream velocity magnitude [m/s]
magUInf 0.5;

// Reference length [m]
lRef 0.25;

// Reference area [m2]
Aref 0.00625;


writeControl adjustableRunTime;
writeInterval 0.001;
}

////////////////////////////

I use a Python script to calculate Lift and Drag forces which simply add pressure and viscous forces in "force.dat" file for horizontal and vertical direction.
Now when I use formula CL=Lift /(0.5*rho*U^2*Aref) or CD=Drag /(0.5*rho*U^2*Aref) using the same Aref I determined in forceoeff function, I get different results. I do not know why? Any one have any idea?

Thank you in advance.
Best

[ansab_sindhu]

Hi,

I am facing a similar issue, even though, my area and length dimensions are also correct. But the results of the force coefficient and forces file don't match.
Any leads will be helpful.

Regards,

Originally Posted by milad1990

Hi all

I am studying wake downstream of a cylinder using LES. I am calculating LIFT and DRAG coefficient using both "forcecoeff" function and "force" function but I can not get similar results. This is my "force" and "force coefficient" functions :


forces_object
{
type forces;
functionObjectLibs ("libforces.so");

enabled true;

//writeControl outputTime;
writeControl adjustableRunTime;
writeInterval 0.001;

patches ("cylinder");

pName p;
Uname U;

////Density only for incompressible flows
rho rhoInf;
rhoInf 1000;

////Centre of rotation
CofR (0 0 0);
}
///////////////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////////////
forceCoeffs1
{
// Mandatory entries
type forceCoeffs;
libs ("libforces.so");
patches ("cylinder");


// Optional entries

// Field names
p p;
U U;
rho rhoInf;
rhoInf 1000;

// Reference pressure [Pa]
pRef 0;

// Include porosity effects?
porosity no;

// Store and write volume field representations of forces and moments
writeFields no;

// Centre of rotation for moment calculations
CofR (0 0 0);

// Lift direction
liftDir (0 1 0);

// Drag direction
dragDir (1 0 0);

// Pitch axis
pitchAxis (0 0 1);

// Freestream velocity magnitude [m/s]
magUInf 0.5;

// Reference length [m]
lRef 0.25;

// Reference area [m2]
Aref 0.00625;


writeControl adjustableRunTime;
writeInterval 0.001;
}

////////////////////////////

I use a Python script to calculate Lift and Drag forces which simply add pressure and viscous forces in "force.dat" file for horizontal and vertical direction.
Now when I use formula CL=Lift /(0.5*rho*U^2*Aref) or CD=Drag /(0.5*rho*U^2*Aref) using the same Aref I determined in forceoeff function, I get different results. I do not know why? Any one have any idea?

Thank you in advance.
Best

[hewei]

Originally Posted by ansab_sindhu

Hi,

I am facing a similar issue, even though, my area and length dimensions are also correct. But the results of the force coefficient and forces file don't match.
Any leads will be helpful.

Regards,

I have met the same problem. I am doing a 2D problem. I found the dimensionless calculation is correct, the dimension varies far away.!!!

[Reptider]

Originally Posted by hewei

I have met the same problem. I am doing a 2D problem. I found the dimensionless calculation is correct, the dimension varies far away.!!!

Hey, could you give me an advise how to make dimensionless calculation to get right values?

[hewei]

Originally Posted by Reptider

Hey, could you give me an advise how to make dimensionless calculation to get right values?

To do dimensionless simulation, you need to make sure the U=1, reference l=1. Then you will get the correct values.