Hot News:

Unser Angebot:

  Foren auf CAD.de (alle Foren)
  OpenFOAM
  Magnus-Effekt

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:  Magnus-Effekt (1632 / mal gelesen)
loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 29. Feb. 2016 10: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


Screenshotfrom2016-02-29095619.png

 
Hallo! Für den Ein oder Anderen mag das Folgende ein leicht zu lösendes Problem sein. Ich halte mich daran allerdings schon länger auf.
Ziel ist es, den Magnus-Effekt zu simulieren. Also eine freie Strömung an einem rotierenden Zylinder.
Ob ich das ganze mit simple, pimple oder IcoFoam mache weiß ich noch nicht.
Hauptproblem ist wie immer die Gitternetzerstellung.
Ich habe mir diverse Tutorials angesehen, komme aber nie zu befriedigenden Ergebnissen, denn
die Simulationen brechen jedes mal mit der Meldung Floating Point Exception ab.
Ich habe es sowohl mit SHM als auch mit GMsh probiert. GMsh ist für mich erstmal gestorben, da es im minutentakt abstürzt und eigentlich noch kein vernünftiges Gitter geliefert hat. Mit SHM komme ich soweit klar, jedoch will es einfach nicht gelingen die Region um den Zylinder zu generieren. Bilder der bisherigen Versuche hängen an. Ich würde mich sehr freuen, wenn es dafür eine Lösung gibt.

Um auf den Hintergrund zurück zukommen. Ich schließe momentan mein Praxissemester ab und möchte meinem Ingenieur an diesem Bsp. die Funktionsweise von openFoam nahebringen bzw. dann später mit ihm auch Kurzanleitungen für Simulationen zu Strömungsvorgängen erarbeiten, da dieses Thema(openFoam) an unserer FH absolutes Neuland ist. Alles was ich bisher weiß, habe ich mir selbst beibringen müssen inkl.der Grundkenntnisse für C++. Also seht es mir nach, wenn ich bei einigen Dingen gerne mal "dumme" Fragen stelle. 

Zum Case:
Der Grundgedanke war, einen Raum (x,y,z) mit z (eine Zelle tiefe) also quasi 2 dimensional zu erstellen.
code für p

Code:
dimensions      [0 2 -2 0 0 0 0];

internalField  uniform 100;

boundaryField
{
    Inlet
    {
        type            zeroGradient;
    }
    Outlet
    {
        type            fixedValue;
        value          $internalField;
    }
    fixedWalls
    {
        type            empty;
    }

    frontAndBack
    {
        type            empty;
    }
    cylinder
    {
        type            fixedValue;
        value          $internalField;
    }
}


code für u

Code:
dimensions      [0 1 -1 0 0 0 0];

internalField  uniform (0 0 0);

boundaryField
{
    Inlet
    {
        type            fixedValue;
        value          uniform (1 0 0);
    }
    Outlet
    {
    type zeroGradient;
    }
    fixedWalls
    {
        type            empty;
    }

    frontAndBack
    {
        type            empty;
    }
   
    cylinder
    {
        type rotatingWallVelocity;
        origin (0 0 0);
        axis (0 0 1);
        omega 10;
    }
}


Winkelgeschwindigkeit soll später variieren, so das die Umfangsgeschw. zwischen 10 - 20 m/s liegt.

hier nun der Code der SHM.dict. Hinweis! Diese Datei habe ich aus einem Tutorial und sie ist in keinem Fall so wie ich es mir vorstelle.
Das angehängte Bild gibt das Ergebnis genau dieser Datei wieder.

Code:
#includeEtc "caseDicts/mesh/generation/snappyHexMeshDict.cfg"

castellatedMesh on;
snap            on;
addLayers      off;

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

castellatedMeshControls
{
    features
    (
      { file "cylinder.eMesh"; level 6; }
     
    );

    refinementSurfaces
    {
       
        cylinder
        {
            level      (6 6);
            cellZone cylinder;
           
        }

    }

    refinementRegions
    {
       
        cylinder
        {
            mode inside;
            levels ((1e-5 5));
        }
    }

    locationInMesh (450 450 0.001); // Offset from (0 0 0) to avoid
                                    // coinciding with face or edge
}

