Ecuación general de transferencia de calor

En dinámica de fluidos, la ecuación general de transferencia de calor es una ecuación diferencial parcial no lineal que describe la producción específica de entropía en un fluido newtoniano sometido a conducción térmica y fuerzas viscosas:

$\rho T\frac{Ds}{Dt}=\mathrm{\nabla }\cdot (k\mathrm{\nabla }T)+\frac{\mu }{2}{(\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}-\frac{2}{3}{\delta }_{ij}\mathrm{\nabla }\cdot 𝐯)}^2+\zeta (\mathrm{\nabla }\cdot 𝐯{\mathbf{)}}^{\mathrm{𝟐}}$

donde $s$ es la entropía específica, $\rho$ es la densidad del fluido, $T$ es la temperatura del fluido, $D/Dt$ es la derivada material, $\kappa$ es la conductividad térmica, $\mu$ es la viscosidad dinámica, $\zeta$ es el segundo parámetro de Lamé, $𝐯$ es la velocidad de flujo, $\mathrm{\nabla }$ es el del utilizado para caracterizar el gradiente y la divergencia, y ${\delta }_{ij}$ es el delta de Kronecker.

Si la velocidad del flujo es despreciable, la ecuación general de transferencia de calor se reduce a la ecuación estándar del calor. También puede aplicarse a flujos estratificados en rotación, como los que se encuentran en la dinámica de fluidos geofísicos.

Derivación

Ampliación de la ecuación de energía del fluido ideal

Para un fluido viscoso newtoniano, las ecuaciones que rigen la conservación de la masa y el momento son la ecuación de continuidad y las ecuaciones de Navier-Stokes:$$\begin{array}{cc}\frac{\partial \rho }{\partial t}& =-\mathrm{\nabla }\cdot (\rho 𝐯)\\ \rho \frac{D𝐯}{Dt}& =-\mathrm{\nabla }p+\mathrm{\nabla }\cdot \sigma \end{array}$$donde $p$ es la presión y $\sigma$ es el tensor de tensión viscoso, con los componentes del tensor de esfuerzo viscoso dados por:$${\sigma }_{ij}=\mu (\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}-\frac{2}{3}{\delta }_{ij}\mathrm{\nabla }\cdot 𝐯)+\zeta {\delta }_{ij}\mathrm{\nabla }\cdot 𝐯$$La energía de una unidad de volumen del fluido es la suma de la energía cinética $\rho v^2/2\equiv \rho k$ y la energía interna $\rho \epsilon$ donde $\epsilon$ es la energía interna específica. En un fluido ideal, como el descrito por las ecuaciones de Euler, la conservación de la energía viene definida por la ecuación:$$\frac{\partial }{\partial t}[\rho (k+\epsilon )]+\mathrm{\nabla }\cdot [\rho 𝐯(k+h)]=0$$donde $h$ es la entalpía específica. Sin embargo, para que la conservación de la energía se cumpla en un fluido viscoso sujeto a conducción térmica, el flujo de energía debido a la advección $\rho 𝐯(k+h)$ debe complementarse con un flujo de calor dado por la ley de Fourier $𝐪=-\kappa \mathrm{\nabla }T$ y un flujo debido a la fricción interna $-\sigma \cdot 𝐯$. Entonces la ecuación general de conservación de la energía es:$$\frac{\partial }{\partial t}[\rho (k+\epsilon )]+\mathrm{\nabla }\cdot [\rho 𝐯(k+h)-\kappa \mathrm{\nabla }T-\sigma \cdot 𝐯]=0$$

Ecuación de la producción de entropía

