Hot News:

Unser Angebot:

  Foren auf CAD.de (alle Foren)
  OpenFOAM
  snappyhexmesh-Problem

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:  snappyhexmesh-Problem (3645 mal gelesen)
bmandt
Mitglied



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

Beiträge: 14
Registriert: 26.10.2009

erstellt am: 30. Mrz. 2010 11:07    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 zusammen,

ich versuche einen Schiffstank mit snappyhexmesh zu vernetzen.

Die boundery box habe ich angebeben, einen Punkt innerhalb der Geometrie auch. Die Vernetzung läuft ohne Fehlermeldung durch. Das Problem ist, dass die Zellen ausserhalb der Geometrie nicht wegfallen. Nach dem snappyhexmesh durchegelaufen ist, habe meine boundery box gleichmäßig vernetzt.

Kennt jemand das Problem?

Als Input dient ein STL, was aus einem 3D-DXF erstellt wurde. Mit Rhino habe ich einmal Flächen aus dem Linien erzeugt, in einem 2. Versuch habe ich noch ein Polygonnetz drüber gelegt. In beiden Fällen tritt das gleiche Problem auf.

Ich hoffe, jemand kennt eine Antwort.

Gruß

Benjamin

===========================================

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

FoamFile
{
    version        2.0;
    format          ascii;

    root            "/home/penfold/mattijs/foam/mattijs2.1/run/icoFoam";
    case            "cavity";
    instance        "system";
    local          "";

    class          dictionary;
    object          autoHexMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

castellatedMesh true;
snap            true;
addLayers      true;

geometry
{
    cross01.stl
    {
        type triSurfaceMesh;
        name cross01;
    }

    refinementBox
    {
        type searchableBox;
        min (126.058 -19.830 0);
        max (133.617 9.830 11.500);
    }
};

castellatedMeshControls
{
    maxLocalCells 10000000;
    maxGlobalCells 50000000;
    minRefinementCells 10;
    nCellsBetweenLevels 2;

    features
    (
    );

    refinementSurfaces
    {
        cross01
        {
            level (2 2);
    faceZone cross01;
            cellZone cross01;
            zoneInside true;
        }
    }

    resolveFeatureAngle 60;

    refinementRegions
    {
        refinementBox
        {
            mode inside;
            levels ((1.0e9 2));
        }
    }

    locationInMesh (127.000.000 -4.970 1.000);
}

snapControls
{
    nSmoothPatch 3;
    tolerance 4.0;
    nSolveIter 30;
    nRelaxIter 5;
}

addLayersControls
{
    relativeSizes true;

    layers
    {
//flowAirDuctReagent_flowAirDuctReagent
        //{
        //    nSurfaceLayers 1;
        //}
//flowAirDuctReagent_inlet1
        //{
        //    nSurfaceLayers 1;
        //}
//flowAirDuctReagent_inlet2
        //{
        //    nSurfaceLayers 1;
        //}
    }

    expansionRatio 1.0;

    finalLayerThickness 0.3;

    minThickness 0.1;

    nGrow 0;

    featureAngle 30;

    nRelaxIter 3;

    nSmoothSurfaceNormals 1;

    nSmoothNormals 3;

    nSmoothThickness 10;

    maxFaceThicknessRatio 0.5;

    maxThicknessToMedialRatio 0.3;

    minMedianAxisAngle 130;

    nBufferCellsNoExtrude 0;
}


meshQualityControls
{
    maxNonOrtho 65;

    maxBoundarySkewness 20;
    maxInternalSkewness 4;

    maxConcave 80;

    minFlatness 0.5;

    minVol 0;

    minArea -1;

    minTwist 0.02;

    minDeterminant 0.001;

    minFaceWeight 0.02;

    minVolRatio 0.01;

    minTriangleTwist -1;

    nSmoothScale 4;
    errorReduction 0.75;
}


// Advanced

// Flags for optional output
// 0 : only write final meshes
// 1 : write intermediate meshes
// 2 : write volScalarField with cellLevel for postprocessing
// 4 : write current intersections as .obj files
debug 0;


// Merge tolerance. Is fraction of overall bounding box of initial mesh.
// Note: the write tolerance needs to be higher than this.
mergeTolerance 1E-6;


// ************************************************************************* //

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

Ocelot
Mitglied
Dipl-Ing. (BA) Maschinenbau


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

Beiträge: 171
Registriert: 29.04.2006

OpenSUSE Leap 42.2 x64
OpenFOAM-plus
OpenFOAM-4.x

erstellt am: 30. Mrz. 2010 20:25    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 bmandt 10 Unities + Antwort hilfreich

Hallo,

prüfe mal die Koordinaten des locationInMesh-Punktes, die x-Koordinate "127.000.000" scheint schon vom Format her nicht ganz korrekt zu sein...

Allgemein lässt sich zu deiner Problemstellung folgendes sagen:

Wenn in der vernetzten Endgeometrie zwar Teile des Umgebungsnetzes, nicht aber die gewünschte Simulationsgeometrie vorhanden sind, stimmt etwas mit dem locationInMesh-Punkt nicht. Manchmal genügt schon als Abhilfe, ungerade Werte für die Punktkoordinaten zu verwenden.

Wenn das fertige snappy-Mesh aus Umgebungsnetz und Simulationsgeometrie besteht, gibt es irgendwo eine Verbindung beider Netze die eigentlich nicht exisitieren dürfte, d.h. Löcher bzw. Fehlstellen in der STL-Geometrie.

Viele Grüße
Johannes

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

bmandt
Mitglied



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

Beiträge: 14
Registriert: 26.10.2009

erstellt am: 31. Mrz. 2010 07:28    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,

das mit dem Zahlenformat ist schonmal ein guter Tipp, habe ich übersehen, danke.

Gibt es eine Möglichkeit ein stl auf Löcher zu checken und ggf. zu schließen?
Dabei muessten es sehr kleine Öffnungen sein und snappyhexmesh sollte doch über sowas hinweggehen, oder sehe ich das falsch?

Viele Grüße

Benjamin

PS: Manchmal wird auch nicht bis zum Ende vernetzt. Fehlermeldung:

Shell refinement iteration 2
----------------------------

Marked for refinement due to refinement shells    : 247104 cells.
Determined cells to refine in = 0.141 s
Selected for internal refinement : 247104 cells (out of 247104)
terminate called after throwing an instance of 'std::bad_alloc'
  what():  std::bad_alloc

This application has requested the Runtime to terminate it in an unusual way.
Please contact the application's support team for more information.

[Diese Nachricht wurde von bmandt am 31. Mrz. 2010 editiert.]

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

TTB
Mitglied
CFD Engineer


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

Beiträge: 353
Registriert: 02.10.2008

BIM HVACTool für Windows OpenFOAM-2.2.x

erstellt am: 31. Mrz. 2010 10:56    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 bmandt 10 Unities + Antwort hilfreich

Hallo Benjamin,

ich gehe auch davon aus, dass deine STL nicht optimal für Snappy ist und deine Geometrie nicht geschlossen ist.

Schaue dir mal deine STL Datei mit der freien Software "admesh" an. Diese Software könnte deine STL auch reparieren.

Weiterhin hat auch OpenFOAM ein paar Tools, wie "surfaceCheck", um deine STL zu prüfen.

Gruß Thomas

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

bmandt
Mitglied



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

Beiträge: 14
Registriert: 26.10.2009

erstellt am: 31. Mrz. 2010 12:19    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 habe die Geometrie mit netfabb geschlossen. Dennoch bricht snappyhexmesh ab.

Surfacecheck meldet:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.6.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
/* OpenFOAM for Windows: www.beuth-hochschule,  Fachbereich 4, Klimalabor      *\

\*---------------------------------------------------------------------------*/
Build  : 1.6.x-b4e37019668d
Exec  : surfacecheck constant\triSurface\cross01.stl
Date  : Mar 31 2010
Time  : 12:16:58
Host  : MANDT
PID    : 2184
Case  : c:/daten/tankKlein
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Reading surface from "constant\triSurface\cross01.stl" ...

Statistics:
Triangles    : 15222
Vertices    : 7603
Bounding Box : (128.515 3.2 0) (131.16 9.83 1.85)

Region  Size
------  ----
kleinerTankPolyHealed(repariert)        15222


Surface has no illegal triangles.

Triangle quality (equilateral=1, collapsed=0):
    0 .. 0.05  : 0.0173433
    0.05 .. 0.1  : 0.00762055
    0.1 .. 0.15  : 0.0129418
    0.15 .. 0.2  : 0.0135987
    0.2 .. 0.25  : 0.0103797
    0.25 .. 0.3  : 0.0124819
    0.3 .. 0.35  : 0.0172776
    0.35 .. 0.4  : 0.01419
    0.4 .. 0.45  : 0.0162265
    0.45 .. 0.5  : 0.029234
    0.5 .. 0.55  : 0.0184601
    0.55 .. 0.6  : 0.0727237
    0.6 .. 0.65  : 0.0995927
    0.65 .. 0.7  : 0.0899356
    0.7 .. 0.75  : 0.0316647
    0.75 .. 0.8  : 0.351991
    0.8 .. 0.85  : 0.0445408
    0.85 .. 0.9  : 0.0243726
    0.9 .. 0.95  : 0.0279858
    0.95 .. 1  : 0.0874392

    min 1e-015 for triangle 137
    max 0.999474 for triangle 2303

Dumping bad quality faces to "badFaces"
Paste this into the input for surfaceSubset


Edges:
    min 0.000991821 for edge 4760 points (131.102 9.83 1.825)(131.103 9.83 1.825
)
    max 0.467405 for edge 21789 points (129.176 8.48 1.85)(129.507 8.81 1.85)

Checking for points less than 1E-6 of bounding box ((2.645 6.63 1.85) meter) apa
rt.
Found 0 nearby points.

Grüße

Benjamin

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

TTB
Mitglied
CFD Engineer


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

Beiträge: 353
Registriert: 02.10.2008

BIM HVACTool für Windows OpenFOAM-2.2.x

erstellt am: 31. Mrz. 2010 17:34    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 bmandt 10 Unities + Antwort hilfreich

Hallo,

ich sehe gerade, dass du mit der Windows-Version arbeitest. Es kann sein, dass da Snappy nicht ganz 100% funktioniert. Hast du die Möglichkeit das in Linux auszuprobieren? Kannst du dein Case veröffenltichen?

Gruß Thomas

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

Ocelot
Mitglied
Dipl-Ing. (BA) Maschinenbau


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

Beiträge: 171
Registriert: 29.04.2006

OpenSUSE Leap 42.2 x64
OpenFOAM-plus
OpenFOAM-4.x

erstellt am: 31. Mrz. 2010 19:54    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 bmandt 10 Unities + Antwort hilfreich

Hallo Benjamin,

ich stimme Thomas zu, zu allererst solltest du prüfen ob snappy überhaupt korrekt funktioniert. Am besten du verwendest dafür eine ganz primitive STL-Geometrie (Quader, Kugel, etc..) von der du weißt, dass sie i.O. ist. Oder einfach das motorBike-Tutorial durchlaufen lassen.

Grüße
Johannes

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

bmandt
Mitglied



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

Beiträge: 14
Registriert: 26.10.2009

erstellt am: 08. Apr. 2010 09:58    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 mal wieder,

ich habe jetzt mal das MotorBike-Tutorial mit der Windows-Version durchlaufen lassen und es hat keine Probleme gegeben.
Ich denke eher, dass das Problem an der Zellenanzahl oder sonstigen Einstellungen liegt.

Ich habe meinen Tank mal verkleinert, aber auch hier werden die Zellen ausserhalb des Tanks nicht entfernt. Laut NETFABB und admesh ist die Geometrie geschlossen.

Viele Grüße

Benjamin

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

t.schumacher
Mitglied
CFD Engineer


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

Beiträge: 184
Registriert: 03.05.2010

erstellt am: 21. Mai. 2010 08:56    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 bmandt 10 Unities + Antwort hilfreich

Hallo,
hast du das Problem mittlerweile gelöst?

ThomasS

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