|
[Sponsors] |
August 16, 2014, 12:59 |
Airfoil Solver Problem
|
#1 |
New Member
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 12 |
Hey I am very new with OpenFoam and I am trying to simmulate a airfoil and show the vortex.
checkMesh: Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 5292 internal points: 0 faces: 9840 internal faces: 4560 cells: 2400 faces per cell: 6 boundary patches: 7 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 2400 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. *Number of regions: 6 The mesh has multiple regions which are not connected by any face. <<Writing region information to "0/cellToRegion" <<Writing region 0 with 400 cells to cellSet region0 <<Writing region 1 with 400 cells to cellSet region1 <<Writing region 2 with 400 cells to cellSet region2 <<Writing region 3 with 400 cells to cellSet region3 <<Writing region 4 with 400 cells to cellSet region4 <<Writing region 5 with 400 cells to cellSet region5 Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology inlet 40 84 ok (non-closed singly connected) outlet 40 84 ok (non-closed singly connected) frontAndBack 4800 5292 ok (non-closed singly connected) top 60 126 ok (non-closed singly connected) bottom 60 126 ok (non-closed singly connected) flügel 40 84 ok (non-closed singly connected) defaultFaces 240 496 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (-40 -40 0) (41 40 0.1) Mesh (non-empty, non-wedge) directions (0 0 0) Mesh (non-empty) directions (0 0 0) ***Number of edges not aligned with or perpendicular to non-empty directions: 3521 <<Writing 2636 points on non-aligned edges to set nonAlignedEdges Boundary openness (7.14637e-20 -3.60722e-19 0) OK. Max cell openness = 1.04083e-15 OK. Max aspect ratio = -1 OK. Minimum face area = 0.00460531. Maximum face area = 82.0002. Face area magnitudes OK. Min volume = 0.00712701. Max volume = 0.4. Total volume = 647.992. Cell volumes OK. Mesh non-orthogonality Max: 89.8603 average: 23.7693 *Number of severely non-orthogonal (> 70 degrees) faces: 380. Non-orthogonality check OK. <<Writing 380 non-orthogonal faces to set nonOrthoFaces Face pyramids OK. ***Max skewness = 20.5, 800 highly skew faces detected which may impair the quality of the results <<Writing 800 skew faces to set skewFaces Coupled point location match (average 0) OK. Failed 2 mesh checks. End ---------------------------------------------------------------------------- And when i try the solver terminal says: Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type LESModel Selecting LES turbulence model oneEqEddy Selecting LES delta type cubeRootVol --> FOAM FATAL ERROR: Case is not 3D or 2D, LES is not applicable From function cubeRootVolDelta::calcDelta() in file cubeRootVolDelta/cubeRootVolDelta.C at line 72. FOAM exiting I would be very thankfoul if someone could help me. |
|
August 17, 2014, 16:26 |
quick reply
|
#2 |
New Member
Jim Conger
Join Date: Jan 2014
Location: California, USA
Posts: 21
Rep Power: 12 |
You don't specify how you created the mesh, so it is difficult to give advise. If you started with an STL model for the airfoil, then take a look at the MeshLab application. It has a number of filters that are very useful for fixing problems with the STL shape. Usually if you get the STL shape as a closed surface, then mesh creation with snappyHexMesh works fine.
|
|
August 18, 2014, 06:57 |
|
#3 |
New Member
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 12 |
i used blockMesh and LES Properties
: /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 1; vertices ( (0 0 0)//0 (1 0 0)//1 (1 40 0)//2 (0 40 0)//3 (0 0 0.1)//4 (1 0 0.1)//5 (1 40 0.1)//6 (0 40 0.1)//7 (0 -40 0)//8 (1 -40 0)//9 (1 0 0)//10 (0 0 0)//11 (0 -40 0.1)//12 (1 -40 0.1)//13 (1 0 0.1)//14 (0 0 0.1)//15 (-40 0 0)//16 (0 0 0)//17 (0 40 0)//18 (-40 40 0)//19 (-40 0 0.1)//20 (0 0 0.1)//21 (0 40 0.1)//22 (-40 40 0.1)//23 (1 0 0)//24 (41 0 0)//25 (41 40 0)//26 (1 40 0)//27 (1 0 0.1)//28 (41 0 0.1)//29 (41 40 0.1)//30 (1 40 0.1)//31 (-40 -40 0)//32 (0 -40 0)//33 (0 0 0)//34 (-40 0 0)//35 (-40 -40 0.1)//36 (0 -40 0.1)//37 (0 0 0.1)//38 (-40 0 0.1)//39 (1 -40 0)//40 (41 -40 0)//41 (41 0 0)//42 (1 0 0)//43 (-40 -40 0.1)//44 (0 -40 0.1)//45 (0 0 0.1)//46 (-40 0 0.1)//47 ); blocks ( hex (0 1 2 3 4 5 6 7) (20 20 1) simpleGrading (1 1 1)//1 hex (8 9 10 11 12 13 14 15) (20 20 1) simpleGrading (1 1 1)//2 hex (16 17 18 19 20 21 22 23) (20 20 1) simpleGrading (1 1 1)//3 hex (24 25 26 27 28 29 30 31) (20 20 1) simpleGrading (1 1 1)//4 hex (32 33 34 35 36 37 38 39) (20 20 1) simpleGrading (1 1 1)//5 hex (40 41 42 43 44 45 46 47) (20 20 1) simpleGrading (1 1 1)//6 ); edges ( polyLine 0 1 ( (0.0000000 0.0000000 0) (0.0005000 0.0023390 0) (0.0010000 0.0037271 0) (0.0020000 0.0058025 0) (0.0040000 0.0089238 0) (0.0080000 0.0137350 0) (0.0120000 0.0178581 0) (0.0200000 0.0253735 0) (0.0300000 0.0330215 0) (0.0400000 0.0391283 0) (0.0500000 0.0442753 0) (0.0600000 0.0487571 0) (0.0800000 0.0564308 0) (0.1000000 0.0629981 0) (0.1200000 0.0686204 0) (0.1400000 0.0734360 0) (0.1600000 0.0775707 0) (0.1800000 0.0810687 0) (0.2000000 0.0839202 0) (0.2200000 0.0861433 0) (0.2400000 0.0878308 0) (0.2600000 0.0890840 0) (0.2800000 0.0900016 0) (0.3000000 0.0906804 0) (0.3200000 0.0911857 0) (0.3400000 0.0915079 0) (0.3600000 0.0916266 0) (0.3800000 0.0915212 0) (0.4000000 0.0911712 0) (0.4200000 0.0905657 0) (0.4400000 0.0897175 0) (0.4600000 0.0886427 0) (0.4800000 0.0873572 0) (0.5000000 0.0858772 0) (0.5200000 0.0842145 0) (0.5400000 0.0823712 0) (0.5600000 0.0803480 0) (0.5800000 0.0781451 0) (0.6000000 0.0757633 0) (0.6200000 0.0732055 0) (0.6400000 0.0704822 0) (0.6600000 0.0676046 0) (0.6800000 0.0645843 0) (0.7000000 0.0614329 0) (0.7200000 0.0581599 0) (0.7400000 0.0547675 0) (0.7600000 0.0512565 0) (0.7800000 0.0476281 0) (0.8000000 0.0438836 0) (0.8200000 0.0400245 0) (0.8400000 0.0360536 0) (0.8600000 0.0319740 0) (0.8800000 0.0277891 0) (0.9000000 0.0235025 0) (0.9200000 0.0191156 0) (0.9400000 0.0146239 0) (0.9600000 0.0100232 0) (0.9700000 0.0076868 0) (0.9800000 0.0053335 0) (0.9900000 0.0029690 0) (1.0000000 0.0005993 0) ) polyLine 4 5 ( ( 0.0000000 0.0000000 0.1) ( 0.0005000 0.0023390 0.1) ( 0.0010000 0.0037271 0.1) ( 0.0020000 0.0058025 0.1) ( 0.0040000 0.0089238 0.1) ( 0.0080000 0.0137350 0.1) ( 0.0120000 0.0178581 0.1) ( 0.0200000 0.0253735 0.1) ( 0.0300000 0.0330215 0.1) ( 0.0400000 0.0391283 0.1) ( 0.0500000 0.0442753 0.1) ( 0.0600000 0.0487571 0.1) ( 0.0800000 0.0564308 0.1) ( 0.1000000 0.0629981 0.1) ( 0.1200000 0.0686204 0.1) ( 0.1400000 0.0734360 0.1) ( 0.1600000 0.0775707 0.1) ( 0.1800000 0.0810687 0.1) ( 0.2000000 0.0839202 0.1) ( 0.2200000 0.0861433 0.1) ( 0.2400000 0.0878308 0.1) ( 0.2600000 0.0890840 0.1) ( 0.2800000 0.0900016 0.1) ( 0.3000000 0.0906804 0.1) ( 0.3200000 0.0911857 0.1) ( 0.3400000 0.0915079 0.1) ( 0.3600000 0.0916266 0.1) ( 0.3800000 0.0915212 0.1) ( 0.4000000 0.0911712 0.1) ( 0.4200000 0.0905657 0.1) ( 0.4400000 0.0897175 0.1) ( 0.4600000 0.0886427 0.1) ( 0.4800000 0.0873572 0.1) ( 0.5000000 0.0858772 0.1) ( 0.5200000 0.0842145 0.1) ( 0.5400000 0.0823712 0.1) ( 0.5600000 0.0803480 0.1) ( 0.5800000 0.0781451 0.1) ( 0.6000000 0.0757633 0.1) ( 0.6200000 0.0732055 0.1) ( 0.6400000 0.0704822 0.1) ( 0.6600000 0.0676046 0.1) ( 0.6800000 0.0645843 0.1) ( 0.7000000 0.0614329 0.1) ( 0.7200000 0.0581599 0.1) ( 0.7400000 0.0547675 0.1) ( 0.7600000 0.0512565 0.1) ( 0.7800000 0.0476281 0.1) ( 0.8000000 0.0438836 0.1) ( 0.8200000 0.0400245 0.1) ( 0.8400000 0.0360536 0.1) ( 0.8600000 0.0319740 0.1) ( 0.8800000 0.0277891 0.1) ( 0.9000000 0.0235025 0.1) ( 0.9200000 0.0191156 0.1) ( 0.9400000 0.0146239 0.1) ( 0.9600000 0.0100232 0.1) ( 0.9700000 0.0076868 0.1) ( 0.9800000 0.0053335 0.1) ( 0.9900000 0.0029690 0.1) ( 1.0000000 0.0005993 0.1) ) polyLine 11 10 ( ( 0.0000000 0 0) ( 0.0005000 -.0046700 0) ( 0.0010000 -.0059418 0) ( 0.0020000 -.0078113 0) ( 0.0040000 -.0105126 0) ( 0.0080000 -.0142862 0) ( 0.0120000 -.0169733 0) ( 0.0200000 -.0202723 0) ( 0.0300000 -.0226056 0) ( 0.0400000 -.0245211 0) ( 0.0500000 -.0260452 0) ( 0.0600000 -.0271277 0) ( 0.0800000 -.0284595 0) ( 0.1000000 -.0293786 0) ( 0.1200000 -.0299633 0) ( 0.1400000 -.0302404 0) ( 0.1600000 -.0302546 0) ( 0.1800000 -.0300490 0) ( 0.2000000 -.0296656 0) ( 0.2200000 -.0291445 0) ( 0.2400000 -.0285181 0) ( 0.2600000 -.0278164 0) ( 0.2800000 -.0270696 0) ( 0.3000000 -.0263079 0) ( 0.3200000 -.0255565 0) ( 0.3400000 -.0248176 0) ( 0.3600000 -.0240870 0) ( 0.3800000 -.0233606 0) ( 0.4000000 -.0226341 0) ( 0.4200000 -.0219042 0) ( 0.4400000 -.0211708 0) ( 0.4600000 -.0204353 0) ( 0.4800000 -.0196986 0) ( 0.5000000 -.0189619 0) ( 0.5200000 -.0182262 0) ( 0.5400000 -.0174914 0) ( 0.5600000 -.0167572 0) ( 0.5800000 -.0160232 0) ( 0.6000000 -.0152893 0) ( 0.6200000 -.0145551 0) ( 0.6400000 -.0138207 0) ( 0.6600000 -.0130862 0) ( 0.6800000 -.0123515 0) ( 0.7000000 -.0116169 0) ( 0.7200000 -.0108823 0) ( 0.7400000 -.0101478 0) ( 0.7600000 -.0094133 0) ( 0.7800000 -.0086788 0) ( 0.8000000 -.0079443 0) ( 0.8200000 -.0072098 0) ( 0.8400000 -.0064753 0) ( 0.8600000 -.0057408 0) ( 0.8800000 -.0050063 0) ( 0.9000000 -.0042718 0) ( 0.9200000 -.0035373 0) ( 0.9400000 -.0028028 0) ( 0.9600000 -.0020683 0) ( 0.9700000 -.0017011 0) ( 0.9800000 -.0013339 0) ( 0.9900000 -.0009666 0) ( 1.0000000 -.0005993 0) ) polyLine 15 14 ( ( 0.0000000 0 0.1) ( 0.0005000 -.0046700 0.1) ( 0.0010000 -.0059418 0.1) ( 0.0020000 -.0078113 0.1) ( 0.0040000 -.0105126 0.1) ( 0.0080000 -.0142862 0.1) ( 0.0120000 -.0169733 0.1) ( 0.0200000 -.0202723 0.1) ( 0.0300000 -.0226056 0.1) ( 0.0400000 -.0245211 0.1) ( 0.0500000 -.0260452 0.1) ( 0.0600000 -.0271277 0.1) ( 0.0800000 -.0284595 0.1) ( 0.1000000 -.0293786 0.1) ( 0.1200000 -.0299633 0.1) ( 0.1400000 -.0302404 0.1) ( 0.1600000 -.0302546 0.1) ( 0.1800000 -.0300490 0.1) ( 0.2000000 -.0296656 0.1) ( 0.2200000 -.0291445 0.1) ( 0.2400000 -.0285181 0.1) ( 0.2600000 -.0278164 0.1) ( 0.2800000 -.0270696 0.1) ( 0.3000000 -.0263079 0.1) ( 0.3200000 -.0255565 0.1) ( 0.3400000 -.0248176 0.1) ( 0.3600000 -.0240870 0.1) ( 0.3800000 -.0233606 0.1) ( 0.4000000 -.0226341 0.1) ( 0.4200000 -.0219042 0.1) ( 0.4400000 -.0211708 0.1) ( 0.4600000 -.0204353 0.1) ( 0.4800000 -.0196986 0.1) ( 0.5000000 -.0189619 0.1) ( 0.5200000 -.0182262 0.1) ( 0.5400000 -.0174914 0.1) ( 0.5600000 -.0167572 0.1) ( 0.5800000 -.0160232 0.1) ( 0.6000000 -.0152893 0.1) ( 0.6200000 -.0145551 0.1) ( 0.6400000 -.0138207 0.1) ( 0.6600000 -.0130862 0.1) ( 0.6800000 -.0123515 0.1) ( 0.7000000 -.0116169 0.1) ( 0.7200000 -.0108823 0.1) ( 0.7400000 -.0101478 0.1) ( 0.7600000 -.0094133 0.1) ( 0.7800000 -.0086788 0.1) ( 0.8000000 -.0079443 0.1) ( 0.8200000 -.0072098 0.1) ( 0.8400000 -.0064753 0.1) ( 0.8600000 -.0057408 0.1) ( 0.8800000 -.0050063 0.1) ( 0.9000000 -.0042718 0.1) ( 0.9200000 -.0035373 0.1) ( 0.9400000 -.0028028 0.1) ( 0.9600000 -.0020683 0.1) ( 0.9700000 -.0017011 0.1) ( 0.9800000 -.0013339 0.1) ( 0.9900000 -.0009666 0.1) ( 1.0000000 -.0005993 0.1) ) ); boundary ( inlet { type patch; faces ( (19 23 20 16)//3 (32 35 39 36)//5 ); } outlet { type patch; faces ( (25 26 30 29)//4 (41 42 46 45)//6 ); } frontAndBack { type empty; faces ( (0 1 2 3)//1 (4 5 6 7)//1 (8 9 10 11)//2 (12 13 14 15)//2 (20 21 22 23)//3 (16 17 18 19)//3 (24 25 26 27)//4 (28 29 30 31)//4 (32 33 34 35)//5 (36 37 38 39)//5 (40 41 42 43)//6 (44 45 46 47)//6 ); } top { type symmetry; faces ( (3 7 6 2)//1 (19 18 22 23)//3 (27 26 30 31)//4 ); } bottom { type symmetry; faces ( (8 9 13 12)//2 (32 33 37 36)//5 (40 41 45 44)//6 ); } flügel { type wall; faces ( (0 4 5 1)//1 (15 14 10 11)//2 ); } ); mergePatchPairs ( ); // ************************************************** *********************** // |
|
August 18, 2014, 06:58 |
|
#4 |
New Member
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 12 |
sorry for double posting ... how can i delete a post ?
|
|
August 18, 2014, 18:02 |
Another approach
|
#5 |
New Member
Jim Conger
Join Date: Jan 2014
Location: California, USA
Posts: 21
Rep Power: 12 |
In my first work with OpenFOAM I attempted to do what you are doing and create the entire problem with blockMesh. This can be done, but is difficult and error prone. The OpenFOAM tools are designed to work with another approach:
1) Use blockMesh only to create the external grid that will contain the problem. 2) Create an STL model of the objects in the grid 3) Use snappyHexMesh to insert the objects into the mesh, including local refinement. It takes a small investment of time to switch approaches, but life will be a lot simpler after that. Look for examples in the distribution with snappyHexMesh running. |
|
August 19, 2014, 08:37 |
|
#6 |
New Member
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 12 |
The Problem is that I have to do it with blockMesh. Actually it is just a studywork with 20 h time expenses. But current i have invested much more time and I only want to finish it.
|
|
August 19, 2014, 16:42 |
|
#7 |
New Member
Simon
Join Date: Jun 2014
Posts: 9
Rep Power: 12 |
Hi,
I don't think you have the mesh that you want - I copied your blockMeshDict, ran it then looked in paraview at the mesh. Its tricky to do this sort of thing in blockMesh - previous suggestion of using snappyHexMesh is a good one, you should find guidelines on generating a 2D mesh from snappy quite easily. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
thobois class engineTopoChangerMesh error | Peter_600 | OpenFOAM | 4 | August 2, 2014 10:52 |
[Commercial meshers] Handling cyclic BC from gambit to openfoam for a cascade airfoil problem - OF 1.6 | maverick | OpenFOAM Meshing & Mesh Conversion | 2 | June 18, 2011 05:36 |
[GAMBIT] Meshing airfoil using .dat file problem | creggie | ANSYS Meshing & Geometry | 10 | June 27, 2010 20:24 |
Plunging Airfoil Problem | jmo | OpenFOAM Running, Solving & CFD | 0 | May 17, 2010 23:36 |
Coupled solver energy equation problem | lucioantonio | FLUENT | 0 | April 3, 2009 11:21 |