| | | Wie Sie mit 3D-Druck glattere Oberflächen erhalten, eine Pressemitteilung
|
Autor
|
Thema: twoPhaseEulerFoam mit kOmegaSSTSato (892 / mal gelesen)
|
Q.E.D. Mitglied
Beiträge: 18 Registriert: 21.12.2015
|
erstellt am: 23. Nov. 2017 07:54 <-- editieren / zitieren --> Unities abgeben:
Hallo alle zusammen! Ich experimentiere zurzeit etwas mit twoPhaseEulerFoam herum. Speziell wollte ich den Einfluss verschiedener Turbulenzmodelle auf das Ergebnis untersuchen. Der Case ist relativ einfach: einer simplen Rohrströmung wird Wasser zugeführt und dabei zerstäubt. Während der Case mit dem mixturekEpsilon-Modell ganz brauchbare Ergebnisse zeigt, konnte ich den Fall mit dem kOmegaSSTSato-Modell bisher nicht nachbilden. Der Fehler, der beim Run entsteht, ist unten aufgeführt. Der Fehler tritt nicht auf, wenn ich bei den turbulenceProperties des Wassers die Turbulenz auf off setze. Turbulenz bei der Luft wird normal berechnet. Offenbar führt der Aufruf der Funktion correct() zu diesem Fehler. In diesem Zusammenhang gibt es bereits eine Bug-Meldung, die sich aber darauf bezieht, dass die Funktion correctNut() nicht aufgerufen wird. Der Fehler wurde bereits behoben. (Siehe BugMeldung.) Leider konnte ich auch die Theorie hinter kOmegaSSTSato nicht beschaffen, um das in OpenFOAM umgesetzte Modell nachzuvollziehen. Meine Frage ist: hat jemand möglicherweise schon Erfahrungen mit kOmegaSSTSato gemacht oder könnte mir eventuell einen Tip geben, woran es eventuell liegen könnte? Beste Grüße Arthur Ps. hier der Fehler: Code:
Starting time loopCourant Number mean: 1.61528e-05 max: 0.021349 Max Ur Courant Number = 0.021349 deltaT = 1.17647e-06 Time = 1.17647e-06 PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 1 Min(alpha.air) = 0 Max(alpha.air) = 1 Constructing momentum equations smoothSolver: Solving for e.air, Initial residual = 1, Final residual = 7.31445e-08, No Iterations 2 smoothSolver: Solving for e.water, Initial residual = 1, Final residual = 2.31345e-13, No Iterations 1 min T.air 300 min T.water 300 GAMG: Solving for p_rgh, Initial residual = 0.62137, Final residual = 7.22273e-05, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 1.75782e-05, Final residual = 3.54548e-08, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 0.485616, Final residual = 0.000412704, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.000143182, Final residual = 2.59894e-08, No Iterations 4 smoothSolver: Solving for omega.air, Initial residual = 9.69396e-05, Final residual = 2.67436e-08, No Iterations 1 smoothSolver: Solving for k.air, Initial residual = 1, Final residual = 2.41552e-08, No Iterations 2 #0 Foam::error: rintStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-dev-debug/src/OSspecific/POSIX/printStack.C:218 #1 Foam::sigFpe::sigHandler(int) at ~/OpenFOAM/OpenFOAM-dev-debug/src/OSspecific/POSIX/signals/sigFpe.C:106 #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ~/OpenFOAM/OpenFOAM-dev-debug/src/OpenFOAM/matrices/lduMatrix/smoothers/GaussSeidel/GaussSeidelSmoother.C:167 #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ~/OpenFOAM/OpenFOAM-dev-debug/src/OpenFOAM/matrices/lduMatrix/smoothers/GaussSeidel/GaussSeidelSmoother.C:208 #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ~/OpenFOAM/OpenFOAM-dev-debug/src/OpenFOAM/matrices/lduMatrix/solvers/smoothSolver/smoothSolver.C:169 #6 Foam::fvMatrix<double>::solveSegregated(Foam::Dictionary const&) at ~/OpenFOAM/OpenFOAM-dev-debug/src/finiteVolume/fvMatrices/fvScalarMatrix/fvScalarMatrix.C:169 (discriminator 1) #7 Foam::fvMatrix<double>::solve(Foam::Dictionary const&) at ~/OpenFOAM/OpenFOAM-dev-debug/src/finiteVolume/lnInclude/fvMatrixSolve.C:82 #8 Foam::fvMatrix<double>::solve() at ~/OpenFOAM/OpenFOAM-dev-debug/src/finiteVolume/lnInclude/fvMatrixSolve.C:326 #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ~/OpenFOAM/OpenFOAM-dev-debug/src/finiteVolume/lnInclude/fvMatrix.C:1373 #10 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::PhaseCompressibleTurbulenceModel<Foam: haseModel> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::PhaseCompressibleTurbulenceModel<Foam: haseModel> > > >::correct() at ~/OpenFOAM/OpenFOAM-dev-debug/src/TurbulenceModels/turbulenceModels/lnInclude/kOmegaSSTBase.C:489 (discriminator 1) #11 Foam::RASModels::kOmegaSSTSato<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::PhaseCompressibleTurbulenceModel<Foam: haseModel> > > >::correct() at ~/OpenFOAM/OpenFOAM-dev-debug/src/TurbulenceModels/phaseCompressible/lnInclude/kOmegaSSTSato.C:164 #12 Foam::twoPhaseSystem::correctTurbulence() at ~/OpenFOAM/OpenFOAM-dev-debug/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C:548 #13 ? at ~/OpenFOAM/OpenFOAM-dev-debug/applications/solvers/multiphase/twoPhaseEulerFoam/contErrs.H:13 #14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #15 ? at ??:? Floating point exception (core dumped)
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
mbuergle Mitglied
Beiträge: 1 Registriert: 07.06.2018
|
erstellt am: 07. Jun. 2018 17:12 <-- editieren / zitieren --> Unities abgeben: Nur für Q.E.D.
|
Q.E.D. Mitglied
Beiträge: 18 Registriert: 21.12.2015
|
erstellt am: 10. Jun. 2018 13:02 <-- editieren / zitieren --> Unities abgeben:
Hallo Matthias, leider nein. Zugegeben, das Problem hatte damals auch nicht die allerhöchste Priorität, sodass ich keine Zeit hatte, das weiter zu verfolgen. Das kOmegaSSTSato ist in der Community auch nicht besonders gefragt, was du sicherlich schon bemerkt hast. Ich hatte es mit OpenFOAM 5.x getestet. Es mag sein, dass es mit einer älteren Version funktioniert. Mein bester Rat an dich wäre zu schauen, ob kOmegaSSTSato in einer älteren OF Version funktioniert. Wenn ja, müsste wahrscheinlich eine Kleinigkeit im Quelltext verändert werden. Dazu einfach schauen, was sich in den gängigeren Turbulenzmodellen verändert hat. Ansonsten kann ich dir zum Modell selbst leider wenig sagen, da ich nach wie vor die Theorie nicht vorliegen habe. Viel Erfolg weiterhin Arthur Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
| Anzeige.:
Anzeige: (Infos zum Werbeplatz >>)
|