Difusión de calor en un horno.

Tengo este problema donde estoy precalentando un horno para alcanzar norte . Entonces, digamos que dentro del horno hay un elemento calefactor en la parte inferior, que a su vez calienta el resto del horno. Suponiendo que el aire en el horno se comporta como un conductor estándar (ignore la convección), ¿cómo puedo medir la temperatura en un punto dado? T y ¿cuál será la temperatura de equilibrio en un punto dado?

Lo que he hecho hasta ahora es hacer algunas suposiciones sobre el sistema.

  • Las paredes están perfectamente aisladas.
  • El horno es unidimensional ( z dirección)
  • El horno tiene longitud L y potencia de salida q
  • El elemento calefactor se coloca en z = 0

A partir de esas suposiciones, he logrado formular la siguiente ecuación usando la ecuación de difusión de calor donde α es la difusividad térmica del aire, C es la capacidad calorífica específica del aire, ρ es la densidad y d es la función delta de dirac.

tu t = α tu z z + q C ρ d ( z )
Y esta ecuación tendría las siguientes condiciones iniciales y de contorno donde T 0 es la temperatura ambiente inicial del horno.
tu ( z , 0 ) = T 0
tu z ( 0 , t ) = tu z ( L , t ) = 0
Este es el punto en el que no estoy seguro de haber realizado el enfoque correcto, ya que se siente incorrecto tener tu z ( 0 , t ) = 0 mientras que también tiene el elemento calefactor en z = 0 . Tal vez podría colocar el elemento calefactor en z = ϵ en cambio, pero no estoy seguro de cómo resolver esto y preferiría tenerlo en z = 0 . Intenté resolver la ecuación usando los métodos que encontré aquí , pero parece que tengo problemas con q ya que no se puede escribir como una suma en X norte ( X ) y q norte ( t ) ya que es solo una constante con la única dependencia espacial que es la función delta de Dirac.

También he intentado encontrar las temperaturas de equilibrio suponiendo un estado estacionario con tu t = 0 lo que da como resultado más problemas ya que no estoy seguro de cómo encontrar la integral indefinida de la función delta de Dirac (no estoy seguro si existe). He intentado usar simplemente la definición de integral definida pero eso da resultados contradictorios.

Si alguien pudiera proporcionar una respuesta detallada sobre cómo resolver el problema, o incluso un punto en la dirección correcta, se lo agradecería, ya que he estado reflexionando sobre esto durante un par de días, pero aún no puedo pensar en una solución.

Tienes un problema con tus suposiciones. Las paredes son perfectamente aislantes, lo que significa que el calor que entra se queda adentro. El elemento calefactor agrega cabeza para siempre. Esto significa que tendrá un horno que se calienta más y más sin alcanzar nunca una temperatura de equilibrio.
Además, no puede tener aire o calor en 1 dimensión.
Usar k tu z = q en z = 0, donde q es el flujo de calor, y pierde el término delta de Dirac en la ecuación diferencial.

Respuestas (1)

Un modelo simple puede adaptarse a su caso. Solucione una ecuación de difusión simple para la densidad de energía térmica tu ( z , t )

(1) tu t = D 2 tu z 2 .
Dónde D la constante de difusión. E imponer las condiciones de contorno:
D [ tu z ] L = 0 ; D [ tu z ] 0 = PAG
dónde PAG es el cable de calefacción de forma de generación de energía. Elaboraré los detalles a continuación.

effects of boundary condition

Para estudiar cómo se relaciona la condición de contorno con la potencia de calentamiento, integremos la ecuación (1) wrt z .

(2) 0 L d z { tu t = D 2 tu z 2 }

Definir tu ( t ) = 0 L tu ( z , t ) d z . Después de la integración, la Ec. (2) rinde:

(3) tu ( t ) t = D [ tu z ] L D [ tu z ] 0
Esta ecuación da un significado dinámico a la derivada de tu en los límites: las derivadas denotan la tasa de cambio de la energía total.

El z = L es un medio aislante para fijar [ tu z ] L = 0 , todas las energías térmicas llegan a z se reflejan: la pendiente desaparece allí y se conserva la energía.

Eliminar el límite en z = L , ecuación (3) se convierte

(4) tu ( t ) t = D [ tu z ] 0
Dónde tu ( t ) es la capacitancia de calor total en el tiempo t . Por lo tanto la pendiente:
(5) tu ( t ) t = D [ tu z ] 0 =  tasa de cambio de la energía térmica total.

En conclusión, lo que debe hacer es fijar las condiciones límite de Neumann:

D [ tu z ] L = 0

La derivada cero significa un límite de reflexión. El límite del aislamiento reflejará todo el flujo de calor hacia él.

y la condición de frontera de Neumann en z = 0 :

D [ tu z ] 0 = PAG

P es la potencia del generador de calor.

Establezca estas dos condiciones, luego resuelva la ecuación de difusión simple. Obtendrás lo que esperas.