Helfen Sie mir, ein Wärmeleitungs-/Emissionsübertragungsproblem zu lösen. Mathematica hat mich enttäuscht

Mein Problem: Ein dünnwandiger Schlauch (Länge L , Durchmesser D und Wandstärke T D ) befindet sich in einem Vakuum. Es wird an einem Ende (at X = 0 ) durch eine Wärmequelle bei konstanter Temperatur T ( 0 ) = T 0 . Die einzige Möglichkeit, Wärme abzuleiten, ist Strahlung. Ich gehe davon aus, dass die Emission nur von der Außenfläche des Rohrs erfolgt. Die Leitfähigkeit des Rohres ist k In [ W / M K ] und der Emissionsgrad ϵ . Was ist das Gleichgewichtstemperaturprofil? T ( X ) im Rohr? (eine numerische Annäherung reicht).

Mein Versuch:

In einem stationären Zustand,

Q ich N = Q Ö u T

Aus dem Fourierschen Gesetz der Wärmeleitung ergibt sich die durch den Endabschnitt eintretende Wärme

Q ich N = k D T D X | X = 0 × π D T

Aus dem Stefan-Boltzmann-Gesetz der Schwarzkörperstrahlung ergibt sich die durch die Außenfläche der Röhre abgeführte Wärme durch

Q Ö u T = 0 L ϵ σ T 4 D X × π D

Wenn man die beiden gleichsetzt, wird das Problem

k T ϵ σ D T D X | X = 0 = 0 L T 4 D X ,       T ( 0 ) = T 0

Der Versuch, dies in Mathematica zu lösen, ist hoffnungslos. Mache ich etwas falsch? Wie finde ich eine lokale Differentialform der Gleichung? Kann ich es weiter vereinfachen?

Vielen Dank für Ihre Hilfe.

Das Problem ist, dass Sie überhaupt keine Wärmebilanz erstellt haben. Siehe: Stealthskater.com/Personal/Thesis.pdf , Seite 3. Obwohl es sich um eine Rute handelt, kann es ganz einfach für ein Rohr angepasst werden.
Danke für den Hinweis, er ist sehr relevant. Aber wie habe ich keine Wärmebilanz erstellt? Wie ist Q ich N = Q Ö u T nicht die Bedingung für einen stationären Zustand?

Antworten (1)

Sie müssen eine unterschiedliche Wärmebilanz auf einem kleinen Segment des Rohrs zwischen x und x + durchführen Δ X .

Bei x = einheizen π D T k ( T X ) X

Bei x + einheizen Δ X = + π D T k ( T X ) X + Δ X

Wärmeverlust durch Strahlung = π D Δ X ϵ σ T 4

Wärmebilanzgleichung:

+ π D T k ( T X ) X + Δ X π D T k ( T X ) X = π D Δ X ϵ σ T 4

Teilen durch Δ X und nehmen die Grenze als Δ X nähert sich Null ergibt:

k T 2 T X 2 = ϵ σ T 4

Was bedeutet eigentlich nur, die ursprüngliche Gleichung zu differenzieren und vorsichtiger mit Vorzeichen umzugehen?
Genaue Formen für die Lösung dieser Gleichung gibt es wahrscheinlich nicht, aber für die physikalische Intuition entspricht dies einem Teilchen, das sich unter dem Einfluss eines Potentials bewegt U ( X ) X 5 .
Danke @Farcher. Ich habe die Wörter „heat out“ bearbeitet, um „heat in“ bei x + einzulesen Δ X . Ich glaube, das war das einzige Zeichenproblem. Ich stehe zu den restlichen Zeichen in den Gleichungen.
@ChesterMiller In keiner Weise habe ich Ihre Ableitung kritisiert. Ich habe lediglich darauf hingewiesen, dass die ursprüngliche Ableitung einen Vorzeichenunterschied zu Ihrer aufweist.
@MichaelSeifert meine Intuition scheint etwas eingerostet zu sein, kannst du das etwas näher erklären (wie die Flugbahn eines Teilchens in U(x) = -ax^5 Potential der Lösung entspricht?) Danke
@uhoh: Wenn sich ein Teilchen in einem solchen Potential bewegen würde, würde es der Bewegungsgleichung gehorchen M X ¨ = D U / D X = 5 A X 4 . Ersetzen X T , T X , k T M , Und 5 A ϵ σ in dieser Gleichung, und Sie erhalten am Ende dieselbe Differentialgleichung.
@Farcher Ich sehe keine Inkonsistenz zwischen den Zeichen in der ursprünglichen Ableitung und in meinem Endergebnis.
@MichaelSeifert ist nicht die Lösung des Problems die zeitliche Entwicklung einer Temperaturverteilung: T (x, t), die hoffentlich eine asymtotisch konstante Gleichgewichtsform hat, wenn t gegen unendlich geht? Wie hilft hier ein einzelnes Teilchen, das schnell nach links x = f(t) beschleunigt?