Hot News:

Unser Angebot:

  Foren auf CAD.de (alle Foren)
  OpenFOAM
  Probleme bei addlayer -snappyHexMesh

Antwort erstellen  Neues Thema erstellen
CAD.de Login | Logout | Profil | Profil bearbeiten | Registrieren | Voreinstellungen | Hilfe | Suchen

Anzeige:

Darstellung des Themas zum Ausdrucken. Bitte dann die Druckfunktion des Browsers verwenden. | Suche nach Beiträgen nächster neuer Beitrag | nächster älterer Beitrag
Autor Thema:  Probleme bei addlayer -snappyHexMesh (3040 mal gelesen)
shake2014
Mitglied


Sehen Sie sich das Profil von shake2014 an!   Senden Sie eine Private Message an shake2014  Schreiben Sie einen Gästebucheintrag für shake2014

Beiträge: 3
Registriert: 07.09.2014

OpenFoam 2.3

erstellt am: 09. Sep. 2014 21:45    Editieren oder löschen Sie diesen Beitrag!  <-- editieren / zitieren -->   Antwort mit Zitat in Fett Antwort mit kursivem Zitat    Unities abgeben: 1 Unity (wenig hilfreich, aber dennoch)2 Unities3 Unities4 Unities5 Unities6 Unities7 Unities8 Unities9 Unities10 Unities

Hallo,

ich versuche eine relativ einfache Geometrie mit shm zu vernetzen. Leider hat es große Probleme beim erstellen von Layern. mein bestes Resultat benötigt 8mio Zellen (s. logfile). trotz dieser großen Zellanzahl ist das Ergebnis unbefriedigend.

Falls jmd. Interesse hat drüber zu schauen, ist das Projekt unter folgendem Link zu erreichen:

https://www.dropbox.com/s/3toe8e4cr5qtaza/dreiD.tar.gz?dl=0

Ich hab fast alle Parameter verändert und ausprobiert. es muss doch möglich sein, dass sHM es schafft, so schwach kann sHM doch nicht sein.

checkMesh meldet keine Fehler.

Code:
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:          9158486
    faces:            26410059
    internal faces:  25855510
    cells:            8627266
    faces per cell:  6.05818
    boundary patches: 9
    point zones:      0
    face zones:      0
    cell zones:      0

Overall number of cells of each type:
    hexahedra:    8409468
    prisms:        29463
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:    188335
    Breakdown of polyhedra by number of faces:
        faces  number of cells
            4  1949
            5  1264
            6  17568
            7  7284
            8  1224
            9  144133
          12  13261
          15  1641
          18  11

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
                  Patch    Faces  Points                  Surface topology
                  front    22880    23817  ok (non-closed singly connected)
                    back    22880    23817  ok (non-closed singly connected)
                  inlet    10400    10857  ok (non-closed singly connected)
                  outlet    10400    10857  ok (non-closed singly connected)
                  ground  352000  353281  ok (non-closed singly connected)
                    top    22000    22321  ok (non-closed singly connected)
                    oben    49532    50016  ok (non-closed singly connected)
                schraege    14924    17460  ok (non-closed singly connected)
                  unten    49533    50018  ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (-90 -50 -16) (130 50 24)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-1.256e-15 -8.33717e-15 -1.218e-15) OK.
    Max cell openness = 2.53573e-15 OK.
    Max aspect ratio = 83.35 OK.
    Minimum face area = 0.000221103. Maximum face area = 1.02978.  Face area magnitudes OK.
    Min volume = 1.46709e-05. Max volume = 1.02174.  Total volume = 879211.  Cell volumes OK.
    Mesh non-orthogonality Max: 66.5288 average: 4.47767
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.07252 OK.
    Coupled point location match (average 0) OK.

Mesh OK.


snappyHexMesh

Code:
Create time

Create mesh for time = 0

Read mesh in = 16.38 s

Overall mesh bounding box  : (-90 -50 -16) (130 50 24)
Relative tolerance        : 1e-06
Absolute matching distance : 0.000244949

Reading refinement surfaces.
Read refinement surfaces in = 1.72 s

Reading refinement shells.
Refinement level 3 for all cells inside refinementBoxGround
Read refinement shells in = 0 s

Setting refinement level of surface to be consistent with shells.
Checked shell refinement in = 2.06 s

Reading features.
Read extendedFeatureEdgeMesh schraege.extendedFeatureEdgeMesh
    points      : 7076
    edges      : 7076
    boundingBox : (-14.5613 -15.075 -4.13438) (14.5613 15.075 4.13438)
    point classification :
        convex feature points          :        0
        concave feature points        :        4
        mixed feature points          :        8
        other (non-feature) points    :    7064
    edge classification :
        external (convex angle) edges  :        0
        internal (concave angle) edges :    2540
        flat region edges              :        0
        open edges                    :    4536
        multiply connected edges      :        0

Refinement level according to distance to "schraege.extendedFeatureEdgeMesh" (7076 points, 7076 edges).
    level 0 for all cells within 0 metre.
Read features in = 0.03 s


Determining initial surface intersections
-----------------------------------------

Edge intersection testing:
    Number of edges            : 12739207
    Number of edges to retest  : 12739207
    Number of intersected edges : 54885
Calculated surface intersections in = 23.09 s

Initial mesh : cells:4076033  faces:12739207  points:4590070
Cells per refinement level:
    0 772082
    1 630770
    2 1686778
    3 986403

Adding patches for surface regions
----------------------------------

Patch Type                Region                       
----- ----                ------                       
membran.stl:

6    wall                oben                         
7    wall                schraege                     
8    wall                unten                       

Added patches in = 0.04 s

Selecting decompositionMethod none
Layer thickness specified as final layer and expansion ratio.

Shrinking and layer addition phase
----------------------------------

Using mesh parameters
{
    maxNonOrtho    75;
    maxBoundarySkewness 20;
    maxInternalSkewness 4;
    maxConcave      80;
    minVol          1e-13;
    minTetQuality  1e-15;
    minArea        -1;
    minTwist        0.02;
    minDeterminant  0.001;
    minFaceWeight  0.02;
    minVolRatio    0.01;
    minTriangleTwist -1;
    nSmoothScale    4;
    errorReduction  0.75;
}

Merging all faces of a cell
---------------------------
    - which are on the same patch
    - which make an angle < 65 degrees
      (cos:0.422618)
    - as long as the resulting face doesn't become concave by more than 90 degrees
      (0=straight, 180=fully concave)

Merging 16 sets of faces.
Setting up information for layer truncation ...

Checking mesh with layer ...
Checking faces in error :
    non-orthogonality > 75  degrees                        : 4
    faces with face pyramid volume < 1e-13                : 0
    faces with face-decomposition tet quality < 1e-15      : 9625
    faces with concavity > 80  degrees                    : 0
    faces with skewness > 4  (internal) or 20  (boundary) : 0
    faces with interpolation weights (0..1)  < 0.02        : 0
    faces with volume ratio of neighbour cells < 0.01      : 0
    faces with face twist < 0.02                          : 0
    faces on cells with determinant < 0.001                : 0
Detected 9629 illegal faces (concave, zero area or negative cell pyramid volume)
Extruding 464168 out of 465989 faces (99.6092%). Removed extrusion at 1435 faces.
Added 4659890 out of 4659890 cells (100%).

patch    faces    layers  overall thickness
                          [m]      [%]
-----    -----    ------  ---      ---
ground  352000  10      0.0609    100   
oben    49532    9.69    0.0477    99.7   
schraege 14924    4.81    0.0258    53.9   
unten    49533    9.69    0.0477    99.7 


selbst wenn ich tet quality auf -1e-30 setze, erhalte ich keine verbesserung.

Spart nicht mit euren Tipps. Vielen Dank für die Mühe.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

Shor-ty
Moderator





Sehen Sie sich das Profil von Shor-ty an!   Senden Sie eine Private Message an Shor-ty  Schreiben Sie einen Gästebucheintrag für Shor-ty

Beiträge: 2463
Registriert: 27.08.2010

OpenFOAM-dev (Foundation)
OpenFOAM-xxxx (ESI)

erstellt am: 10. Sep. 2014 16:10    Editieren oder löschen Sie diesen Beitrag!  <-- editieren / zitieren -->   Antwort mit Zitat in Fett Antwort mit kursivem Zitat    Unities abgeben: 1 Unity (wenig hilfreich, aber dennoch)2 Unities3 Unities4 Unities5 Unities6 Unities7 Unities8 Unities9 Unities10 Unities Nur für shake2014 10 Unities + Antwort hilfreich

Hi,

das kenn ich doch irgendwoher 

Code:
patch    faces    layers  overall thickness
                          [m]      [%]
-----    -----    ------  ---      ---
ground  352000  10      0.0609    100  
oben    49532    9.69    0.0477    99.7  
schraege 14924    4.81    0.0258    53.9  
unten    49533    9.69    0.0477    99.7 

Sieht gar nicht sooo schlecht aus.
Tipp:

Code:

blockMesh
surfaceTransformPoints -scale "(10000 10000 10000)"  bla.stl bla.stl
transformPoints -scale "(10000 10000 10000)"
snappyHexMesh -overwrite

surfaceTransformPoints -scale "(0.00001 0.00001 0.00001)" bla.stl bla.stl
transformPoints -scale "(0.00001 0.00001 0.00001)"


Viel Glück,

wenn ich mal Zeit hab schau ichs mir an.


------------------
Best regards,
Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

shake2014
Mitglied


Sehen Sie sich das Profil von shake2014 an!   Senden Sie eine Private Message an shake2014  Schreiben Sie einen Gästebucheintrag für shake2014

Beiträge: 3
Registriert: 07.09.2014

OpenFoam 2.3

erstellt am: 11. Sep. 2014 07:55    Editieren oder löschen Sie diesen Beitrag!  <-- editieren / zitieren -->   Antwort mit Zitat in Fett Antwort mit kursivem Zitat    Unities abgeben: 1 Unity (wenig hilfreich, aber dennoch)2 Unities3 Unities4 Unities5 Unities6 Unities7 Unities8 Unities9 Unities10 Unities

Guten Morgen Tobi,

