Fenómenos de transporte

Conducción del calor

Simulación de la
conducción

Difusión unidimensional

Simulación de la difusión

Movimiento browniano

Sedimentación

Solución analítica

Solución numérica

Ley de Fourier

Sea J la densidad de corriente de energía (energía por unidad de área y por unidad de tiempo), que se establece en la barra debido a la diferencia de temperaturas entre dos puntos de la misma. La ley de Fourier afirma que hay una proporcionalidad entre el flujo de energía J y el gradiente de temperatura.

Siendo K una constante característica del material denominada conductividad térmica.

Consideremos un elemento de la barra de longitud dx y sección S. La energía que entra en el elemento de volumen en la unidad de tiempo es JS, y la que sale es J’S. La energía del elemento cambia, en la unidad de tiempo, en una cantidad igual a la diferencia entre el flujo entrante y el flujo saliente.

Esta energía, se emplea en cambiar la temperatura del elemento. La cantidad de energía absorbida o cedida (en la unidad de tiempo) por el elemento es igual al producto de la masa de dicho elemento por el calor específico y por la variación de temperatura.

Igualando ambas expresiones, y teniendo en cuenta la ley de Fourier, se obtiene la ecuación diferencial que describe la conducción térmica

Solución analítica

Supongamos una barra metálica de longitud L, conectada por sus extremos a dos focos de calor a temperaturas T_a y T_b respectivamente. Sea T₀ la temperatura inicial de la barra cuando se conectan los focos a los extremos de la barra.

Al cabo de cierto tiempo, teóricamente infinito, que en la práctica depende del tipo de material que empleamos, se establece un estado estacionario en el que la temperatura de cada punto de la barra no varía con el tiempo. Dicho estado está caracterizado por un flujo J constante de energía. La ley de Fourier establece que la temperatura variará linealmente con la distancia x al origen de la barra.

La temperatura en cualquier punto x a lo largo de la barra, en un instante determinado, T(x, t) es la solución de la ecuación diferencial, que es una combinación de dos términos, la que corresponde al régimen permanente más la del régimen transitorio.

Las condiciones de contorno, es decir, la temperatura T₀ en el instante inicial (t=0), y las temperaturas en los extremos T_a (para x=0) y T_b (para x=L) que permanecen invariables, nos permiten obtener los valores de los coeficientes k_n

Para n par

Para n impar

Así, la temperatura en cualquier punto de la barra x, en un instante t, se compone de la suma de un término proporcional a x, y de una serie rápidamente convergente que describe el estado transitorio.

Actividades

En este programa estudiaremos la conducción del calor a lo largo de una barra metálica cuyos extremos están conectados a dos focos de calor que tienen distintas temperaturas. Observaremos la evolución de la temperatura de cada punto de la barra a medida que transcurre el tiempo.

Examinaremos los factores que determinan la conducción del calor a lo largo de una barra metálica, probando barras de distintos materiales, con distintas temperaturas fijas de los extremos e inicial de la barra.

• Seleccionar en el control selección el Metal de la barra. Las unidades de las magnitudes están expresadas en el Sistema Internacional de Unidades de Medida.

• Se introduce la temperatura fija en el extremo izquierdo de la barra T_a.
• La temperatura fija en el extremo derecho de la barra T_b.
• La temperatura inicial de la barra T₀.
• El instante t en el que queremos representar la distribución de temperaturas a lo largo de la barra, y a continuación, se pulsa en el botón titulado Gráfica.
• Se introduce otro instante t en el que queremos observar la distribución de temperaturas, y a continuación, se pulsa en el botón titulado Gráfica, y así sucesivamente.

Nota: la longitud de la barra metálica se ha fijado en el valor de 0.5 m.

Cuestiones:

• Examinar la evolución de la distribución de temperaturas con el tiempo. Comprobar que el régimen estacionario es independiente de la temperatura inicial, solamente depende de la temperatura de los focos frío y caliente.

• Examinar el comportamiento de barras hechas de distintos materiales, con la misma temperatura inicial y fijas en los extremos.

Solución numérica

Para simplificar, tomemos una escala de tiempos tal que x=a²t, la ecuación diferencial que describe la conducción térmica se transforma en otra más sencilla.

La solución de la ecuación diferencial en derivadas parciales puede obtenerse utilizando el procedimiento de redes del siguiente modo

Consideremos un sistema de coordenadas posición – tiempo (x en el eje horizontal y x en el vertical). Construyamos un retículo trazando rectas paralelas al eje x equidistantes un intervalo fijo h y tal que h=L/n, siendo L la longitud de la barra y n el número de intervalos en los que se ha dividido la barra. Tracemos rectas paralelas al eje X, equidistantes una cantidad k.

Podemos calcular la temperatura en los puntos de la barra x=ih (i=1, 2, 3, 4...n) y en el instante x=(j+1)k, a partir de los datos de la temperatura de la barra en los puntos x de la barra en el instante anterior x=jk (j=0, 1, 2, 3...) sin más que aplicar el procedimiento de recurrencia esquematizado en la figura, y cuya fórmula es

La distribución inicial de partida (j=0) está dada por la temperatura inicial de la barra T₀, y las temperaturas fijas en los extremos T_a y T_b a partir de las cuales y aplicando la fórmula de recurrencia, pude calcularse, sucesivamente, las temperaturas de cada una de los puntos de la malla (i, j).

Para practicar este método con una calculadora o con un pequeño programa de ordenador, se puede seguir el siguiente procedimiento:

1. Constrúyase una tabla como la indicada y rellénese la columna izquierda, la columna derecha y la fila inferior con las temperaturas fijas en el extremo izquierdo de la barra T_a, del extremo derecho de la barra T_b, y con la temperatura inicial T₀.

2. Completar la primera fila vacía aplicando la fórmula de recurrencia. Las cifras obtenidas corresponden a la distribución de temperatura en un instante posterior al inicial.

3. A partir de los datos de la primera fila, completar la segunda fila vacía, y así sucesivamente.