Hallo,
ich habe hier ein Problem mit einer Berechnung in OpenFoam. Ich habe nun schon ein paar Tutorials hinter mir, aber es zeigt sich, dass eine selbst konstruierte Aufgabe so ihre Tücken hat.
Ich möchte den Druckverlust in einem Hosenrohr wie in der angehängten Abb. berechnen. Ich gehe von einer turbulenten Stömung aus. Die Reynoldszahl liegt bei 1.5e6. Das Modell konnte ich nach längerem Probieren mit snappyHexMesh vernetzen. Am den beiden Inlets sei die Stömungsgeschwindigleit U=53m/s der Eigangsdruck würde 5 bar betragen.
Wenn ich das ganze rechnen lasse, dann endet die Rechnung mit dem k-eps-Modell nach ca. 16 Rechenschritten. Auch das RNG k-eps-Modell liefert nach dem 16 Rechenschritt einen Abbruch:
Meldung: time step continuity errors : sum local = 4.769675e+51, global = 1.944662e+48, cumulative = 1.867398e+48
Wenn ich mir die k-Werte anschaue, werden diese immer größer. In ParaView die Vektoren anzeigen lasse, dann sieht es auch so aus, als würde am Auslass nichts herauskommen.
Hier noch meine Randbedingungen:
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 1.618;
boundaryField
(
hosenrohr_wall
{
type zeroGradient;
}
hosenrohr_inlet1
{
type fixedValue;
value uniform 1.618;
}
hosenrohr_inlet2
{
type fixedValue;
value uniform 1.618;
}
hosenrohr_outlet
{
type zeroGradient;
}
}
object epsilon;
}
dimensions [0 2 -3 0 0 0 0];
internalField uniform 12.074;
boundaryField
{
hosenrohr_wall
{
type zeroGradient;
}
hosenrohr_inlet1
{
type fixedValue;
value uniform 12.074;
}
hosenrohr_inlet2
{
type fixedValue;
value uniform 12.074;
}
hosenrohr_outlet
{
type zeroGradient;
}
}
object U;
}
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
hosenrohr_wall
{
type fixedValue;
value uniform (0 0 0);
}
hosenrohr_inlet1
{
type surfaceNormalFixedValue;
refValue uniform -41.778;
}
hosenrohr_inlet2
{
type surfaceNormalFixedValue;
refValue uniform -41.778;
}
hosenrohr_outlet
{
type zeroGradient;
}
}
object p;
}
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
hosenrohr_wall
{
type zeroGradient;
}
hosenrohr_inlet1
{
type fixedValue;
value uniform 0; // 2.7e5;
}
hosenrohr_inlet2
{
type fixedValue;
value uniform 0; // 2.7e5;
}
hosenrohr_outlet
{
type fixedValue;
value uniform 0;
}
}
Gibt es eventuell besser geeignete Boundary-Bedingung?
Besten Dank für die Unterstützung,
Andreas
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP