|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
A P
Join Date: Feb 2020
Location: Indiana, USA
Posts: 10
Rep Power: 7 ![]() |
I am running simulation for 2D airfoil NACA0015 using simpleFoam. I am getting approximate half the values of Cl for various angle of attack. When I halve the value of area I get the correct result of Cl. I am confused what I am doing wrong. I am attaching condition for angle of attack 6 deg and Re10^5. Please guide me through it.
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application simpleFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 4500; deltaT 1; writeControl timeStep; writeInterval 50; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable true; functions { forces1 { type forceCoeffs; libs ("libforces.so"); writeControl timeStep; writeInterval 1; patches ( "AIRFOIL"); rho rhoInf log true; rhoInf 1.225; liftDir (-0.104 0.994 0); dragDir (0.994 0.104 0); CofR (0 0 0); pitchAxis (0 0 1); magUInf 1.48; lRef 1; Aref 0.5; } } // ************************************************** *********************** // /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (1.47 0.155 0); boundaryField { AIRFOIL { type noSlip; } INLET { type fixedValue; value uniform (1.47 0.155 0); } OUTLET { type zeroGradient; } TOP { type fixedValue; value uniform (1.47 0.155 0); } BOTTOM { type fixedValue; value uniform (1.47 0.155 0); } frontAndBackPlanes { type empty; } } // ************************************************** *********************** // |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
You Wu
Join Date: Nov 2018
Posts: 16
Rep Power: 8 ![]() |
Hi prakashPra,
Your problem may be caused by your grid. I think there may be something wrong with your force calculation; Force coefficients are calculated by: force/(0.5*rho*Aref*U^2) However, the "force" is NOT calculated according to "Aref" in controlDict. It is decided according to your mesh. So even if the Aref is set to be the same as the experimental wing area, wrong forces will result in incorrect force coefficients. You should pay attention to the thickness of your mesh. In 2D simulation, Your 2D mesh will be extruded into 3D mesh (the division in 3rd direction is 1). You should set the extruding length equal to the wing span, then you could get right forces and force coefficients. |
|
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Local spikes in force coefficients | sjsivabharathy | OpenFOAM Running, Solving & CFD | 0 | May 1, 2019 15:14 |
using chemkin | JMDag2004 | OpenFOAM Pre-Processing | 2 | March 8, 2016 23:38 |
Sectional force coefficients | tfuwa | OpenFOAM Post-Processing | 1 | June 27, 2013 07:49 |
[OpenFOAM] spatial distribution of force coefficients in paraview&openfoam | kkpal | ParaView | 1 | April 21, 2013 08:15 |
Force can not converge | colopolo | CFX | 13 | October 4, 2011 23:03 |