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

Lift and drag coefficient with strange values for NACA airfoil

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 15, 2008, 07:07
Default Hi everybody. I'm trying to
  #1
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
Hi everybody.

I'm trying to simulate the flow around a NACA airfoil with the simpleFoam solver using the Spalart Allmaras model with Re=1e6. Strangely, the boundary layer is very thick resembling flow separation. Also, the drag coefficient is 0.2 which, for an airfoil at 0 AOA, is quite high and the lift coefficient is 0.07 while xfoil gives me 0.5. I can understand that the drag should be higher using a turbulence model, but why the lift coeff is so low?.
I'm posting here all the files in the system directory just in case

--------------------------------------------------
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5 |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application simpleFoam;

startFrom startTime;

startTime 4000;

stopAt endTime;

endTime 8000;

deltaT 1;

writeControl timeStep;

writeInterval 100;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

graphFormat raw;

runTimeModifiable yes;

functions
(
forces
{
type forces;
functionObjectLibs ("libforces.so"); //Lib to load -> dylib on Mac and so on Linux
patches (profile1 profile2);//Name of patche to integrate forces
rhoInf 1.0; //Reference density for fluid - can be changed later ...
CofR ( 0 0 0);
}
forceCoeffs
{
type forceCoeffs;
functionObjectLibs ("libforces.so");
patches (profile1 profile2);
rhoInf 1.0;
CofR (0 0 0);
liftDir (0 1 0);
dragDir (1 0 0);
pitchAxis (0 0 0);
magUInf 100.0;
lRef 0.1;
Aref 0.01;
} );

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

----------------------------------------------

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5 |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear limited 0.7;
laplacian(nu,U) Gauss linear limited 0.7;
laplacian((1|A(U)),p) Gauss linear limited 1;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}

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

-----------------------------------------------

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5 |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p PCG
{
preconditioner DIC;
tolerance 1e-08;
relTol 0;
};

U PBiCG
{
preconditioner DILU;
tolerance 1e-08;
relTol 0;
};
k PBiCG
{
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
};
epsilon PBiCG
{
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
};
R PBiCG
{
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
};
nuTilda PBiCG
{
preconditioner DILU;
tolerance 1e-08;
relTol 0.1;
};
}
/*
k BICCG 1e-06 0;
epsilon BICCG 1e-06 0;
R BICCG 1e-06 0;
nuTilda BICCG 1e-06 0;
*/
SIMPLE
{
nNonOrthogonalCorrectors 2;
}


PISO
{
nCorrectors 1;
nNonOrthogonalCorrectors 1;
/* pRefCell 0;
pRefValue 0;*/
}

relaxationFactors
{
p 0.3;
U 0.7;
k 0.5;
epsilon 0.5;
/* R 0.7;*/
nuTilda 0.7;
}

// ************************************************** *********************** //
antonio_ing is offline   Reply With Quote

Old   December 15, 2008, 10:21
Default I found the porblem. Close to
  #2
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
I found the porblem. Close to the Trailing edge the mesh was not so strong. I have increased there the number of cells there and now everything is better...!
antonio_ing is offline   Reply With Quote

Old   December 16, 2008, 05:51
Default what did you use to mesh it?
  #3
Member
 
Leonardo Honfi Camilo
Join Date: Mar 2009
Location: Delft, Zuid Holland, The Netherlands
Posts: 47
Rep Power: 8
lhcamilo is on a distinguished road
what did you use to mesh it?
lhcamilo is offline   Reply With Quote

Old   December 17, 2008, 09:00
Default gmsh. It is nice and powerful.
  #4
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
gmsh. It is nice and powerful. I wrote the mesh file through matlab and then compiled it with gmsh.
antonio_ing is offline   Reply With Quote

Old   December 19, 2008, 05:49
Default Hi Antonio, did you use tet m
  #5
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Hi Antonio,
did you use tet mesh, or a transfinite one?
How did you specifie the b.l. stretching?
I'm managing to use gmsh for airfoil meshing, but I have problems with wall resolving..
ivan_cozza is offline   Reply With Quote

Old   December 29, 2008, 04:29
Default you are right ivan. Actually i
  #6
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
you are right ivan. Actually i have used the transfinite algorithm close to the profile surface and in the first part of the wake. With this method you can easily define the vertical or horizontal stretching of the cells close to the profile. The gmsh tutorials t3.geo and t6.geo helped me a lot in doing that
antonio_ing is offline   Reply With Quote

Old   December 29, 2008, 11:17
Default Antonio, I tryed to do the sa
  #7
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Antonio,
I tryed to do the same, and my guess was to use the tet meshing far from the foil in order to reduce the number of elements. Did you do something similar?
With which kind of results?
ivan_cozza is offline   Reply With Quote

Old   December 29, 2008, 16:23
Default actualy I did the opposite sin
  #8
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
actualy I did the opposite since it is more important to have a cartesian grid close to the airfoil in order to simulate correctely the boundary layer evolution.

