Lista Primeira Prova

Lista para a primeira prova de IT – 385.
UNIVERSIDADE FEDERAL RURAL DO RIO DE JANEIRO

INSTITUTO DE TECNOLOGIA
DEPARTAMENTO DE ENGENHARIA QUÍMICA
Métodos Computacionais Aplicado II
Prof. Luís Américo Calçada
Lista Geral Primeira Prova.
_________________________________
Questão 1) Considere a equação de transferência de calor,
𝜕2 𝑇 1 𝜕𝑇
= . , sujeita as condições de fronteira:
𝜕𝑥 2 𝛼 𝜕𝜏
T(t, 0) = 0 e T(t, L) = 0.
Resolva o problema de valor inicial se a temperatura inicial é dada por,

πx 2πx
𝑇(0, 𝑥 ) = 3 sin − sin
𝐿 𝐿
__________________________________
Questão 2) Considere a equação de transferência de calor. Determine o perfil
de temperatura variando com a posição e o tempo.
𝜕2 𝑇 1 𝜕𝑇
= . , sujeita as condições de fronteira/condições de contorno,
𝜕𝑥 2 𝛼 𝜕𝜏
∂u/∂x (t, 0) = 0, u(t, L) = 0 e u(0, x) = f(x) como condição inicial.
_________________________________
Questão 3) A Figura 1 representa a condução de calor em um sistema semi-
infinito.
1
A transferência de calor é dada pela equação,
𝜕2𝑇 1 𝜕𝑇
= .
𝜕𝑥 2 𝛼 𝜕𝜏
Condição inicial:
T(x,0) = 25℃
T(0,𝜏) = 100℃ para 𝜏 > 0
T(∞,𝜏) = 25℃ para 𝜏 > 0
Esse problema pode ser resolvido pela técnica da combinação de variáveis

sendo dada por
𝑇(𝑥, 𝜏) − 𝑇0 𝑥
= erf ( )
𝑇𝑖 − 𝑇0 2√𝛼𝜏
Pede-se para determinar, considerando o bloco de alumínio:

3.1 O perfil de temperatura no regime permanente e o fluxo de calor em
x=0,01 m.
3.2 O perfil de temperatura no regime transiente para os tempos de 10,
15 e 20 segundos.
3.3 Considerando regime transiente, o valor do fluxo de calor em função
do tempo para x=0,01 e em x=0,1 cm.
3.4 Um gráfico de posição em função do tempo do primeiro ponto em que
o fluxo de calor é nulo.
Dados o valor de  de 8,418. 10-5 m2/s, quando necessário use o Octave para
facilitar os cálculos.
2
Questão 4). Calcule os produtos internos PI e PII entre as funções e discuta os

resultados
1 +𝜋
𝑃𝐼 = ∫ cos 𝑘𝑥. cos 𝑙𝑥. 𝑑𝑥
𝜋 −𝜋
1 +𝜋
𝑃𝐼𝐼 = ∫ sin 𝑘𝑥. sin 𝑙𝑥. 𝑑𝑥
𝜋 −𝜋
________________________________________
𝐾𝑔 𝐾𝑔
Questão 5) Água com viscosidade 0,001 e densidade de 1000 é
𝑚 𝑠𝑒𝑔 𝑚3
alimentada entre duas placas a uma vazão Q sujeito a um diferencial de pressão

P1-P2 de 0,01 atm. O fluido escoa na direção axial x entre as duas placas
distanciadas em 2D=1,0 cm e com 0,5 m de comprimento. Para o sistema
descrito, desconsidere os efeitos de entrada e saída, para o escoamento de
fluido incompressível, escoamento laminar.
Considere somente o escoamento na direção x de fluido Newtoniano, observe

