| |
| Materialien für den 3D-Druck im Hochtemperaturbereich, eine Pressemitteilung
|
Autor
|
Thema: pimpleFoam stoppt nach erstem Zeitschritt (1796 / mal gelesen)
|
agm Mitglied Automatisierungsingenieur
Beiträge: 12 Registriert: 21.01.2021
|
erstellt am: 13. Feb. 2021 15:35 <-- editieren / zitieren --> Unities abgeben:
Hallo ich versuche mich gerade daran einen Rührer in einem Behälter zu simulieren. Nachdem ich nun die Geometrie fertig habe und die ersten Probleme mit der Konfiguration bewältigen konnte, scheint pimpleFoam zu rechnen. Leider beendet sich die Simulation nach dem ersten Zeitschritt ohne eine Fehlermeldung (denke ich).
Code: /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _7bdb509494-20201222 OPENFOAM=2012 Arch : "LSB;label=32;scalar=64" Exec : pimpleFoam Date : Feb 13 2021 Time : 14:50:04 Host : Gauss PID : 462 I/O : uncollated Case : /home/amai/of/cases/amai/Mixer nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: solidBody Applying solid body motion to cellZone rotatingZone Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone rotatingZone PIMPLE: no residual control data found. Calculations will employ 3 corrector loops Reading field p Reading field U Reading/calculating face flux field phi AMI: Creating addressing and weights between 758424 source faces and 760892 target faces AMI: Patch source sum(weights) min:0.7520201 max:17.15187 average:1.001055 AMI: Patch target sum(weights) min:0 max:1.979961 average:1.00046 Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave RAS { RASModel kOmegaSST; turbulence on; printCoeffs on; alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.5555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; decayControl false; kInf 0; omegaInf 0; } No MRF models present No finite volume options present Constructing face velocity Uf Courant Number mean: 6.88888e-11 max: 0.003455231 Starting time loop Courant Number mean: 6.88888e-11 max: 0.003455231 deltaT = 0.025 Time = 0.025 PIMPLE: iteration 1 AMI: Creating addressing and weights between 758424 source faces and 760892 target faces AMI: Patch source sum(weights) min:0.6784376 max:17.15187 average:1.00108 AMI: Patch target sum(weights) min:0 max:1.981575 average:1.000505 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 7.653184e-07, No Iterations 27 time step continuity errors : sum local = 1.49118e-09, global = -7.6762e-13, cumulative = -7.6762e-13 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0919353, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.09198667, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.09107858, No Iterations 3 amai@Gauss:~/of/cases/amai/Mixer$
Entsprechend dem controlDict sollte die Simulation jedoch deutlich länger laufen (10s)
Code: /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //application pimpleFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 10; deltaT 0.025; writeControl adjustable; writeInterval 0.05; purgeWrite 0; writeFormat ascii; writePrecision 7; writeCompression no; timeFormat general; timePrecision 6; runTimeModifiable yes; adjustTimeStep yes; maxCo 1; maxDeltaT 1; // ************************************************************************* //
Habt Ihr eine Idee wo mein Fehler liegt, im Moment bin ich etwas ratlos. Gruß Andreas Georg
[Diese Nachricht wurde von agm am 13. Feb. 2021 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 15. Feb. 2021 15:19 <-- editieren / zitieren --> Unities abgeben: Nur für agm
hey um das ganze weiterrechnen zu lassen kannst Du mal maxCo 1; im controlDict rauswerfen ... und den timeStep kleiner wählen ... ggf adjustable auch raus... dann kannst Du das wirte intervall klein setzen damit er wenigstens einen zeitschrit schreibt ... dann kannst Du Dir das mal in paraFoam ansehen und so eine idee bekommen warum die kiste explodiert Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
agm Mitglied Automatisierungsingenieur
Beiträge: 12 Registriert: 21.01.2021
|
erstellt am: 16. Feb. 2021 17:35 <-- editieren / zitieren --> Unities abgeben:
Hallo @hulli1 danke ich habe beides (maxCo und TimeStep) probiert hat leider nicht geholfen. Jetzt habe ich nochmal checkMesh laufen lassen. Da kommen dann ein paar Fehler die nicht verstehe (Anfänger / Hobby)
Code: amai@Gauss:~/of/cases/amai/Mixer$ checkMesh -allGeometry -allTopology -writeAllFields /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : _7bdb509494-20201222 OPENFOAM=2012 Arch : "LSB;label=32;scalar=64" Exec : checkMesh -allGeometry -allTopology -writeAllFields Date : Feb 16 2021 Time : 16:42:18 Host : Gauss PID : 503 I/O : uncollated Case : /home/amai/of/cases/amai/Mixer nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Enabling all (cell, face, edge, point) topology checks. Enabling all geometry checks. Writing mesh quality as fields 15 ( faceZone cellDeterminant nonOrthoAngle cellVolumeRatio skewness minTetVolume faceWeight aspectRatio cellRegion cellAspectRatio minPyrVolume wallDistance cellVolume cellZone cellShapes ) Time = 0
Mesh stats points: 64240792 faces: 185933019 internal faces: 181875355 cells: 60901619 faces per cell: 6.039386 boundary patches: 8 point zones: 0 face zones: 1 cell zones: 1 Overall number of cells of each type: hexahedra: 59063867 prisms: 609322 wedges: 0 pyramids: 0 tet wedges: 3009 tetrahedra: 56 polyhedra: 1225365 Breakdown of polyhedra by number of faces: faces number of cells 4 90531 5 67642 6 240220 7 1131 8 980 9 609721 10 258 11 48 12 168908 13 86 14 32 15 44902 17 6 18 900 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Topological cell zip-up check OK. <<Number of duplicate (not baffle) faces found: 9. This might indicate a problem. <<Number of faces with non-consecutive shared points: 12. This might indicate a problem. <<Writing 33 faces with non-standard edge connectivity to set edgeFaces <<Writing 1 cells with zero or one non-boundary face to set oneInternalFaceCells <<Writing 286 cells with two non-boundary faces to set twoInternalFacesCells *Number of regions: 3 The mesh has multiple regions which are not connected by any face. <<Writing region information to "0/cellToRegion" <<Writing region 0 (fully disconnected) with 1 cells to cellSet region0 <<Writing region 1 (fully disconnected) with 53253910 cells to cellSet region1 <<Writing region 2 (fully disconnected) with 7647708 cells to cellSet region2 Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology Bounding box airInlet 296 325 ok (non-closed singly connected) (-0.01249262 -0.01249511 -0.235) (0.01249973 0.01249507 -0.2345703) waterInlet 290 320 ok (non-closed singly connected) (-0.4124962 0.3995686 -0.06248615) (-0.3875003 0.4 -0.03752626) waterOutlet 296 325 ok (non-closed singly connected) (0.3875076 -0.4 0.5375047) (0.4124998 -0.3995705 0.5624949) tankTop 32907 33350 ok (non-closed singly connected) (-0.4999788 -0.499985 0.7977873) (0.4999999 0.4999851 0.8) tankWalls 149496 152519 ok (non-closed singly connected) (-0.4999994 -0.4999994 -0.235) (0.5 0.4999994 0.8) Mixer 2355063 2458743 ok (closed singly connected) (-0.15 -0.1499999 -0.02000003) (0.15 0.1499999 0.075) AMI1 758424 760923 ok (non-closed singly connected) (-0.1699989 -0.1699993 -0.03000076) (0.17 0.1699996 0.0900088) AMI2 760892 762619 ok (non-closed singly connected) (-0.1699989 -0.1699999 -0.03000091) (0.17 0.1699996 0.09000063) Checking faceZone topology for multiply connected surfaces... FaceZone Faces Points Surface topology Bounding box rotatingZone 1519316 1522713 ok (closed singly connected) (-0.1699989 -0.1699999 -0.03000091) (0.17 0.1699996 0.0900088) Checking basic cellZone addressing... CellZone Cells Points VolumeBoundingBox rotatingZone 53253911 55606033 0.01057497 (-0.1699989 -0.1699993 -0.03000076) (0.17 0.1699996 0.0900088) Checking geometry... Overall domain bounding box (-0.4999994 -0.4999994 -0.235) (0.5 0.4999994 0.8) Mesh has 3 geometric (non-empty/wedge) directions (1 1 1) Mesh has 3 solution (non-empty) directions (1 1 1) Boundary openness (1.137245e-15 -2.181296e-16 2.356427e-14) OK. Max cell openness = 6.023385e-16 OK. Max aspect ratio = 23.28093 OK. Minimum face area = 7.158894e-09. Maximum face area = 0.0001037044. Face area magnitudes OK. Min volume = 6.058035e-13. Max volume = 1.040871e-06. Total volume = 0.7833632. Cell volumes OK. Mesh non-orthogonality Max: 65.00749 average: 4.565471 Non-orthogonality check OK. Face pyramids OK. Max skewness = 2.871118 OK. Coupled point location match (average 0) OK. ***Error in face tets: 806 faces with low quality or negative volume decomposition tets. <<Writing 804 faces with low quality or negative volume decomposition tets to set lowQualityTetFaces Min/max edge length = 1.321443e-05 0.01020106 OK. <<Writing 860938 near (closer than 1.75249e-06 apart) points to set nearPoints *There are 4613 faces with concave angles between consecutive edges. Max concave angle = 79.97198 degrees. <<Writing 4613 faces with concave angles to set concaveFaces Face flatness (1 = flat, 0 = butterfly) : min = 0.5475574 average = 0.9998917 *There are 499 faces with ratio between projected and actual area < 0.8 Minimum ratio (minimum flatness, maximum warpage) = 0.5475574 <<Writing 499 warped faces to set warpedFaces Cell determinant (wellposedness) : minimum: 0 average: 1.138269 ***Cells with small determinant (< 0.001) found, number of cells: 333 <<Writing 333 under-determined cells to set underdeterminedCells ***Concave cells (using face planes) found, number of cells: 362891 <<Writing 362891 concave cells to set concaveCells Face interpolation weight : minimum: 0.05663294 average: 0.4920535 Face interpolation weight check OK. Face volume ratio : minimum: 0.01339107 average: 0.9619096 Face volume ratio check OK. Calculating AMI weights between owner patch: AMI1 and neighbour patch: AMI2 AMI: Creating addressing and weights between 758424 source faces and 760892 target faces AMI: Patch source sum(weights) min:0.7520201 max:17.15187 average:1.001055 AMI: Patch target sum(weights) min:0 max:1.979961 average:1.00046 Failed 3 mesh checks. Writing fields with mesh quality parameters Writing non-orthogonality (angle) to nonOrthoAngle Writing face interpolation weights (0..0.5) to faceWeight Writing face skewness to skewness Writing cell determinant to cellDeterminant Writing aspect ratio to aspectRatio Writing approximate aspect ratio to cellAspectRatio Writing cell shape (hex, tet etc.) to cellShapes Writing cell volume to cellVolume Writing cell volume ratio to cellVolumeRatio Writing minTetVolume to minTetVolume Writing minPyrVolume to minPyrVolume Writing cell region to cellRegion Selecting patchDistMethod meshWave Writing wall distance to wallDistance Writing cell zoning to cellZone Writing face zoning to faceZone End amai@Gauss:~/of/cases/amai/Mixer$
Was habe ich da verbrochen? Mein Workflow war
Code:
rm -r 0blockMesh > blockMesh.log surfaceFeatureExtract > surfaceFeatureExtract.log decomposePar > decomposePar.log mpirun -np 16 snappyHexMesh -parallel > snappyHexMesh.log reconstructParMesh -latestTime > reconstructPar.log cp -r 0.05/polyMesh constant/polyMesh renumberMesh -latestTime > renumberMesh.log createPatch -overwrite pimpleFoam
Ist an dem Workflow vieleicht etwas falsch? Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2466 Registriert: 27.08.2010 ESI-OpenCFD OpenFOAM v2312
|
erstellt am: 01. Mrz. 2021 08:05 <-- editieren / zitieren --> Unities abgeben: Nur für agm
Ich frag mich inständig, wieso ich keine Emails mehr von CAD.de bekomme. Das ist schon nervig Interessante Ausgabe. Der Solver sollte eigentlich weiterlaufen. Anscheinend (nach dem Log-File) wird vor der Druck-Korrektur einfach aufgehört - die Frage ist warum. FOAM gibt normalerweise einen Fehler aus wenn was nicht korrekt ist. Auch das Netz ist per se mal in Ordnung. Ich hab hier leider gerade nur meinen Laptop, auf dem ich ein Netz mit 12 Kernen nicht erstellen kann und ich daher davon ausgehe, dass es auch etwas größer ist (schau es mir mal heute Abend an). Ansonsten ist es nicht möglich eine konkrete Aussage anhand Deiner Daten zu geben. Was ich jedoch empfehle ist, mit moveDynamicMesh zu prüfen, ob das AMI etc. korrekt funktioniert.
------------------ Glück Auf, Tobi OpenFOAM® Community - Knowledge Base Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
agm Mitglied Automatisierungsingenieur
Beiträge: 12 Registriert: 21.01.2021
|
erstellt am: 02. Mrz. 2021 20:44 <-- editieren / zitieren --> Unities abgeben:
|
Shor-ty Moderator
Beiträge: 2466 Registriert: 27.08.2010 ESI-OpenCFD OpenFOAM v2312
|
erstellt am: 02. Mrz. 2021 20:56 <-- editieren / zitieren --> Unities abgeben: Nur für agm
|
agm Mitglied Automatisierungsingenieur
Beiträge: 12 Registriert: 21.01.2021
|
erstellt am: 02. Mrz. 2021 21:40 <-- editieren / zitieren --> Unities abgeben:
Das dauert noch, die 0,2s haben 9h gebraucht Wenn die Strömung etwas weiter entwickelt ist, stelle ich nochmal Bilder ein, wo dann auch der Behälter in dem sich der Rührer befindet, zu erkennen ist. Der ist aber im Moment sooooo laaaangweilig. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2466 Registriert: 27.08.2010 ESI-OpenCFD OpenFOAM v2312
|
erstellt am: 02. Mrz. 2021 21:43 <-- editieren / zitieren --> Unities abgeben: Nur für agm
|
agm Mitglied Automatisierungsingenieur
Beiträge: 12 Registriert: 21.01.2021
|
erstellt am: 04. Mrz. 2021 21:37 <-- editieren / zitieren --> Unities abgeben:
Danke, gibt es ein guten Test Case (am besten mehrphasig) dazu? Zum Hintergrund, bei der Geometrie handelt es sich um einem Rührer mit kleinen Bohrungen in der Bodenplatte durch die Luft in die Wasserströmung einströmt. Ziel ist eine gleichmäßige Gasverteilung in der Biologie einer Kläranlage. Das habe ich vor 25 Jahren mal als 2D Problem einphasig als Diplomarbeit gerechnet - jetzt möchte ich das ganz langsam zu Ende bringen (mehr als 2D war damals nicht Aufgabenstellung ;-) ). Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2466 Registriert: 27.08.2010 ESI-OpenCFD OpenFOAM v2312
|
erstellt am: 06. Mrz. 2021 12:48 <-- editieren / zitieren --> Unities abgeben: Nur für agm
|