
[Sponsors] 
February 9, 2011, 08:55 
airfoil study

#1 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
Hi,
I am studying a turbulent flow that involves a laminar profile NLF_0414f. I have imposed conditions written at the end of the post, but I got to the end of the calculations of the Cl and Cd absurd values, over 10 ^ 10. Could you help me to understand where I'm wrong? nut dimensions [0 2 1 0 0 0 0]; internalField uniform 1.79e05; boundaryField { inlet { type freestream; freestreamValue uniform 1.79e05; } outlet { type freestream; freestreamValue uniform 1.79e05; } wall { type nutSpalartAllmarasWallFunction; value uniform 0; } frontAndBack { type empty; } } nutilda dimensions [0 2 1 0 0 0 0]; internalField uniform 0.1; boundaryField { inlet { type freestream; freestreamValue uniform 0.1; } outlet { type freestream; freestreamValue uniform 0.1; } wall { type fixedValue; value uniform 0; } frontAndBack { type empty; } } U dimensions [0 1 1 0 0 0 0]; internalField uniform (116 2 0); boundaryField { inlet { type freestream; freestreamValue uniform (116 2 0); } outlet { type freestream; freestreamValue uniform (116 2 0); } wall { type fixedValue; value uniform (0 0 0); } frontAndBack { type empty; } } P dimensions [0 2 2 0 0 0 0]; internalField uniform 110983; boundaryField { inlet { type freestreamPressure; } outlet { type freestreamPressure; } wall { type zeroGradient; } frontAndBack { type empty; } } RAS PROPERTIES RASModel SpalartAllmaras; turbulence on; printCoeffs on; TRANSPORT PROPERTIES transportModel Newtonian; rho rho [ 1 3 0 0 0 0 0 ] 1.288; nu nu [0 2 1 0 0 0 0] 1.79e05; _____________CONTROL DICT_______________ application simpleFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 0.035; deltaT 5e06; writeControl runTime; writeInterval 0.005; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression compressed; timeFormat general; timePrecision 6; graphFormat raw; runTimeModifiable yes; adjustTimeStep no; maxCo 0.5; maxDeltaT 1e02; functions { /*forces { type forces; functionObjectLibs ( "libforces.so" ); // lib to load outputControl timeStep; outputInterval 1; patches ( wall // change to your patch name ); // name of fields //pName p; //UName U; log true; // dump to file rhoName rhoInf; rhoInf 1.225; CofR ( 0 0 1 ); } */ forcesCoeffs { type forceCoeffs; functionObjectLibs ( "libforces.so" ); // lib to load outputControl timeStep; outputInterval 1; patches ( wall// change to your patch name ); // name of fields //pName p; //UName U; log true; // dump to file rhoName rhoInf; rhoInf 1.288; CofR ( 0 0 0 ); liftDir ( 0.087 0.996 0 ); dragDir ( 0.996 0.897 0 ); pitchAxis ( 0 0 1); magUInf 138.835; lRef 1; Aref 1; } } // Definition of terms: // rhoInf  reference density // CofR  Centre of rotation // dragDir  Direction of drag coefficient // liftDir  Direction of lift coefficient // pitchAxis  Pitching moment axis // magUinf  free stream velocity magnitude // lRef  reference length // Aref  reference area // ************************************************** *********************** // 

February 9, 2011, 09:38 

#2 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
Hello, salvoblack,
first thing that catches my eye: the freestream values of nutilda are much too high. Values of 3 to 5 times the molecular viscosity are appropriate for external aerodynamic cases (see http://turbmodels.larc.nasa.gov/spalart.html ). Greetings, Felix. 

February 9, 2011, 10:31 

#3 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
thanks felixL, i tried to put nutilda=0,0000716 but i have also in this case values of Cl and Cd too high.
what do you think about about the other boundary conditions??? 

February 9, 2011, 10:43 

#4 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
The other BCs are looking alright. Are the reference values for calculating the aerodynamic coefficients correct? (ref. area = 1mē)
What does checkMesh say? Greetings, Felix. 

February 9, 2011, 10:51 

#5 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
this is a problem!
i use also fluent (with the same mesh, obviously) and in the window "reference values" i put area=1 and length=1 and i have no problems. how can i see if these values are correct for openfoam?? p.s checkMesh says OK 

February 9, 2011, 10:59 

#6 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
Hello,
reference Area = chord length * wing span For 2D cases the wing span can be set arbitrarily but it has to be accounted for when calculating the ref area. Use the checkMesh output to evaluate the spanwise dimension of your domain. Greetings, Felix. 

February 9, 2011, 11:04 

#7 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
sorry felixl but i don't understand, i'm really new with OF
This is what CheckMesh says in the terminal: Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 148992 internal points: 0 faces: 295680 internal faces: 146688 cells: 73728 boundary patches: 4 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 73728 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface topology outlet 192 386 ok (nonclosed singly connected) inlet 768 1538 ok (nonclosed singly connected) wall 576 1152 ok (nonclosed singly connected) frontAndBack 147456 148992 ok (nonclosed singly connected) Checking geometry... Overall domain bounding box (26.0276 31.0934 0.823828) (30 29.304 0.823828) Mesh (nonempty, nonwedge) directions (1 1 0) Mesh (nonempty) directions (1 1 0) All edges aligned with or perpendicular to nonempty directions. Boundary openness (1.12341e17 1.82382e18 1.02153e19) OK. Max cell openness = 5.36428e13 OK. Max aspect ratio = 818.404 OK. Minumum face area = 4.62385e09. Maximum face area = 4.28893. Face area magnitudes OK. Min volume = 7.61852e09. Max volume = 3.86456. Total volume = 4779.57. Cell volumes OK. Mesh nonorthogonality Max: 89.9384 average: 7.88302 *Number of severely nonorthogonal faces: 1021. Nonorthogonality check OK. <<Writing 1021 nonorthogonal faces to set nonOrthoFaces Face pyramids OK. Max skewness = 0.28683 OK. Mesh OK. End Thanks a lot 

February 9, 2011, 11:08 

#8 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
Hello,
This line tells you the overall dimension of your computational domain: Overall domain bounding box (26.0276 31.0934 0.823828) (30 29.304 0.823828) The bold numbers are the minimum and maximum z coordinates. I.e. your 2D domain is 1.6477 metres wide which is the same as your wing span. Now calculating the reference area should be straight forward: ref area = 1.6477 m * chord length Greetings, Felix. 

February 9, 2011, 11:27 

#9 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
this is absurd.
1)i've taken a new case from this site : Negative drag coefficient... 2)i just replaced the mesh, and put as you said aref=1.6477. in the case of naca 0006 the results are good. in my case, with my own mesh, the results are wrong. what's the problem??? it is the mesh!!!! i'm really desperate 

