Airfoil Solver Problem
 

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.

quick reply
 

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.

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
(
);

// ************************************************** *********************** //

sorry for double posting ... how can i delete a post ?

Another approach
 

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.

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.

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.