Mathematik online lernen im Mathe-Forum. Nachhilfe online
Startseite » Forum » Gleichung Umformen (Prandtl-Colebrook Fließformel)

Gleichung Umformen (Prandtl-Colebrook Fließformel)

Universität / Fachhochschule

Tags: Fließgeschwindigkeit, Gleichungen, Gleichungssystem, Prandtl Colebrook, Umformen, umstellen, Wolfram Alpha

 
Antworten Neue Frage stellen Im Forum suchen
Neue Frage
WolframA

WolframA aktiv_icon

18:07 Uhr, 10.08.2021

Antworten
Guten Tag!

Ich versuche die Prandtl Colebrook Fließformel nun schon seit einiger Zeit nach „IE“ umzuformen (für einen Nachweis).

Kollegen konnten mir bis jetzt leider nicht weiterhelfen und Wolfram Alpha auch nicht.


Die Formel lautet ( soll nach „IE“ umgeformt sein):

vm = -2*lg((2,51*v)/(d*wurzel(2*g*d*IE))+k/(3,71*d))*wurzel(2*g*IE*d)

Im Anhang findet man ebenfalls ein Bild der Formel.
Ich hoffe jemand kann mir weiterhelfen.
Schonmal danke im Voraus!

Screenshot_20210810-180609_Gallery

Für alle, die mir helfen möchten (automatisch von OnlineMathe generiert):
"Ich benötige bitte nur das Ergebnis und keinen längeren Lösungsweg."
Online-Nachhilfe in Mathematik
Antwort
supporter

supporter aktiv_icon

18:14 Uhr, 10.08.2021

Antworten
Diese Gleichung kann man algebraisch nicht nach I_E umstellen.
Verwende ein Näherungsverfahren, falls die anderen Variablen bekannt sind.
Antwort
pivot

pivot aktiv_icon

21:43 Uhr, 10.08.2021

Antworten
Hallo,

ich definiere x=2gIEd. Dann ist strukturell die Gleichung

a=b(cx+d)x

Im Prinzip ist das eine quadratische Gleichung die erst einmal nach x aufgelöst werden kann.

Gruß
pivot
WolframA

WolframA aktiv_icon

00:34 Uhr, 11.08.2021

Antworten
Danke für deine schnelle Antwort!

Ich habe einen Lösungsansatz nach deinem Vorschlag angefangen(im Anhang als Bild). Ich wüsste aber nicht wie ich jetzt auf x kommen kann da ich ja auch noch den Log drin habe den ich auflösen muss.


Screenshot_20210811-002600_WhatsApp
Antwort
pivot

pivot aktiv_icon

02:24 Uhr, 11.08.2021

Antworten
Ach, ich habe leider den Logarithmus übersehen. Da stimme ich supporter zu, dass die Gleichung nicht algebraisch lösbar ist. Dann doch eine Approximationsmethode wie das Newtonverfahren verwenden. Dazu würde ich auch die Gleichung strukuturell darstellen. Dafür erst einmal alles auf die linke Seite.

f(x)=2lg(Ax-0,5+B)Cx-0,5+E=0

Ich hoffe du kannst die Konstanten selbst identifizieren. So ist z.B. E=vM.
Wenn lg(x) für den 10er-Logarithmus steht, dann würde ich verwenden, dass lg(x)=ln(x)ln(10). Das macht die Ableitung leichter.

f(x)=2ln(Ax-0,5+B)ln(10)Cx-0,5+E=0
Antwort
Roman-22

Roman-22

04:35 Uhr, 11.08.2021

Antworten
Die Gleichung lässt sich nach IE auflösen, wenn man die Lambertsche W-Funktion verwendet. Allerdings ist diese keine elementare Funktion und in aller Regel auch nur näherungsweise bestimmbar. de.wikipedia.org/wiki/Lambertsche_W-Funktion
Wenn du allerdings ein Mathe-Programm oder eine Programmier-Toolbox verwendest, in der diese Funktion bereits implementiert ist, kann die folgende Darstellung Sinn machen:

