April 23, 2022, 05:22
|
wrong forceCoeffs for standart airFoil2D simpleFoam case
|
#1
|
|
Member
Join Date: Sep 2020
Posts: 30
Rep Power: 7
|
Hello, FOAMers. I try to calculate Cd и Cl coeffs for airFoil2D using simpleFoam (standard case) OpenFOAM v1912. I'm based on this thread: Problem in simpleFoam airFoil2D forceCoeffs.
controlDict data:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1912 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application simpleFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 500;
deltaT 1;
writeControl timeStep;
writeInterval 50;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
functions
{
forces
{
type forceCoeffs;
libs ( "libforces.so" );
// Store and write volume field representations of forces and moments
writeFields yes;
// Centre of rotation for moment calculations
CofR (0 0 0);
patches ( walls );
p p;
U U;
rho rhoInf;
log true;
rhoInf 1.225;
liftDir (-0.139173 0.990268 0);
dragDir (0.990268 0.139173 0);
pitchAxis (0 0 1);
magUInf 26;
lRef 1;
Aref 1;
}
}
// ************************************************************************* //
So I get other drag and lift coeffs (Cd = 0514466, Cl = 1.7015) as opposed to the above thread.
Code:
SIMPLE solution converged in 313 iterations
forceCoeffs forces execute:
Coefficients
Cd : 0.0514466 (pressure: 0.0433982 viscous: 0.00804843)
Cs : 5.33394e-17 (pressure: 5.38832e-17 viscous: -5.43772e-19)
Cl : 1.7015 (pressure: 1.70157 viscous: -6.34881e-05)
CmRoll : -0.0425376 (pressure: -0.0425392 viscous: 1.5872e-06)
CmPitch : 98.4881 (pressure: 98.4116 viscous: 0.0764274)
CmYaw : 0.00128616 (pressure: 0.00108495 viscous: 0.000201211)
Cd(f) : -0.0168143
Cd(r) : 0.0682609
Cs(f) : 0.00128616
Cs(r) : -0.00128616
Cl(f) : 99.3388
Cl(r) : -97.6373
End
I understand that a lot of time has passed from 2012 year and methods of calcutating could be change. I have tried this standard case with help OpenFOAMv9 and got similar results.
Could you give me advise? Why do results differ from old ones? Tnx in advanсe.
Last edited by Reptider; April 23, 2022 at 08:54.
|
|
|