Hot News:

Unser Angebot:

  Foren auf CAD.de (alle Foren)
  OpenFOAM
  reactingFoam

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:  reactingFoam (3247 mal gelesen)
rauhpeha
Mitglied


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

Beiträge: 4
Registriert: 07.06.2011

erstellt am: 07. Jun. 2011 09:09    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

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


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

Beiträge: 353
Registriert: 02.10.2008

BIM HVACTool für Windows OpenFOAM-2.2.x

erstellt am: 07. Jun. 2011 19:40    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 rauhpeha 10 Unities + Antwort hilfreich

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


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

Beiträge: 4
Registriert: 07.06.2011

erstellt am: 09. Jun. 2011 11: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

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


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

Beiträge: 353
Registriert: 02.10.2008

BIM HVACTool für Windows OpenFOAM-2.2.x

erstellt am: 09. Jun. 2011 13:21    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 rauhpeha 10 Unities + Antwort hilfreich

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


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

Beiträge: 4
Registriert: 07.06.2011

erstellt am: 10. Jun. 2011 11:14    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


0.zip

 
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



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

Beiträge: 57
Registriert: 07.02.2006

erstellt am: 28. Jun. 2011 21:39    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 rauhpeha 10 Unities + Antwort hilfreich

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


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

Beiträge: 4
Registriert: 07.06.2011

erstellt am: 04. Jul. 2011 18:02    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

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



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

Beiträge: 57
Registriert: 07.02.2006

erstellt am: 03. Aug. 2011 16:13    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 rauhpeha 10 Unities + Antwort hilfreich

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

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