IE=23,712Dlg(e)gkPr2W(10vm22,51kPrν23,71D2lg(e))

EDIT: Diese Umformung ist grundfalsch. Ich bin da offenbar von einer wesentlich anderen Grundformel ausgegangen. Die richtige Umstellung nach IE findet sich weiter unten in einer weiteren Antwort von mir.
WolframA

WolframA aktiv_icon

16:00 Uhr, 11.08.2021

Antworten
Danke für deine Antwort!

Ich müsste mir das mal genauer ansehen,da mir das noch nicht ganz klar ist.

@pivot
WolframA

WolframA aktiv_icon

16:07 Uhr, 11.08.2021

Antworten
@ Roman-22

Danke dafür das ist schonmal eine große Hilfe!

Ich hab im WolframAlpha Rechner meine Werte für W versucht zu rechnen leider scheint das aber für Werte wie zb. 0,0000036 nicht zu funktionieren.

Ich habe jetzt denke ich trotzdem eine funktionierende Methode mit der ich auf IE relativ genau kommen kann.

Trotzdem wäre eine "richtige" Lösung besser.
Kann man irgendwo irgendwie W für Werte um die 0,00000036 rechnen?


Antwort
pivot

pivot aktiv_icon

17:37 Uhr, 11.08.2021

Antworten
Ich habe für W(0.00000036)=0.00000036
www.wolframalpha.com/input/?i=lambertW%280.00000036%29

Das heißt 0.00000036exp(0.00000036) müsste 0.00000036 ergeben. Das schein so zu sein:
www.wolframalpha.com/input/?i=0.00000036*exp%280.00000036%29

Antwort
pivot

pivot aktiv_icon

19:02 Uhr, 11.08.2021

Antworten
Das liegt wohl daran, dass e0,000000036
Antwort
pivot

pivot aktiv_icon

19:04 Uhr, 11.08.2021

Antworten
Das liegt wohl daran, dass e0,00000036/approx1 ist.
Antwort
pivot

pivot aktiv_icon

19:05 Uhr, 11.08.2021

Antworten
Das liegt wohl daran, dass e0,000000361 ist.
WolframA

WolframA aktiv_icon

10:08 Uhr, 12.08.2021

Antworten
Ok anscheinend macht es einen Unterschied ob man nur W() eingibt oder lambertW().

Das hat jetzt geklappt!

Wie Roman-22 schon gesagt hat ist die Formel nur ungefähr richtig.
Da es in meinem Fall aber einen großen Unterschied macht ob da jetzt 0,007 oder 0,008 rauskommt ist die Formel leider ohne irgendwelche Faktoren "dazuzubasteln" nur bedingt nutzbar.

Das hat mir aber sehr weitergeholfen!
Falls jemand trotzdem eine noch genauere Formel hat wäre das super!
Antwort
Roman-22

Roman-22

13:19 Uhr, 12.08.2021

Antworten
> Falls jemand trotzdem eine noch genauere Formel hat wäre das super!
Die angegeben Formel ist schon ganz genau! Nur die Berechnung der konkreten Werte der Lambert-W Funktion ist halt nicht elementar mit ein paar Additionen oder Multiplikationen zu bewerkstelligen. Aber Rechenknechte wie Wolfram Alpha sollten dir die Werte dieser Funktion mit "beliebiger" Genauigkeit liefern können und sicher zwischen 0,007 und 0,008 unterscheiden können.
Und wenn du ohnedies WA nutzt, musst du ja nicht die angegebene Umformung nutzen, sondern kannst deine Gleichung direkt eingeben und numerisch lösen lassen. Wenn du die Sache selbst programmieren möchtest, musst du entweder einen einigermaßen stabilen Algorithmus zur Lösung nichtlinearer Gleichungen implementieren (besser: eine fertige Toolbox verwenden) oder einen entsprechenden Algorithmus zur Bestimmung der Werte der Lambert-W Funktion einbauen.

P.S.: Mein Rechenknecht liefert mir übrigens
B1
und ist sich damit mit WA ziemlich einig
B2


Antwort
HAL9000