que o fluido segue a direção das setas segundo a representação. A simplificação
da equação do movimento para a componente x resultando em
Regime transiente
𝜕𝑣 𝜕𝑃 𝜕2 𝑣𝑥
𝜌 ( 𝜕𝑡𝑥 ) = − 𝜕𝑥 + 𝜇 .
𝜕𝑦 2
Regime Permanente
𝜕𝑃 𝜕2 𝑣𝑥
0 = − 𝜕𝑥 + 𝜇 .
𝜕𝑦 2
a) Para o regime permanente, demonstre que a equação pode ser

colocada na forma abaixo e obtenha o perfil de velocidade e calcule a
vazão Q.
3
∆𝑃 𝑑2 𝑣𝑥
− =𝜇 , sujeita as condições de contorno,
𝐿 𝑑𝑦 2
𝑑𝑣𝑥
= 0 𝑒𝑚 𝑦 = 0 𝑒 𝑣𝑥 = 0 𝑒𝑚 𝑦 =0,5 cm.
𝑑𝑦
Observe que P=P1-P2
b) Para o regime transiente, discretizando a derivada segunda no espaço

por diferenças finitas na direção y gera o sistema de EDO’s abaixo
descrito. Obtenha o perfil de velocidade em função do espaço e tempo.
Faça gráficos para apresentar os resultados e compare com o regime
permanente calculado analiticamente e o obtido pelo método numérico.
Em quanto tempo o sistema entre em regime permanente?
𝑑𝑣𝑥 (1) 1 ∆𝑃 𝜇 𝑣𝑥 (2)−𝑣𝑥 (1)

= − 𝜌 ∆𝐿 + 𝜌 ( ).
𝑑𝑡 ∆𝑦 2
𝑑𝑣𝑥 (2) 1 ∆𝑃 𝜇 𝑣𝑥 (3)−2𝑣𝑥 (2)+𝑣𝑥 (1)

= − 𝜌 ∆𝐿 + 𝜌 ( ).
𝑑𝑡 ∆𝑦 2
𝑑𝑣𝑥 (3) 1 ∆𝑃 𝜇 𝑣𝑥 (4)−2𝑣𝑥 (3)+𝑣𝑥 (2)

= − 𝜌 ∆𝐿 + 𝜌 ( ).
𝑑𝑡 ∆𝑦 2
𝑑𝑣𝑥 (4) 1 ∆𝑃 𝜇 𝑣𝑥 (5)−2𝑣𝑥 (4)+𝑣𝑥 (3)

= − 𝜌 ∆𝐿 + 𝜌 ( ).
𝑑𝑡 ∆𝑦 2
𝑑𝑣𝑥 (5) 1 ∆𝑃 𝜇 0−2𝑣𝑥 (5)+𝑣𝑥 (4)

= − 𝜌 ∆𝐿 + 𝜌 ( ).
𝑑𝑡 ∆𝑦 2
Observe que em t=0, vx=0 para todos os pontos em y e y=D/5.
Questão 6) As informações abaixo foram obtidas na página 531 do livro Gordon J. Van
Wylen & Richard E. Sonntag, segunda edição referente a propriedades do Nitrogênio
Superaquecido. De posse destes dados calcule:
a) O coeficiente de fugacidade de gases utilizando a equação abaixo para a pressão

variando de 0,1 a 4,0 Mpa e temperatura de 125 K.
4,0 𝑑𝑃
ln(∅) = ∫0,1 (𝑍 − 1) .
𝑃
b) O calor específico a pressão constante pela equação

𝜕ℎ
C𝑃 = ( )
𝜕𝑇 𝑃
4
𝐽
Dados 𝑅 = 8,314 𝑚𝑜𝑙.𝐾 , 𝑚𝑎𝑠𝑠𝑎 𝑚𝑜𝑙𝑒𝑐𝑢𝑙𝑎𝑟 𝑑𝑜 𝑁2 = 28 𝑔𝑚𝑜𝑙.
Tabela com os dados necessários para a solução do problema.
a) Teste se os dados satisfazem as relações de Maxwell. Escolha uma

