Autor
|
Thema: Felder plotten (1035 mal gelesen)
|
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 04. Feb. 2020 11:54 <-- editieren / zitieren --> Unities abgeben:
Hi Leute, ich habe mal eine best-practice Frage. Ich habe einen Fall bei dem eine Kugel umströmt wird (laminarer case - creeping flow). Mit snappyHeyMesh habe ich refinment regions um die kugel drumherum gesetzt und dann habe ich das ganze rechnen lassen. Wenn ich mir das Ganze in paraview ansehe (z.B. U cell data das orangene symbol) erhalte ich ein pixeliges Abbild, von dem was ich erwarten würde. Was mich aber daran nervt, sind stufenartige Übergange zwischen den refinement regions. Jetzt habe ich aber gelernt, dass paraview das Ganze in Abhängigkeit zur Schnittlage interpoliert und diese Stufen recht normal wären. Dummerweise werden diese noch deutlicher, wenn ich mir die interpolierten pointata Felder ansehe. Wie plottet Ihr denn Eure Felder (z. B. U) mit paraview um diese Stufen wegzubekommen? Für Eure Tipps bin ich wie immer dankbar VLG H Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 04. Feb. 2020 13:17 <-- editieren / zitieren --> Unities abgeben: Nur für hulli1
|
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 04. Feb. 2020 13:33 <-- editieren / zitieren --> Unities abgeben:
|
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 04. Feb. 2020 13:34 <-- editieren / zitieren --> Unities abgeben:
|
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 04. Feb. 2020 13:39 <-- editieren / zitieren --> Unities abgeben: Nur für hulli1
|
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 04. Feb. 2020 14:04 <-- editieren / zitieren --> Unities abgeben:
Hi Tobi, das ganze macht mich schon auch stutzig. Anbei die ux Felder, da sieh es nicht janz sooo schlimm aus ... Oder meinst Du es liegt an der Aspect-ratio sprich kleine Object riesen Domäne ??? LG H [Diese Nachricht wurde von hulli1 am 04. Feb. 2020 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 04. Feb. 2020 14:08 <-- editieren / zitieren --> Unities abgeben:
hier nich meine schemes und die solution /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.1.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //<P>ddtSchemes { default backward; }<P>gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear; }<P>divSchemes { default none; div(phi,U) Gauss filteredLinear2V 0.2 0; div(phi,k) Gauss limitedLinear 1; div(phi,B) Gauss limitedLinear 1; div(phi,nuTilda) Gauss limitedLinear 1; div(B) Gauss linear; div((nuEff*dev2(T(grad(U))))) Gauss linear; div(((nu*rho)*dev(grad(U).T()))) Gauss linear; //div((viscousTerm*dev(grad(U).T()))) Gauss linear; div((nuEff*dev2(T(grad(U))))) Gauss linear 1; div((viscousTerm*dev2(grad(U).T()))) Gauss linear 1; div(phi,T) Gauss limitedLinear 1; }<P>laplacianSchemes { default Gauss linear corrected; // default none; // laplacian(nuEff,U) Gauss linear corrected; // laplacian((1|A(U)),p) Gauss linear corrected; //laplacian(DkEff,k) Gauss linear corrected; //laplacian(DBEff,B) Gauss linear corrected; //laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; //laplacian((nu*rho),U) Gauss linear corrected; //laplacian((voidfraction2|A(U)),p) Gauss linear corrected; //laplacian(viscousTerm,U) Gauss linear corrected; }<P>interpolationSchemes { default linear; interpolate(U) linear; }<P>snGradSchemes { default corrected; }<P>fluxRequired { default no; p ; }<P> // ************************************************************************* //<P> fvSolutions<P>/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.1.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //<P>solvers {<P> p { solver GAMG; tolerance 1e-6; relTol 0.1;<P> smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2;<P> cacheAgglomeration true;<P> nCellsInCoarsestLevel 50;//10; agglomerator faceAreaPair; mergeLevels 1; };<P> pFinal { solver GAMG; tolerance 1e-6; relTol 0;<P> smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2;<P> cacheAgglomeration true;<P> nCellsInCoarsestLevel 50;//10; agglomerator faceAreaPair; mergeLevels 1; };<P> /* p { solver PCG; preconditioner DIC; tolerance 1e-06; relTol 0.05; }<P> pFinal { solver PCG; preconditioner DIC; tolerance 1e-06; relTol 0; } */<P> U { solver PBiCGStab; preconditioner DILU; tolerance 1e-05; relTol 0; } "(U|k|epsilon)Final" { $U; tolerance 1e-05; relTol 0; } k { solver PBiCGStab; preconditioner DILU; tolerance 1e-05; relTol 0; }<P> B { solver PBiCGStab; preconditioner DILU; tolerance 1e-05; relTol 0; }<P> nuTilda { solver PBiCGStab; preconditioner DILU; tolerance 1e-05; relTol 0; } T { solver PBiCGStab; preconditioner DILU; tolerance 1e-05; relTol 0; } }<P>PISO { nCorrectors 2; nNonOrthogonalCorrectors 4; pRefPoint (0.002 0.002 0.039); pRefValue 0; } PIMPLE { nOuterCorrectors 1; nCorrectors 2; nNonOrthogonalCorrectors 2; pRefPoint (0.002 0.002 0.039); pRefValue 0; }<P> relaxationFactors // <- more stable / -> more convergent { fields { p 0.3; pFinal 1; //Underrelaxation of the final iteration: < 1 only if alot of nOuterCorrectors (~ 100) are done } equations { U 0.9; UFinal 1; //Underrelaxation of the final iteration: < 1 only if alot of nOuterCorrectors (~ 100) are done "(nut|k|epsilon)" 0.9; } }<P> // ************************************************************************* //
[Diese Nachricht wurde von hulli1 am 04. Feb. 2020 editiert.] [Diese Nachricht wurde von hulli1 am 04. Feb. 2020 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 04. Feb. 2020 16:35 <-- editieren / zitieren --> Unities abgeben: Nur für hulli1
Hi, es gibt im Forum die code Tags. Bitte das nächste mal verwenden. Also Dein Profil kann prinzipiell (zeit-punktuell) passen. Du bekommst ja Ablösungen. Was mich aber gerade stutzig macht ist das hier:
Code:
div(phi,U) Gauss filteredLinear2V 0.2 0;
Ich hab dieses numerische Schema noch nie verwendet. Ich gehe mal davon aus, das hier Dein Problem her rührt. Übrigens, sag bitte nicht Du verwendest OpenFOAM 2.1.x. Das war ja fast die Version mit der ich angefangen hab (vor 9 Jahren). ------------------ Viele Grüße, Tobi OpenFOAM® Tutorials | Publikationen | Für Anfänger wiki.openfoam.com Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 04. Feb. 2020 16:36 <-- editieren / zitieren --> Unities abgeben: Nur für hulli1
|
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 05. Feb. 2020 10:41 <-- editieren / zitieren --> Unities abgeben:
Guten Morgen !!! Vielen Vielen Dank !!! .. Haha nein nein keine Sorge im Moment habe ich 5x laufen ich weiß ich sollte mal updaten gerade gerade lauft alles. Ich copy paste die scrpite nur immer mal ein bisschen herum Top der case ist der Hammer. Aber was viel interessanter ist das adaptive Refinement, was Du da eingebaut hast. Das ist ja der Hammer. Es ist mir ein bisschen peinlich aber Du mir sagen wie Du den case startest. Ich finde kein blockmesh ;((( normalerweise gehe ich so vor blockMesh surfaceFeatureExtract snappyHexmesh pisoFoam bzw in parallel Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 05. Feb. 2020 14:45 <-- editieren / zitieren --> Unities abgeben:
Hi Tobi, dein Case habe ich zum laufen geracht Danke dafür !!! Nur eins noch wie hast Du das mit dem adapitve Mesh Refinement rechnen lassen. Wenn ich pimpleMyDyFoam nehme bekomme ich die Meldung Unknown dynamicFvMesh type dynamicRefine2DFvMesh. Muss ich das mit den Versionen von openfoam.com rechnen lassen ??? sprich v1812 ??? LG H Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 05. Feb. 2020 17:04 <-- editieren / zitieren --> Unities abgeben: Nur für hulli1
|
hulli1 Mitglied
Beiträge: 61 Registriert: 23.01.2020 --
|
erstellt am: 06. Feb. 2020 12:06 <-- editieren / zitieren --> Unities abgeben:
|