Hi All,
I need some help wi
Hi All,
I need some help with simpleFoam. Here is my problem... I am simulating airflow through the following structure. http://www.cfdonline.com/OpenFOAM_D...ges/1/9078.jpg 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. http://www.cfdonline.com/OpenFOAM_D...ges/1/9079.jpg 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 
Hi All,
Below is the output
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 zipup check OK. Face vertices OK. Faceface 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.71605e17 2.01321e17 6.00016e17) OK. Max cell openness = 6.12672e16 OK. Max aspect ratio = 17.4782 OK. Minumum face area = 3.1167e10. Maximum face area = 1.90368e06. Face area magnitudes OK. Min volume = 7.24431e15. Max volume = 5.97908e10. Total volume = 4.3833e05. Cell volumes OK. Mesh nonorthogonality Max: 84.0588 average: 33.8839 *Number of severely nonorthogonal faces: 5045. Nonorthogonality check OK. <<Writing 5045 nonorthogonal faces to set nonOrthoFaces Face pyramids OK. Max skewness = 2.12424 OK. Min/max edge length = 2.42777e05 0.00284903 OK. All angles in faces OK. All face flatness OK. Mesh OK. End 
Hi All,
Is the Nonorthogon
Hi All,
Is the Nonorthogonality in the checkMesh the cause for the convergence problem? Thanks in Advance Senthil 
Senthil: that is a possibility
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 Gagnon,
Thanks for the s
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? http://www.cfdonline.com/OpenFOAM_D...ges/1/9141.jpg Thanks in advance Senthil 
I am not sure about this. Howe
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 
Quote:
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? 
All times are GMT 4. The time now is 23:56. 