região da tabela para fazer o teste.
b) Demonstre que os dados satisfazem a igualdade abaixo.
𝜕ℎ 𝜕𝑠
C𝑃 = ( ) = 𝑇 ( )
𝜕𝑇 𝑃 𝜕𝑇 𝑃
Questão 7) Vamos considerar uma partícula de catalisador de espessura 1

aonde ocorre uma reação de ordem 1 (Rice e Do, 2012). Vamos abordar o
problema unidimensional.
𝑑2𝑦
− 𝜙2 𝑦 = 0
𝑑𝑥 2
𝑑𝑦
C.C.1, 𝑥 = 0, =0
𝑑𝑥
C.C.2, 𝑥 = 1, 𝑦 = 1
A EDO tem solução analítica na forma,
5
cosh (𝜙𝑥)
𝑦=
cosh (𝜙)
onde  é um parâmetro relacionado ao sistema reacional e a transferência de
massa. Este problema está apresentado no Rice e Do (2012) no tópico de
resíduos ponderados. Responda ou calcule:
a) A derivada analítica dy/dx e numérica para =1 nos pontos x = 0,0, 0,25,
0,75 e 1,0.
b) Calcule o fator de efetividade analítico e numérico dado pela equação abaixo
1
𝐸𝑡𝑎 = ∫ 𝑦𝑑𝑥
0
Questão 8) Resolver de forma numérica por diferenças finitas no espaço e

aplicando a LSODE para integrar no tempo para o problema de transferência de
calor, Holman na sua décima edição. O sistema representa a condução de calor
em regime transiente em uma placa de largura 2L.
𝜕 2 𝜃 1 𝜕𝜃
=
𝜕𝑥 2 𝛼 𝜕𝜏
Condição inicial e condições de contorno:

a) ϴ = ϴi = Ti – T1, 𝜏=0, 0 ≤ x ≤ 2L
b) ϴ = 0 em x = 0, 𝜏 > 0
c) ϴ = 0 em x = 2L, 𝜏 > 0.
Dados: considerando o material alumínio, 2L = 10 cm, a temperatura

inicial de 25℃ e T1 = 500 ℃
Questão 9) Em Mecânica dos Fluidos a avaliação da turbulência é feita a partir de estudos

da velocidade instantânea em um dado tempo. Em termos de velocidade num dado ponto
em regime permanente, as figuras abaixo ilustram o escoamento laminar e turbulento
dentro de um tubo vz = vz(t). A figura abaixo ilustra como a velocidade varia no tempo para
os regimes de escoamento laminar e turbulento.
6
A velocidade no escoamento turbulento é definida por,
_
v z = v z + vz
,
_
onde Vz é a velocidade em um dado instante, v z velocidade média temporal
vz
e flutuação da velocidade.
A velocidade média temporal é calculada por
𝜏 𝜏
∫ 𝑉𝑧 𝑑𝑡 ∫ 𝑉𝑧 𝑑𝑡
𝑉̅𝑧 = 0 𝜏 = 0𝜏 ,
∫0 𝑑𝑡
onde Vz é a velocidade média em dado tempo.
a) Para os dados da tabela abaixo de velocidade obtida em períodos de 10

segundos com intervalo de 1s, construa um programa geral que calcule
as derivadas de velocidade no tempo e a velocidade média temporal.
t(seg) 0 1 2 3 4 5 6 7 8 9 10
Vz(m/s) 0 4,5 3,9 4,6 3,0 3,3 4,9 4,0 4,1 3,4 4,6
Questão 10) Considere o problema a seguir e dado em aula do modelo presa

predador proposto por Lotka-Volterra.
𝑑𝑥
=𝐴∗𝑥−𝐵∗𝑥∗𝑦 𝑡 = 0, 𝑥 = 4
𝑑𝑡
7
𝑑𝑦
= −𝐶 ∗ 𝑦 + 0,5 ∗ 𝑥 ∗ 𝑦, 𝑡 = 0, 𝑦 = 3
𝑑𝑡
onde, y é o número de indivíduos de algum predador (exemplo: lobo), x é o
número da indivíduos da sua presa (exemplo: coelho) e t representa o tempo.
Os parâmetros A=1.00; B=.5; C=0.75; D =0.5.
a) Integre o sistema de EDO usando Euler Simples avançando 0,1 segundos,