snapControls
{
    explicitFeatureSnap    true;
}

addLayersControls
{
    layers
    {
    }

    relativeSizes      true; // false, usually with firstLayerThickness
    expansionRatio      1.2;
    finalLayerThickness 0.05;
    minThickness        1e-3;
//  firstLayerThickness 0.01;

//  maxThicknessToMedialRatio 0.6;
}

meshQualityControls
{
//    minTetQuality -1e+30;
}

writeFlags
(
    scalarLevels
    layerSets
    layerFields
);

mergeTolerance 1e-7;



------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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: 29. Feb. 2016 10:39    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 loopster 10 Unities + Antwort hilfreich

Hey,..

ich hab deine RB nicht angeschaut aber dein Netz ist wirklich furchtbar (:
Ist das nur ein Zylinder? Wenn ja dann schau dir doch das hier an: http://holzmann-cfd.de/index.php/de/tutorial/kugelumstroemung
Zudem lieftert dir sHM immer ein 3d Netz. Heißt also das du dein Netz erstmal 3d aufbaust und ein patch (front / back) dann auf eine
Zelle extrudierst. Dazu wäre dann das hier für dich gut: http://holzmann-cfd.de/index.php/de/tutorial/hydraulischer-sprung

Sollte es damit noch nicht gehen, einfach nochmals meldung machen. Achja, wahrscheinlich wäre es am geschicktesten wenn du gleich mit dem Tutorial für den hydraulischen Sprung anfängst.

PS: Danke für die Nennung des Effekts ... (:

------------------
Viele Grüße,
Tobias Holzmann

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: 29. Feb. 2016 14: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 loopster 10 Unities + Antwort hilfreich

Hey,

nur zur Info:

Code:

    cylinder
    {
        type rotatingWallVelocity;
        origin (0 0 0);
        axis (0 0 1);
        omega 10;
    }

Es wäre schöner wenn du das mit 6DoF machst (: Also nen Rigid-Körper nimmst.

------------------
Viele Grüße,
Tobias Holzmann

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: 29. Feb. 2016 15: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 loopster 10 Unities + Antwort hilfreich


MagnusEffect.png

 
Tada 

Das gefällt mir, da kann man seine Karman-Straße schön iwo hin lenken.

------------------
Viele Grüße,
Tobias Holzmann

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 29. Feb. 2016 16: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

Ok, ich bin verwirrt. Erstmal danke für die Infos.
Frage Nr. 1  Bei dem Kugel Case habe benutzt Du sowohl .stl files in triSurface als auch die ideasToFoam mit vorheriger Nutzung von Salome. Das Zusammenspiel der einzelnen Gitter-Sektionen habe ich noch nicht ganz verstanden. Ich habe mir alle .stl Files in MeshLab angesehen und blickte da nicht ganz durch.

Und die 6DoF Sache ist mir auch noch ein Rätsel.(lese mich grad ein)

Zu guter Letzt! Genau sowas will ich haben!  Nur wie??? 

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 29. Feb. 2016 17:39    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

Ich habe versucht den hydraulischen Sprung laufen zu lassen. Dort bekomme ich jedoch die Fehlermeldung, das GroovyBC nicht bekannt ist. Nach einer Änderung wie z.Bsp. zu fixedValue folgt das nicht erkennen der Turbulenzproperties "RASModels"  usw. Ich will darin nicht weiter rumwerken. Lange Rede....das mit den Baffles?! Tatsächlich lauter Dinge die mir neu sind. 
Je mehr ich in die Dateistrukturen von openFoam vordringe desto mehr und undurchsichtiger wird es. Ich muss gestehen, dass das Frustrationspotenzial sehr hoch ist, da es an wirklichen Erfolgserlebnissen mangelt. 
Zu den Salome files. Wurde dort mit STEP files gearbeitet und dann ähnlich den bekannten YouTube-tutorials gearbeitet?https://www.youtube.com/watch?v=1zQbU-E4k1U

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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: 29. Feb. 2016 20: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 Nur für loopster 10 Unities + Antwort hilfreich

Hi,

um dich mal etwas zu entwirren.

Salome nutze ich zur CAD Erstellung und oft zur Erstellung des *.unv Netzes. Das ist nicht anderes als ein Hexnetz das man mit "ideasUnvToFoam" konvertieren kann. Ist also identisch zu blockMesh. Das erstellt dir nur ein Hintergrund-Netz. Mit Baffles und allem anderen brauchst du dich nicht rumärgern. Es ging eigentlich nur um die Vernetzung mittels sHM und dann dem Extrudieren auf 2D.

a) Hintergrundnetz,
b) STL (Cylinder)
c) Vernetzen
d) Setup
e) Fertig.