If you give me your mail I can send you some screenshots of my mesh or the code I have used in gmsh. Hopefully this will help you
antonio_ing is offline   Reply With Quote

Old   December 30, 2008, 01:48
Default hi everyone i could not
  #9
Senior Member
 
NAVEEN.K.M
Join Date: Mar 2009
Location: Bangalore, Karnataka, india
Posts: 114
Rep Power: 8
naveen is on a distinguished road
hi everyone

i could not solve naca0012 case because my vertices are not working correctly...could u please send me naca 0012 vertices..
naveen is offline   Reply With Quote

Old   January 20, 2009, 04:52
Default hi everybody i am trying
  #10
Senior Member
 
NAVEEN.K.M
Join Date: Mar 2009
Location: Bangalore, Karnataka, india
Posts: 114
Rep Power: 8
naveen is on a distinguished road
hi everybody

i am trying to work on C type domain of naca 0012 airfoil case in openfoam 1.4.1,but i dont hav the vertices and edges for naca 0012 airfoil,can u please send me the vertices and edges for naca 0012 airfiol..
naveen is offline   Reply With Quote

Old   January 20, 2009, 05:22
Default hi everybody i am trying
  #11
Senior Member
 
NAVEEN.K.M
Join Date: Mar 2009
Location: Bangalore, Karnataka, india
Posts: 114
Rep Power: 8
naveen is on a distinguished road
hi everybody

i am trying to work on C type domain of naca 0012 airfoil case in openfoam 1.4.1,but i dont hav the vertices and edges for naca 0012 airfoil,can u please send me the vertices and edges for naca 0012 airfiol..
naveen is offline   Reply With Quote

Old   January 20, 2009, 05:38
Default Naveen, http://www.aerospac
  #12
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Naveen,

http://www.aerospaceweb.org/question...ls/q0100.shtml

here you can find the naca 4 digit equation, but did you ever try to google "naca 0012"?
ivan_cozza is offline   Reply With Quote

Old   January 20, 2009, 05:44
Default I have left the code, case and
  #13
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,753
Rep Power: 21
hjasak will become famous soon enough
I have left the code, case and results for you on:

http://powerlab.fsb.hr/ped/kturbo/Op...P/naca0012.tgz

It is unbelievable that all the work you did on this since December 2008 is to ask the same question over and over again. Which school did you go to?

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   January 20, 2009, 05:49
Default hi yes i tried that webs
  #14
Senior Member
 
NAVEEN.K.M
Join Date: Mar 2009
Location: Bangalore, Karnataka, india
Posts: 114
Rep Power: 8
naveen is on a distinguished road
hi

yes i tried that website but i need the edges of naca 0012 airfoil to work in openfoam 1.4.1...can u send me the tutorial of naca 0012 airfoil in openfoam 1.4.1...including how to solve and postprocessing..
naveen is offline   Reply With Quote

Old   January 20, 2009, 05:58
Default Dear Naveen, do you really ex
  #15
caw
Member
 
Christian Winkler
Join Date: Mar 2009
Location: Mannheim, Germany
Posts: 63
Rep Power: 8
caw is on a distinguished road
Dear Naveen,
do you really expect other people to do all the work for you?

Have a look here (google is your friend)

http://www.basiliscus.com/ProaSectio.../AppendixD.pdf

http://www.ppart.de/aerodynamics/profiles/NACA4.html

Christian
caw is offline   Reply With Quote

Old   January 20, 2009, 06:03
Default hi thanks for replying..
  #16
Senior Member
 
NAVEEN.K.M
Join Date: Mar 2009
Location: Bangalore, Karnataka, india
Posts: 114
Rep Power: 8
naveen is on a distinguished road
hi

thanks for replying....
naveen is offline   Reply With Quote

Old   September 13, 2012, 12:21
Default
  #17
Member
 
R. P.
Join Date: Jul 2010
Location: Brazil
Posts: 63
Rep Power: 6
Rophys is on a distinguished road
Hi all,

I'm doing some simulations using the OpenFoam 2.1 and I'd like to measure the following aerodynamic forces: Drag Coefficient, Lift Coefficient, Axial-force coefficient, Normal-force coefficient, Pitching-moment coefficient and center of pressure.

How can I state this in the controlDict ?

Many thanks.
Rophys 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
Fluent Good Lift coefficient BAD drag coefficient Rif Main CFD Forum 4 March 9, 2010 11:52
lift & drag coefficient on airfoil n. natik FLUENT 7 December 17, 2009 09:30
NACA 23020 airfoil drag and lift calculation. Zmur CFX 2 December 23, 2008 17:35
Naca 0012 lift/drag values Raj FLUENT 5 August 9, 2006 16:27
Naca airfoil with to much drag Andreas CFX 6 March 17, 2006 07:13


All times are GMT -4. The time now is 21:13.