|
[Sponsors] |
Foam exciting while Force calculation in OF 1.6 |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Member
Devesh Baghel
Join Date: Mar 2009
Posts: 84
Rep Power: 18 ![]() |
Hi Foamers,
I am trying to simulate flow over airfoil by using Spalat-Allmaras model. I want to calculate Pressure & Viscous forces to find out Lift & Drag coefficients. Below are the description of "controlDict file" contains...... /************************************************** ****************** application simpleFoam; startFrom startTime; startTime 0; stopAt endTime; endTime 10; deltaT 5.0e-2; writeControl timeStep; writeInterval 100; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression uncompressed; timeFormat general; timePrecision 6; runTimeModifiable yes; functions ( forces { type forces; functionObjectLibs ("libforces.so"); //Lib to load outputControl outputTime; patches (wall); //change to your patch name pName p; // name of the field UName U; // name of the field log true; //Dump to file rhoInf 1.225; //Reference density for fluid CofR (0 0 0); //Origin for moment calculations } forceCoeffs { type forceCoeffs; functionObjectLibs ("libforces.so"); outputControl outputTime; patches (wall); //change to your patch name pName p; // name of the field UName U; // name of the field log true; //Dump to file rhoInf 1.225; //density CofR (0 0 0); liftDir (0 1 0); //Direction for lift dragDir (1 0 0); //Direction for Drag pitchAxis (0 0 0); //Pitching axis magUInf 26.0; //Free stream velocity magnitude lRef 0.6; //reference length "frontal length at velocity inlet" Aref 1.632653; //reference area "frontal area at velocity inlet" } ); // ************************************************** *********************** // while running the simpleFoam.....,after some time I get the message shown below " A component used as a directory in "/data/OpenFOAM/user1-1.6/Experiments_by_ourself/Tutorial_solved/Airfoil/SimpleFoam/forceCoeffs/5" is not, in fact, a directory. From function Foam::mkDir(const fileName&, mode_t) in file POSIX.C at line 372. FOAM exiting...................................... " Now, I got POSIX.C file, where at line 372, below codes were written..... case ENOTDIR: { FatalErrorIn("Foam::mkDir(const fileName&, mode_t)") << "A component used as a directory in " << pathName << " is not, in fact, a directory." << exit(FatalError); return false; } I din't find any solution/hint to resolve it ![]() Please help on this problem...... Thanks alot |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Devesh Baghel
Join Date: Mar 2009
Posts: 84
Rep Power: 18 ![]() |
Hi all,
I got the solution for the problem, I pasted there. ![]() thanks |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Tamas Regert
Join Date: Jan 2010
Posts: 3
Rep Power: 17 ![]() |
Hello!
I'm a newcommer to openFOAM and I would like to ask how did you manage to apply the forces function that you used here for a "type wall"? For me it works only for "type patch" boundaries. I have also an airfoil and using also Spalart-Allmaras model, thus wall BC was inevitable. Sorry if I ask some too obvious question and thank you for your help in advance. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Member
Devesh Baghel
Join Date: Mar 2009
Posts: 84
Rep Power: 18 ![]() |
Hi,
Have you defined the exact name of "wall" (as in Boundary File in constant/polymesh folder) in the above function? Other than that, I did not change anything else. |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Tamas Regert
Join Date: Jan 2010
Posts: 3
Rep Power: 17 ![]() |
Hello,
thank you for your answer. Finally I realised that the type 'wall' is accepted by the forces function as a patch, just I made a mistake when I defined the patches. As I have multiple wall faces (airfoil pressure and suction side), I had to write: patches(Airfoil_pressureside Airfoil_Suctionside); but in the beginning I made a mistake because I put a comma between the individual surfaces. This way here above works fine. Best regards, Tamas |
|
![]() |
![]() |
![]() |
![]() |
#6 | |
Senior Member
Sandy Lee
Join Date: Mar 2009
Posts: 213
Rep Power: 19 ![]() |
Quote:
Hi Devesh, Could you help me to solve the similar problem? http://www.cfd-online.com/Forums/ope...e-posix-c.html Thanks. Sandy |
||
![]() |
![]() |
![]() |
![]() |
#7 |
Member
Join Date: Nov 2009
Location: Germany
Posts: 96
Rep Power: 17 ![]() |
||
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[blockMesh] BlockMesh FOAM warning | gaottino | OpenFOAM Meshing & Mesh Conversion | 7 | July 19, 2010 14:11 |
Problem with rhoSimpleFoam | matteo_gautero | OpenFOAM Running, Solving & CFD | 0 | February 28, 2008 06:51 |
[blockMesh] Axisymmetrical mesh | Rasmus Gjesing (Gjesing) | OpenFOAM Meshing & Mesh Conversion | 10 | April 2, 2007 14:00 |
[Gmsh] Import gmsh msh to Foam | adorean | OpenFOAM Meshing & Mesh Conversion | 24 | April 27, 2005 08:19 |
Warning 097- | AB | Siemens | 6 | November 15, 2004 04:41 |