GroovyBC ist übrigens von Bernhard Gscheider entwickelt und ist im swak4Foam enthalten. Ein MUSS für Foamer 

------------------
Viele Grüße,
Tobias Holzmann

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 01. Mrz. 2016 11: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


Screenshotfrom2016-03-01105943.png

 
Hallo und nochmals Danke für die Hilfe.
Im Bild sieht man was ich bisher erreichen konnte.
Zum Werdegang des Ganzen.
Ich habe mit SolidWorks einen Zylinder erstellt, diesen dann mit Salome in Gruppen eingeteilt, um die
Zylinderwand implementieren zu können.
Mit BlockMesh habe ich dann ein Hintergrundgitter erstellt.
inlet
outlet
UpDownWalls (oben und unten)
fixedWalls (vorne/hinten) da ich, sofern richtig verstanden die extrudeDict an die UpDownWalls anlehnen möchte.
Dies hat bisher noch nicht funktioniert.
Ohne extrude komme ich zu diesem Ergebnis (Bild)
Leider habe ich es bisher nicht geschafft die Kanten an den Zylindergrenzen vernünftig hin zubekommen.
Code:
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      snappyHexMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Which of the steps to run
castellatedMesh true;
snap            true;
addLayers      true;


// Geometry. Definition of all surfaces. All surfaces are of class
// searchableSurface.
// Surfaces are used
// - to specify refinement for any mesh cell intersecting it
// - to specify refinement for any mesh cell inside/outside/near
// - to 'snap' the mesh boundary to the surface
geometry
{
    cylinderWall.stl
    {
        type triSurfaceMesh;
        name cylinderWall;
    }
    cylinderUpDown.stl
    {
        type triSurfaceMesh;
        name cylinderUpDown;
    }
    refinementBox
    {
        type searchableBox;
        min (-0.15 -0.15 -0.025);
        max ( 0.15 0.15 0.025);
    }
};

// Settings for the castellatedMesh generation.
castellatedMeshControls
{
    maxLocalCells 20000000;
    maxGlobalCells 10000000;
    minRefinementCells 2;
    maxLoadUnbalance 0.00;
    nCellsBetweenLevels 10;

    // Explicit feature edge refinement
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    // Specifies a level for any cell intersected by its edges.
    // This is a featureEdgeMesh, read from constant/triSurface for now.
    features
    (
        {file "cylinderWall.eMesh"; level 3;}
        {file "cylinderUpDown.eMesh"; level 0;}
       
    );

    // Surface based refinement
    // ~~~~~~~~~~~~~~~~~~~~~~~~

    // Specifies two levels for every surface. The first is the minimum level,
    // every cell intersecting a surface gets refined up to the minimum level.
    // The second level is the maximum level. Cells that 'see' multiple
    // intersections where the intersections make an
    // angle > resolveFeatureAngle get refined up to the maximum level.

    refinementSurfaces
    {
        cylinderWall
        {
            level (3 3);
         
        }
       
       
       
    }

    // Resolve sharp angles
    resolveFeatureAngle 90;


    // Region-wise refinement
    // ~~~~~~~~~~~~~~~~~~~~~~

    // Specifies refinement level for cells in relation to a surface. One of
    // three modes
    // - distance. 'levels' specifies per distance to the surface the
    //  wanted refinement level. The distances need to be specified in
    //  descending order.
    // - inside. 'levels' is only one entry and only the level is used. All
    //  cells inside the surface get refined up to the level. The surface
    //  needs to be closed for this to be possible.
    // - outside. Same but cells outside.

    refinementRegions
    {
    refinementBox
    {
    levels ((1 3));
    mode inside;
    }
    }


    // Mesh selection
    // ~~~~~~~~~~~~~~

    // After refinement patches get added for all refinementSurfaces and
    // all cells intersecting the surfaces get put into these patches. The
    // section reachable from the locationInMesh is kept.
    // NOTE: This point should never be on a face, always inside a cell, even
    // after refinement.
    locationInMesh (0.1 0.1 0.01);


    // Whether any faceZones (as specified in the refinementSurfaces)
    // are only on the boundary of corresponding cellZones or also allow
    // free-standing zone faces. Not used if there are no faceZones.
    allowFreeStandingZoneFaces true;
}

