Autor
|
Thema: Vorgabewerte via SimpleFoam (2037 mal gelesen)
|
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 29. Mai. 2011 20:37 <-- editieren / zitieren --> Unities abgeben:
Hallo zusammen, hab da gerade ein kleines Problem. Versuch einfach aus Interesse die Karmansche Wirbelstraße nachzubilden. Die Umströmung an einem Zylinder hat schon ganz gut funktioniert. --> http://www.youtube.com/watch?v=oxyvV61zbxk Das gleiche versuche ich derzeit mit zwei umströmten Zylindern (Bild)... Ich habe mir gedacht, zu Beginn solve ich mit dem SimpleFOAM um nach ein paar Iterationsschritten, diese weiter mit dem PimpleFoam weiterrechnen zu lassen. Jedoch bricht mir der pimpleFoam immer wieder ab - egal wie ich die Co -Zahl und meine dT einstell.
Code:
DILUPBiCG: Solving for k, Initial residual = 0.605555, Final residual = 8.00151e-06, No Iterations 40 DILUPBiCG: Solving for Ux, Initial residual = 0.752058, Final residual = 9.07024e-06, No Iterations 72 DILUPBiCG: Solving for Uy, Initial residual = 0.696289, Final residual = 9.13927e-06, No Iterations 78 GAMG: Solving for p, Initial residual = 0.772051, Final residual = 0.00619517, No Iterations 15 time step continuity errors : sum local = 5.87453e-05, global = 1.89729e-06, cumulative = 6.9077e-07 GAMG: Solving for p, Initial residual = 0.0174493, Final residual = 5.43005e-07, No Iterations 30 time step continuity errors : sum local = 3.79406e-07, global = -3.0262e-08, cumulative = 6.60509e-07 DILUPBiCG: Solving for epsilon, Initial residual = 0.868715, Final residual = 7.56692e-06, No Iterations 89 bounding epsilon, min: -4664.21 max: 1.81625e+06 average: 37844.9 DILUPBiCG: Solving for k, Initial residual = 0.810382, Final residual = 4.79514e-06, No Iterations 98 bounding k, min: -461.751 max: 101135 average: 2273.11 ExecutionTime = 31.68 s ClockTime = 32 sCourant Number mean: 0.110016 max: 25.3441 deltaT = 1.64736e-06 --> FOAM Warning : From function Time: perator++() in file db/Time/Time.C at line 688 Increased the timePrecision from 9 to 10 to distinguish between timeNames at time 2250 Time = 2250.000293 DILUPBiCG: Solving for Ux, Initial residual = 0.983884, Final residual = 0.0430883, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.92528, Final residual = 0.021659, No Iterations 3 GAMG: Solving for p, Initial residual = 0.597389, Final residual = 0.00283059, No Iterations 12 time step continuity errors : sum local = 9.57969e-07, global = -1.61293e-08, cumulative = 6.44379e-07 GAMG: Solving for p, Initial residual = 1.93431e-09, Final residual = 1.93431e-09, No Iterations 0 time step continuity errors : sum local = 0.000247586, global = -1.61293e-08, cumulative = 6.2825e-07 DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 7.793e-07, No Iterations 13 bounding epsilon, min: -193920 max: 8.74493e+06 average: 216491 DILUPBiCG: Solving for k, Initial residual = 0.732745, Final residual = 12978.9, No Iterations 1001 bounding k, min: -6.71451e+15 max: 1.20902e+17 average: 2.19224e+13 DILUPBiCG: Solving for Ux, Initial residual = 0.614391, Final residual = 6.89226e-06, No Iterations 8 DILUPBiCG: Solving for Uy, Initial residual = 0.713802, Final residual = 5.92312e-06, No Iterations 13 GAMG: Solving for p, Initial residual = 6.10787e-06, Final residual = 8.24413e-07, No Iterations 4 time step continuity errors : sum local = 0.0294446, global = -2.75492e-08, cumulative = 6.00701e-07 GAMG: Solving for p, Initial residual = 1.79165e-15, Final residual = 1.79165e-15, No Iterations 0 time step continuity errors : sum local = 1.79146e+18, global = 23.4613, cumulative = 23.4613 DILUPBiCG: Solving for epsilon, Initial residual = 0.0716592, Final residual = 9.37785, No Iterations 1001 bounding epsilon, min: -9.99488e+30 max: 1.12314e+31 average: 4.7151e+25 #0 Foam::error: rintStack(Foam::Ostream&) in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so" #2 in "/lib/libc.so.6" #3 Foam::LimitedScheme<double, Foam::limitedLinearLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libfiniteVolume.so" #4 Foam::limitedSurfaceInterpolationScheme<double>::weights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libfiniteVolume.so" #5 Foam::fv::gaussConvectionScheme<double>::fvmDiv(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) const in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libfiniteVolume.so" #6 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::Div<double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libincompressibleRASModels.so" #7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::Div<double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libincompressibleRASModels.so" #8 Foam::incompressible::RASModels::kEpsilon::correct() in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libincompressibleRASModels.so" #9 in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/applications/bin/linux64GccDPOpt/pimpleFoam" #10 __libc_start_main in "/lib/libc.so.6" #11 in "/home/tobi/OpenFOAM/OpenFOAM-1.7.x/applications/bin/linux64GccDPOpt/pimpleFoam" Gleitkomma-Ausnahme
Ist mein Netz vllt zu grob etc? Wenn ich direkt mit dem pimpleFoam anfang zu rechnen, dann klappt das, jedoch bin ich nach 1 Tag rechnen bei einer Zeit von 0.3 Sek. gewesen und am Strömungsprofil hat sich noch nicht viel getan.
hmmm ... :/
Falsche RB ? Grüße Tobi Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
tH3f0rC3 Mitglied
Beiträge: 47 Registriert: 12.03.2011
|
erstellt am: 05. Jun. 2011 22:27 <-- editieren / zitieren --> Unities abgeben: Nur für Shor-ty
1. Netz zu grob... -> Hast du mal checkMesh 2. Welches Turbulenzmodell verwendest du denn? Hast du da den dimensionslosen Wandabstand berücksichtigt, bzw., dir mit yPlusRAS angeschaut? Grüße tH3ff0rC3 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. Jun. 2011 22:37 <-- editieren / zitieren --> Unities abgeben:
Hey dirk, Problem schon lang gelöst Mein Netz war sehr fein... Hät das am nachlauf etwas gröber auführen dürfen. Aber es hat am Schluss funktioniert. Ich hatte ne falsche RB formuliert das war alles. Aber danke für die Info
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
| Anzeige.:
Anzeige: (Infos zum Werbeplatz >>)
|