Hot News:

Unser Angebot:

  Foren auf CAD.de (alle Foren)
  OpenFOAM
  Vorgabewerte via SimpleFoam

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:  Vorgabewerte via SimpleFoam (2037 mal gelesen)
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. Mai. 2011 20:37    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


netz.png

 
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 s

Courant 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



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

Beiträge: 47
Registriert: 12.03.2011

erstellt am: 05. Jun. 2011 22:27    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 Shor-ty 10 Unities + Antwort hilfreich

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





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: 05. Jun. 2011 22:37    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

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 >>)

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