// Settings for the snapping.
snapControls
{
    //- Number of patch smoothing iterations before finding correspondence
    //  to surface
    nSmoothPatch 4;

    //- Relative distance for points to be attracted by surface feature point
    //  or edge. True distance is this factor times local
    //  maximum edge length.
    tolerance 1.2;

    //- Number of mesh displacement relaxation iterations.
    nSolveIter 150;

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

    // Feature snapping

        //- Number of feature edge snapping iterations.
        //  Leave out altogether to disable.
        nFeatureSnapIter 6;

        //- Detect (geometric only) features by sampling the surface
        //  (default=false).
        implicitFeatureSnap true;

        //- Use castellatedMeshControls::features (default = true)
        explicitFeatureSnap true;

        //- Detect points on multiple surfaces (only for explicitFeatureSnap)
        multiRegionFeatureSnap false;
}

// Settings for the layer addition.
addLayersControls
{
    // Are the thickness parameters below relative to the undistorted
    // size of the refined cell outside layer (true) or absolute sizes (false).
    relativeSizes false;

    // Per final patch (so not geometry!) the layer information
    layers
    {
        cylinderWall
        {
            nSurfaceLayers 6;
            finalLayerThickness 0.0014;
        }
    }

    // Expansion factor for layer mesh
    expansionRatio 1.3;

    // Wanted thickness of final added cell layer. If multiple layers
    // is the thickness of the layer furthest away from the wall.
    // Relative to undistorted size of cell outside layer.
    // See relativeSizes parameter.
    finalLayerThickness 0.0013;

    // 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.0001;

    // 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.
    // Note: changed(corrected) w.r.t 17x! (didn't do anything in 17x)
    nGrow 0;

    // Advanced settings

    // When not to extrude surface. 0 is flat surface, 90 is when two faces
    // are perpendicular
    featureAngle 0;

    // At non-patched sides allow mesh to slip if extrusion direction makes
    // angle larger than slipFeatureAngle.
    slipFeatureAngle 90;

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

    // Number of smoothing iterations of surface normals
    nSmoothSurfaceNormals 3;

    // Number of smoothing iterations of interior mesh movement direction
    nSmoothNormals 8;

    // Smooth layer thickness over surface patches
    nSmoothThickness 10;

    // Stop layer growth on highly warped cells
    maxFaceThicknessRatio 0.15;

    // Reduce layer growth where ratio thickness to medial
    // distance is large
    maxThicknessToMedialRatio 0.3;

    // Angle used to pick up medial axis points
    // Note: changed(corrected) w.r.t 17x! 90 degrees corresponds to 130 in 17x.
    minMedianAxisAngle 90;


    // Create buffer region for new layer terminations
    nBufferCellsNoExtrude 0;


    // Overall max number of layer addition iterations. The mesher will exit
    // if it reaches this number of iterations; possibly with an illegal
    // mesh.
    nLayerIter 50;

    nRelaxedIter 20;
}

