Hot News:

Unser Angebot:

  Foren auf CAD.de (alle Foren)
  OpenFOAM
  Problem Wandschichten-SnappyHM

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:  Problem Wandschichten-SnappyHM (5862 mal gelesen)
Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 19. Dez. 2011 19:14    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 möchte einen Spalt von 1mm mit snappyHexMesh vernetzen, und möchte 2-3 Wandschichten aufbauen. Nun dachte ich letzte Woche ich hatte Erfolg, doch ich sah nur das refinement Gebiet von snappy. Ich kann nach 20 Versuchen immernoch keine layers sehen, und dies ist 4Wochen vor Abgabe der Arbeit ein mieses Gefühl.

Mein log.snappy sieht wie folgt aus:

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.6                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                       |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.6-f802ff2d6c5a
Exec   : snappyHexMesh -overwrite
Date   : Dec 19 2011
Time   : 18:07:57
Host   : c026-12
PID    : 15297
Case   : /home/ft/ft8864/Documents/SkiSim_Case2/SKISim_Case2a/SKISim_Konvergenz
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

--> FOAM Warning :
    From function dlLibraryTable:  pen(const fileName& functionLibName)
    in file db/dlLibraryTable/dlLibraryTable.C at line 79
    could not load libparMetisDecomp.so: cannot open shared object file: No such file or directory
Create mesh for time = 0

Read mesh in = 0.66 s

Overall mesh bounding box  : (0 -0.005 -0.0025) (0.04 0.005 0.0025)
Relative tolerance         : 1e-06
Absolute matching distance : 4.15331e-08

Reading refinement surfaces.
Read refinement surfaces in = 0 s

Reading refinement shells.
Refinement level 2 for all cells outside refinementBox
Read refinement shells in = 0 s

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


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

Edge intersection testing:
    Number of edges             : 312500
    Number of edges to retest   : 312500
    Number of intersected edges : 40231
Calculated surface intersections in = 0.55 s

Initial mesh : cells:100000  faces:312500  points:112761
Cells per refinement level:
    0 100000

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

Patch Region
----- ------
motorBike:

5 motorBike_patch0

Added patches in = 0 s

Selecting decompositionMethod hierarchical

Refinement phase
----------------

Found point (0.01 0 0) in cell 44850 on processor 0
Reading external feature lines.
Read feature lines in = 0.06 s


Surface refinement iteration 0
------------------------------

Marked for refinement due to surface intersection : 80077 cells.
Marked for refinement due to curvature/regions    : 0 cells.
Determined cells to refine in = 0.02 s
Selected for refinement : 80077 cells (out of 100000)
Edge intersection testing:
    Number of edges             : 2090171
    Number of edges to retest   : 2090171
    Number of intersected edges : 160997
Refined mesh in = 6.28 s
After refinement surface refinement iteration 0 : cells:660539  faces:2090171  points:770035
Cells per refinement level:
    0 19923
    1 640616

Surface refinement iteration 1
------------------------------

Marked for refinement due to surface intersection : 320646 cells.
Marked for refinement due to curvature/regions    : 0 cells.
Determined cells to refine in = 0.11 s
Selected for refinement : 320649 cells (out of 660539)
Edge intersection testing:
    Number of edges             : 9317618
    Number of edges to retest   : 9277461
    Number of intersected edges : 643783
Refined mesh in = 67.48 s
After refinement surface refinement iteration 1 : cells:2905082  faces:9317618  points:3510328
Cells per refinement level:
    0 19900
    1 320174
    2 2565008

Surface refinement iteration 2
------------------------------

No cells marked for refinement since reached limit 2000000.
Determined cells to refine in = 0.11 s
Selected for refinement : 0 cells (out of 2905082)
Stopping refining since too few cells selected.


Removing mesh beyond surface intersections
------------------------------------------

Found point (0.01 0 0) in cell 344214 in global region 2 out of 5 regions.
Keeping all cells in region 2 containing point (0.01 0 0)
Selected for keeping : 1683074 cells.
Edge intersection testing:
    Number of edges             : 5664444
    Number of edges to retest   : 964314
    Number of intersected edges : 643783

Shell refinement iteration 0
----------------------------

Marked for refinement due to refinement shells    : 0 cells.
Determined cells to refine in = 11.14 s
Selected for internal refinement : 152 cells (out of 1683074)
Edge intersection testing:
    Number of edges             : 5667648
    Number of edges to retest   : 7308
    Number of intersected edges : 644084
