Termocupla expuesta a un pulso rectangular de temperatura

Por: publicado el Oct 14, 2022

Problema 5.5 (Process Dynamics and Control - Seborg, Edgar, M, D)

diagrama p5.5

Una termocupla tiene las siguientes características cuando es sumergido en un tanque agitado:

Masa termocupla = 1 gCapacidad calorıˊfica = 0.25 cal/(cm2 °C)Coeficiente de transferencia de calor = 20 cal/(cm2 h °C)Aˊrea de superficie de la termocupla = 3 cm2\begin{array}{rl} Masa\ termocupla\ & =\ 1\ g\\ Capacidad\ calorífica\ & =\ 0.25\ cal/(cm^2\ °C)\\ Coeficiente\ de\ transferencia\ de\ calor\ & =\ 20\ cal/(cm^2\ h\ °C)\\ Área\ de\ superficie\ de\ la\ termocupla\ & =\ 3\ cm^2\\ \end{array}

Halle la función transferencia de la termocupla que relaciona la temperatura de la termocupla TT y la de la de los alrededores TAT_A.

Si la termocupla esta inicialmente fuera tanque, con una temperatura ambiente de 23 °C. ¿Cuál será la máxima temperatura alcanzada si es puesta en el tanque a temperatura de 80 °C, y se la saca después de 20s?

Grafique las la variación de temperatura de la termocupla.

Hallando la función transferencia

Escribiendo el balance de energía en estado transitorio y en estado estacionario.

mCpdTdt=UA(TA−T)mC_p\frac{dT}{dt}=UA(T_A-T)

0=UA(TAs−Ts)0=UA(T_{As}-T_s)

Restando ambas ecuaciones y pasando a variables desviación

mCpd(T−Ts)dt=UA[(TA−TAs)−(T−Ts)]m C_p \frac{d(T-T_s)}{dt}=UA\big[(T_A-T_{As})-(T-T_s)\big]

mCpdT′dt=UA(TA′−T′)mC_p\frac{dT'}{dt}=UA(T'_A-T')

Aplicando la transformada de Laplace, sabiendo que para el estado estacionario T’(t=0)=0

mCps T(s)=UA(TA′(s)−T′(s))mC_ps\space T(s)=UA(T'_A(s)-T'(s))

Despejando la función tranferencia y haciendo τ=(mCp)/(UA)\tau = (mC_p)/(UA) en segundos

T′(s)TA′(s)=1τs+1\frac{T'(s)}{T'_A(s)}=\frac{1}{\tau s+1}

τ=1 g ×0.25 calg °C20calcm2 h °C×1 h360 s3 cm2\tau = \frac{1\ g\ \times 0.25\ \frac{cal}{g\ °C}}{20\frac{cal}{cm^2\ h\ °C}\times \frac{1\ h}{360\ s}3\ cm^2}

τ=15\tau = 15

T′(s)TA′(s)=115s+1    ... (1)\mathbf{\frac{T'(s)}{T'_A(s)}=\frac{1}{15 s+1}}\space\space\space\space\textbf{... (1)}

Describiendo la perturbación

TA′(t)=TA(t)−TAs{0si t<080−23 °Csi 0<t<200si t>20T'_A(t)= T_A(t)-T_{As} \begin{cases} 0 &\text{si } t < 0 \\ 80-23\ °C &\text{si } 0< t < 20\\ 0 &\text{si } t>20\\ \end{cases}

TA′(t)=57u(t)−57u(t−20)T'_A(t)=57u(t)-57u(t-20)

Aplicando la antitransformada

TA′(s)=57s−57e−20ssT'_A(s)= \frac{57}{s}-57\frac{e^{-20s}}{s}

Reemplanzando en la ecuación (1)

T′(s)=57s(15s+1)−57e−20ss(15s+1)T'(s)=\frac{57}{s(15 s+1)}-\frac{57e^{-20s}}{s(15 s+1)}

Realizando la antitransformada (de tablas). Recuerde que e−ase^{-as} crea un desfase en la antitransformada igual a t=t−at= t-a

T´(t)=57(1−e−t/15)u(t)−57(1−e−t/15)u(t)∣t=t−20T´(t)=57(1-e^{-t/15})u(t)-57(1-e^{-t/15})u(t)|_{t=t-20}

T´(t)=57(1−e−t/15)u(t)−57(1−e−(t−20)/15)u(t−20)T´(t)=57(1-e^{-t/15})u(t)-57(1-e^{-(t-20)/15})u(t-20)

T(t)=T′(t)+TsT(t)=T'(t)+T_s

T(t)=57(1−e−t/15)u(t)−57(1−e−(t−20)/15)u(t−20)+23T(t)=57(1-e^{-t/15})u(t)-57(1-e^{-(t-20)/15})u(t-20)+23

Expresando esta función de manera más entendible

T(t)={57(1−e−t/15)+23si 0<t<2057(−e−t/15+e−(t−20)/15)+23si t>20T(t)= \begin{cases} 57(1-e^{-t/15})+23 &\text{si } 0<t< 20 \\ 57(-e^{-t/15}+e^{-(t-20)/15})+23 &\text{si } t > 20\\ \end{cases}

Las funciones son monotonicas en ambos rangos, por lo que para el máximo bastará con analizar el punto anguloso, la temperatura máxima se da un instante antes de sacar la termocupla del tanque a 80 °C, es decir a los 20 s.

T(t=20)=57(1−e−20/15)+23T(t=20)=57(1-e^{-20/15})+23

Tmax=64.974 °C\mathbf{T_{max}=64.974\ °C}

Gráfica de la función

Respuesta temperatura

Referencias

  • Seborg, D. E.; Edgar, T. F.; Mellichamp, D. A.; Doyle, F. J. (2016). Process Dynamics and Control (4th edition). John Wiley & Sons, Inc. ISBN 978-1-119-28591-5.