// Generic mesh quality settings. At any undoable phase these determine
// where to undo.
meshQualityControls
{
    #include "meshQualityDict"

relaxed
    {
        //- Maximum non-orthogonality allowed. Set to 180 to disable.
        maxNonOrtho 180;
    }
    // Advanced

    //- Number of error distribution iterations
    nSmoothScale 4;
    //- Amount to scale back displacement at error points
    errorReduction 0.75;
}
writeFlags
(
    scalarLevels    // write volScalarField with cellLevel for postprocessing
    layerSets      // write cellSets, faceSets of faces in layer
    layerFields    // write volScalarField for layer coverage
);

// Advanced
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-7;


Code:
convertToMeters 0.1;
a 20;
b 20;
c 1;
xa -4;
xb  4;
ya -4;
yb  4;
za -0.25;
zb  0.25;


vertices
(
    ($xa $ya $za)
    ($xb $ya $za)
    ($xb $yb $za)
    ($xa $yb $za)
    ($xa $ya $zb)
    ($xb $ya $zb)
    ($xb $yb $zb)
    ($xa $yb $zb)
);

blocks
(
    hex (0 1 2 3 4 5 6 7) ($a $b $c) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
   
    inlet
    {
    type patch;
    faces
    (
    (0 4 7 3)
    );
    }
    outlet
    {
        type patch;
        faces
        (
            (2 6 5 1)
        );
    }
    fixedWalls
    {
        type empty;
        faces
        (
           
           
            (1 5 4 0)
            (3 7 6 2)
        );
    }
    UpDownWalls
    {
    type empty;
    faces
    (
    (0 3 2 1)
    (4 5 6 7)
    );
    }
);

mergePatchPairs
(
);



------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 01. Mrz. 2016 11:24    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


Screenshotfrom2016-03-01112115.png


Screenshotfrom2016-03-01112121.png

 
Zusatz. ExtrudeMesh liefert mir so etwas....ich glaube zu wissen wo hier das Problem liegt, habe aber keine Rechte Idee wie ich das lösen kann.

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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: 01. Mrz. 2016 11:26    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 loopster 10 Unities + Antwort hilfreich

Hey,


  • Ich rate dir dazu ab featureEdge refinement zu verwenden wenn du damit keine Erfahrung hast, siehe auch http://holzmann-cfd.de/index.php/de/tutorial/featureedge-refinement
  • Deine Files bin ich jetzt nicht durchgegangen (anhang wäre übrigens besser)
  • Dein Netz sieht wirklich schlecht aus, da kannst auch nix extrudieren (möglich vllt noch mit flattenMesh das gerade zu biegen)
  • Es sieht so aus als würde dein Hintergrundnetz höher als der Zylinder sein

Sofern du es heute Abend noch nicht geschaft hast, wird das Tutorial auf meiner Seite für diesen Fall zu finden sein.

------------------
Viele Grüße,
Tobias Holzmann

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 01. Mrz. 2016 11:47    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


Screenshotfrom2016-03-01114303.png


cylinder1.tar.gz

 
Vielen Dank erstmal. Ich hänge jetzt einfach mal den ganzen Case an.

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 01. Mrz. 2016 13:21    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


Screenshotfrom2016-03-01132020.png

 
soho......ich hoffe dass ich jetzt bereit bin für weitere Schritte.

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 01. Mrz. 2016 14:33    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


Screenshotfrom2016-03-01143243.png


Screenshotfrom2016-03-01143246.png

 
TADA!!!    

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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: 01. Mrz. 2016 14: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 Nur für loopster 10 Unities + Antwort hilfreich

Gratuliere. Wahrscheinlich denkst du dir  gerade das es wirklich einfach war, oder?

Es wäre dennoch interessant zu wissen, was dein Problem war.

------------------
Viele Grüße,
Tobias Holzmann

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 01. Mrz. 2016 14: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

