Enunciado

Processos em batelada são frequentemente usados em operações químicas e farmacêuticas para atingir uma desejada composição química do produto final. Os processos de transferência de calor presentes são tipicamente transientes, envolvendo um líquido com volume fixo que pode ser aquecido da temperatura ambiente até uma temperatura especificada pelo processo ou resfriado da temperatura do processo para a temperatura ambiente. Considere um processo em batelada no qual um fármaco (o fluido frio, f ) é despejado em um vaso isolado termicamente e fortemente agitado (um reator de mistura ideal) e aquecido pela passagem de um fluido quente ( q ) através de uma serpentina submersa, de parede delgada e área superficial A s . A vazão mássica, m ˙ q , a temperatura de entrada média, T q ,   e n t , e o calor específico, c p ,   q , do fluido quente são conhecidos, como o são a temperatura inicial, T f , i < T q ,   e n t , o volume, V f , a massa específica, ρ f , e o calor específico, c p ,   f , do fármaco. A transferência de calor do fluido quente para o fármaco é governada por um coeficiente global de transferência de calor U .

(a) Partindo de princípios básicos, deduza expressões que possam ser usadas para determinar as variações de T f e T q ,   s a i com o tempo durante o processo de aquecimento.

(b) Seja um volume de fármaco V f = 1   m 3 , com uma massa específica ρ f = 1000 k g m 3 , calor específico c p ,   f = 2000 J k g K e uma temperatura inicial de T f , i = 25 ° C . Uma serpentina com comprimento L = 40   m , diâmetro do tubo D = 50   m m e diâmetro da serpentina C = 500   m m está submersa no vaso. O fluido quente é alimentado na serpentina a T q ,   e n t = 200 ° C e auma vazão de m ˙ q = 2,4   k g / s . O coeficiente convectivo na superfície externa da serpentina pode ser aproximado por h e = 1000   W / m 2 K e as propriedades do fluido quente são c p ,   q = 2500   J / k g K , μ q = 0,002   N s / m 2 , k q = 0,260   W / m K e P r q = 20 . Para as condições anteriores, calcule e represente graficamente a temperatura do fármaco T f e a temperatura de saída do fluido quente T q ,   s a i como uma função do tempo para o intervalo 0 ≤ t ≤ 3600   s . Quanto tempo demora para se atingir uma temperatura da batelada de T f = 160 ° C ? O operador do processo pode controlar o tempo de aquecimento variando m ˙ q . Para 1 ≤ m ˙ q ≤ 5   k g / s , explore o efeito da vazão mássica no tempo t f necessário para atingir um valor de T f = 160 ° C .

Passo 1

Temos a seguinte situação: um fluido quente escoa através de uma serpentina que está imersa no fármaco. O fármaco, que tinha uma temperatura inicial T f , i é então aquecido ao longo do tempo, devido ao processo de transferência de calor.

No item (a) devemos encontrar relações que expressem a temperatura do fármaco ao longo do tempo T f ( t ) e também para a temperatura de saída do fluido quente da serpentina T q ,   s a i   ( t ) .

Vamos começar considerando que o processo ocorre sem perdas. Ou seja, todo o calor perdido pelo fluido quente é recebido pelo fármaco. Também vamos considerar que o processo de agitação é ideal e por isso a temperatura do fármaco é constante em todo o volume.

Além disso, também temos a informação de que a taxa de transferência de calor é controlada por um coeficiente global U . Então, podemos escrever:

q = U A s ∆ T m l

Abrindo o termo ∆ T m l teremos:

q = U A s ∆ T e n t - ∆ T s a i ln ⁡ ∆ T e n t ∆ T s a i

q = U A s T q ,   e n t - T f t - T q ,   s a i t - T f t ln ⁡ T q ,   e n t - T f t T q ,   s a i t - T f t

q = U A s T q ,   e n t - T q ,   s a i t ln ⁡ T q ,   e n t - T f t T q ,   s a i t - T f t

Passo 2

Como dissemos anteriormente, a taxa de transferência de calor também pode ser expressa como a energia perdida pelo fluido quente ao longo do escoamento:

q = m ˙ q c p q ∆ T q

q = m ˙ q c p q T q ,   e n t - T q ,   s a i t

Igualando esta expressão com a anterior:

m ˙ q c p q T q ,   e n t - T q ,   s a i t = U A s T q ,   e n t - T q ,   s a i t ln ⁡ T q ,   e n t - T f t T q ,   s a i t - T f t

m ˙ q c p q = U A s ln ⁡ T q ,   e n t - T f t T q ,   s a i t - T f t

ln ⁡ T q ,   e n t - T f t T q ,   s a i t - T f t = U A s m ˙ q c p q

T q ,   e n t - T f t T q ,   s a i t - T f t = exp ⁡ U A s m ˙ q c p q

T q ,   s a i t - T f t = T q ,   e n t - T f t exp ⁡ U A s m ˙ q c p q

T q ,   s a i t = T f t + T q ,   e n t - T f t exp ⁡ - U A s m ˙ q c p q

E temos uma expressão para a temperatura de saída do fluido quente. Vamos agora obter uma expressão para a temperatura do fármaco.

Passo 3

Assim como expressamos a taxa de transferência de calor em termos do calor perdido pelo fluido quente, também podemos escrevê-la como a energia recebida pelo fármaco. Então:

q = m f c p f d T f d t

q = ρ f V f c p f d T f d t

Igualando com a expressão em função do fluido quente:

ρ f V f c p f d T f d t = m ˙ q c p q T q ,   e n t - T q ,   s a i t

Substituindo a relação de T q ,   s a i t :

ρ f V f c p f d T f d t = m ˙ q c p q T q ,   e n t - T f t + T q ,   e n t - T f t exp ⁡ - U A s m ˙ q c p q

ρ f V f c p f d T f d t = m ˙ q c p q T q ,   e n t - T f t 1 - exp ⁡ - U A s m ˙ q c p q

1 T q ,   e n t - T f t d T f d t = m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q

d T f T q ,   e n t - T f t = m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q d t

Integrando dos dois lados:

∫ T f ,   i T f d T f T q ,   e n t - T f t = ∫ 0 t m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q d t

- ln ⁡ T q ,   e n t - T f + ln ⁡ T q ,   e n t - T f ,   i = m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

ln ⁡ T q ,   e n t - T f ,   i T q ,   e n t - T f = m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

T q ,   e n t - T f ,   i T q ,   e n t - T f = exp ⁡ m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

T q ,   e n t - T f = ( T q ,   e n t - T f ,   i ) exp ⁡ - m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

T f = T q ,   e n t - ( T q ,   e n t - T f ,   i ) exp ⁡ - m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

E chegamos em uma função da temperatura do fármaco que depende do tempo de processo.

Passo 4

Vamos agora para o item (b). Neste item recebemos alguns dados de propriedades do fluido quente e do fármaco, bem como das dimensões da serpentina. Devemos agora utilizar estas informações nas relações que encontramos no item (a) para graficar as temperaturas do fármaco e da saída de fluido quente, ao longo do intervalo 0 ≤ t ≤ 3600   s . Vamos lá?

Como primeiro passo antes de utilizar as nossas expressões ainda temos uma variável para obter: o coeficiente global. Como a serpentina tem parede delgada, podemos considerar que as áreas de troca térmica interna e externa são praticamente iguais. Neste caso, o coeficiente global se reduz à:

U = 1 1 h i + 1 h e

Sabemos que h e = 1000   W / m 2 K pelo enunciado mas ainda precisamos calcular o valor do coeficiente de convecção externo.

Para isso, começamos calculando o número de Reynolds:

R e D = 4 m ˙ q π D μ q

R e D = 4 × 2,4 π × 50 × 10 - 3 × 0,002

R e D = 30557,75

Como R e D > 2300 , o regime de escoamento é turbulento.

Para regime turbulento no interior de dutos cilíndricos:

N u D = h i D k q = 0,023 R e D 0,8 P r q 0,3

h i = k q D × 0,023 R e D 0,8 P r q 0,3

h i = 0,26 50 × 10 - 3 × 0,023 × 30557,75 0,8 × 20 0,3

h i = 1138   W / m 2 K

E o coeficiente global será:

U = 1 1 1138 + 1 1000

U = 532,3   W / m 2 K

Passo 5

Com o auxílio de Excel, ou outra ferramenta gráfica, iremos aplicar este valor nas expressões de temperatura do fármaco e temperatura de saída do fluido quente e, então fazer o gráfico destas funções para valores de tempo no intervalo 0 ≤ t ≤ 3600   s :

T f = T q ,   e n t - ( T q ,   e n t - T f ,   i ) exp ⁡ - m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

T f = 200 - ( 200 - 25 ) × exp ⁡ - 2,4 × 2500 1000 × 1 × 2000 × 1 - exp ⁡ - 532,3 × π × 50 × 10 - 3 × 40 2,4 × 2500 t

T f = 200 - ( 175 ) × exp ⁡ - 0,003 × 1 - exp ⁡ - 0,557 t

E

T q ,   s a i t = T f t + T q ,   e n t - T f t exp ⁡ - U A s m ˙ q c p q

T q ,   s a i t = T f t + 200 - T f t exp ⁡ - 532,3 × π × 50 × 10 - 3 × 40 2,4 × 2500

T q ,   s a i t = T f t + 200 - T f t exp ⁡ - 0,557

Aplicando estas funções para diversos valores de tempo no intervalo 0 ≤ t ≤ 3600   s devemos obter curvas parecidas com:

Ainda no item (b), devemos verificar em quanto tempo a temperatura do fármaco atinge T f = 160 ° C . Substituindo este valor na expressão:

T f = 200 - ( 175 ) × exp ⁡ - 0,003 × 1 - exp ⁡ - 0,557 t

160 = 200 - ( 175 ) × exp ⁡ - 0,003 × 1 - exp ⁡ - 0,557 t

- 40 - 175 = exp ⁡ - 0,003 × 1 - exp ⁡ - 0,557 t

0,229 = exp ⁡ - 0,003 × 1 - exp ⁡ - 0,557 t

ln ⁡ 0,229 = - 0,003 × 1 - exp ⁡ - 0,557 t

t = - ln ⁡ 0,229 0,003 × 1 - exp ⁡ - 0,557

t = 1152   s

Passo 6

Achou que tinha acabado? Então respira fundo porque ainda tem mais uma coisa a se fazer. Recebemos a informação de que é possível alterar a vazão mássica de fluido quente. Sabendo disso, devemos verificar como a vazão influencia no tempo necessário para atingir T f = 160 ° C . Vamos lá?

Se alteramos a vazão isto irá influenciar no número de Reynolds e, por consequência, no valor de h i e do coeficiente global. Para ver como isto influencia no tempo de processo, vamos fazer os cálculos para os dois extremos de vazão: m ˙ q = 1   k g / s e m ˙ q = 5   k g / s .

Começando com a vazão mínima, o número de Reynolds será:

R e D = 4 m ˙ q π D μ q

R e D = 4 × 1 π × 50 × 10 - 3 × 0,002

R e D = 12732,4

Como R e D > 2300 , o regime de escoamento é turbulento.

Para regime turbulento no interior de dutos cilíndricos:

N u D = h i D k q = 0,023 R e D 0,8 P r q 0,3

h i = k q D × 0,023 R e D 0,8 P r q 0,3

h i = 0,26 50 × 10 - 3 × 0,023 × 12732,4 0,8 × 20 0,3

h i = 564,9   W / m 2 K

E o coeficiente global será:

U = 1 1 564,9 + 1 1000

U = 361   W / m 2 K

Aplicando este valor na expressão de temperatura do fármaco:

T f = T q ,   e n t - ( T q ,   e n t - T f ,   i ) exp ⁡ - m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

T f = 200 - ( 200 - 25 ) × exp ⁡ - 1 × 2500 1000 × 1 × 2000 × 1 - exp ⁡ - 361 × π × 50 × 10 - 3 × 40 1 × 2500 t

160 = 200 - ( 175 ) × exp ⁡ - 0,00 125 × 1 - exp ⁡ - 0 , 90 7 t

- 40 - 175 = exp ⁡ - 0,00 125 × 1 - exp ⁡ - 0 , 90 7 t