Refined mesh in = 8.2 s
After refinement shell refinement iteration 0 : cells:1684138  faces:5667648  points:2302709
Cells per refinement level:
    0 0
    1 80000
    2 1604138

Shell refinement iteration 1
----------------------------

Marked for refinement due to refinement shells    : 0 cells.
Determined cells to refine in = 1.71 s
Selected for internal refinement : 0 cells (out of 1684138)
Stopping refining since too few cells selected.


Splitting mesh at surface intersections
---------------------------------------

Introducing baffles for 644084 faces that are intersected by the surface.

Edge intersection testing:
    Number of edges             : 6307648
    Number of edges to retest   : 5141672
    Number of intersected edges : 1284130
Created baffles in = 13.14 s


After introducing baffles : cells:1684138  faces:6307648  points:2302709
Cells per refinement level:
    0 0
    1 80000
    2 1604138

Introducing baffles to block off problem cells
----------------------------------------------

markFacesOnProblemCells : marked 1920631 additional internal faces to be converted into baffles.
Analyzed problem cells in = 5.33 s


Introducing baffles to delete problem cells.

Edge intersection testing:
    Number of edges             : 8228279
    Number of edges to retest   : 5781742
    Number of intersected edges : 1284135
Created baffles in = 27.79 s


After introducing baffles : cells:1684138  faces:8228279  points:2302709
Cells per refinement level:
    0 0
    1 80000
    2 1604138

Remove unreachable sections of mesh
-----------------------------------

Keeping all cells in region 10088 containing point (0.01 0 0)
Selected for keeping : 720000 cells.
Edge intersection testing:
    Number of edges             : 2445000
    Number of edges to retest   : 0
    Number of intersected edges : 320000
Split mesh in = 4.13 s


After subsetting : cells:720000  faces:2445000  points:1006507
Cells per refinement level:
    0 0
    1 80000
    2 640000

Handling cells with snap problems
---------------------------------

Introducing baffles for 320000 faces that are intersected by the surface.

Edge intersection testing:
    Number of edges             : 2445000
    Number of edges to retest   : 1282000
    Number of intersected edges : 320000
Created baffles in = 2.63 s


After introducing baffles : cells:720000  faces:2445000  points:1006507
Cells per refinement level:
    0 0
    1 80000
    2 640000

Introducing baffles to block off problem cells
----------------------------------------------

markFacesOnProblemCells : marked 0 additional internal faces to be converted into baffles.
Analyzed problem cells in = 1.76 s


Introducing baffles to delete problem cells.

Edge intersection testing:
    Number of edges             : 2445000
    Number of edges to retest   : 0
    Number of intersected edges : 320000
Created baffles in = 1.71 s


After introducing baffles : cells:720000  faces:2445000  points:1006507
Cells per refinement level:
    0 0
    1 80000
    2 640000

Remove unreachable sections of mesh
-----------------------------------

Keeping all cells in region 0 containing point (0.01 0 0)
Selected for keeping : 720000 cells.
Edge intersection testing:
    Number of edges             : 2445000
    Number of edges to retest   : 0
    Number of intersected edges : 320000
Split mesh in = 1.88 s


After subsetting : cells:720000  faces:2445000  points:1006507
Cells per refinement level:
    0 0
    1 80000
    2 640000
dupNonManifoldPoints : Found : 0 non-manifold points (out of 1006507)
Edge intersection testing:
    Number of edges             : 2445000
    Number of edges to retest   : 0
    Number of intersected edges : 320000
getDuplicateFaces : found 0 pairs of duplicate faces.

Detected unsplittable baffles : 0

Merge refined boundary faces
----------------------------

mergePatchFaces : Merging 0 sets of faces.
Removing 0 straight edge points.
Refined mesh : cells:720000  faces:2445000  points:1006507
Cells per refinement level:
    0 0
    1 80000
    2 640000
Writing mesh to time constant
Written mesh in = 19.77 s.

Morphing phase
--------------

Constructing mesh displacer ...
Using mesh parameters
{
    maxNonOrtho     65;
    maxBoundarySkewness 20;
    maxInternalSkewness 4;
    maxConcave      80;
    minFlatness     0.5;
    minVol          1e-13;
    minArea         -1;
    minTwist        0.02;
    minDeterminant  0.001;
    minFaceWeight   0.02;
    minVolRatio     0.01;
    minTriangleTwist -1;
    nSmoothScale    4;
    errorReduction  0.75;
}