February 9, 2011, 11:51 

#10 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
Care to post the log file of your case? It seems to me that the simulation crashed or something, leading to unphysical results and hence the lift and drag coefficients in the order of 10^10.
Please provide the log, fvSchemes and fvSolution files, otherwise I can't really tell why the simulation could've crashed. Greetings, Felix. 

February 9, 2011, 12:05 

#11 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
what do you mean for log file??
fvscheme ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear; } divSchemes { default none; div(phi,U) Gauss linearUpwind Gauss linear; div(phi,nuTilda) Gauss linearUpwind Gauss linear; div((nuEff*dev(grad(U).T()))) Gauss linear; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss linear corrected; laplacian((1A(U)),p) Gauss linear corrected; laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; laplacian(1,p) Gauss linear corrected; } interpolationSchemes { default linear; interpolate(U) linear; } snGradSchemes { default corrected; } fluxRequired { default no; p ; } fvsolution solvers { p { solver GAMG; tolerance 1e06; relTol 0.1; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } U { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e08; relTol 0.1; } nuTilda { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e08; relTol 0.1; } } SIMPLE { nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 0; } relaxationFactors { default 0; p 0.3; U 0.7; nuTilda 0.7; } 

February 9, 2011, 12:12 

#12 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
To create a logfile, run simpleFoam with a command like this:
simpleFoam > log & This will create a file called log where all the residual information will be stored. Also simpleFoam will run in the background then. You can observe the logfile by using the command tail f log Greetings, Felix. 

February 9, 2011, 12:18 

#13 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
ok
this is the file 

February 9, 2011, 12:35 

#14 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
There is a problem with nuTilda, the transport equation isn't being solved. I'm guessing it's got something to do with the boundary setup, though I don't have any recommendations right now. My suggestion would be to rebuild the case from scratch and change the freestream BCs at in and outlet to something more straightforward (i.e. fixedValue at inlet for everything but pressure, zeroGradient at outlet for everything but pressure. Pressure is vice versa).
Greetings, Felix. 

February 9, 2011, 12:48 

#15 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
but the strange thing is that if i use the exemplum of the naca 0006 i have no problems, but if i replace the mesh i have these absurd values.
i don't change anything about the equation and the BCs!!!!! if you note in the checkMesh there are nonorthogonal faces. Is this a problem?? Have you got a good case for an airfoil so i can make some trials???? Thank you very very much. 

February 9, 2011, 13:14 

#16 
Senior Member
Felix L.
Join Date: Feb 2010
Location: Hamburg
Posts: 165
Rep Power: 10 
The nonorthogonal faces shouldn't be that much of a problem. I had worse meshes and I was able to get decent results. Could you post some screenshots of your mesh? Like an overview and some details of the vicinity of the airfoil.
Greetings, Felix. 

February 9, 2011, 13:37 

#17 
Member
Join Date: Oct 2010
Location: Naples
Posts: 50
Rep Power: 8 
felix, i sent you a private message


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Low Speed Airfoil  Mancusi  FLUENT  7  April 3, 2014 06:11 
[GAMBIT] Meshing airfoil using .dat file problem  creggie  ANSYS Meshing & Geometry  10  June 27, 2010 19:24 
ICEM CFD Airfoil Tutorial Input  Tom  ANSYS Meshing & Geometry  1  September 25, 2009 06:50 
Modeling Backflow for a 3D Airfoil (Wing of Finite Span)  Josh  CFX  9  August 18, 2009 11:31 
Airfoil boundary condition  Frank  Main CFD Forum  1  April 21, 2008 18:36 