0,229 = exp ⁡ - 0,00 125 × 1 - exp ⁡ - 0 , 90 7 t

ln ⁡ 0,229 = - 0,00 125 × 1 - exp ⁡ - 0 , 90 7 t

t = - ln ⁡ 0,229 0,0 0125 × 1 - exp ⁡ - 0,90 7

t = 1980   s

Passo 7

Fazendo agora para a vazão máxima, o número de Reynolds será:

R e D = 4 m ˙ q π D μ q

R e D = 4 × 5 π × 50 × 10 - 3 × 0,002

R e D = 63662

Como R e D > 2300 , o regime de escoamento é turbulento.

Para regime turbulento no interior de dutos cilíndricos:

N u D = h i D k q = 0,023 R e D 0,8 P r q 0,3

h i = k q D × 0,023 R e D 0,8 P r q 0,3

h i = 0,26 50 × 10 - 3 × 0,023 × 63662 0,8 × 20 0,3

h i = 2047,1   W / m 2 K

E o coeficiente global será:

U = 1 1 2047,1 + 1 1000

U = 671,8   W / m 2 K

Aplicando este valor na expressão de temperatura do fármaco:

T f = T q ,   e n t - ( T q ,   e n t - T f ,   i ) exp ⁡ - m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

T f = 200 - ( 200 - 25 ) × exp ⁡ - 5 × 2500 1000 × 1 × 2000 × 1 - exp ⁡ - 671,8 × π × 50 × 10 - 3 × 40 5 × 2500 t

160 = 200 - ( 175 ) × exp ⁡ - 0,00 6 25 × 1 - exp ⁡ - 0 , 338 t

- 40 - 175 = exp ⁡ - 0,00 6 25 × 1 - exp ⁡ - 0 , 338 t

0,229 = exp ⁡ - 0,00 6 25 × 1 - exp ⁡ - 0 , 338 t

ln ⁡ 0,229 = - 0,00 6 25 × 1 - exp ⁡ - 0 , 338 t

t = - ln ⁡ 0,229 0,0 06 25 × 1 - exp ⁡ - 0 , 338

t = 823   s

Então, como esperado, um aumento na vazão de fluido quente aumenta a turbulência do escoamento, gerando um coeficiente global de maior valor e uma troca térmica mais eficiente. Portanto, aumentando a vazão de fluido quente é possível diminuir o tempo de processo.

Resposta

(a) T f = T q ,   e n t - ( T q ,   e n t - T f ,   i ) exp ⁡ - m ˙ q c p q ρ f V f c p f 1 - exp ⁡ - U A s m ˙ q c p q t

T q ,   s a i t = T f t + T q ,   e n t - T f t exp ⁡ - U A s m ˙ q c p q .

(b) t f = 1152   s ; um aumento na vazão de fluido quente aumenta a turbulência do escoamento, gerando um coeficiente global de maior valor e uma troca térmica mais eficiente. Portanto, aumentando a vazão de fluido quente é possível diminuir o tempo de processo.

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Água c p = 4180 J k g K entra em um tubo de 12 cm de diâmetro e 8,5 m de comprimento, a uma temperatura de 75°C e uma taxa de 0,35 k g s , sendo resfriada por um refrigerante evaporando a -10°C. Se o
A água entra em um tubo de 5 mm de diâmetro e 13 m de comprimento a 45°C com velocidade de 0,3 m/s. O tubo é mantido a uma temperatura constante de 5°C. O comprimento necessário do tubo para que a águ
Ar a 110°C entra em um duto de 18 cm de diâmetro e 9 m de comprimento com velocidade de 4,5m/s. O duto é quase isotérmico a 85°C. A taxa de perda de calor do ar no duto é: 760W ; 890W ; 1210W ; 1370W
A água entra em um tubo de 5 mm de diâmetro e 13 m de comprimento a 45°C, com uma velocidade de 0,3 m/s. O tubo é mantido a uma temperatura constante de 8°C. A temperatura de saída da água é: 4 , 4 °
Qual a propriedade de um fluido é responsável pelo desenvolvimento da camada limite hidrodinâmica? Para que tipo de fluido não haverá nenhuma camada limite hidrodinâmica do tubo?