Obsérvese que las relaciones termodinámicas para la energía interna y la entalpía vienen dadas por:$$\begin{array}{cc}\rho d\epsilon & =\rho Tds+\frac{p}{\rho }d\rho \\ \rho dh& =\rho Tds+dp\end{array}$$También podemos obtener una ecuación para la energía cinética tomando el producto punto de la ecuación de Navier-Stokes con la velocidad del flujo $𝐯$ para obtener:$$\rho \frac{Dk}{Dt}=-𝐯\cdot \mathrm{\nabla }p+v_i\frac{\partial {\sigma }_{ij}}{\partial x_j}$$El segundo término del lado derecho puede ampliarse de la siguiente manera:$$\begin{array}{cc}{v}_i\frac{\partial {\sigma }_{ij}}{\partial x_j}& =\frac{\partial }{\partial x_j}\left({\sigma }_{ij}{v}_i\right)-{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}\\ & \equiv \mathrm{\nabla }\cdot (\sigma \cdot 𝐯)-{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}\end{array}$$Con la ayuda de la relación termodinámica para la entalpía y el último resultado, podemos entonces poner la ecuación de la energía cinética en la forma:$$\rho \frac{Dk}{Dt}=-\rho 𝐯\cdot \mathrm{\nabla }h+\rho T𝐯\cdot \mathrm{\nabla }s+\mathrm{\nabla }\cdot (\sigma \cdot 𝐯)-{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}$$Ahora expandiendo la derivada temporal de la energía total, tenemos:$$\frac{\partial }{\partial t}[\rho (k+\epsilon )]=\rho \frac{\partial k}{\partial t}+\rho \frac{\partial \epsilon }{\partial t}+(k+\epsilon )\frac{\partial \rho }{\partial t}$$Entonces expandiendo cada uno de estos términos, encontramos que:$$\begin{array}{cc}\rho \frac{\partial k}{\partial t}& =-\rho 𝐯\cdot \mathrm{\nabla }k-\rho 𝐯\cdot \mathrm{\nabla }h+\rho T𝐯\cdot \mathrm{\nabla }s+\mathrm{\nabla }\cdot (\sigma \cdot 𝐯)-{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}\\ \rho \frac{\partial \epsilon }{\partial t}& =\rho T\frac{\partial s}{\partial t}-\frac{p}{\rho }\mathrm{\nabla }\cdot (\rho 𝐯)\\ (k+\epsilon )\frac{\partial \rho }{\partial t}& =-(k+\epsilon )\mathrm{\nabla }\cdot (\rho 𝐯)\end{array}$$Y recogiendo términos, nos quedamos con:$$\frac{\partial }{\partial t}[\rho (k+\epsilon )]+\mathrm{\nabla }\cdot [\rho 𝐯(k+h)-\sigma \cdot 𝐯]=\rho T\frac{Ds}{Dt}-{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}$$Añadiendo ahora la divergencia del flujo de calor debido a la conducción térmica a cada lado, tenemos que:$$\frac{\partial }{\partial t}[\rho (k+\epsilon )]+\mathrm{\nabla }\cdot [\rho 𝐯(k+h)-\kappa \mathrm{\nabla }T-\sigma \cdot 𝐯]=\rho T\frac{Ds}{Dt}-\mathrm{\nabla }\cdot (\kappa \mathrm{\nabla }T)-{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}$$Sin embargo, sabemos que por la conservación de la energía en el lado izquierdo es igual a cero, lo que nos deja con:$$\rho T\frac{Ds}{Dt}=\mathrm{\nabla }\cdot (\kappa \mathrm{\nabla }T)+{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}$$El producto del tensor de esfuerzo viscoso y el gradiente de velocidad puede expandirse como:$$\begin{array}{cc}{\sigma }_{ij}\frac{\partial v_i}{\partial x_j}& =\mu (\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}-\frac{2}{3}{\delta }_{ij}\mathrm{\nabla }\cdot 𝐯)\frac{\partial v_i}{\partial x_j}+\zeta {\delta }_{ij}\frac{\partial v_i}{\partial x_j}\mathrm{\nabla }\cdot 𝐯\\ & =\frac{\mu }{2}{(\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}-\frac{2}{3}{\delta }_{ij}\mathrm{\nabla }\cdot 𝐯)}^2+\zeta (\mathrm{\nabla }\cdot 𝐯{)}^2\end{array}$$Así se llega a la forma final de la ecuación para la producción de entropía específica:$$\rho T\frac{Ds}{Dt}=\mathrm{\nabla }\cdot (\kappa \mathrm{\nabla }T)+\frac{\mu }{2}{(\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}-\frac{2}{3}{\delta }_{ij}\mathrm{\nabla }\cdot 𝐯)}^2+\zeta (\mathrm{\nabla }\cdot 𝐯{)}^2$$En el caso en que la conducción térmica y las fuerzas viscosas están ausentes, la ecuación para la producción de entropía se colapsa a $Ds/Dt=0$ - mostrando que el flujo de fluido ideal es isentrópico.

Aplicación

Esta ecuación se deduce en la sección 49, al comienzo del capítulo sobre "Conducción térmica en fluidos" del sexto volumen del Curso de Física Teórica (Course of Theoretical Physics) de L.D. Landau y E.M. Lifshitz. Puede utilizarse para medir la transferencia de calor y el flujo de aire en un frigorífico doméstico, hacer un análisis armónico de los regeneradores, o para comprender la física de los glaciares.

Véase también

• Disipación
• Ecuación del calor

Enlaces externos

• Plantilla:Cite book
• Plantilla:Cite journal