es freut mich das du dich erinnerst  . Ich habe es als 2D Modell simuliert. Die Ergebnisse sind/waren auch gut. Nun, wo ich etwas mehr Zeit hab, will ich es als 3D Modell simulieren. Was ja der eigentliche Zweck war.

Das Ergebnis wäre gut, wenn es nicht 8Mio Zellen wären. Die zu berechnen ist eine große Herausforderung für mein PC 

Deinen Tipp habe ich umgesetzt. Auf Anhieb leider erfolglos. Es sind überhaupt keine Layer eingepflegt worden. Aber der Gedanke dahinter gefällt mir. Ich versuche es mit einer etwas kleineren Dimensionierung. Hoffentlich klappt es.

Gruß Shake 

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

Shor-ty
Moderator





Sehen Sie sich das Profil von Shor-ty an!   Senden Sie eine Private Message an Shor-ty  Schreiben Sie einen Gästebucheintrag für Shor-ty

Beiträge: 2463
Registriert: 27.08.2010

OpenFOAM-dev (Foundation)
OpenFOAM-xxxx (ESI)

erstellt am: 11. Sep. 2014 08:17    Editieren oder löschen Sie diesen Beitrag!  <-- editieren / zitieren -->   Antwort mit Zitat in Fett Antwort mit kursivem Zitat    Unities abgeben: 1 Unity (wenig hilfreich, aber dennoch)2 Unities3 Unities4 Unities5 Unities6 Unities7 Unities8 Unities9 Unities10 Unities Nur für shake2014 10 Unities + Antwort hilfreich

Hallo,

du musst in sHM deinen locationPoint natürlich mit 1000 multiplizieren, sowie die Layer definition wie:

Code:

finalyLayerThickness
minLayerThickness

------------------
Best regards,
Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

shake2014
Mitglied


Sehen Sie sich das Profil von shake2014 an!   Senden Sie eine Private Message an shake2014  Schreiben Sie einen Gästebucheintrag für shake2014

Beiträge: 3
Registriert: 07.09.2014

OpenFoam 2.3

erstellt am: 11. Sep. 2014 15:18    Editieren oder löschen Sie diesen Beitrag!  <-- editieren / zitieren -->   Antwort mit Zitat in Fett Antwort mit kursivem Zitat    Unities abgeben: 1 Unity (wenig hilfreich, aber dennoch)2 Unities3 Unities4 Unities5 Unities6 Unities7 Unities8 Unities9 Unities10 Unities

Hi,

ja die Parameter habe ich berücksichtigt gehabt. Es liegt mehr an diesen Fehler:

Code:
Checking mesh with layer ...
Checking faces in error :
    non-orthogonality > 65  degrees                                        : 27
    faces with face pyramid volume < 1e-13                          : 40
    faces with face-decomposition tet quality < 1e-30            : 47808
    faces with concavity > 80  degrees                                    : 0
    faces with skewness > 4  (internal) or 20  (boundary)      : 0
    faces with interpolation weights (0..1)  < 0.02                  : 135677
    faces with volume ratio of neighbour cells < 0.01            : 0
    faces with face twist < 0.02                                                : 56
    faces on cells with determinant < 0.001                            : 4048947

Detected 4232555 illegal faces (concave, zero area or negative cell pyramid volume)
Extruding 0 out of 135991 faces (0%). Removed extrusion at 84519 faces.
Added 1359910 out of 1359910 cells (100%).


Ich versuche es jetzt mal mit einem kleineren Raum und einer groberen Darstellung.

Vielen Dank Tobi.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

silencebreak
Mitglied



Sehen Sie sich das Profil von silencebreak an!   Senden Sie eine Private Message an silencebreak  Schreiben Sie einen Gästebucheintrag für silencebreak

Beiträge: 12
Registriert: 22.05.2017

OpenFoam version 2.4.0
CentOS 6.9

erstellt am: 21. Jun. 2017 15:09    Editieren oder löschen Sie diesen Beitrag!  <-- editieren / zitieren -->   Antwort mit Zitat in Fett Antwort mit kursivem Zitat    Unities abgeben: 1 Unity (wenig hilfreich, aber dennoch)2 Unities3 Unities4 Unities5 Unities6 Unities7 Unities8 Unities9 Unities10 Unities Nur für shake2014 10 Unities + Antwort hilfreich

Hey,

Ich stehe vor einem ähnlichen Problem. Konntest du es damals lösen?

Viele Grüße
Oli

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

Anzeige.:

Anzeige: (Infos zum Werbeplatz >>)

Darstellung des Themas zum Ausdrucken. Bitte dann die Druckfunktion des Browsers verwenden. | Suche nach Beiträgen

nächster neuerer Beitrag | nächster älterer Beitrag
Antwort erstellen


Diesen Beitrag mit Lesezeichen versehen ... | Nach anderen Beiträgen suchen | CAD.de-Newsletter

Administrative Optionen: Beitrag schliessen | Archivieren/Bewegen | Beitrag melden!

Fragen und Anregungen: Kritik-Forum | Neues aus der Community: Community-Forum

(c)2023 CAD.de | Impressum | Datenschutz