CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Airfoil Solver Problem

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 16, 2014, 11:59
Default Airfoil Solver Problem
  #1
New Member
 
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 11
boundary is on a distinguished road
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.
boundary is offline   Reply With Quote

Old   August 17, 2014, 15:26
Default quick reply
  #2
New Member
 
Jim Conger
Join Date: Jan 2014
Location: California, USA
Posts: 21
Rep Power: 12
james.conger is on a distinguished road
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.
james.conger is offline   Reply With Quote

Old   August 18, 2014, 05:57
Default
  #3
New Member
 
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 11
boundary is on a distinguished road
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
(
);

// ************************************************** *********************** //
Attached Files
File Type: gz 0.tar.gz (1.1 KB, 0 views)
File Type: gz system.tar.gz (1.4 KB, 0 views)
boundary is offline   Reply With Quote

Old   August 18, 2014, 05:58
Default
  #4
New Member
 
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 11
boundary is on a distinguished road
sorry for double posting ... how can i delete a post ?
boundary is offline   Reply With Quote

Old   August 18, 2014, 17:02
Default Another approach
  #5
New Member
 
Jim Conger
Join Date: Jan 2014
Location: California, USA
Posts: 21
Rep Power: 12
james.conger is on a distinguished road
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.
james.conger is offline   Reply With Quote

Old   August 19, 2014, 07:37
Default
  #6
New Member
 
Adam Apfel
Join Date: Aug 2014
Posts: 4
Rep Power: 11
boundary is on a distinguished road
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.
boundary is offline   Reply With Quote

Old   August 19, 2014, 15:42
Default
  #7
New Member
 
Simon
Join Date: Jun 2014
Posts: 9
Rep Power: 11
momentum-analysis is on a distinguished road
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.
momentum-analysis is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09: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 04:36
[GAMBIT] Meshing airfoil using .dat file problem creggie ANSYS Meshing & Geometry 10 June 27, 2010 19:24
Plunging Airfoil Problem jmo OpenFOAM Running, Solving & CFD 0 May 17, 2010 22:36
Coupled solver energy equation problem lucioantonio FLUENT 0 April 3, 2009 10:21


All times are GMT -4. The time now is 12:47.