Die Stefan-Gleichung beschreibt die Beziehung der Schmierungstheorie zwischen Quetschkraf, Viskosität, Scheibengröße, Spalt und Annäherungsgeschwindigkeit. Für zwei parallele kreisförmige Scheiben mit dem Radius (R), die sich einander nähern, lautet die Gleichung unter Verwendung der Halbschicht-Konvention (H(t)):

[ F(t)=\frac{3\pi \mu U(t)R^4}{8H(t)^3}, ]

wobei (F) die Normalkraft, (\mu) die dynamische Viskosität, (U=-dH/dt) die Annäherungsgeschwindigkeit jeder Scheibe und (H) der halbe Abstand zwischen den Scheiben ist. In einer pädagogischen Pilotanlage ermöglicht die Messung von (F), (H) und (U) den Studierenden, die Schmierungstheorie zu testen und die Viskosität anhand von Kraft-Weg-Daten zu bestimmen.

Das zentrale Ergebnis ist die starke Abhängigkeit von der dritten Potenz des Kehrwerts des Spalts: Eine Verringerung des halben Spalts erhöht die erforderliche Quetschkraft um den Faktor (H^{-3}), während eine Vergrößerung des Scheibenradius die Kraft um den Faktor (R^4) erhöht. Die Gleichung ist zuverlässig, wenn der Spalt schmal ist, die Strömung langsam ist und viskose Effekte die Trägheit dominieren.

Was die Stefan-Gleichung beschreibt

Das physikalische Experiment

Das Experiment verwendet zwei parallele kreisförmige Scheiben, die durch eine dünne Flüssigkeitsschicht getrennt sind. Die Scheiben bewegen sich aufeinander zu und zwingen die Flüssigkeit, radial nach außen durch den enger werdenden Spalt zu fließen.

Der Versuchsaufbau stellt einen vereinfachten Quetschströmungsprozess dar, der für Schmierung, Polymerverarbeitung, Beschichtung, Umformung und rheologische Messungen relevant ist.

Die Schlüsselvariablen

Die Variablen sind:

(R): Scheibenradius
(H(t)): Halber Abstand zwischen den Scheiben
(U(t)): Annäherungsgeschwindigkeit jeder Scheibe
(F(t)): Angewandte Normalkraft
(\mu): Dynamische Viskosität
(\rho): Flüssigkeitsdichte
(r): Radialer Abstand vom Scheibenzentrum

Die gesamte Dicke der Flüssigkeitsschicht beträgt (2H). Diese Konvention ist wichtig, da alternative Referenzen den gesamten Spalt oft als (h=2H) definieren, was zu einem anders aussehenden numerischen Koeffizienten führt.

Die zugrunde liegenden Annahmen

Die Herleitung setzt Folgendes voraus:

Kriechströmung mit vernachlässigbarer Trägheit
Ein schmaler Spalt, (H_0/R \ll 1)
Newtonsches, inkompressibles Flüssigkeitsverhalten
Kein Schlupf an den Scheibenoberflächen
Parallele Scheiben und achsensymmetrische Strömung
Gleichmäßige Annäherungsgeschwindigkeit über jede Scheibe
Vernachlässigbare Randeffekte und Flüssigkeitsleckagen

Die relevante Reynolds-Zahl ist

[ Re=\frac{H_0U_0}{\nu}\ll 1, ]

wobei (\nu=\mu/\rho) die kinematische Viskosität ist.

Wie die Gleichung hergeleitet wird

Beginn mit der radialen Quetschströmung

Während sich die Scheiben nähern, kann sich die Flüssigkeit nicht durch die festen Scheiben bewegen. Sie bewegt sich daher radial nach außen.

Da der Spalt im Vergleich zum Scheibenradius klein ist, dominieren die Geschwindigkeitsgradienten über den Spalt die Geschwindigkeitsgradienten in radialer Richtung. Dies ist die zentrale Näherung der Schmierungstheorie.

Reduktion der Impulsgleichung

Unter Kriechströmungsbedingungen sind die Trägheitsterme in den Navier-Stokes-Gleichungen vernachlässigbar. Das radiale Impulsgleichgewicht reduziert sich auf

[ 0=-\frac{\partial P}{\partial r} +\mu\frac{\partial^2 u_r}{\partial z^2}, ]

wobei (u_r) die radiale Geschwindigkeit und (z) die Koordinate durch die Flüssigkeitsschicht ist.

Der Druck variiert hauptsächlich mit dem Radius, während die radiale Geschwindigkeit stark über den schmalen Spalt variiert.

Anwendung der Randbedingungen an den Scheiben

Die Haftbedingung (No-Slip) erfordert, dass die radiale Geschwindigkeit an beiden Scheibenoberflächen verschwindet:

