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

SimpleFoam Convergence problem

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

Like Tree1Likes
  • 1 Post By louisgag

Reply
 
LinkBack Thread Tools Display Modes
Old   September 8, 2008, 01:27
Default Hi All, I need some help wi
  #1
Senior Member
 
Senthil Kabilan
Join Date: Mar 2009
Posts: 113
Rep Power: 8
skabilan is on a distinguished road
Hi All,

I need some help with simpleFoam. Here is my problem... I am simulating airflow through the following structure.




The boundary conditions include zero pressure at the outlet (largest cylinder) and the eight inlets are prescribed a constant pressure of 13 pa. The problem converges with out a problem.

I have another mesh which consists of only three cylinders (subset of the above problem) as shown below.



Once again the outlet (largest cylinder) is maintained at zero and the two inlets are prescribed a pressure of 6 pa (this was the average pressure at the the same locations obtained from the first simulation). The problem does not converge. The solver setting remain the same as the first simulation.

Thanks all in advance!
Senthil
skabilan is offline   Reply With Quote

Old   September 8, 2008, 12:53
Default Hi All, Below is the output
  #2
Senior Member
 
Senthil Kabilan
Join Date: Mar 2009
Posts: 113
Rep Power: 8
skabilan is on a distinguished road
Hi All,

Below is the output from checkMesh for the case that does not converge.

bigbox76% checkMesh . weibel_2gen_ss_pos6
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : checkMesh . weibel_2gen_ss_pos6
Date : Sep 08 2008
Time : 09:50:59
Host : bigbox
PID : 5555
Root : /files0/skabilan/uw_workdir/openfoam/weibel/weibel_pressure_simulation
Case : weibel_2gen_ss_pos6
Nprocs : 1
Create time

Create polyMesh for time = constant

Time = constant

Mesh stats
points: 135357
edges: 934993
faces: 1589196
internal faces: 1569040
cells: 789559
boundary patches: 4
point zones: 0
face zones: 0
cell zones: 0

Number of cells of each type:
hexahedra: 0
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 789559
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 1548 799 ok (not multiply connected)
out2 538 281 ok (not multiply connected)
out3 663 346 ok (not multiply connected)
w1 17407 8751 ok (not multiply connected)

Checking geometry...
Domain bounding box: (-0.0341182 -0.0678843 -0.00901863) (0.0348 0.106513 0.00904051)
Boundary openness (-4.71605e-17 2.01321e-17 -6.00016e-17) OK.
Max cell openness = 6.12672e-16 OK.
Max aspect ratio = 17.4782 OK.
Minumum face area = 3.1167e-10. Maximum face area = 1.90368e-06. Face area magnitudes OK.
Min volume = 7.24431e-15. Max volume = 5.97908e-10. Total volume = 4.3833e-05. Cell volumes OK.
Mesh non-orthogonality Max: 84.0588 average: 33.8839
*Number of severely non-orthogonal faces: 5045.
Non-orthogonality check OK.
<<Writing 5045 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 2.12424 OK.
Min/max edge length = 2.42777e-05 0.00284903 OK.
All angles in faces OK.
All face flatness OK.

Mesh OK.

End
skabilan is offline   Reply With Quote

Old   September 10, 2008, 03:40
Default Hi All, Is the Non-orthogon
  #3
Senior Member
 
Senthil Kabilan
Join Date: Mar 2009
Posts: 113
Rep Power: 8
skabilan is on a distinguished road
Hi All,

Is the Non-orthogonality in the checkMesh the cause for the convergence problem?

Thanks in Advance
Senthil
skabilan is offline   Reply With Quote

Old   September 10, 2008, 17:03
Default Senthil: that is a possibility
  #4
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Senthil: that is a possibility (84 is pretty high). Try increasing the nNonOrthogonalCorrectors value in you system/fvSolution dictionary. Or even better, try to make the mesh more orthogonal.
s.m likes this.
louisgag is offline   Reply With Quote

Old   September 11, 2008, 01:32
Default Hi Gagnon, Thanks for the s
  #5
Senior Member
 
Senthil Kabilan
Join Date: Mar 2009
Posts: 113
Rep Power: 8
skabilan is on a distinguished road
Hi Gagnon,

Thanks for the suggestion. The mesh did produce good results for a simple positive pressure simulation (i.e, positive pressure at the inlet and 0 pressures at the outlets). It looks like it might be a problem with the boundary type that I am specifying in /0/p and /0/U file for the following steadystate case. What boundary conditions needs to be specified for the following case?



Thanks in advance
Senthil
skabilan is offline   Reply With Quote

Old   September 12, 2008, 15:10
Default I am not sure about this. Howe
  #6
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
I am not sure about this. However, when I run my simulations I always prescribe a pressure on the outlet, a velocity on the inlet and have velocity on outet as zerogradient and pressure at inlet as zeroGradient.

Looking at your problem, I could suggest that you set either the inlet or outlet velocity at fixedValue and leave the other one as zeroGradient, but that's just me using my hunch.

good luck,

-Louis
louisgag is offline   Reply With Quote

Old   May 31, 2013, 03:21
Default
  #7
s.m
Senior Member
 
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 5
s.m is on a distinguished road
Quote:
Originally Posted by louisgag View Post
Senthil: that is a possibility (84 is pretty high). Try increasing the nNonOrthogonalCorrectors value in you system/fvSolution dictionary. Or even better, try to make the mesh more orthogonal.
hi louis
what do you mean "try to make the mesh more orthogonal" ?
how should we make the mesh more orthogonal, if we use SnappHexMesh for doing th mesh, we should increase the "maxNonOrtho 45;" e.g to 65 or else?
s.m 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
SimpleFoam convergence problems brahim OpenFOAM Running, Solving & CFD 20 June 9, 2015 09:09
Convergence problem using simpleFoam steady state vvqf OpenFOAM Running, Solving & CFD 12 May 18, 2011 07:51
SimpleFoam solution convergence pattern philippose OpenFOAM Running, Solving & CFD 0 June 26, 2008 14:18
SimpleFoam convergence problems schnitzlein OpenFOAM Running, Solving & CFD 6 June 24, 2005 09:51
Axial 2D calculation in simpleFoam prblem with convergence rybakov OpenFOAM Running, Solving & CFD 3 May 16, 2005 02:00


All times are GMT -4. The time now is 15:54.