CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

getting wrong value of drag

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 29, 2013, 20:08
Default getting wrong value of drag
  #1
Member
 
shiv
Join Date: Jun 2012
Location: Lucknow, IN
Posts: 51
Rep Power: 5
shash is on a distinguished road
Hi openfoamers,

I am simulating a 3d flow around a sphere enclosed in a large cube, but problem is that i am getting my drag =1.04 instead of expt. drag =0.66 for reynolds no. of 300, parameters: fluid:air with rho=1.2, nu=1e-04, dia. of sphere=1m, velocity at inlet=0.03m/s ; Can someone please suggest what i am doing wrong.

my controldict where drag lib parameters given :
Quote:
application icoFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 50;

deltaT 0.005;

writeControl timeStep;

writeInterval 200;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;


// ************************************************** *********************** //

functions
{

forces
{
type forceCoeffs;
functionObjectLibs ("libforces.so");
enabled true;

outputControl timeStep;
outputInterval 1;

patches (sphere);
rhoName rhoInf;
rhoInf 1.2;
log true;

CofR (0 0 0);
liftDir (0 1 0);
dragDir (1 0 0);
pitchAxis (0 0 0);

magUInf 0.03;
lRef 1;
Aref 0.785;
}
}

//************************************************** ***************************
shash is offline   Reply With Quote

Old   June 29, 2013, 20:45
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,205
Rep Power: 17
immortality is on a distinguished road
what does vortices look like?is there any?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 30, 2013, 16:28
Default
  #3
Member
 
shiv
Join Date: Jun 2012
Location: Lucknow, IN
Posts: 51
Rep Power: 5
shash is on a distinguished road
Hi Ehsan,
you pointed out correctly, i donot know why i am not getting the double threaded vortices for Re =300 thou with similar settings i have got as expected karman vortex street in 2d case. Got any ideas whats wrong when i am using 3d geometry??
shash is offline   Reply With Quote

Old   June 30, 2013, 16:43
Default
  #4
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,205
Rep Power: 17
immortality is on a distinguished road
Hi Shiv
are you using a turbulent model?which?how-Re or low-Re?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 30, 2013, 17:22
Red face
  #5
Member
 
shiv
Join Date: Jun 2012
Location: Lucknow, IN
Posts: 51
Rep Power: 5
shash is on a distinguished road
no, i am using icofoam solver so its laminar flow.

Ehsan, i noted something that when i use [nu=1e-04, U=0.03, Re=300] i am not getting vortices though when i use [nu=0.01, U=3, Re=300] i am getting vortices , both the cases have same Re why variation in results
shash is offline   Reply With Quote

Old   July 1, 2013, 11:26
Default
  #6
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,205
Rep Power: 17
immortality is on a distinguished road
its unsteady if wakes there are.maybe deltaT isn't low enough.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How can I calculate drag coefficients when there is a blocking object in front hjcaroline FLUENT 0 June 10, 2012 22:49
CFD-Post: Pressure and Viscous Drag on Isosurface Shirock CFX 4 June 4, 2012 01:02
Calculation of Drag Coefficient, Help Please teek22 CFX 1 April 26, 2012 18:41
Airfoil negative drag coefficent James Forrest Main CFD Forum 7 April 9, 2011 07:06
anyone has got this paper? Al Main CFD Forum 3 October 8, 2007 02:04


All times are GMT -4. The time now is 18:04.