Checking initial mesh ...
Checking faces in error :
    non-orthogonality >  65 degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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 2244000 illegal faces (concave, zero area or negative cell pyramid volume)
Checked initial mesh in = 4.67 s


Smoothing patch points ...
Smoothing iteration 0
Found 0 non-mainfold point(s).

Scaling iteration 0
Moving mesh using diplacement scaling : min:1  max:1
Checking faces in error :
    non-orthogonality >  65 degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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
Successfully moved mesh

Smoothing iteration 1
Found 0 non-mainfold point(s).

Scaling iteration 0
Moving mesh using diplacement scaling : min:1  max:1
Checking faces in error :
    non-orthogonality >  65 degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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
Successfully moved mesh

Smoothing iteration 2
Found 0 non-mainfold point(s).

Scaling iteration 0
Moving mesh using diplacement scaling : min:1  max:1
Checking faces in error :
    non-orthogonality >  65 degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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
Successfully moved mesh

Patch points smoothed in = 13.19 s


Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:1.60516e-05 min:5.95374e-06 max:1.6127e-05
Calculated surface displacement in = 0.11 s


Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Displacement smoothed in = 1.73 s


Moving mesh ...

Iteration 0
Moving mesh using diplacement scaling : min:1  max:1
Checking faces in error :
    non-orthogonality >  65 degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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
Successfully moved mesh
Moved mesh in = 4.05 s


Repatching faces according to nearest surface ...
Repatched 0 faces in = 1.66 s


Edge intersection testing:
    Number of edges             : 2445000
    Number of edges to retest   : 1308916
    Number of intersected edges : 320000
Snapped mesh : cells:720000  faces:2445000  points:1006507
Cells per refinement level:
    0 0
    1 80000
    2 640000
Writing mesh to time constant
Written mesh in = 17.19 s.

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