b) Escreva um programa em Octave para solução usando LSODE que construa
faça o gráfico.
Questão 10 ) Com relação ao polinômio de Jacobi e quadratura de Gauss.

a) Construa os polinômios de Jacobi de primeira, segunda e terceira ordem
para alfa e beta iguais a zero.
b) Determine as raízes dos polinômios.
c) Determine os pesos da quadratura de Gauss para os casos abaixo e
compare com os tabelados.
c.1) um ponto interno
c.2) dois pontos internos
c.3) três pontos internos.

Lista Primeira Prova

Enviado por

Direitos autorais:

Formatos disponíveis

Lista Primeira Prova

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Lista Primeira Prova

Enviado por

Direitos autorais:

Formatos disponíveis

Lista para a primeira prova de IT – 385.

UNIVERSIDADE FEDERAL RURAL DO RIO DE JANEIRO

Resolva o problema de valor inicial se a temperatura inicial é dada por,

A transferência de calor é dada pela equação,

Esse problema pode ser resolvido pela técnica da combinação de variáveis

Pede-se para determinar, considerando o bloco de alumínio:

Questão 4). Calcule os produtos internos PI e PII entre as funções e discuta os

alimentada entre duas placas a uma vazão Q sujeito a um diferencial de pressão

Considere somente o escoamento na direção x de fluido Newtoniano, observe

a) Para o regime permanente, demonstre que a equação pode ser

Observe que P=P1-P2

b) Para o regime transiente, discretizando a derivada segunda no espaço

𝑑𝑣𝑥 (1) 1 ∆𝑃 𝜇 𝑣𝑥 (2)−𝑣𝑥 (1)

𝑑𝑣𝑥 (2) 1 ∆𝑃 𝜇 𝑣𝑥 (3)−2𝑣𝑥 (2)+𝑣𝑥 (1)

𝑑𝑣𝑥 (3) 1 ∆𝑃 𝜇 𝑣𝑥 (4)−2𝑣𝑥 (3)+𝑣𝑥 (2)

𝑑𝑣𝑥 (4) 1 ∆𝑃 𝜇 𝑣𝑥 (5)−2𝑣𝑥 (4)+𝑣𝑥 (3)

𝑑𝑣𝑥 (5) 1 ∆𝑃 𝜇 0−2𝑣𝑥 (5)+𝑣𝑥 (4)

Observe que em t=0, vx=0 para todos os pontos em y e y=D/5.

a) O coeficiente de fugacidade de gases utilizando a equação abaixo para a pressão

b) O calor específico a pressão constante pela equação

Tabela com os dados necessários para a solução do problema.

a) Teste se os dados satisfazem as relações de Maxwell. Escolha uma

b) Demonstre que os dados satisfazem a igualdade abaixo.

Questão 7) Vamos considerar uma partícula de catalisador de espessura 1

A EDO tem solução analítica na forma,

Questão 8) Resolver de forma numérica por diferenças finitas no espaço e

Condição inicial e condições de contorno:

Dados: considerando o material alumínio, 2L = 10 cm, a temperatura

Questão 9) Em Mecânica dos Fluidos a avaliação da turbulência é feita a partir de estudos

A velocidade no escoamento turbulento é definida por,

onde Vz é a velocidade média em dado tempo.

a) Para os dados da tabela abaixo de velocidade obtida em períodos de 10

Questão 10) Considere o problema a seguir e dado em aula do modelo presa

a) Integre o sistema de EDO usando Euler Simples avançando 0,1 segundos,

Questão 10 ) Com relação ao polinômio de Jacobi e quadratura de Gauss.

Você também pode gostar