---

Hallo zusammen,

ich habe eine einfache Kanalströmung, bei der ich in der Mitte mit createBaffles und der fan BC einen Drucksprung vorgebe. Dies klappt auch wunderbar bei z.B. simpleFoam. Verwende ich allerdings einen kompressiblen Solver wie rhoSimpleFoam oder rhoCentralFoam (soweit ich es nach meinen Tests festellen konnte liegt es daran), bricht die Simulation mit folgender Ausgabe sofort ab:

Code:

---

#0  Foam::error:  rintStack(Foam::Ostream&) in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::heThermo<Foam:  siThermo, Foam:  ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam:  erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::he(Foam::Field<double> const&, Foam::Field<double> const&, int) const in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4  Foam::heThermo<Foam:  siThermo, Foam:  ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam:  erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::init() in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5  Foam::heThermo<Foam:  siThermo, Foam:  ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam:  erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#6  Foam::hePsiThermo<Foam:  siThermo, Foam:  ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam:  erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::hePsiThermo(Foam::fvMesh const&, Foam::word const&) in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#7  Foam:  siThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam:  siThermo, Foam:  ureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam:  erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#8  Foam::autoPtr<Foam:  siThermo> Foam::basicThermo::New<Foam:  siThermo>(Foam::fvMesh const&, Foam::word const&) in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#9  Foam:  siThermo::New(Foam::fvMesh const&, Foam::word const&) in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#10 
in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12 
in "/home/lhoang/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
^C
[1]+  Gleitkomma-Ausnahme     rhoCentralFoam > log

---

Hat dies etwas mit den thermophysical properties zu tun?

Code:

---

thermoType
{
    type                hePsiThermo;
    mixture             pureMixture;
    transport           sutherland;//const;
    thermo              hConst;
    energy              sensibleInternalEnergy;//sensibleEnthalpy;
    equationOfState     perfectGas;
    specie              specie;
}

mixture
{
    specie
    {
        nMoles          1;
        molWeight       28.9;
    }
    thermodynamics
    {
        Cp              1006.4;
        Hf              0;
    }
    transport
    {
As 1.4792e-06;
Ts 116;
    }
}

---

Oder ist die fan BC nicht mit einem kompressiblen Solver kompatibel (aufgrund der Kompressibilität bzw. Dichteabhängigkeit)?

Viele grüße
Christof

[Diese Nachricht wurde von cristopf am 04. Okt. 2016 editiert.]

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hi Christof,

die Randbedingungen funktionieren für beide Varianten, wäre nicht gut wenns nur für Inkompressible Löser ging. Wie du sicher schon weißt, hast du eine Division mit Null. Aufgrund deiner Fehlermeldung (leider keine Angaben wann dies passiert - siehe auch die Forumsregeln) kann man dazu jetzt nicht viel sagen. Allerdings hoffe ich das du dein Druck auch entsprechend um 1bar erhöht hast, da ein p0 uniform 0 im kompressiblen Fall sicher nicht funktioniert. Bezüglich rhoCentralFoam kann ich kein Statement geben, da selber noch nie verwendet. Scheint aber eher ein Beispiel darzustellen, wie man einen eigenen Solver schreibt (siehe *.C File). Hat ja ein ganz spezielles Schemata. Sollte aber auch funktionieren, denke ich.

------------------
Viele Grüße,
Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hi Tobias,

achje danke für den Tipp. Bei Value bei der fan BC hatte ich tatsächlich noch 0 für den Druck stehen. Das habe ich leider komplett übersehen. Habe den Wert angepasst und es hat angefangen zu rechnen. Denke mal, dass dann auch die richtigen Ergebnisse geliefert werden. Das Thema hat sich damit also erledigt.

Vielen Dank!

Gruß Christof

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Ja wunderbar, dann bitte das Thema noch als gelöst markieren. Für diejenigen die es interessiert. Christof verwendet     equationOfState    perfectGas; entsprechend wird diese Gleichung verwendet »klick. Demnach ist die Dichte überall Null und dann muss man eben nur noch schauen wo die Dichte im Nenner steht und woala hat man seine Division durch Null.

------------------
Viele Grüße,
Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP