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

PotentialFoam few problems

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 16, 2008, 06:54
Default Hallo, I'm trying to run si
  #1
New Member
 
Chris
Join Date: Mar 2009
Location: Europe
Posts: 19
Rep Power: 8
chris_sev is on a distinguished road
Hallo,

I'm trying to run simulation in openFoam. It's an external flow around particular shape. I used fluentMeshToFoam to convert mesh, which has been done in HyperMesh. Firstly I dealed with potentialFoam solver (according to other users advices which I found in this forum).
I performed simulation but results are wrong. When I run ParaView to check velocity distribution nothing seems to happened there. Whole model is blue which means velocity is zero everywhere. I'm wondering what I did wrong. I went through all dictionary files once again but In my opinion everything is fine. I applied boundary conditions correctly. I'm not pretty sure how solution schemes work (I've just began playing with OF) so I have considered different scenarios but always with the same results.

Thought there might be something wrong with the mesh. I have run checkMesh but I'm not sure what's the explanation of those results. Could you please have a look on it:

--------------------------------------------------------------------
Create polyMesh for time = constant

Time = constant

Mesh stats
points: 373149
edges: 2289739
faces: 3709523
internal faces: 3590683
cells: 1792916
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0

Number of cells of each type:
hexahedra: 0
prisms: 124346
wedges: 0
pyramids: 4196
tet wedges: 0
tetrahedra: 1664374
polyhedra: 0

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface
inlet 260 154 ok (not multiply connected)
outlet 260 154 ok (not multiply connected)
wall_side 720 418 ok (not multiply connected)
wall_water 20152 10270 ok (not multiply connected)
wall_top 1156 610 ok (not multiply connected)
wall_ship 96292 49038 ok (not multiply connected)

Checking geometry...
Domain bounding box: (-1.82397e+06 -1.32002e+06 8294.26) (1.80603e+06 1.3199
8e+06 593733)
Boundary openness (9.88848e-18 9.87072e-18 1.38914e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 110560, number of cells 1
<<Writing 1 cells with high aspect ratio to set highAspectRatioCells
Minumum face area = 361.351. Maximum face area = 4.56545e+10. Face area mag nitudes OK.
Min volume = 212.664. Max volume = 3.14113e+15. Total volume = 5.60991e+18. Cell volumes OK.
Mesh non-orthogonality Max: 87.7938 average: 21.4714
*Number of severely non-orthogonal faces: 9750.
Non-orthogonality check OK.
<<Writing 9750 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
***Max skewness = 10.2348, 2 highly skew faces detected which may impair the qu ality of the results
<<Writing 2 skew faces to set skewFaces
Min/max edge length = 6.34815 377938 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : average = 0.997814 min = 0.809597
All face flatness OK.

Failed 2 mesh checks.

End
-------------------------------------------------------------------------------- ---


It says 'Failed 2 mesh checks.' It founds 1 cell with high aspect ratio 110560. This value is definitely to high. I've checked this mesh in HyperMesh and for this particular cell value of max aspect ratio is 1.8, so what's the reason of the discrepancy between OF and HM?

There are also two faces with highly skew. And again there's nothing wrong with those faces in HM.

What do you think about this? Do you think the fact that 'Failed 2 mesh checks.' could be the reason of my problems with results? In my opinion I should obtain some velocity distribution, probably wrong in areas which were pointed by checkMesh, but I shouldn't get zero everywhere.

thanks in advance!

Chris
chris_sev is offline   Reply With Quote

Old   July 22, 2008, 10:15
Default hi, I messed with initial v
  #2
New Member
 
Chris
Join Date: Mar 2009
Location: Europe
Posts: 19
Rep Power: 8
chris_sev is on a distinguished road
hi,

I messed with initial values. I found a problem and performed few more simulations. Now, velocity distribution looks much more better but values of maximum velocity fluctuates between 569 and 1260 m/s (depending on pressure solver type), when initial velocity is 19.97 m/s. Lowest values (569) I have obtained for smoothSolver with GaussSeidel smoother. But honestly I have no idea what else could be improved. I really must to work with this mesh. I didn't met any problems dealing with this in CFX. What do you think about this, is there any way in OpenFOAM to handle with my mesh and my case? Maybe I should exploit the potential of smoothSolver, but how?
I would be very appreciate for any hints.!!!
thank you!!!

chris
chris_sev is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
POTENTIALFOAM ERROR for a naca profile dinonettis OpenFOAM Running, Solving & CFD 7 September 9, 2010 10:59
PotentialFOAM for channel flow albcem OpenFOAM Running, Solving & CFD 0 October 16, 2008 12:18
POTENTIALFOAM ERROR for a naca profile dinonettis OpenFOAM Running, Solving & CFD 2 April 10, 2008 12:17
Using potentialFoam with simpleBuoyantFoam lasb OpenFOAM Running, Solving & CFD 4 August 10, 2007 08:31
A fundamental problem about Pressure equation of the potentialFoam solver dbxmcf OpenFOAM Running, Solving & CFD 0 October 6, 2006 11:32


All times are GMT -4. The time now is 17:37.