[ u_r=0 \quad \text{bei} \quad z=-H \text{ und } z=H. ]

Die zweimalige Integration der vereinfachten Impulsgleichung über den Spalt ergibt ein parabolisches radiales Geschwindigkeitsprofil. Die mit diesem Profil verbundene Volumenstromrate muss der Flüssigkeitsmenge entsprechen, die durch die sich nähernden Scheiben verdrängt wird.

Verwendung der Kontinuitätsgleichung

Für zwei sich symmetrisch nähernde Scheiben nimmt der halbe Spalt mit der Geschwindigkeit (U) ab:

[ \frac{dH}{dt}=-U. ]

Der sich ändernde Spalt verdrängt Flüssigkeit mit einer Rate, die proportional zur lokalen Scheibenfläche ist. Die achsensymmetrische Kontinuitätsgleichung setzt dann die radiale Strömung mit der Rate der Spaltverringerung in Beziehung.

Die Kombination der Kontinuitätsgleichung mit dem parabolischen Geschwindigkeitsprofil ergibt den radialen Druckgradienten

[ \frac{\partial P}{\partial r}

-\frac{3\mu U r}{2H^3}. ]

Das negative Vorzeichen zeigt an, dass der Druck von der Mitte zum äußeren Rand hin abnimmt.

Integration des Druckfeldes

Unter der Annahme des Drucks am Scheibenrand als Referenzdruck, (P(R)=0), ergibt die Integration von (r) bis (R)

[ P(r)=\frac{3\mu U}{4H^3}\left(R^2-r^2\right). ]

Der Druck ist daher in der Mitte am höchsten und nimmt quadratisch zum Rand hin ab.

Integration des Drucks zur Ermittlung der Kraft

Die gesamte Normalkraft ist der über die kreisförmige Scheibenfläche integrierte Druck:

[ F=2\pi\int_0^R P(r),r,dr. ]

Das Einsetzen der Druckverteilung ergibt

[ F

2\pi\int_0^R \frac{3\mu U}{4H^3} \left(R^2-r^2\right)r,dr, ]

und somit

[ \boxed{ F=\frac{3\pi\mu U R^4}{8H^3} }. ]

Dies ist die Stefan-Gleichung unter der angegebenen Halbspalt- und symmetrischen Annäherungskonvention.

Warum die Skalierung wichtig ist

Die Kraft nimmt bei schließendem Spalt rapide zu

Die Gleichung enthält den Faktor

[ F\propto \frac{1}{H^3}. ]

Eine Halbierung des halben Spalts erhöht die vorhergesagte Kraft um den Faktor acht, wenn (U), (R) und (\mu) konstant bleiben.

Diese starke Abhängigkeit ist der Grund, warum kleine Fehler bei der Spaltmessung zu großen Fehlern bei der berechneten Viskosität oder der vorhergesagten Kraft führen können.

Der Scheibenradius hat einen noch stärkeren Effekt

Die Kraft skaliert mit der vierten Potenz des Radius:

[ F\propto R^4. ]

Eine größere Scheibe erzeugt einen wesentlich höheren Widerstand, da sie mehr Flüssigkeit verdrängt und den radialen Abstand vergrößert, über den der Druck wirkt.

Viskosität und Annäherungsgeschwindigkeit haben lineare Effekte

Die Kraft ist direkt proportional sowohl zur Viskosität als auch zur Annäherungsgeschwindigkeit:

[ F\propto \mu U. ]

Eine viskosere Flüssigkeit oder eine schnellere Kompression erfordert eine größere Normalkraft, vorausgesetzt, die Flüssigkeit bleibt Newtonsch und die Kriechströmungsannahmen bleiben gültig.

Anwendung der Gleichung in einer Pilotanlage

Vorbereitung und Inspektion der Apparatur

Installieren Sie die kreisförmigen Scheiben, stellen Sie sicher, dass sie parallel sind, und messen Sie den Scheibenradius. Füllen Sie den Spalt vorsichtig, um eingeschlossene Luft zu vermeiden, da Luftkompressibilität die Kraftmessungen verfälschen kann.

Die Temperatur sollte überwacht werden, da die Viskosität signifikant mit der Temperatur variieren kann.

Steuerung der Bewegung

Legen Sie ein bekanntes Annäherungsgeschwindigkeitsprofil (U(t)) fest oder erzwingen Sie eine kontrollierte Weg-Zeit-Historie (H(t)). Die Steuerung sollte die tatsächliche Scheibenposition aufzeichnen, anstatt sich nur auf die Soll-Position zu verlassen.

Der halbe Spalt wird aus dem gemessenen Scheibenabstand berechnet:

[ H(t)=\frac{\text{gesamter Scheibenabstand}}{2}. ]

Aufzeichnung von Kraft und Weg

Zeichnen Sie während des Quetschvorgangs Folgendes auf:

Normalkraft (F(t))
Scheibenposition oder Spalt (H(t))
Annäherungsgeschwindigkeit (U(t))
Flüssigkeitstemperatur
Zeit

Eine Kraft-Weg-Kurve ist besonders nützlich, da sie zeigt, wie der Widerstand steigt, wenn sich der Spalt verengt.

Schätzung der Viskosität

Durch Umstellen der Stefan-Gleichung erhält man

[ \boxed{ \mu

\frac{8F H^3}{3\pi U R^4} }. ]

Studierende können eine momentane scheinbare Viskosität aus jedem gemessenen Datenpunkt berechnen oder das Modell an den vollständigen Kraft-Zeit- oder Kraft-Spalt-Datensatz anpassen.

Bei konstanter Viskosität sollten die berechneten Werte über den gültigen Betriebsbereich hinweg annähernd konsistent bleiben.

Validierung der vorhergesagten Skalierung

Das Experiment kann das Modell testen, indem ein Parameter nach dem anderen variiert wird:

Tragen Sie (F) gegen (U) auf, um die lineare Geschwindigkeitsabhängigkeit zu testen.
Tragen Sie (F) gegen (\mu) auf, indem Sie Flüssigkeiten unterschiedlicher Viskosität verwenden.
Tragen Sie (F) gegen (H^{-3}) auf, um die Spaltabhängigkeit zu testen.
Vergleichen Sie Scheiben unterschiedlicher Radien, um die (R^4)-Abhängigkeit zu testen.

Ein doppelt-logarithmisches Diagramm kann die vorhergesagten Potenzgesetz-Exponenten experimentell leichter identifizierbar machen.

Verständnis der Kompromisse

Die Schmalspalt-Näherung kann versagen

Das Modell erfordert

[ \frac{H_0}{R}\ll 1. ]

Wenn der Spalt im Verhältnis zum Scheibenradius nicht ausreichend klein ist, werden radiale Effekte und Randeffekte wichtig, und die eindimensionale Schmierungstheorie verliert an Genauigkeit.

Trägheit kann Hochgeschwindigkeitstests verfälschen

Die Kriechströmungsbedingung erfordert

[ Re=\frac{H_0U_0}{\nu}\ll 1. ]

Bei hohen Annäherungsgeschwindigkeiten oder Flüssigkeiten mit niedriger Viskosität können Trägheitskräfte messbar werden. Die angewandte Quetschkraft sollte die angegebene praktische Bedingung erfüllen

[ F\ll \frac{\mu^2}{\rho} \left(\frac{R}{H_0}\right)^4, ]

was dazu beiträgt, Trägheitsstörungen im Pilotanlagenexperiment gering zu halten.

Randleckagen reduzieren die gemessene Kraft

Die ideale Herleitung setzt eine klar definierte kreisförmige Grenze und einen Druck gleich dem Referenzdruck am Scheibenrand voraus. Wenn Flüssigkeit ungleichmäßig entweicht, sind das Druckfeld und die gemessene Kraft niedriger als vorhergesagt.

Mangelnde Parallelität kann die gleiche Art von Fehler erzeugen, indem ein ungleichmäßiger Spalt entsteht.

Newtonsches Verhalten trifft möglicherweise nicht zu

Die Stefan-Gleichung verwendet eine konstante dynamische Viskosität. Polymerschmelzen, konzentrierte Suspensionen, Gele und viele strukturierte Flüssigkeiten können scherverdünnend, scherverdickend, viskoelastisch sein oder eine Fließgrenze aufweisen.

Für diese Materialien kann die gemessene Kraft eine scheinbare oder ratenabhängige Viskosität darstellen und keine einzelne Materialkonstante.

Die Spaltmessung ist besonders empfindlich

Da die Kraft mit (H^{-3}) variiert, kann ein kleiner Nullpunktfehler oder eine Unsicherheit bei der Scheibentrennung einen großen Viskositätsfehler verursachen.

Die Apparatur sollte vor der Datenerfassung auf Weg, Kraft, Scheibenparallelität und den anfänglichen Spalt kalibriert werden.

Transiente Effekte können von Bedeutung sein

Die Gleichung beschreibt die viskose Quetschantwort unter quasi-stationären Annahmen. Zu Beginn und am Ende eines Bewegungszyklus können Stellantriebsbeschleunigung, Ansprechzeit des Instruments, eingeschlossene Luft und Flüssigkeitsträgheit transiente Abweichungen erzeugen.

Daten sollten in einem Bereich interpretiert werden, in dem die befohlene Bewegung und die gemessene Bewegung übereinstimmen.

Anwendung auf Ihr Projekt

Das zuverlässigste pädagogische Experiment kombiniert kontrollierte Bewegung, genaue Spaltmessung, Temperaturkontrolle und dimensionslose Überprüfungen, bevor die Kraftdaten interpretiert werden.

Wenn Ihr Hauptfokus auf der Validierung der Schmierungstheorie liegt: Verwenden Sie eine niedrigviskose oder gut charakterisierte Newtonsche Flüssigkeit, ein kleines Spalt-Radius-Verhältnis und mehrere Annäherungsgeschwindigkeiten, um (F\propto U) und (F\propto H^{-3}) zu testen.
Wenn Ihr Hauptfokus auf der Viskositätsmessung liegt: Stellen Sie die Gleichung um, um (\mu=8FH^3/(3\pi U R^4)) zu berechnen, und vergleichen Sie das Ergebnis mit einer unabhängigen Viskosimeter-Messung.
Wenn Ihr Hauptfokus auf der Untersuchung von Skaleneffekten liegt: Ändern Sie den Scheibenradius oder den anfänglichen Spalt systematisch und vergleichen Sie die gemessene Kraft mit den vorhergesagten (R^4)- und (H^{-3})-Abhängigkeiten.
Wenn Ihr Hauptfokus auf der Polymer- oder Nicht-Newtonschen Verarbeitung liegt: Behandeln Sie das Ergebnis als Messung der scheinbaren Viskosität und untersuchen Sie, wie sich die abgeleitete Viskosität mit der Annäherungsgeschwindigkeit und dem Spalt ändert.
Wenn Ihr Hauptfokus auf der Gewinnung belastbarer Daten liegt: Bestätigen Sie (Re\ll1), (H_0/R\ll1), Scheibenparallelität, Flüssigkeitseinschluss und Temperaturstabilität, bevor Sie das Stefan-Modell anpassen.

Wenn die Geometrie und die Betriebsbedingungen innerhalb ihrer Annahmen gehalten werden, bietet die Stefan-Gleichung eine direkte Brücke zwischen der gemessenen Quetschkraft und der viskositätsgetriebenen Physik der Schmalspaltströmung.

Zusammenfassungstabelle:

Thema	Wichtige Informationen
Gleichung	F = 3πμUR⁴/(8H³)
Hauptvariablen	F: Kraft; μ: Viskosität; U: Annäherungsgeschwindigkeit; R: Scheibenradius; H: halber Spalt
Kernskalierung	Kraft nimmt mit μ, U und R⁴ zu und steigt stark als H⁻³
Hauptannahmen	Kriechströmung, schmaler Spalt, Newtonsche Flüssigkeit, kein Schlupf, parallele Scheiben, vernachlässigbare Leckage
Experimentelle Messungen	Normalkraft, Scheibenabstand, Annäherungsgeschwindigkeit, Temperatur und Zeit
Viskositätsberechnung	μ = 8FH³/(3πUR⁴)
Validierungstests	Variieren Sie Geschwindigkeit, Viskosität, Spalt oder Scheibenradius und vergleichen Sie die gemessenen Potenzgesetz-Trends
Häufige Fehlerquellen	Eingeschlossene Luft, mangelnde Parallelität, Spaltmessfehler, Randleckagen, Trägheit und Temperaturänderungen
Pädagogische Anwendungen	Schmierungstheorie, Rheologie, Polymerverarbeitung, Beschichtung, Umformung und Ausbildung in Strömungsmechanik

Anwendung der Quetschströmungstheorie mit LABPARK

Bringen Sie die Stefan-Gleichung aus dem Lehrbuch in ein kontrolliertes, messbares Experiment. LABPARK bietet pädagogische und berufliche Pilotanlagen für Verfahrenstechnik für die chemische Verfahrenstechnik, Bioprozess- und Biotechnologie sowie für die Ausbildung in Umwelt- und Wasseraufbereitung.

Unsere Pilotanlagen helfen Universitäten, Forschungsinstituten und Unternehmen, Strömungsmechanik durch kontrollierte Bewegung, Kraft- und Wegmessung, Temperaturüberwachung und wiederholbare Betriebsbedingungen zu lehren und zu bewerten. Studierende und Forscher können die Schmierungstheorie validieren, Viskosität schätzen, Newtonsches und Nicht-Newtonsches Verhalten untersuchen und verstehen, wie Spalt, Scheibenradius und Annäherungsgeschwindigkeit die Prozessleistung beeinflussen.

Kontaktieren Sie LABPARK, um eine Pilotanlagenlösung zu besprechen, die auf Ihre Lehr-, Forschungs- oder industriellen Ausbildungsziele zugeschnitten ist.