Einfach...naja. 
Es waren mehrere Dinge. Zum einen hattest Du Recht. Das Hintergrundgitter war zu hoch. Was mir noch nicht ganz klar ist, da ich beim Erstellen der .stl immer auf die Skalierung achte. Ich habe also das Hintergrundgitter sehr genau anpassen müssen, sowie die refinemnet regions, searchableBox usw. Ganz flach ist das Gitter an der Zylinderkante immernoch nicht, aber simpleFoam scheint das nicht zu stören.
für extrudeMesh habe ich die Ober und Unterseite des Hintergrundgitters nochmal aufgeteilt und mich einfach für eine seite entschieden.
Ich habe übrigens die
Code:

    {
        type rotatingWallVelocity;
        origin (0 0 0);
        axis (0 0 1);
        omega 5000;
    }

Variante gewählt.

Das von Dir erwähnte 6DoF hat etwas mit Freiheitsgraden zu tun?
Da bin ich in der Kürze der Zeit nicht hintergestiegen.

Für GroovyBC reichen meine c++ Kenntnisse momentan leider nicht aus. Bin ja schon froh, dass ich "halbwegs" verstehe, was welcher code bei openFoam macht.

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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: 01. Mrz. 2016 16: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 Nur für loopster 10 Unities + Antwort hilfreich

Hey,

6DoF sind 6 Freiheitsgrade für einen Körper (Rigid) bspw. Du könntest somit den Zylinder als Festkörper simulieren und diesen dann in Rotation bringen. Wäre ne aufwendigere aber auch schönere Simulation. Ich habe noch ein Skalar mit genommen und dann sieht das alles auchschon anders aus. Zudem hab ich mittels Groovy noch die Rotationsgeschwindigkeit variiert 

Trotzdem, coole Sache das es selber hingebracht hast.

------------------
Viele Grüße,
Tobias Holzmann

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 02. Mrz. 2016 11: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 und dankeschön!
Als nächstes steht an, die Ergebnisse über der Zeit für U p und auch die Kräfte die auf den Zylinder wirken darzustellen.
Das ist ja das, worum beim Magnus-Effekt geht. Ziel soll sein, das Verhältnis zwischen Anström- und Rotationsgeschwindigkeit so einzustellen, das ein optimaler Auf-, bzw Abtrieb entsteht. Mein Ing und ich haben uns erstmal dafür entschieden weiter mit einem inkompressiblem Medium zu arbeiten, da die Anströmgeschwindigkeit im normal Fall recht gering ist.
Dazu bräuchte ich auch noch ein paar Infos, bezüglich der Zusatzfunktionen. Ich hatte etwas gelesen über Gnuplot oder ähnliches?

Ich wäre auch hier dankbar für Infos.

Das nutzen eines Festkörpers klingt sehr interessant.

P.S.: die Rotationsgeschwindigkeit habe ich über die Winkelgeschwindigkeit omega variiert.

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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: 02. Mrz. 2016 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 loopster 10 Unities + Antwort hilfreich

Hey,


  • Drag & Drop Kräfte -> schau ins motorBike totorial da sind die Funktionen glaub gegeben (controlDict)
  • Optimierung würde ich mittels Dakota durchführen (sofern Min- oder Maximierungsproblem)
  • Winkelgeschwindigkeit ist aber bei die ein fixer Wert oder?

Meine Interpretation des Magnus-Effekts ist nun auch verfügbar: Magnus Effekt
Die Umsetzung ist hier zu finden http://www.holzmann-cfd.de/index.php/de/openfoam

------------------
Viele Grüße,
Tobias Holzmann

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

loopster
Mitglied
Student


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

Beiträge: 35
Registriert: 29.09.2015

ACER 7750G,
Intel Core i5 2450M,
2.5-3.1 GHz,
AMDRadeon HD7670M 2GB,
Intel® Sandybridge Mobile ,
DDR3 8GB,
Windows 7 Home/Prem,
Linux/Ubuntu 14.04LTS,

erstellt am: 03. Mrz. 2016 11:01    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!
Ja die Winkelgeschwindigkeit ist Fix.
Habe den Case einfach mit drei verschiedenen Geschwindigkeiten laufen lassen um erstmal einen Vergleich zu haben.

Deine Version schaue ich mir gleich mal an.

Ich melde mich bestimmt wieder, mit weiteren Fragen.  

------------------
"Die ersten Menschen waren nicht die letzten Affen." (Kästner)

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