HAL9000

16:43 Uhr, 12.08.2021

Antworten
Aus der definierenden Eigenschaft W(x)eW(x)=x der LambertW-Funktion kann man ja implizit die Ableitungen berechnen und damit auch eine Taylorentwicklung im Nullpunkt, deren Anfang lautet so

W(x)=x-x2+32x3+O(x4)

Reicht sicher aus für das vorliegende Zahlenbeispiel. ;-)

WolframA

WolframA aktiv_icon

10:06 Uhr, 13.08.2021

Antworten
Danke an alle die sich die Mühe gemacht haben mir zu Antworten!

@ Roman-22
Hm ich glaube irgendwas passt da trotzdem nicht ganz.
Ich habe mit der Prandtl Colebrook Formel mir einen Wert ausgerechnet und IE mit 0,005 angenommen, dann habe ich die Werte die ich verwendet habe in die Umgeformte Formel eingesetzt.
Das Problem ist das für IE nicht 0,005 rauskommt (nur mit "dazugebastelten Faktor").

Ich bin mir ziemlich sicher,dass die Formel richtig eingegeben wurde da es für die Prandtl Colebrook Formel eine Tabelle für vm gibt und die alle mit meinen Weren übereinstimmen.

Die Umgeformte Formel habe ich auch sicher richtig eingegeben.

Im Bild (im Anhang) sieht man meine Exceltabelle mit den Formeln.

Kann das sein das ich was übersehe?

Prandtl Colebrook nach IE mit Faktor
Antwort
Enano

Enano

13:30 Uhr, 13.08.2021

Antworten
Mit welchen Werten bist du denn auf vm=0,353ms gekommen, doch sicher nicht mit den darunter in der Tabelle aufgeführten oder hast du 0,00335(=0,0050,67) anstatt 0,005 für I_E eingesetzt, wie man gem. deiner künstlerisch wertvollen Zeichnung vermuten könnte?

Wie hast du den Füllungsgrad von 0,5 berücksichtigt?

Wo hast du denn den Faktor von 0,67 her bzw. wie lautet deine Bastelanweisung dafür?
Antwort
Roman-22

Roman-22

16:52 Uhr, 13.08.2021

Antworten
Wie du mit deinen im Bild ersichtlichen Werten auf die 0,353ms gekommen bist, ist mir auch nicht erklärlich. Wobei ich allerdings deinen künstlerischen Pfeilen nicht weiter gefolgt bin, sondern nur mit deinen Werten den zugehörigen vm Wert ausgerechnet habe:
B1

Zu deinem Wert vm=0,353ms würde der Wert IE=0,01735 passen.

Meine ursprünglich angegebene Formel mit der Lambert-W Funktion für IE ist allerdings leider tatsächlich falsch. Ich kann es nicht mehr genau nachvollziehen, aber ich bin damals offenbar von einer anderen Grundfunktion ausgegangen. Dass meine Formel falsch ist erkennt man u.a. ja auch daran, dass sie nicht mehr Einheitenkonsistent ist. Allein der falsche 10er-Exponent vm2 (müsste dimensionslos sein) hätte mich von einem Absenden des Beitrags abhalten müssen. Sorry für die gestiftete Verwirrung.

Aber wenn du ohnedies mit Excel arbeitest, kannst du doch einfach die Grundformel verwenden und dir IE mit dem Solver berechnen lassen.
Antwort
Roman-22

Roman-22

00:39 Uhr, 14.08.2021

Antworten
Der Vollständigkeit halber und als kleine Buße gebe ich hier die richtige Formel für IE an (im Bild unten rot hinterlegt) und deren Überprüfung anhand zweier konkreter numerischer Beispiele.
Natürlich spielt die Lambert-Funktion wieder eine zentrale Rolle.

Ich denke aber nach wie vor, dass in Excel die Verwendung der Ausgangsformel und des Solvers vermutlich der sinnvollste und auch einfachste Ansatz ist.

B
Diese Frage wurde automatisch geschlossen, da der Fragesteller kein Interesse mehr an der Frage gezeigt hat.