Ayúdame a resolver un problema de transferencia de emisión/conducción de calor. Mathematica me ha fallado

Mi problema: un tubo de pared delgada (longitud L , diámetro D y espesor de pared t D ) está en el vacío. Se sujeta por un extremo (en X = 0 ) por una fuente de calor a temperatura constante T ( 0 ) = T 0 . La única forma en que puede disipar el calor es por radiación. Supongo que la emisión solo ocurre desde la superficie exterior del tubo. La conductividad del tubo es k en [ W / metro k ] y la emisividad ϵ . ¿Cuál es el perfil de temperatura de equilibrio? T ( X ) en el tubo? (una aproximación numérica servirá).

Mi intento:

En un estado estable,

q i norte = q o tu t

A partir de la ley de conducción térmica de Fourier, el calor que entra por la sección final es

q i norte = k d T d X | X = 0 × π D t

A partir de la ley de Stefan-Boltzmann de radiación de cuerpo negro, el calor disipado a través de la superficie exterior del tubo viene dado por

q o tu t = 0 L ϵ σ T 4 d X × π D

Igualando los dos, el problema se convierte en

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

Tratar de resolver esto en Mathematica es inútil. ¿Estoy haciendo algo mal? ¿Cómo puedo encontrar una forma diferencial local de la ecuación? ¿Puedo simplificarlo más?

Gracias por tu ayuda.

El problema es que realmente no has construido un balance de calor en absoluto. Ver: stealthskater.com/Personal/Thesis.pdf , página 3. Aunque para una caña, se puede adaptar para un tubo con bastante facilidad.
Gracias por la referencia, es muy relevante. Sin embargo, ¿cómo no he construido un balance de calor? Cómo es q i norte = q o tu t ¿No es la condición para un estado estacionario?

Respuestas (1)

Debe hacer un balance de calor diferencial en un pequeño segmento del tubo entre x y x + Δ X .

Calentar en x = π D t k ( T X ) X

Calentar en x + Δ X = + π D t k ( T X ) X + Δ X

Calor perdido por radiación = π D Δ X ϵ σ T 4

Ecuación de balance de calor:

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

dividiendo por Δ X y tomando el límite como Δ X se aproxima a cero da:

k t 2 T X 2 = ϵ σ T 4

¿Qué en efecto es solo diferenciar la ecuación original y ser más cuidadoso con los signos?
Probablemente no existan formas exactas para la solución de esta ecuación, pero para la intuición física, esto será equivalente a una partícula que se mueve bajo la influencia de un potencial. tu ( X ) X 5 .
Gracias @Farcher. Edité las palabras calentar para leer calor en x + Δ X . Creo que este fue el único problema de señal. Me quedo con el resto de los signos en las ecuaciones.
@ChesterMiller De ninguna manera estaba criticando su derivación. Todo lo que estaba señalando era que la derivación original tiene una diferencia de signo con la tuya.
@MichaelSeifert mi intuición parece estar un poco oxidada, ¿puedes explicar eso un poco más (¿cómo la trayectoria de una partícula en U(x) = -ax^5 potencial es equivalente a la solución?) Gracias
@uhoh: si una partícula se moviera en tal potencial, entonces obedecería la ecuación de movimiento metro X ¨ = d tu / d X = 5 a X 4 . Reemplazar X T , t X , k t metro , y 5 a ϵ σ en esta ecuación, y terminas con la misma ecuación diferencial.
@Farcher No veo ninguna inconstancia entre los signos en la derivación original y en mi resultado final.
@MichaelSeifert ¿no es la solución al problema la evolución temporal de una distribución de temperatura: T (x, t) que, con suerte, tendrá una forma de equilibrio asintóticamente constante a medida que t tiende al infinito? ¿Cómo ayuda aquí una sola partícula que acelera rápidamente hacia la izquierda x = f(t)?