Using mesh parameters
{
    maxNonOrtho     65;
    maxBoundarySkewness 20;
    maxInternalSkewness 4;
    maxConcave      80;
    minFlatness     0.5;
    minVol          1e-13;
    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 < 30 degrees
      (cos:0.866025)
    - as long as the resulting face doesn't become concave by more than 90 degrees
      (0=straight, 180=fully concave)

Merging 0 sets of faces.

No faces merged ...

Merging all points on surface that
- are used by only two boundary faces and
- make an angle with a cosine of more than 0.866025.

No straight edges simplified and no points removed ...
Constructing mesh displacer ...

Checking mesh manifoldness ...

Checking initial mesh ...
Checking faces in error :
    non-orthogonality >  65 degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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 2244000 illegal faces (concave, zero area or negative cell pyramid volume)

Handling points with inconsistent layer specification ...

Handling non-manifold points ...

Checking patch manifoldness ...
Outside of local patch is multiply connected across edges or points at 0 points.
Set displacement to zero for all 0 non-manifold points

Handling feature edges ...
Set displacement to zero for points on 0 feature edges

Handling cells with warped patch faces ...
Set displacement to zero on 0 warped faces since layer would be > 0.5 of the size of the bounding box.

Growing non-extrusion points by one layer ...
Set displacement to zero for an additional 0 points.

patch               faces    layers avg thickness[m]
                                    near-wall overall
-----               -----    ------ --------- -------
motorBike_patch0    320000   2      1.36e-12  2.86e-12

medialAxisSmoothingInfo : Calculate distance to Medial Axis ...
shrinkMeshDistance : Smoothing normals ...
    Iteration 0   residual 1.64e-18
Inserting points on patch ski
Inserting points on patch snow
Inserting points on patch sides
Inserting points on patch inlet
Inserting points on patch outlet
shrinkMeshDistance : Smoothing normals ...
shrinkMeshDistance : Smoothing normals in interior ...
    Iteration 0   residual 0.0273

Layer addition iteration 0
--------------------------
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

Determining displacement for added points according to pointNormal ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

shrinkMeshMedialDistance : Smoothing using Medial Axis ...
shrinkMeshMedialDistance : Reduce layer thickness at 4000 nodes where thickness to medial axis distance is large
shrinkMeshDistance : Removing isolated regions ...
Number isolated points extrusion stopped : 0
shrinkMeshDistance : Smoothing field ...
    Iteration 0   residual 2.44e-15
shrinkMeshMedialDistance : Moving mesh ...
Iteration 0
Moving mesh using diplacement scaling : min:1  max:1
Checking faces in error :
    non-orthogonality > 65  degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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
shrinkMeshMedialDistance : Successfully moved mesh
shrinkMeshMedialDistance : Finished moving mesh ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.
truncateDisplacement : Unextruded 0 faces due to non-consecutive vertices being extruded.

Setting up information for layer truncation ...

Checking mesh with layer ...
Checking faces in error :
    non-orthogonality > 65  degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 4152016
    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 4152016 illegal faces (concave, zero area or negative cell pyramid volume)
Extruding 0 out of 320000 faces. Removed extrusion at 180000 faces.


Layer addition iteration 1
--------------------------
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

Determining displacement for added points according to pointNormal ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

shrinkMeshMedialDistance : Smoothing using Medial Axis ...
shrinkMeshMedialDistance : Reduce layer thickness at 0 nodes where thickness to medial axis distance is large
shrinkMeshDistance : Removing isolated regions ...
Number isolated points extrusion stopped : 0
shrinkMeshDistance : Smoothing field ...
    Iteration 0   residual 0
shrinkMeshMedialDistance : Moving mesh ...
Iteration 0
Moving mesh using diplacement scaling : min:1  max:1
Checking faces in error :
    non-orthogonality > 65  degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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
shrinkMeshMedialDistance : Successfully moved mesh
shrinkMeshMedialDistance : Finished moving mesh ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.
truncateDisplacement : Unextruded 0 faces due to non-consecutive vertices being extruded.

Setting up information for layer truncation ...

Checking mesh with layer ...
Checking faces in error :
    non-orthogonality > 65  degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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 2244000 illegal faces (concave, zero area or negative cell pyramid volume)
Extruding 0 out of 320000 faces. Removed extrusion at 0 faces.
Edge intersection testing:
    Number of edges             : 2445000
    Number of edges to retest   : 0
    Number of intersected edges : 320000
Writing 0 added cells to cellSet addedCells
Writing 0 faces inside added layer to faceSet layerFaces
Layer mesh : cells:720000  faces:2445000  points:1006507
Cells per refinement level:
    0 0
    1 80000
    2 640000
Writing mesh to time constant
Written mesh in = 62.4 s.
Finished meshing in = 279 s.
End



Ich vermute es liegt an diesem kleinen Abstand von 1mm.
Code:
Checking faces in error :
    non-orthogonality > 65  degrees                        : 0
    faces with face pyramid volume < 1e-13                 : 2244000
    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

Was sagt mir dieser Auszug genau? Kann es sein, dass die Vernetzung dort im Spalt "kein Platz findet"?

Bilder werd ich noch anfügen.
@Shorty: Ich weiß du empfiehlst Salome. Doch bis ich mich da einarbeite vergehen nochmals 2 Wochen? Es muss doch auch mit snappyHM möglich sein.

Vielleicht hat ja jemand den entscheidenden Tip.

Hier noch meine Einstellungen:

Code:
layers
    {
     
        motorBike_patch0
        {
            nSurfaceLayers 2;
        }

    }
   
    // Expansion factor for layer mesh
    expansionRatio 1.1;//1.3

    //- Wanted thickness of final added cell layer. If multiple layers
    //  is the thickness of the layer furthest away from the wall.
    //  See relativeSizes parameter.
    finalLayerThickness 0.3;//0.3

    //- Minimum thickness of cell layer. If for any reason layer
    //  cannot be above minThickness do not add layer.
    //  Relative to undistorted size of cell outside layer.
    minThickness 0.1; //0.1

    //- If points get not extruded do nGrow layers of connected faces that are
    //  also not grown. This helps convergence of the layer addition process
    //  close to features.
    nGrow 1;


    // Advanced settings

    //- When not to extrude surface. 0 is flat surface, 90 is when two faces
    //  make straight angle.
    featureAngle 30;

    //- Maximum number of snapping relaxation iterations. Should stop
    //  before upon reaching a correct mesh.
    nRelaxIter 3;


Gruß

[Diese Nachricht wurde von Burnquist am 19. Dez. 2011 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: 20. Dez. 2011 12:30    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 Burnquist 10 Unities + Antwort hilfreich

Hallo,

hast ja noch vier Wochen  Mein Vorschlag: Wäre es nicht einfach möglich nur BlockMesh zu verwenden? Warum Snappy?

Sieht deine Geometrie sehr kompliziert aus?

Gruß Thomas

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

Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 20. Dez. 2011 12:40    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


bl1.png

 
Die layer Bildung ist doch nur mit blockMesh nicht möglich, oder irre ich mich da?
Ich hab die Vorgabe, die Netzfeinheitstudie mit den Wandschichten durchzuführen.

Im Bild sind die zwei grauen Blöcke meine stl (Schnittdarstellung). Die Vernetzung im Inneren weißt nur das refinment Gebiet auf.

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: 20. Dez. 2011 13:57    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 Burnquist 10 Unities + Antwort hilfreich

Klar kannst du mit blockMesh eine Verfeinerung der Wandregionen vornehmen. Das geht sogar sehr einfach.

Ich teste das mal eben mit Snappy, wenn mein Laptop will. Ansonsten heut Abend.

[Diese Nachricht wurde von Shor-ty am 20. Dez. 2011 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: 20. Dez. 2011 16:13    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 Burnquist 10 Unities + Antwort hilfreich

Hallo,

definitiv eine Aufgabe für blockMesh. Für deine Geometrie ist ein Block ja vollkommen ausreichend.

Ausnahme: nicht simpleGrading sondern eher mit:

hex (XX XX XX XX XX XX XX XX) (XX XX XX) edgeGrading (1 1 1 1 2 2 2 2 3 3 3 3)

http://www.openfoam.org/docs/user/blockMesh.php

Gruß Thomas

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

Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 20. Dez. 2011 16:57    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

Ok, danke, wieder was gelernt 

Allerdings muss ich dieses Flächen später mit Rauigkeiten vesehen (unregelm. und regelm.) In diesem Fall fällt bM weg.

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: 21. Dez. 2011 12:40    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 Burnquist 10 Unities + Antwort hilfreich

Ich hab mich mal nochmals zurückerinnert und ich hab dir im anderen Betrag schon geschrieben, dass dies mit sHM funktioniert.
http://ww3.cad.de/foren/ubb/uploads/Shor-ty/sHM_layer.png

Grüße Tobi

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

Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 21. Dez. 2011 13:02    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

Bist du dir sicher, dass dies layers sind? Für mich sieht das wie das nächste refinementLvl aus. Bei dir sind die Zellen vielleicht so schmal, weil du kein aspectratio von 1:1:1 hast.

Nochmal kurz zum Verständnis:
Ich definiere all meine 6 Flächen im blockMesh mit patches. Die boundary-file gibt mir den Namen der stl aus (hier: motorBike_patch0), den ich im SHMdict unter addlayers einfüge.
nCellsBetweenLevels = 3
locationInMesh liegt bei mir mittig im Spalt (also zwischen den stls)
mode: inside --> wobei ich auch schon auf outside gestellt habe; ohne Erfolg

expansionRatio 1.1;//1.3

  finalLayerThickness 0.3;//0.3

  minThickness 0.1; //0.1

  nGrow 1;

  featureAngle 30;

Hast du etwas anderes definiert?

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: 21. Dez. 2011 13: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 Burnquist 10 Unities + Antwort hilfreich

Falsch.

Das sind meine Layer. Das Refinement ist bei mir auf 0 gesetzt damit ich nur die Layer erhalte.
Würd dir mal ne andere Version empfehlen. Wie ich gerade gesehen habe verwendest du die 1.6.

Vllt liegts da dran. Ich kann meine Layer so generieren wie ich will. In dem Beispiel is vllt das
expansionRatio etwas schlecht gewählt aber joa ich kann ja heut Abendn ochmals machen (hats ja gleich).


BlockMesh:
ich würde deinen Flächen schon die entsprechenden Typen zuordnen. Also nicht nur Patches sondern auch Wall. Die Boundary gibt dir keine STL aus. Sie gibt dir die Flächennamen aus, die du vorgeben musst. Das hat nichts mit deinen STL`s zu tun. Lediglich erhältst du über sHM ggf. in der Boundarydatei einen neuen Namen in den neu erstellten Zeitordnern.

Wenn bei dir diese Fläche mit "motorBike_patch0" benannt wird, dann is es richtig, dass du diese Angabe in sHM in die Layergeneration einfügen bzw. definieren musst.

du kannst auch mal mit dem folgenden Ausdruck arbeiten

Code:

"*."

und schauen was rauskommt.
Mode Inside - Outside hat damit nichts zu tun.

Inside --> dort wo deine LocationPoint steht wird vernetz (die von der STL eingschlossen wird)
Outside --> das negativ.

Bsp. Kugel. Du hast ne Kugel als STL vorliegen. Setzt du den locationPoint in die Kugel hinein, willst aber eine Umströmung simulieren, verwendest du "mode outside".
Oder du nimmst "mode inside" dabei muss aber der location point dann außerhalb deiner Kugel sein.

Ich hab nicht deine expliziten Angaben genommen, sondern eigene.

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

Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 21. Dez. 2011 19:00    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

Danke für die Antwort.

Habe heute einen Spalt mit 1cm höhe vernetzt. Klappt wunderbar, ich sehe 6 gleichmäßige Schichten.

Der 1mm Spalt geht nach wie vor nicht. Darf ich dir mal meine stl schicken und du schaust kurz ob das mit Salome geht? Dann verwende ich Salome

[Diese Nachricht wurde von Burnquist am 21. Dez. 2011 editiert.]

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: 21. Dez. 2011 19:16    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 Burnquist 10 Unities + Antwort hilfreich

Hey ho...

joa gerne aber ich teste es mit sHM. Auf meinem Profil ist meine Adresse.

Salome kann keine STL einlesen; nur ausgeben.
In Salome machst du die Geometrie erneut mit dem Geomodul und kannst die direkt vernetzen.

Aber schick erstmal deine STL dnan kuck ich mal was sich machen lässt. Salome ist nämlich gewöhnungsbedürftig << zumindest zu Beginn.


Tobi

[Diese Nachricht wurde von Shor-ty am 21. Dez. 2011 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: 22. Dez. 2011 09:36    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 Burnquist 10 Unities + Antwort hilfreich

Hallo,

ich verstehe jetzt aber immer noch nicht, warum blockMesh ausgeschlossen ist? Möchstest du später die Rauhigkeit auf den STLs bauen, sprich kleine Berge erzeugen und deswegen nur Snappy verwenden?

Gruß Thomas

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

Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 22. Dez. 2011 10:23    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

Richtig.
Meine untere Platte weißt später Rauigkeiten quer zur Strömung auf, meine Obere feine Rinnen in Längsrichtung.

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

Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 22. Dez. 2011 15:06    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


b.png

 
Hab die Vernetzung eines 2mm hohen Spaltes geschafft.
Wie im Bild zu sehen hat Snappy anscheinend ab diesem Schritt Probleme. Die oberen Schichten werden verzerrt. 1mm klappt dann nichtmehr.

Werde weiterhin an den Parametern experimentieren.

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: 22. Dez. 2011 16:29    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 Burnquist 10 Unities + Antwort hilfreich

Hallo,

und wenn du erst einmal in Meter arbeitest und nach dem Vernetzen runterskalierst?

Befehl: transformPoints

Gruß Thomas

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

Burnquist
Mitglied



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

Beiträge: 27
Registriert: 12.12.2011

erstellt am: 04. Jan. 2012 23:00    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


NeuesBild.PNG

 
Poste mal den aktuellen Stand (rauigkeitsbehafteter Spalt mit 1mm Dicke)

Die Layer Bildung funktioniert nur, wenn folgendes definiert ist:

Im blockMeshDict:
blocks
(
    hex (0 1 2 3 4 5 6 7) (30 30 20)    simpleGrading (1  1  1)  // In z-Richtung müssen 20 Zellen definiert sein. Bei 45,45,30 oder andere Kombinationen sehe ich keine layer.
);

Im snappyDict:
rau1mmzu
        {
            // Surface-wise min and max refinement level
            level (2 2);                                                                    //Bei 2,3 bzw. 3,4  funktioniert es nicht.
        }
refinementBox

        {
            mode inside;
            levels ((1E15 1));                                                        //Bei 2 bzw. 3 funktioniert es nicht.
        }

Somit komme ich laut log.file auf maximal 77150 Zellen.
Ziel sind immernoch 1Mio Zellen für die Netzunabhängigkeitsstudie.

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: 05. Jan. 2012 17:43    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 Burnquist 10 Unities + Antwort hilfreich

Kann sein, dass durch die Einstellungen der Feinheit deine Layer falsch definiert sind und deswegen keine Layer entstehen.
Eine mögliche Lösung ist der vom Thomas genannte Weg. Erstelle deine Gemetrie in Meter, vernetze diese und reskaliere dann wieder.
Nen Versuch is es Wert!

Grüße
PS: Sorry ich komm grad nich dazu irgendwie was zu testen.

[Diese Nachricht wurde von Shor-ty am 05. Jan. 2012 editiert.]

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