Autor
|
Thema: reactingFoam (3247 mal gelesen)
|
rauhpeha Mitglied
Beiträge: 4 Registriert: 07.06.2011
|
erstellt am: 07. Jun. 2011 09:09 <-- editieren / zitieren --> Unities abgeben:
Hallo zusammen, das Thema meiner Diplomarbeit ist eine Verbennungssimulation, wozu reactingFoam zum Einsatz kommt. Aum Tutorial habe ich mich entlang gehangelt, und meine Simulation so ähnlich wie möglich zur Tutorial-Simulation aufgebaut. Dennoch geht es nicht. Kann mir vielleicht jemand helfen? Leider kann ich nicht einmal die unten stehende Fehlermeldung interpretieren. So einen Fehler hatte ich bisher noch nicht. Create time Create mesh for time = 0 Reading chemistry properties Reading g
Reading thermophysicalProperties Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics> Selecting thermodynamics package hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>> Selecting chemistryReader foamChemistryReader Selecting chemistrySolver ode Selecting ODE solver SIBS ODEChemistryModel: Number of species = 5 and reactions = 1 Reading field U Reading/calculating face flux field phi Creating turbulence model. Selecting turbulence model type RASModel Selecting RAS turbulence model kEpsilon #0 Foam::error: rintStack(Foam::Ostream&) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so" #2 in "/lib/libc.so.6" #3 Foam::compressible::RASModels::mutWallFunctionFvPatchScalarField::calcMut() const in "/opt/openfoam170/lib/linux64GccDPOpt/libcompressibleRASModels.so" #4 Foam::compressible::RASModels::mutWallFunctionFvPatchScalarField::updateCoeffs() in "/opt/openfoam170/lib/linux64GccDPOpt/libcompressibleRASModels.so" #5 Foam::fvPatchField<double>::evaluate(Foam::Pstream::commsTypes) in "/opt/openfoam170/applications/bin/linux64GccDPOpt/reactingFoam" #6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/opt/openfoam170/applications/bin/linux64GccDPOpt/reactingFoam" #7 Foam::compressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libcompressibleRASModels.so" #8 Foam::compressible::RASModel::adddictionaryConstructorToTable<Foam::compressible::RASModels::kEpsilon>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libcompressibleRASModels.so" #9 Foam::compressible::RASModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libcompressibleRASModels.so" #10 Foam::compressible::turbulenceModel::addturbulenceModelConstructorToTable<Foam::compressible::RASModel>::NewturbulenceModel(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libcompressibleRASModels.so" #11 Foam::compressible::turbulenceModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libcompressibleTurbulenceModel.so" #12 in "/opt/openfoam170/applications/bin/linux64GccDPOpt/reactingFoam" #13 __libc_start_main in "/lib/libc.so.6" #14 in "/opt/openfoam170/applications/bin/linux64GccDPOpt/reactingFoam" Floating point exception
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
TTB Mitglied CFD Engineer
Beiträge: 353 Registriert: 02.10.2008 BIM HVACTool für Windows OpenFOAM-2.2.x
|
erstellt am: 07. Jun. 2011 19:40 <-- editieren / zitieren --> Unities abgeben: Nur für rauhpeha
Hallo rauhpeha, du hast irgendwo eine Division durch Null: "Floating point exception". Vergleiche mal ganz genau deine Randbedingungen im Zeitordner Null mit eines aus dem Tutorial, welches reactingFoam verwendet. Evtl. hast du irgendwo als Beispiel eine Null für den Druck im "internalField", obwohl ein Wert ungleich Null gefordert wird. Des Weiteren verwende auch mal für k und Epsilon kleine Werte als "internalField", sprich keine Null als Initialisierungen. Gruß Thomas Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
rauhpeha Mitglied
Beiträge: 4 Registriert: 07.06.2011
|
erstellt am: 09. Jun. 2011 11:27 <-- editieren / zitieren --> Unities abgeben:
Danke, der Druck stand noch auf 0. Das war der Fehler. Kannst du mir vielleicht noch etwas weiterhelfen? reactingFoam läuft jetzt, zumindest bricht es (zum jetzigen Zeoitpunkt) nicht ab. Aber es findet bisher scheinbar keine Zündung statt, denn die Temperaturen verändern sich nicht, und das rho-Residuum ist und bleibt 0. Was ist hier zu tun? mfg Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
TTB Mitglied CFD Engineer
Beiträge: 353 Registriert: 02.10.2008 BIM HVACTool für Windows OpenFOAM-2.2.x
|
erstellt am: 09. Jun. 2011 13:21 <-- editieren / zitieren --> Unities abgeben: Nur für rauhpeha
Hallo, ohne dein Case zu kennen, ist eine Ferndiagnose immer schwierig. @Rho-Residuum: Meinst du dieses:
Code: "diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0"
Bei manchen Solver ist das Normal. Daran sollte es nicht liegen, ich vermute eher, dass du noch irgendwo bei den Reaktionen oder im Zeitordner "0" ein Fehler hast. Gruß Thomas Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
rauhpeha Mitglied
Beiträge: 4 Registriert: 07.06.2011
|
erstellt am: 10. Jun. 2011 11:14 <-- editieren / zitieren --> Unities abgeben:
Hallo nochmal, hoffentlich hat der upload funktioniert. Wenn nicht, kann ich dir alle Textdateien gern als mail schicken; die Gitterdateien sind 200MB groß; das wäre zu viel. Kann man irgendwo einstellen, dass es ab einer bestimmten Temperatur eine Zündung gibt? Weiß openFoam, dass Methan bei 600°C zündet? Kann ich im Volumen eine "Zündstelle" vorgeben? Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
pajofego Mitglied
Beiträge: 57 Registriert: 07.02.2006
|
erstellt am: 28. Jun. 2011 21:39 <-- editieren / zitieren --> Unities abgeben: Nur für rauhpeha
Stimmen die Arrhenius Koeffizienten und deren Einheit? Schon mal ein chemkin file mit folgender Reaktionskinetik ausprobiert? ELEMENTS H O C N END SPECIES CH4 O2 N2 CO2 H2O END REACTIONS CH4 + 2O2 => CO2 + 2H2O 6.70091E+12 0.0 4.84149E+04! A, b, E_a FORD / CH4 0.2 / FORD / O2 1.3 / END FORD / CH4 0.2 / FORD
[Diese Nachricht wurde von pajofego am 28. Jun. 2011 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
rauhpeha Mitglied
Beiträge: 4 Registriert: 07.06.2011
|
erstellt am: 04. Jul. 2011 18:02 <-- editieren / zitieren --> Unities abgeben:
Hallo, verzeih die späte Meldung. Ich hab einiges ausprobiert und viel gelesen.... Trotz einiger Anstrengungen zündet die Simulation leider nicht. Am Rand lass ich die Temperatur schon auf 1.500K ansteigen, um über den Flammpunkt zu kommen. Doch das ist auch noch nicht wirklich erfolgsträchtig. Meine chemkin-Datei habe ich aus den Tutorials übernommen. Du findest sie nachfolgend. species ( O2 H2O CH4 CO2 N2 ); reactions ( irreversibleArrheniusReaction CH4 + 2O2 = CO2 + 2H2O (5.2e16 0 14906) ); Sind die unteschiedlichen Faktoren wegen der unterschiedlichen Einheiten (SI- gegen amerikanische Einheiten)? Ich hab auch eine combustion-Datei eingefügt. Doch ist mir bei einigen Einträgen nicht klar, wofür sie stehen. Ist die combustionProperties-Datei überhaupt der richtige Ansatz? Danke für deine Antworten Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
pajofego Mitglied
Beiträge: 57 Registriert: 07.02.2006
|
erstellt am: 03. Aug. 2011 16:13 <-- editieren / zitieren --> Unities abgeben: Nur für rauhpeha
Versuch doch mal bitte den von mir geposteten Reaktionsmechanismus und schau doch mal bitte ob du eine Reaktion festellen kannst. Darüber hinaus ist die Simulation sehr rechenintensiv, d.h. lange Rechenzeiten. Bei mir läuft mein Fall seit 4 Tagen ohne Reaktion und ist erst bei einem Zeitpunkt von 0.18 sec, dabei hat sich noch keine quasistationäre Strömung ausgebildet. Verwende auch die simplefunctions um Werte an den Patches über die Zeit zu beobachten. Wann stellt sich eine stationäre Temperatur am Auslass bei dir ein und etc.? Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
| Anzeige.:
Anzeige: (Infos zum Werbeplatz >>)
|