CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO
description
Transcript of CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO
![Page 1: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/1.jpg)
CAPÍTULO 3CAPÍTULO 3
ESTRATÉGIAS DE CÁLCULOESTRATÉGIAS DE CÁLCULO
11 de agosto de 2014
![Page 2: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/2.jpg)
REVISÃO DOS CAPÍTULOS ANTERIORES
![Page 3: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/3.jpg)
1.1. INTRODUÇÃO GERALINTRODUÇÃO GERAL
![Page 4: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/4.jpg)
É o campo da Engenharia que congrega os conceitos e os métodos destinados à
concepção, ao projeto e à operação de processos químicos
ENGENHARIA DE PROCESSOS
Implícitos na sua prática encontram-se o emprego intensivo de recursos computacionais e o atendimento a requisitos de
natureza econômica, material, energética, de preservação ambiental e de segurança.
em que se encontram integrados equipamentos de reação, separação, integração material e
energética e controle.
![Page 5: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/5.jpg)
H2O
2 NH4Cl
1CaCO3
CaO
Ca(OH)2
CaCl2
2 H2O
2 NH3
2 NH4OH2 NaCl
2 NaHCO3 Na2CO3
CO2
CO2
2 CO2
2
3
4
5
6
PROCESSO SOLVAY
Etapa 1: Calcinação do calcáreo CaCO3 (matéria prima na reação originalmente desejada) que se decompõe na cal CaO para a Etapa 2 e no CO2 para a Etapa 5.
CaCO3 + 2 NaCl Na2CO3 + CaCl2
![Page 6: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/6.jpg)
A Engenharia de Processos surgiu com a “Fertilização” da Eng. Química tradicional com elementos de:
Resultando:
Utilização mais organizada e mais eficiente dos conhecimento específicos da Engenharia Química no Projeto de Processos:
- Projeto mais rápido e mais eficiente.
- Processos mais econômicos, seguros e limpos.
CIÊNCIAS BÁSICAS
FUNDAMENTOS
ENG. DE EQUIPAMENTOS
ENG. DE PROCESSOS
Engenharia de Sistemas:No tratamento de conjuntos complexos de elementos interdependentes
Inteligência Artificial:Na resolução de problemas combinatórios
![Page 7: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/7.jpg)
Teoria de Projeto
Eng. Naval
Eng. Elétrica
Eng.Química
Eng. Mecânica
Conhecimento específico
de cada área utilizado intuitivamente
Aplicável a todas as áreas
As engenharias
experimentaram um ganho expressivo
Engenharia de Processos
![Page 8: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/8.jpg)
O Projeto resulta de um conjunto de ações desenvolvidas
DesdeA decisão de se produzir um determinado produto
AtéÀ conclusão do Projeto
As ações são numerosas e diversificadas !!!
PROJETO DE PROCESSOS
É o conjunto de documentos elaborados por uma equipe de engenheiros com detalhes suficientes para a construção e a
operação de uma planta industrial
![Page 9: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/9.jpg)
Investigar mercado para o produto
Investigar disponibilidade de matéria prima
Estabelecer as condições da reação e sub-produtos
Estabelecer o número e o tipo dos reatores
Definir o número e o tipo dos separadores
Definir o número e o tipo de trocadores de calor
Estabelecer malhas de controle
Definir o fluxogramado processo
Calcular as dimensõesdos equipamentos
Calcular o consumo de matéria prima
Calcular o consumo de utilidades
Calcular o consumo de insumos
Calcular a vazão dascorrentes intermediárias
Investigar reagentesplausíveis Avaliar a
lucratividadedo processo
![Page 10: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/10.jpg)
À luz da Engenharia de Processos elas são organizadas da seguinte forma
quanto à sequência no Projeto
![Page 11: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/11.jpg)
SELEÇÃO DAROTA QUÍMICA
Investigar mercado para o produto
Investigar reagentesplausíveis
Investigar a disponibilidade
das matérias primas
Definir as condições das reações e
identificar os sub-produtos gerados
SÍNTESE
Estabelecer o número e o tipo
dos reatores
Definir o número e o tipo dos
separadores
Definir o número e o tipo de
trocadores de calor
Estabelecer malhas de controle
Definir o fluxograma do
processo
ANÁLISE
Calcular o consumo de utilidades
Calcular a vazão das correntes
intermediárias
Calcular as dimensões dos equipamentos
Calcular o consumo dos insumos
Calcular o consumo de matéria prima
Avaliar a lucratividadedo processo
![Page 12: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/12.jpg)
INTRODUÇÃO GERAL
1
INTRODUÇÃO À
SÍNTESE DE PROCESSOS
8
6
SÍNTESE DESISTEMAS DE SEPARAÇÃO
7
SÍNTESE
SÍNTESE DE
SISTEMAS DE
INTEGRAÇÃO ENERGÉTICA
INTRODUÇÃO À
ANÁLISE DE PROCESSOS
2
ESTRATÉGIAS
DE CÁLCULO
3
OTIMIZAÇÃOAVALIAÇÃO
ECONÔMICA
4 5
ANÁLISE
![Page 13: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/13.jpg)
Projeto: primeiro passo
DIFICULDADE:
MULTIPLICIDADE DE SOLUÇÕES
Seleção da Rota Química
![Page 14: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/14.jpg)
Rotas para a produção de fenol
![Page 15: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/15.jpg)
DIFICULDADE NA SÍNTESE
MULTIPLICIDADE DE SOLUÇÕES
Síntese de um Fluxograma para a Rota Química
Projeto: segundo passo
![Page 16: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/16.jpg)
Equipamentos disponíveis para a geração de um fluxograma
RM
Reator demistura
RT
Reator tubular
DS
Coluna de destilaçãosimples
DE
Coluna de destilaçãoextrativa
A
Aquecedor
R
Resfriador
T
Trocador deIntegração
Este problema é simples
O espaço das soluções é constituido apenas de 8 fluxogramas
EXEMPLO
![Page 17: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/17.jpg)
DS
RM
R
A
A,B
P,A
P
A
(7)
RM
A,B
P,A
DS
P
A
T
(8)
RM
R
A
A,B
P,A
P
A
DE
(9)
DSRT RAA,B A,P
P
A
(11)
RM
A,B
P,A
P
A
T DE
(10)
DSRT A,P
P
A
T
A,B
(12)
RT RAA,B A,P
P
A
DE
(13)
RT A,P
P
A
T
A,B
DE
(14)
![Page 18: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/18.jpg)
UM RISCO INERENTE À SÍNTESE . . .
![Page 19: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/19.jpg)
EXPLOSÃO COMBINATÓRIA !!!
Muitas soluções para analisar
![Page 20: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/20.jpg)
DIFICULDADE NA ANÁLISE
MULTIPLICIDADE DE SOLUÇÕES
Análise dos fluxogramas gerados na Síntese
Projeto: terceiro passo
![Page 21: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/21.jpg)
1 2
Q = 10.000 kgA/h
x = 0,02 kgAB/kgAo
W1
kgB/hW2
kgB/h
y1
kgAB/kgBy2
kgAB/kgB
x1
x2
kgAB/kgAkgAB/kgA
Modelo Matemático1. Q(xo - x1) - W1 y1 = 02. y1 - k x1 = 03. Q(x1 -x2) - W2 y2 = 04. y2 - k x2 = 0
Avaliação EconômicaL = R - CR = pAB (W1 y1 + W2 y2 )C = pB (W1 + W2)pAB = 0,4 $/kgAB : pB = 0,01 $/kgB
Para cada par de valores x1,x2 resultam valores de W1, W2, y1, y2 e Lucro
EXEMPLO
dimensionamento de 2 extratores em série
![Page 22: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/22.jpg)
0,005
0,010
0,015
0,020
0,025
0,030
0,035
02468
101214161820
0,002
0,004
0,006
0,008
0,010
0,012
0,014
0,016
0,018
0,020
MULTIPLICIDADE NA ANÁLISE
Dificuldade: infinidade de soluções viáveis
A cada par (x1,x2) corresponde uma solução viável
![Page 23: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/23.jpg)
Todo problema com Multiplicidade de Soluções
exige a busca da sua
OTIMIZAÇÃO
Solução Ótima
através de
![Page 24: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/24.jpg)
O Projeto de Processos é um problema complexo de otimização.
Constata-se, assim, que ...
![Page 25: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/25.jpg)
Primeiro fator de complexidade
multiplicidade de soluções nos três níveis
Nível Tecnológico: determinar a rota química ótima.
Nível Paramétrico (Análise): determinar as dimensões ótimas de equipamentos e correntes.
Nível Estrutural (Síntese): determinar a estrutura ótima.
![Page 26: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/26.jpg)
Segundo fator de complexidade
Os 3 problemas são interdependentes.
A solução ótima de um está atrelada à solução ótima dos outros dois.
![Page 27: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/27.jpg)
Super-estrutura
Uma primeira abordagem para o Projeto
![Page 28: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/28.jpg)
SUPER - ESTRUTURA
Uma estrutura que abriga qualquer uma das estruturas alternativas para um sistema.
ExemploSuper-estrutura para algarismos
![Page 29: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/29.jpg)
Equipamentos disponíveis para a geração de um fluxograma
RM
Reator demistura
RT
Reator tubular
DS
Coluna de destilaçãosimples
DE
Coluna de destilaçãoextrativa
A
Aquecedor
R
Resfriador
T
Trocador deIntegração
Este problema é simples
O espaço das soluções é constituido apenas de 8 fluxogramas
EXEMPLO
![Page 30: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/30.jpg)
DS
RM
R
A
A,B
P,A
P
A
(7)
RM
A,B
P,A
DS
P
A
T
(8)
RM
R
A
A,B
P,A
P
A
DE
(9)
DSRT RAA,B A,P
P
A
(11)
RM
A,B
P,A
P
A
T DE
(10)
DSRT A,P
P
A
T
A,B
(12)
RT RAA,B A,P
P
A
DE
(13)
RT A,P
P
A
T
A,B
DE
(14)
![Page 31: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/31.jpg)
DE
DS
RT
RM
T
R
A
Superestrutura para o exemplo ilustrativo
Contém todos os equipamentos e todas as conexões lógicas.
Abriga todos os fluxogramas possíveis do exemplo.
![Page 32: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/32.jpg)
DE
DS
RT
RM
T
R
A
Monta-se o modelo matemático da superestrutura
Recai-se em problema complexo de otimização
Solução
O resultado é o fluxograma ótimo
![Page 33: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/33.jpg)
Super-estrutura do Problema evidenciando o Fluxograma 7
DE
DS
RT
RM
T
R
A
DS
RM
R
A
A,B
P,A
P
A
(7)
![Page 34: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/34.jpg)
Busca Orientada por Árvore de Estados
Uma abordagem alternativa
![Page 35: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/35.jpg)
Nível TecnológicoSeleção de uma Rota
Fluxograma ?Dimensões ?
Nível EstruturalSíntese de um
FluxogramaDimensões ? Lucro?
Nível ParamétricoAnálise do Fluxograma
Dimensionamentodos Equipamentos
e das Correntes. Lucro.
Solução Ótima: Reagentes = D,E; Fluxograma = 3; x = 4 demais dimensões.
RaizRota Química ?Fluxograma ?Dimensões ?
Busca Orientada por Árvore de Estados
P?? ?
D+E P+FD,E P,F
??
A+B P+CA,B P,C
??
1 PAB Cx
?
T D
2PA
B Cx
?T A
P3DE Fx
?
DM
PF
4DE x
?
M E
L
x
6
x o = 3x*
8
L
xx o = 4x*
L
10
xx o = 6x*
L
x
7
x o = 5x*
![Page 36: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/36.jpg)
P?? ?
D+E P+FD,E P,F
??
L
x4
10
?
P3DE Fx
Nível TecnológicoSeleção de uma Rota
Fluxograma ?Dimensões ?
Nível EstruturalSíntese de um
FluxogramaDimensões ? Lucro?
Nível ParamétricoAnálise do Fluxograma
Dimensionamentodos Equipamentos
e das Correntes. Lucro.Solução Ótima: Reagentes = D,E; Fluxograma = 3; x = 4 demais dimensões.
RaizRota Química ?Fluxograma ?Dimensões ?
Solução do Problema de Projeto por Busca Orientada
Vantagem
Varre todas as soluções sem
repetiçõessem omitir a ótima
Desvantagem
Explosão Combinatória
(outros métodos)
![Page 37: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/37.jpg)
2. INTRODUÇÃO À ANÁLISE DE PROCESSOS2. INTRODUÇÃO À ANÁLISE DE PROCESSOS
![Page 38: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/38.jpg)
Nível TecnológicoSeleção de uma Rota
Fluxograma ?Dimensões ?
Nível EstruturalSíntese de um
FluxogramaDimensões ? Lucro?
Nível ParamétricoAnálise do Fluxograma
Dimensionamentodos Equipamentos
e das Correntes. Lucro.
Solução Ótima: Reagentes = D,E; Fluxograma = 3; x = 4 demais dimensões.
RaizRota Química ?Fluxograma ?Dimensões ?
Busca Orientada por Árvore de Estados
P?? ?
D+E P+FD,E P,F
??
A+B P+CA,B P,C
??
1 PAB Cx
?
T D
2PA
B Cx
?T A
P3DE Fx
?
DM
PF
4DE x
?
M E
L
x
6
x o = 3x*
8
L
xx o = 4x*
L
10
xx o = 6x*
L
x
7
x o = 5x*
![Page 39: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/39.jpg)
INTRODUÇÃO GERAL
1
INTRODUÇÃO À
SÍNTESE DE PROCESSOS
8
6
SÍNTESE DESISTEMAS DE SEPARAÇÃO
7
SÍNTESE
SÍNTESE DE
SISTEMAS DE
INTEGRAÇÃO ENERGÉTICA
INTRODUÇÃO À
ANÁLISE DE PROCESSOS
2
ESTRATÉGIAS
DE CÁLCULO
3
OTIMIZAÇÃOAVALIAÇÃO
ECONÔMICA
4 5
ANÁLISE
![Page 40: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/40.jpg)
“Bola de Cristal”
OBJETIVO DA ANÁLISE
Prever e avaliar
o desempenho físico e econômico
ou ainda inexistente (em fase de projeto)
de um processo já existente (em operação)
![Page 41: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/41.jpg)
Consiste em
(a) prever as dimensões dos principais equipamentos e as condições das correntes, necessárias para atender às especificações técnicas estabelecidas para o projeto.
BaseModelo Matemático
Prever e avaliar o desempenho FÍSICO
(b) prever o comportamento do processo em condições diferentes daquelas para qual foi dimensionado.
![Page 42: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/42.jpg)
Consiste em Verificar se o processo atende aos critérios econômicos de lucratividade de forma a justificar a sua montagem e a sua operação.
BaseCritério Econômico
Prever e avaliar o desempenho ECONÔMICO
![Page 43: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/43.jpg)
ESTRATÉGIAS
DE CÁLCULO
3
AVALIAÇÃO
ECONÔMICA
4
INTRODUÇÃO À
ANÁLISE DE PROCESSOS
2
OTIMIZAÇÃO
5
Resumo da Análise de ProcessosCada Capítulo gera subsídios para os Módulos Computacionais
MODELOFÍSICO
MODELOECONÔMICO OTIMIZAÇÃO
Variáveis Especificadas
Variáveis de Projeto
Parâmetros Econômicos
ParâmetrosFísicos Dimensões Calculadas Lucro
![Page 44: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/44.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 45: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/45.jpg)
ROTEIRO PARA A RESOLUÇÃO DE PROBLEMASNA ANÁLISE DE PROCESSOS
2. Escrever o modelo matemático.
1. Reconhecer ou desenhar o fluxograma: equipamentos, correntes, variáveis do processo.
7. Avaliar criticamente o resultado.
6. Resolver o problema.
5. Estabelecer uma estratégia de cálculo.
4. Efetuar o Balanço de Informação.
3. Identificar as variáveis conhecidas e as metas de projeto.
fundamental
mais importante
![Page 46: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/46.jpg)
FINAL DA REVISÃO
![Page 47: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/47.jpg)
MOTIVAÇÃO PARA O CAPÍTULO 3
![Page 48: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/48.jpg)
Forma Geral dos Modelos Matemáticos de Processos
f1(x1, x2, ..., xi ,..., xM) = 0f2(x1, x2, ..., xi ,..., xM) = 0 . . . . . .fN(x1, x2, ..., xi,..., xM) = 0
N equaçõesM incógnitas
constituído do conjunto dos modelos dos equipamentos.
Sistema de equações algébricas
![Page 49: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/49.jpg)
Exemplo: Modelo do Processo Ilustrativo
![Page 50: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/50.jpg)
Partindo dos modelos dos equipamentos
01. f11 - f12 - f13 = 002. W15 - f23 = 003. f31 - f32 = 004. k – (3 + 0,04 Td) = 005. k – x13 / x12 = 006. (f11 Cp1 + f31 Cp3) (T1 - Td) + W15 Cp2l (T15 - Td) = 007. Vd - (f11 /1 + W15/2 + f31/3) = 008. r - f13/f11 = 009. T2 – Td = 010. T3 – Td = 0
11. f13 - f14 = 012. f23 - f24 - W5 = 013. W6 - W7 = 014. W6 [3 + Cpv (T6 – T7)] - Qe = 015. Qe – [(f13Cp1 + f23Cp2l)(Te - T3) + W5 2] = 016. Qe - Ue Ae e = 017. e - (T6- Te) = 018. T4 – Te = 019. T5 – Te = 0
extrator
evaporador
![Page 51: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/51.jpg)
20. W8 - W9 = 021. W5 - W10 = 022. Qc - W8 Cp3 (T9 - T8) = 023. W5 [2 + Cp2g (T5 – T10)] - Qc = 024. Qc - Uc Ac c = 025. c - [(T5 - T9) - (T10 - T8)]/ln[(T5 - T9)/(T10 - T8)] = 0
26. W11 - W12 = 027. W10 - W13 = 028. Qr - W11 Cp3 (T12 - T11) = 029. Qr - W10 Cp2l (T10 - T13) = 030. Qr - Ur Ar r = 031. r - [(T10 - T12) - (T13 - T11)]/ln[(T10 - T12)/(T13 - T11)] = 0
32. W13 + W14 - W15 = 033. W13 (T15 - T13) + W14 (T15 - T14) = 0
34. f11 + f31 - W1 = 035. x11 - f11 / W1 = 036. f12 + f22 – W2 = 037. x12 - f12/ W2 = 038. f13 + f23 – W3 = 039. x13 - f13 / W3 = 040. f14 + f24 - W4 = 041. x14 - f14/ W4 = 0
condensador
resfriador
misturador
correntes multicomponentes
![Page 52: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/52.jpg)
Modelo Completo
01. f11 - f12 - f13 = 002. W15 - f23 = 003. f31 - f32 = 004. k – (3 + 0,04 Td) = 005. k – x13 / x12 = 006. (f11 Cp1 + f31 Cp3) (T1 - Td) + W15 Cp2l (T15 - Td) = 007. Vd - (f11 /1 + W15/2 + f31/3) = 008. r - f13/f11 = 009. T2 – Td = 010. T3 – Td = 011. f13 - f14 = 012. f23 - f24 - W5 = 013. W6 - W7 = 014. W6 [3 + Cpv (T6 – T7)] - Qe = 015. Qe – [(f13Cp1 + f23Cp2l)(Te - T3) + W5 2] = 016. Qe - Ue Ae e = 017. e - (T6- Te) = 018. T4 – Te = 019. T5 – Te = 0
20. W8 - W9 = 021. W5 - W10 = 022. Qc - W8 Cp3 (T9 - T8) = 023. W5 [2 + Cp2g (T5 – T10)] - Qc = 024. Qc - Uc Ac c = 025. c - [(T5 - T9) - (T10 - T8)]/ln[(T5 - T9)/(T10 - T8)] = 026. W11 - W12 = 027. W10 - W13 = 028. Qr - W11 Cp3 (T12 - T11) = 029. Qr - W10 Cp2l (T10 - T13) = 030. Qr - Ur Ar r = 031. r - [(T10 - T12) - (T13 - T11)]/ln[(T10 - T12)/(T13 - T11)] = 032. W13 + W14 - W15 = 033. W13 (T15 - T13) + W14 (T15 - T14) = 0
34. f11 + f31 - W1 = 035. x11 - f11 / W1 = 036. f12 + f22 – W2 = 037. x12 - f12/ W2 = 038. f13 + f23 – W3 = 039. x13 - f13 / W3 = 040. f14 + f24 - W4 = 041. x14 - f14/ W4 = 0
![Page 53: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/53.jpg)
Consiste em representar o processo matematicamente utilizando os conhecimentos relativos aos Fundamentos e Equipamentos
Consiste em utilizar técnicas de processamento de informação na resolução de problemas.
Competem ao Engenheiro Químico
(a) Formulação (Modelagem Matemática):
(b) Resolução :
É um pré-requisito para esta Disciplina
Palavras-chave : Formulação e Resolução !!!
Formulação e Resolução do Modelo CIÊNCIAS BÁSICAS
FUNDAMENTOS
ENG. DE EQUIPAMENTOS
Objeto deste Capítulo
![Page 54: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/54.jpg)
A complexidade dos modelos exige o estabelecimento prévio de uma
Estratégia de CálculoTema deste Capítulo
Fontes de complexidade:
Em geral, os modelos de processos são complexos.
(c) presença de reciclos
(b) não-linearidade de equações
(a) grande número de equações e de variáveis
Desafio: como viabilizar a resolução de modelos tão complexos, e como faze-lo da forma mais eficiente possível ???
COMPLEXIDADE DOS MODELOS
![Page 55: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/55.jpg)
MODELOFÍSICO
MODELOECONÔMICO OTIMIZAÇÃO
Variáveis Especificadas
Variáveis de Projeto
Parâmetros Econômicos
ParâmetrosFísicos Dimensões Calculadas Lucro
Objetivo de uma Estratégia de Cálculo
Minimizar o esforço computacional envolvido na resolução dos modelos (problemas de dimensionamento, simulação e otimização de processos).
![Page 56: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/56.jpg)
3.1.1 Representação 3.1.2 Métodos Numéricos3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e representação 3.2.2 Resolução: partição, abertura, ordenação de equações3.3 Dimensionamento e Simulação de Equipamentos3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global 3.4.2 Estratégia Modular3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Projeto 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.1 Equações Não - Lineares
![Page 57: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/57.jpg)
3.1 EQUAÇÕES NÃO – LINEARES
Motivação para o estudo de equações não-lineares
No dimensionamento e na simulação de equipamentos e de processos ocorrem sistemas de equações que só podem ser
resolvidos por métodos iterativos de tentativas
empregados na resolução deequações
Assim, este Capítulo começa com métodos de resolução de equações.
![Page 58: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/58.jpg)
x
x*1 2 3 4 5 6 7 8
x1
x2
x3
x4
x5
x6
x7
x8
6
o
Consiste em decompor o sistema em sub-sistemas
PARTIÇÃO"partitioning"
1. f1(xo,x1) = 02. f2(x1,x2) = 03. f3(x2,x3,x6) = 04. f4(x3,x4) = 05. f5(x4,x5) = 06. f6(x5,x6) = 07. f7(x6,x7) = 08. f8(x7,x8) = 0
Resolvem-se os sub-sistemas sequencialmente
1, 2[ ]
Parte AcíclicaCálculo Direto
xo* x2[ 3, 4 , 5 ,6 ]
Parte CíclicaCálculo Iterativo
x6[7, 8]
x8
Parte AcíclicaCálculo Direto
![Page 59: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/59.jpg)
3.1 Equações Não-Lineares
3.1.2 Métodos Numéricos3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e representação 3.2.2 Resolução: partição, abertura, ordenação de equações3.3 Dimensionamento e Simulação de Equipamentos3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global 3.4.2 Estratégia Modular3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Projeto 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.1.1 Representação
![Page 60: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/60.jpg)
Na abordagem aqui adotada, uma vez formulada representando um fenômeno físico, a equação
f (x1*, ..., xi - 1
*, xi, xi + 1*,…, xM
*) = 0
passa ser vista como um “processador de informação” assim representada :
3.1 EQUAÇÕES NÃO – LINEARES
3.1.1 Representação
fj
. . . . . . . . .
. . . . . . . . .
x1
x2 x i - 1
x i + 1xM
x iincógnita
variáveis conhecidas
![Page 61: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/61.jpg)
A dificuldade da resolução de
f (x1*, ..., xi - 1
*, xi , xi + 1*,…, xM
*) = 0
depende da sua forma funcional.
incógnita x2: x1* x2 + ln x1
* = 0
incógnita x1: x1 x2* + ln x1 = 0
Solução analítica simples: x2 = - (ln x1
*) / x1*
Solução numérica por tentativas
Exemplo
x1 x2 + ln x1 = 0
![Page 62: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/62.jpg)
3.1 Equações Não-Lineares 3.1.1 Representação
3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e representação 3.2.2 Resolução: partição, abertura, ordenação de equações3.3 Dimensionamento e Simulação de Equipamentos3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global 3.4.2 Estratégia Modular3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Dimensionamento 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.1.2 Métodos Numéricos
![Page 63: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/63.jpg)
Métodos de Aproximações Sucessivas
Há duas famílias importantes de métodos numéricos para a resolução de equações não-lineares.
Métodos deRedução de Intervalos
Por diferentes raciocínios lógicos, promovem a redução do intervalo até que se torne menor do que uma tolerância pré-estabelecida.
Por diferentes raciocínios lógicos, testam valores sucessivos até que a diferença relativa entre 2 valores sucessivos se torne menor do que uma tolerância pré-estabelecida.
Partem de um intervalo inicial.(limites inferior e superior)
Partem de um valor inicial.
3.1.2 Métodos Numéricos
![Page 64: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/64.jpg)
Dados os limites superior xs e inferior xi , define-se o intervalo de incerteza xs - xi .
Qualquer valor no interior ou na fronteira do intervalo serve como solução.
xs
xi
(a) Métodos de Redução de Intervalos
Este é reduzido sucessivamente até se tornar menor do que uma tolerância pré-estabelecida: xs - xi .
xi
xs
xi
xs
f (x)
![Page 65: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/65.jpg)
Um método típico de Redução de Intervalos
Método da Bisseção ou Busca Binária
A cada iteração, o intervalo de incerteza é reduzido à metade.
![Page 66: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/66.jpg)
x
ALGORITMO
SE ABS(fi) < ABS(fs) ENTÃO Solução = xi SENÃO Solução = xs
f(x)
x
xi
fi
xs
fs
x
f
xs
fs
xi
fi
x
f
f(x)
Se Sinal (f) = Sinal (fi): Então atualizar : xi = x : fi = f
Estabelecer xi, xs, (tolerância)Calcular fi em xi
Calcular fs em xs
REPETIR
x = (xi + xs)/2Calcular f em x
Senão atualizar : xs = x : fs = f
ATÉ xs - xi
x fBISS
f (x)
![Page 67: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/67.jpg)
Exemplo: x1 x2* + ln x1 = 0
Solução para = 0,1 : x = 0,4375 f = 0,048
xi fi x f xs fs
0,00005 -11,51 1 2 1
0,00005 -11,51
0,5 0,307
0,5 0,307
0,375 -0,231
0,5 0,307
0,25 -0,88
0,375 -0,231
0,4375 0,048
0,5 0,307
0,25 -0,88
0,375 -0,231
0,4375 0,048
0,5
0,25
0,125
0,0625
f = x1 x2* + ln x1
x2* = 2 : xi = 0 : xs = 1: = 0,1
Se Sinal (f) = Sinal (fi): Então atualizar : xi = x : fi = f
Estabelecer xi, xs, (tolerância)Calcular fi em xi
Calcular fs em xs
REPETIR
x = (xi + xs)/2Calcular f em x
Senão atualizar : xs = x : fs = f
ATÉ xs - xi
![Page 68: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/68.jpg)
EFICIÊNCIA DO MÉTODO
Nt : número total de cálculos da função para alcançar o intervalo
Nm : número de cálculos da função, no meio do intervalo, necessário para alcançar o intervalo
Como o método se inicia com o cálculo da função nos limites inferior e superior, então:
Nt = Nm + 2
![Page 69: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/69.jpg)
Solução para = 0,1 : x = 0,4375 f = 0,048
xi fi x f xs fs
0,00005 -11,51 1 2 1
0,00005 -11,51
0,5 0,307
0,5 0,307
0,375 -0,231
0,5 0,307
0,25 -0,88
0,375 -0,231
0,4375 0,048
0,5 0,307
0,25 -0,88
0,375 -0,231
0,4375 0,048
0,5
0,25
0,125
0,0625
= 0,5Nm
ln = Nm ln 0,5
Nt = 2 + ln / ln 0,5
10% : = 0,1 N = 5,3 Nt = 6
1% : = 0,01 N = 8,6 Nt = 9
Nt = 2 – 1,4 ln
![Page 70: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/70.jpg)
Atribui-se um valor inicial para a incógnita.
(b) Métodos de Aproximações Sucessivas
xi xs
x1 x2 x3
Esse valor é atualizado sucessivamente até que o erro relativo entre duas aproximações sucessivas, abs [(xk - xk-1) / xk], seja menor do que uma tolerância pré-estabelecida.
x4 x5
![Page 71: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/71.jpg)
Um método típicoMétodo da Substituição Direta
Se a incógnita aparecer em mais de termo da equação, ela é explicitada parcialmente:
f(xi ) = 0 xi = F(xi)
Exemplo
x1 = e - x1
x2
*
x1 = - (1/ x2*) ln x1
F(x1)
F(x1)
x1 x2* + ln x1 = 0
f(x1)
![Page 72: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/72.jpg)
xi = F(xi)
F(x)
x
A solução é o valor de xi em que F(xi) = xi .
0,2
0,2
![Page 73: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/73.jpg)
ALGORITMO
Convergir = |(F-x)/x| < (erro relativo)
Estabelecer xinicial, (tolerância)
REPETIR x = F Calcular a Função F em xATÉ Convergir
xsolução = F
F = xinicial
Como dar a partida ?
Em cada iteração, o valor arbitrado para xi é o valor de F(xi - 1) obtido na iteração anterior.
![Page 74: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/74.jpg)
F(x)
x
Convergir = |(F-x)/x| < (erro relativo)
Estabelecer xinicial, (tolerância)F = xinicial
REPETIR x = F Calcular a Função F em xATÉ Convergirxsolução = F
x1x2x3
Executando o Algoritmo
Em cada iteração, o valor arbitrado para xi é o valor de F(xi - 1) obtido na iteração anterior.
![Page 75: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/75.jpg)
Condição para Convergência |F´(x)| < 1
F(x)
xx1
x2x3
convergência monotônicaderivada positiva
convergência oscilatória derivada negativa
x1
x3
x2
F(x)
x
Na direção da Solução
A convergência pode ser monotônica ou oscilatória
![Page 76: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/76.jpg)
Condição para Divergência |F´(x)| > 1
F(x)
xx1x3
x2
F(x)
xx1x2x3
divergência monotônicaderivada positiva
divergência oscilatóriaderivada negativa
Afastamento da Solução
A divergência pode ser monotônica ou oscilatória
![Page 77: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/77.jpg)
Exemplo: x1 x2* + ln x1 = 0
x1 = F(x1)(x2
* = 2 : x1 inicial = 0,5)
F(x1) = - (1/ x2*) ln x1 F(x1) = e - x
1 x
2*
Divergência Oscilatória F’(x1) = - 1,17
Convergência Oscilatória F’(x1) = - 0,85
Solução: x = 0,4263
F(x)
x1
x2
x3
x
x1
x3
x2
F(x)
x
x F 0,5 0,346 0,3080,346 0,529 0,529
0,529 0,317 0,400
0,317 0,573 0,806
0,573 0,278 0,515
x F 0,5 0,367 0,2640,367 0,479 0,302
0,479 0,383 0,199
0,383 0,464 0,210
0,464 0,395 0,149
Exemplo: x1 x2* + ln x1 = 0
oscilando
para maior
oscilando
para menor
![Page 78: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/78.jpg)
F(x)
xx1
x2x3
F(x)
x
x1x2x3
Equações que representam processos de operação estável, também apresentam comportamento estável, ou seja,
convergem.
![Page 79: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/79.jpg)
Em resumo
Equações Não-Lineares podem ser resolvidas por métodos:
- redução de intervalos (ex.: bisseção)
- aproximações sucessivas (ex.: substituição direta)
Esses métodos serão evocados a seguir em
Sistemas de Equações.
![Page 80: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/80.jpg)
3.1 Equações Não-Lineares 3.1.1 Representação 3.1.2 Métodos Numéricos
3.2.2 Resolução: partição, abertura, ordenação de equações3.3 Dimensionamento e Simulação de Equipamentos3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global 3.4.2 Estratégia Modular3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Dimensionamento 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.2 Sistemas de Equações Não-Lineares3.2.1 Estrutura e Representação
![Page 81: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/81.jpg)
A equação
f (x1, ..., xi-1, xi, xi+1,…, xM) = 0
pode ser representada como um “processador de informação”
3.2 SISTEMAS DE EQUAÇÕES NÃO – LINEARES 3.2.1 Estrutura e Representação
fj
. . . . . . . . .
. . . . . . . . .
x1
x2 x i - 1
x i + 1xM
x i
![Page 82: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/82.jpg)
Um sistema de equações pode ser representado por um
sistema de processadores
Os elementos desse sistema são as equações.
As conexões são as variáveis comuns.
Durante a resolução de um problema, os processadores transmitem informação de uns para os outros.
f1(xo,x1) = 0f2(x1,x2) = 0f3(x2,x3) = 0
1 2 3x x
1x
2x
30
![Page 83: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/83.jpg)
f1(xo,x1) = 0f2(x1,x2) = 0f3(x2,x3) = 0
1 2 3x x
1x
2x
30
Estrutura Acíclica
f1(xo,x1,x3) = 0f2(x1,x2) = 0f3(x2,x3) = 0
1 2 3x0
x1
x2
x3
x3
Estrutura Cíclica
Estruturas Básicas
Os sistemas de equações podem assumir as mais variadas estruturas.
![Page 84: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/84.jpg)
1 2 3x x
1x
2x
30
Estrutura Acíclica
1 2 3x0
x1
x2
x3
x3
Estrutura Cíclica
Estrutura acíclica: resolução trivial por encadeamento sucessivo a partir que qualquer variável conhecida (xo, por exemplo).
Estrutura cíclica: solução somente por tentativas (exemplo: conhecida xo, o cálculo de x1 depende de x3 ainda não calculada).
?
f1(xo, x1) = 0f2(x1, x2) = 0f3(x2, x3) = 0
f1(xo, x1, x3 ) = 0f2(x1, x2) = 0f3(x2, x3) = 0
![Page 85: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/85.jpg)
Quanto mais complexa a estrutura, mais difícil é a resolução do sistema.
![Page 86: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/86.jpg)
Um Sistema de Equações Típico de um Modelo de Processo
1. f1(xo*, x1) = 02. f2(x1, x2) = 03. f3(x2, x3, x6) = 04. f4(x3, x4) = 05. f5(x4, x5) = 06. f6(x5, x6) = 07. f7(x6, x7) = 08. f8(x7, x8) = 0
![Page 87: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/87.jpg)
Representações
Através da representação é possível enxergar como os processadores trocam informações.
E, assim, conceber métodos eficientes de resolução.
![Page 88: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/88.jpg)
X0 X1 X2 X3 X4 X5 X6 X7 X8
1 1 1 0 0 0 0 0 0 02 0 1 1 0 0 0 0 0 03 0 0 1 1 0 0 1 0 04 0 0 0 1 1 0 0 0 05 0 0 0 0 1 1 0 0 06 0 0 0 0 0 1 1 0 07 0 0 0 0 0 0 1 1 08 0 0 0 0 0 0 0 1 1
Matriz Incidência (Numérica)
X0 X1 X2 X3 X4 X5 X6 X7 X8
1 * *2 * *3 * * *4 * *5 * *6 * *7 * *8 * *
Matriz Incidência (Gráfica)
Matrizes Esparsas !
1. f1(xo*,x1) = 02. f2(x1,x2) = 03. f3(x2,x3,x6) = 04. f4(x3,x4) = 05. f5(x4,x5) = 06. f6(x5,x6) = 07. f7(x6,x7) = 08. f8(x7,x8) = 0
Representação Matricial
Característica em Processos
o número de variáveis em cada equação é pequeno: nem todas as variáveis figuram em todas as equações.
![Page 89: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/89.jpg)
1. f1(xo*,x1) = 0
2. f2(x1,x2) = 03. f3(x2,x3,x6) = 04. f4(x3,x4) = 05. f5(x4,x5) = 06. f6(x5,x6) = 07. f7(x6,x7) = 08. f8(x7,x8) = 0
Representação Gráfica (Grafo)
Ciclo !
Proporciona uma visão mais clara da estrutura do sistema
x6
1 2 3 4 5 6 7 8x1 x2 x3 x4 x5 x6
x7 x8xo
![Page 90: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/90.jpg)
3.1 Equações Não-Lineares 3.1.1 Representação 3.1.2 Métodos Numéricos3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e Representação
3.3 Dimensionamento e Simulação de Equipamentos3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global 3.4.2 Estratégia Modular3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Dimensionamento 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.2.2 Resolução: partição, abertura, ordenação de equações
![Page 91: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/91.jpg)
3.2.2 Resolução: partição, abertura, ordenação de equações
Os sistemas de equações podem ser resolvidos por - métodos simultâneos- método seqüencial.
![Page 92: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/92.jpg)
Métodos Simultâneos
Calcular F1
x1(k+1) = F1
Calcular F2
x2(k+1) = F2
TESTE
TESTE
x1 = x1(k+1)
x1k
x2k
x1(k+1)
x2(k+1)
x2 = x2(k+1)
Diversos métodos são descritos em livros texto e abordados em disciplinas de Métodos Numéricos.
Exemplo: Newton-Raphson, Wegstein, ...
Todas as variáveis são alteradas simultaneamente.
![Page 93: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/93.jpg)
Método Sequencial
É um método alternativo em que as equações são acionadas uma-a-uma, passando informação de uma para a outra, numa sequência lógica previamente estabelecida.
Este método é chamado de Resolução por Equações ("equation oriented").
![Page 94: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/94.jpg)
Este método pode ser implementado através do
ALGORITMO DE ORDENAÇÃO DE EQUAÇÕES
A.O.E.
![Page 95: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/95.jpg)
É um algoritmo de atribuição de tarefas
Algoritmo de Ordenação de Equações (AOE)
1. Atribui a cada equação a tarefa de calcular uma das incógnitas do sistema.
2. Ao mesmo tempo, organiza as equações segundo uma
Sequencia de Cálculo
que minimiza o esforço computacional.
![Page 96: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/96.jpg)
Outros resultados
4. Em problemas com graus de liberdade, indica as variáveis de projeto compatíveis com o esforço computacional mínimo.
3. Efetua automaticamente a partição do sistema em conjuntos cíclicos e acíclicos de equações, minimizando o número de equações envolvidas em cálculos iterativos.
5. Em problemas com ciclos, indica as variáveis de abertura(a definir).
Partição ???
![Page 97: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/97.jpg)
x
x*1 2 3 4 5 6 7 8
x1
x2
x3
x4
x5
x6
x7
x8
6
o
Consiste em decompor o sistema em sub-sistemas
PARTIÇÃO"partitioning"
1. f1(xo,x1) = 02. f2(x1,x2) = 03. f3(x2,x3,x6) = 04. f4(x3,x4) = 05. f5(x4,x5) = 06. f6(x5,x6) = 07. f7(x6,x7) = 08. f8(x7,x8) = 0
Resolvem-se os sub-sistemas sequencialmente
1, 2[ ]
Parte AcíclicaCálculo Direto
xo* x2[ 3, 4 , 5 ,6 ]
Parte CíclicaCálculo Iterativo
x6[7, 8]
x8
Parte AcíclicaCálculo Direto
![Page 98: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/98.jpg)
O Algoritmo simplesmente formaliza
ações intuitivas e óbvias
utilizando os seguintes termos básicos
Equações de Incógnita Única
Variáveis de Frequência Unitária
Ciclos
![Page 99: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/99.jpg)
x
x*1 2 3 4 5 6 7 8
x1
x2
x3
x4
x5
x6
x7
x8
6
o
1. f1(xo,x1) = 02. f2(x1,x2) = 03. f3(x2,x3,x6) = 04. f4(x3,x4) = 05. f5(x4,x5) = 06. f6(x5,x6) = 07. f7(x6,x7) = 08. f8(x7,x8) = 0
Equações de Incógnita Única
São equações em que todas as variáveis têm os seus valores conhecidos, menos uma!
Exemplo: equação 1
Uma vez resolvida para x1
a equação 2 fica com incógnita única
podendo ser resolvida para x2
Não há mais equações de incógnita única
![Page 100: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/100.jpg)
Enquanto houver equações com incógnita única
(a) atribuir (vincular) essa incógnita à respectiva equação
O Algoritmo pode começar assim:
(c) remover a variável da lista das incógnitas
(b) colocar a equação no primeira posição disponível na Sequencia de Cálculo
x
x*1 2 3 4 5 6 7 8
x1
x2
x3
x4
x5
x6
x7
x8
6
o
2. x2
1. x1Não há mais equações de incógnita única
Sequênciade Cálculo
![Page 101: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/101.jpg)
Varáveis de Frequência Unitária
São variáveis que pertencem a uma só equação
Só pode ser calculada pela Eq. 8 depois de x7
x
x*1 2 3 4 5 6 7 8
x1
x2
x3
x4
x5
x6
x7
x8
6
o
Exemplo: x8
Então x7 só pode ser calculada pela Eq.7 ficando com frequência unitária
Só pode ser calculada pela Eq. 7 depois de x6
Não há mais variáveis de frequência unitária
![Page 102: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/102.jpg)
Enquanto houver variáveis de frequência unitária
(a) atribuir (vincular) essa incógnita à respectiva equação
O Algoritmo pode prosseguir assim:
(c) remover a variável da lista das incógnitas
(b) colocar a equação na última posição disponível na Sequencia de Cálculo
1. x1
2. x2
-
-
-
-x
x*1 2 3 4 5 6 7 8
x1
x2
x3
x4
x5
x6
x7
x8
6
o
8. x8
7. x7
Não há mais variáveis de frequência unitária
Sequênciade Cálculo
![Page 103: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/103.jpg)
Ciclos
x
3 4 5 6x
3x
4x
5
6
x2 x6
x6 = f6(x5) = f6(f5(x4)) = f6(f5(f4(x3))) = f6(f5(f4(f3(x2,x6)))) = F(x6)
São conjuntos cíclicos de equações em que cada variável vem a ser função dela mesma.
Solução exclusivamente por métodos iterativos
![Page 104: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/104.jpg)
Preparação do sub-sistema cíclico para resolução por tentativas
(d) Estabelecer o esquema de convergência
(a) Selecionar uma Equação Final
(b) Retornar à Etapa 2 (VFU)
(c) Identificar a Variável de Abertura (não atribuída a qualquer equação)
3 4 5
X3
X4
X5
3 4 5 6
X3
X4
X5
Equação Final
3 4 5
X3
X4
X5
3 4 5 6
X3
X4
X5
X6Variável de Abertura 1. x1
2. x2
8. x8
7. x7
5. x5
4. x4
3. x3
6. final
x6
Sequênciade Cálculo
![Page 105: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/105.jpg)
1 2
Xo
*X
1X
27 8
X6
X7
X8
3 4 5
X3
X4
X5
3 4 5 6
X3
X4
X5
Equação Final
3 4 5
X3
X4
X5
3 4 5 6
X3
X4
X5
EQUAÇÃO VARIÁVEL
1 x12 x2
7 x7
6 final
8 x8
EQUAÇÃO VARIÁVEL
1 x12 x2
7 x7
3 x34 x45 x56 final
8 x8
X6Variável de Abertura
x6
META DO ALGORITMO
Sequencia de Cálculo Resultante
![Page 106: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/106.jpg)
ALGORITMO PRONTO
![Page 107: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/107.jpg)
Algoritmo de Ordenação de Equações
Enquanto houver equações
Enquanto houver equações com incógnita única
(a) atribuir (vincular) essa incógnita à respectiva equação.(b) colocar a equação no primeira posição disponível na Sequencia de Cálculo.(c) remover a variável (X na vertical).
Enquanto houver variáveis de frequência unitária
(a) atribuir (vincular) essa variável à respectiva equação.(b) colocar a equação no última posição disponível na Seqüência de Cálculo.(c) remover a equação (X na horizontal).Se ainda houver equações (ciclo!)
(a) selecionar uma equação que contenha pelo menos uma variável de freqüência igual à menor freqüência dentre todas as variáveis (Final).(b) colocar essa equação na última posição disponível na Seqüência de Cálculo.(c ) remover equação (X na horizontal).
![Page 108: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/108.jpg)
ESTABELECIMENTO DO ESQUEMA DE CONVERGÊNCIA
Insere-se um Promotor de Convergência
São apresentados dois Promotores de Convergência baseados nos métodos para equações não-lineares:
(a) Bisseção
(b) Substituição Direta
![Page 109: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/109.jpg)
Uma vez ordenadas, as equações podem ser resolvidas na seqüência estabelecida, com o mínimo de esforço computacional.
![Page 110: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/110.jpg)
x
ALGORITMO
SE ABS(fi) < ABS(fs) ENTÃO Solução = xi SENÃO Solução = xs
f(x)
x
xi
fi
xs
fs
x
f
xs
fs
xi
fi
x
f
f(x)
Se Sinal (f) = Sinal (fi): Então atualizar : xi = x : fi = f
Estabelecer xi, xs, (tolerância)Calcular fi em xi
Calcular fs em xs
REPETIR
x = (xi + xs)/2Calcular f em x
Senão atualizar : xs = x : fs = f
ATÉ xs - xi
x fBISS
f (x)
Relembrando o Método da Bisseção
![Page 111: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/111.jpg)
A cada iteração:- arbitra-se x6a .- resolve-se sucessivamente as equações 3, 4 e 5.- pela equação 6 calcula-se f6 (x5, x6).- avalia-se a convergência pelo critério do método da bisseção.
x
xi
fi
xs
fs
x
f
f(x)
(a) BISSSEÇÃOx6a
3 4 5 6x2 x3 x4 x5 x6
BISSf6 (x5, x6)
f3 (x2,x3,x4) = 0 f4 (x3,x4) = 0 f5 (x4,x5) = 0 f6 (x5,x6) = 0
![Page 112: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/112.jpg)
(c)
F'(x) < 0 |F'(x)| < 1
convergência oscilatória
F(x)
xx1x2 x3
(a)
F'(x) > 0 |F'(x) < 1
convergência monotonica
x1x2
F(x)
x
x3
ALGORITMO
Estabelecer xinicial, (tolerância)
F = xinicial
xsolução = F
Convergir = |(F - x)/x| <
REPETIR
x = F
Calcular a Função F em x
ATÉ Convergir
RELEMBRANDO O MÉTODO DA SUBSTITUIÇÃO DIRETA
x = F ( x)f (x) = 0 explicitando x
A solução é a interseção de F(x) com a reta de 45º onde ela é igual ao próprio x.
![Page 113: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/113.jpg)
(b) SUBSTITUIÇÃO DIRETA
Arbitra-se x6c inicial.A cada iteração:- toma-se x6a = x6c . - resolve-se sucessivamente as equações 3, 4, 5 e 6, que calcula x6c.- avalia-se a convergência através do erro relativo ABS (x6c – x6a) / x6a
x1
x2
x3
x6c
x6a
x6c
3 4 5 6x2 x3 x4 x5 x6
SDx6a
![Page 114: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/114.jpg)
COMPARAÇÃO DOS PROMOTORES DE CONVERGÊNCIA
(b) Substituição Direta
Arbitra-se x6a . A cada iteração, a eq.6 calcula x6c = f6(x5) : x6a = x6c (até convergir).
f6 (x5, x6)(a) Bisseção
x6a
3 4 5 6x2 x3 x4 x5 x6
BISS
f3 (x2,x3,x6) = 0 f4 (x3,x4) = 0 f5 (x4,x5) = 0 f6 (x5,x6) = 0
Arbitra-se x6a. A cada iteração, a eq.6 calcula f6 (x5, x6) (até convergir)
x6c
3 4 5 6x2 x3 x4 x5 x6
SDx6a
f3 (x2,x3,x4) = 0 f4 (x3,x4) = 0 f5 (x4,x5) = 0 f6 (x5,x6) = 0
![Page 115: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/115.jpg)
APLICAÇÃO AO SISTEMA ILUSTRATIVO
![Page 116: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/116.jpg)
1. f1(xo*, x1) = 02. f2(x1, x2) = 03. f3(x2, x3, x6) = 04. f4(x3, x4) = 05. f5(x4, x5) = 06. f6(x5, x6) = 07. f7(x6, x7) = 08. f8(x7, x8) = 0
Do ponto de vista prático, o Algoritmo é aplicado sobre a Matriz Incidência
![Page 117: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/117.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X *
2 * *
3 * * *
4 * *
5 * *
6 * *
7 * *
8 * *
1 -
2 -
3 -
4 -
5 -
6 -
7 -
8 -
Seqüência
Equações de Incógnita Única (EIU)
Círculo na variável inscrição no último lugar x na horizontal
![Page 118: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/118.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X *
3 * * *
4 * *
5 * *
6 * *
7 * *
8 * *
Seqüência
1 - x1
2 -
3 -
4 -
5 -
6 -
7 -
8 -
Equações de Incógnita Única (EIU)
Círculo na variável inscrição no primeiro lugar x na vertical
![Page 119: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/119.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 * *
7 * *
8 * *
1 - x1
2 - x2
3 -
4 -
5 -
6 -
7 -
8 -
Seqüência
Equações de Incógnita Única (EIU)
Círculo na variável inscrição no primeiro lugar x na vertical
Não há mais Equações de Incógnita Única (EIU)
![Page 120: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/120.jpg)
Estado atual da Sequência de Cálculo
![Page 121: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/121.jpg)
EQUAÇÃO VARIÁVEL1 x12 x2
1 2
XO
*X
1X
2
![Page 122: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/122.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 * *
7 * *
8 * *
Seqüência
1 - x1
2 - x2
3 -
4 -
5 -
6 -
7 -
8 -
Variáveis de Frequência Unitária (VFU)
Círculo na variável inscrição no último lugar x na horizontal
![Page 123: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/123.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 * *
7 * *
8 X O
1 - x1
2 - x2
3 -
4 -
5 -
6 -
7 -
8 - x8
Seqüência
Variáveis de Frequência Unitária (VFU)
Círculo na variável inscrição no último lugar x na horizontal
![Page 124: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/124.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 * *
7 X O
8 X O
1 - x1
2 - x2
3 -
4 -
5 -
6 -
7 - x7
8 - x8
Seqüência
Variáveis de Frequência Unitária (VFU)
Círculo na variável inscrição no último lugar x na horizontal
![Page 125: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/125.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 * *
7 X O
8 X O
1 - x1
2 - x2
3 -
4 -
5 -
6 -
7 - x7
8 - x8
Seqüência
Não há mais Variáveis de Frequência Unitária (VFU)
![Page 126: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/126.jpg)
Estado atual da Sequência de Cálculo
![Page 127: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/127.jpg)
EQUAÇÃO VARIÁVEL1 x12 x2
7 x78 x8
1 2
XO
*X
1X
27 8
X7
X8
X6
![Page 128: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/128.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 * *
7 X O
8 X O
1 - x1
2 - x2
3 -
4 -
5 -
6 -
7 - x7
8 - x8
Seqüência
Ciclo!
Equação Final ?
x
3 4 5 6x
3x
4x
5
6
x2 x6
![Page 129: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/129.jpg)
x
3 4 5 6x
3x
4x
5
6
x2 x6
Em princípio, qualquer uma do ciclo
A figura motiva a 6 !
![Page 130: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/130.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 X X
7 X O
8 X O
Equação Final: 6
1 - x1
2 - x2
3 -
4 -
5 -
6 final
7 - x7
8 - x8
Seqüência
![Page 131: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/131.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 X X
7 X O
8 X O
1 - x1
2 - x2
3 -
4 -
5 -
6 - final
7 - x7
8 - x8
Seqüência
Variáveis de Frequência Unitária (VFU)
Volta-se a buscar...
![Page 132: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/132.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 * *
6 X X
7 X O
8 X O
1 - x1
2 - x2
3 -
4 -
5 -
6 - final
7 - x7
8 - x8
Seqüência
Aqui há 2 VFU: X5 e X6.
A ordem em que são escolhidas não importa. Uma não afeta a outra.
Ela poderiam ser até resolvidas em paralelo.
Na verdade, elas poderiam ser calculadas em paralelo porque
![Page 133: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/133.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 * *
5 X O
6 X X
7 X O
8 X O
1 - x1
2 - x2
3 -
4 -
5 - x5
6 - final
7 - x7
8 - x8
Seqüência
Variáveis de Frequência Unitária (VFU)
Círculo na variável inscrição no último lugar x na horizontal
![Page 134: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/134.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X * *
4 X O
5 X O
6 X X
7 X O
8 X O
1 - x1
2 - x2
3 -
4 - x4
5 - x5
6 - final
7 - x7
8 - x8
Seqüência
Variáveis de Frequência Unitária (VFU)
Círculo na variável inscrição no último lugar x na horizontal
![Page 135: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/135.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X O X
4 X O
5 X O
6 X X
7 X O
8 X O
1 - x1
2 - x2
3 - x3
4 - x4
5 - x5
6 - final
7 - x7
8 - x8
Seqüência
Variáveis de Frequência Unitária (VFU)
Círculo na variável inscrição no último lugar x na horizontal
![Page 136: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/136.jpg)
X0* X1 X2 X3 X4 X5 X6 X7 X8
1 X O
2 X O
3 X O X
4 X O
5 X O
6 X X
7 X O
8 X O
Variável de Abertura: x6
1 - x1
2 - x2
3 - x3
4 - x4
5 - x5
6 - final
7 - x7
8 - x8
x6
Seqüência
![Page 137: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/137.jpg)
SEQUÊNCIA DE CÁLCULO FINAL
EQUAÇÃO VARIÁVEL1 x12 x23 x34 x45 x56 final7 x78 x8
1 2 3 4 5 6 7 8
XO
*X
1X
2X
3X
4X
5X
6X
7X
8
X6
X6
variável de abertura
equação final
x6
![Page 138: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/138.jpg)
Mostrar o Programa AOE.xls
![Page 139: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/139.jpg)
Estratégia de Cálculo para
4 situações típicas em Engenharia de Processos.
Antecipando...
![Page 140: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/140.jpg)
PROCESSO
OTIMIZAÇÃO
*LEE x
13 2 1
x4
x3
x2
x1
x2
x3
x5 5
AVALIAÇÃO
ECONÔMICA
PROCESSO
OTIMIZAÇÃO
*LEE x
13 21x
4x
32x x
2x
3x
x4
5 5AVALIAÇÃO
ECONÔMICA
PROCESSO*
LEE x14 3 2 1
x4
x3
x2
x1
x2
x3
x4
AVALIAÇÃO
ECONÔMICASol.únicasem ciclo
Otimizaçãocom ciclo
Sol.únicacom ciclo
Otimizaçãosem ciclo
PROCESSO*LEE x
14 3 21x
4x
3x
2x 2 x 3 x
4AVALIAÇÃO
ECONÔMICA
x4
x1
x1
![Page 141: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/141.jpg)
1 f1(x1, x2)2 f2(x2, x3, x4) = 03 f3(x3, x4) = 04 f4(x4) = 0
Sistema 1
G = 0 : solução única, sem variável de projetoCiclo potencial: a confirmar, haverá uma variável de abertura
1 * *2 * * *3 * *4 *
x1 x2 x3 x4
Matriz Incidência
1 2 3 4x1 x2 x3 x4
x4Grafo
![Page 142: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/142.jpg)
Algoritmo de Ordenação de Equações
Enquanto houver equações com incógnita única
(c) remover a variável.
4 x4
3 x3
2 x2
1 x1
Seqüência de CálculoEquação Variável
Matriz Incidência x1 x2 x3 x4
1 * *
2 * * *
3 * *
4 *
X
X
X X
(a) atribuir (vincular) essa incógnita à respectiva equação.(b) colocar a equação na primeira posição disponível na Seqüência de Cálculo.
![Page 143: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/143.jpg)
G = 0 : Solução Única, sem variável de projetoProcesso : sequência direta (sem ciclos)
PROCESSO
LEE*4 3 2 1
x4 x3 x2 x1 AVALIAÇÃO
ECONÔMICA
x4 x3 x2 x1
1 2 3 4x1 x2 x3 x4
x4Grafo
Como a eq. 4 é de incógnita única, x4 é a primeira e o seu valor é transmitido para a eq.2 desfazendo o ciclo em potencial
![Page 144: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/144.jpg)
Sistema 2
1 f1(x1,x2)2 f2(x2,x3,x4) = 03 f3(x3,x4) = 04 f4(x4,x5) = 0
G = 1 : problema de otimização, haverá uma variável de projeto.Ciclo potencial: a confirmar, haverá uma variável de abertura.
1 * *x1 x2 x3 x4 x5
2 * * *3 * *4 * *
Matriz Incidência
1 2 3 4x1 x2 x3 x4
x4
x5
Grafo
![Page 145: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/145.jpg)
Algoritmo de Ordenação de Equações
Enquanto houver equações com incógnita única
(a) atribuir (vincular) essa incógnita à respectiva equação.(b) colocar a equação na primeira posição disponível na Seqüência de Cálculo.(c) remover a variável.
4 x5
3 x3
2 x2
1 x1
Seqüência de Cálculo
Equação VariávelMatriz Incidência
x1 x2 x3 x4 x5
1 * *
2 * * *
3 * *
4 * *
X
X
X X
X
Enquanto houver equações Enquanto houver variáveis de freqüência unitária
(c) remover a equação.
x4 variável de
projeto
(a) atribuir (vincular) essa variável à respectiva equação.
(b) colocar a equação na última posição disponível na Seqüência de Cálculo.
x4
![Page 146: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/146.jpg)
1 2 3 4x1 x2 x3 x4
x4
x5
Grafo
PROCESSO
OTIMIZAÇÃOLEE*
3 2 1x3
4
x2 x1
x5
x4
AVALIAÇÃO
ECONÔMICA
x1 x2 x3 x5
Como x4 resultou como Variável de Projeto o ciclo desapareceu
G = 1 : Otimização, uma variável de projetoProcesso : sequência direta (sem ciclos)
![Page 147: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/147.jpg)
Sistema 3
1 f1(x1,x2)2 f2(x1,x2,x3,x4) = 03 f3(x3,x4) = 04 f4(x4) = 0
G = 0: solução única, sem variável de projetoCiclos potenciais: a confirmar haverá variáveis de abertura
x1 x2 x3 x4
1 * *2 * * * *3 * *4 *
Matriz Incidência
1 2 3 4x2 x3 x4
x4
x1
Grafo
![Page 148: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/148.jpg)
Algoritmo de Ordenação de Equações
Enquanto houver equações com incógnita única
(a) atribuir (vincular) essa incógnita à respectiva equação.(b) colocar a equação na primeira posição disponível na Seqüência de Cálculo.(c) remover a variável.
4 x4
3 x3
1 x2
2 final
Seqüência de CálculoEquação Variável
Matriz Incidência x1 x2 x3 x4
1 * *
2 * * * *
3 * *
4 *
X
X XXX
X
x1: Variável de Abertura
x1
![Page 149: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/149.jpg)
x1 : variável de abertura
PROCESSO
LEE*4 3 21
x4 x3 x2
x1
AVALIAÇÃO
ECONÔMICA
x1 x2 x3 x4
1 2 3 4x2 x3 x4
x4
x1
Grafo
Aqui, o ciclo persistiu e teve que ser aberto
G = 0 : Solução Única, sem variável de projetoProcesso : sequência com ciclo e uma variável de abertura
![Page 150: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/150.jpg)
Sistema 4
1 f1(x1,x2)2 f2(x1,x2,x3,x4) = 03 f3(x3,x4) = 04 f4(x4,x5) = 0
G = 1: problema de otimização com uma variável de projetoCiclos potenciais: poderá haver variáveis de abertura
x1 x2 x3 x4 x5
1 * *
2 * * * *
3 * *
4 * *
Matriz Incidência
1 2 3 4x2 x3 x4
x4
x1
Grafo
x5
![Page 151: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/151.jpg)
Matriz Incidência x1 x2 x3 x4 x5
1 * *
2 * * * *
3 * *
4 * *
X
X X X
X 4 x5
3 x3
1 x2
2 final
Seqüência de Cálculo
Equação Variável
X
X
Enquanto houver equações
Enquanto houver variáveis de freqüência unitária
(a) atribuir (vincular) essa variável à respectiva equação.(b) colocar a equação na última posição disponível na Seqüência de Cálculo.(c) remover a equação.
Algoritmo de Ordenação de Equações
Enquanto houver equações com incógnita única
(a) atribuir (vincular) essa incógnita à respectiva equação.(b) colocar a equação na primeira posição disponível na Seqüência de Cálculo.(c) remover a variável.
Sobram x1 e x4
Uma de projeto
Outra de abertura
Eqs 1 e 3 independentes: qualquer uma
![Page 152: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/152.jpg)
x4: variável de aberturax1 : variável de projeto
1 2 3 4x2 x3
x4
x1
Grafo
x5
E*PROCESSO
OTIMIZAÇÃOLE
2 31x
4x
3
x4
2x
x1
5
x2
x3
x4
x5 AVALIAÇÃO
ECONÔMICA
G = 1: Otimização, com uma variável de projeto Processo: Sequência Cíclica com uma variável de abertura
Opção 1:
![Page 153: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/153.jpg)
x4: variável de projetox1 : variável de projeto
1 2 3 4x2 x3
x4
x1
Grafo
x5
E*PROCESSO
OTIMIZAÇÃOLE
2 31x
4x
3
x1
2x
x4
5
x2
x3
x4
x5 AVALIAÇÃO
ECONÔMICA
Opção 2:
G = 1: Otimização, com uma variável de projeto Processo: Sequência Cíclica com uma variável de abertura
![Page 154: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/154.jpg)
RESUMO DOS 4 TIPOS DE PROBLEMAS
![Page 155: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/155.jpg)
PROCESSO
OTIMIZAÇÃO
*LEE x
13 2 1
x4
x3
x2
x1
x2
x3
x5 5
AVALIAÇÃO
ECONÔMICA
PROCESSO
OTIMIZAÇÃO
*LEE x
13 21x
4x
32x x
2x
3x
x4
5 5AVALIAÇÃO
ECONÔMICA
PROCESSO*
LEE x14 3 2 1
x4
x3
x2
x1
x2
x3
x4
AVALIAÇÃO
ECONÔMICASol.únicasem ciclo
Otimizaçãocom ciclo
Sol.únicacom ciclo
Otimizaçãosem ciclo
PROCESSO*LEE x
14 3 21x
4x
3x
2x 2 x 3 x
4AVALIAÇÃO
ECONÔMICA
x4
x1
x1
![Page 156: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/156.jpg)
REGRAS COMPLEMENTARES NA APLICAÇÃO DOALGORITMO DE ORDENAÇÃO DE EQUAÇÕES
- Variáveis discretas- Variáveis de cálculo direto e iterativo- Variáveis limitadas- Ciclos múltiplos- Variáveis de abertura e de projeto- Eliminação de ciclos.
![Page 157: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/157.jpg)
Variáveis Discretas
Seus valores são limitados a um conjunto finito.
Exemplos: - tipos de insumos: utilidades, solventes, catalisadores.- diâmetros comerciais de tubos.- número de estágios.
Em problemas com G > 0 elas têm preferência como Variáveis de Projeto.Assim:
- assumem apenas os valores viáveis atribuídos pelo otimizador.- não sendo calculadas, não há risco de assumirem valores inviáveis.
![Page 158: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/158.jpg)
Para x = 1 : y = 3 a = 3 (não existe !)
Para: a = 0,5 : x = 1 y = 0,5Para: a = 1 : y = 1 x = 1
x
ya = 1
a = 0,5
V = 3 : N = 1 : G = 2 (duas variáveis de projeto)
Exemplo: y = a xO parâmetro só pode assumir 2 valores: 1 ou 0,5
Logo: a tem que ser uma das duas variáveis de projeto
![Page 159: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/159.jpg)
Variáveis de Cálculo Direto ou Iterativo
Exemplo
43
21
4321
TT
TTln
)TT(TT
Nesta equação:
- é uma variável de cálculo direto (dadas as temperaturas)- qualquer T é de cálculo iterativo (dado e as demais T’s)
As variáveis de cálculo direto têm preferência para a condição de calculadas.
Variáveis de cálculo direto são aquelas que podem ser facilmente explicitadas numa equação e calculadas sem necessidade de iterações.
![Page 160: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/160.jpg)
Varáveis Limitadas
Os seus valores variam entre limites bem definidos.
Exemplos:- frações mássicas ou molares- temperaturas em trocadores de calor
Variáveis limitadas devem ter preferência para atuar como variáveis de abertura e de
projeto.
Durante a execução do Algoritmo, a atribuição deve ser postergada ao máximo
para que essa preferência seja concretizada.
![Page 161: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/161.jpg)
Ciclos Múltiplos
f1(xo,x1,x3) 0
f2(x1,x2) 0
f3(x2,x3) 0
f4(x3,x4) 0
f5(x4,x5,x7) 0
f6(x5,x6) 0
f7(x6,x7) 0
=
=
=
=
=
=
=
1
2
4
5
6
1. x
2. x
3. final
4. x
5. x
6. x
7. final
x3
x7
Ciclos em Sequência
Primeira entrada de x7: eq. 5
Primeira entrada de x3: eq. 1
Fechar o ciclo com a final mais próxima
Um sistema de equações pode exibir diversos ciclos.
![Page 162: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/162.jpg)
Ciclos Aninhados (“nested”)
0
0
0
0
f x x x
f x x x
f x x x
f x x
f x x x
f x x x
f x x
1 o 1 7
2 1 2 6
3 2 3 5
4 3 4
5 3 4 5
6 5 6 7
7 6 7
0
0
0
( , , )
( , , )
( , , )
( , )
( , , )
( , , )
( , )
=
=
=
=
=
=
=
X4
X7
1. x1
4. x3
6. x5
3. x2
5. final
7. x6
2. final
Ciclos Múltiplos
Primeira entrada de x7: eq. 7
Primeira entrada de x4: eq. 4
Fechar o ciclo com a final mais próxima
![Page 163: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/163.jpg)
Variáveis de Abertura e de Projeto Simultâneas
(a)
1. x1
2. x2
4. x4
6. x 6
3. x 3
x7
x5
5. final
7. x8
Escolha ConvenienteCiclo com 3 equações
1. x1
2. x2
4. x4
6. x 6
3. x 3
x7
x5
5. final
(b)
7. x8
Escolha InconvenienteCiclo com 4 equações
Em problemas com G > 0 e com ciclo, a variável de abertura deve ser aquela que fecha, com a Equação Final, um ciclo com o
menor número de equações.
![Page 164: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/164.jpg)
Eliminação de Ciclos
31. x31 = 1 – x11* 32. x32 = 1 – x12
04. x13 = k x12 / [1 + (k – 1) x12]07. W3 = W1* x11* r / x13
01. W2 = W1* x31 / x32
02. W1* x11* – W2 x12 – W3 x13 = 0
x12
Equação Final 02. W1* x11* – W2 x12 – W3 x13 = 0 W2 da eq.01:02’. W1* x11* - [W1* x31 / x32] x12 – W3 x13 = 0 W3 da eq. 07:02’. W1* x11* - [W1* x31 / x32] x12 – [W1* x11* r / x13] x13 = 0 x13 da eq. 04 e x32 da eq.32:02’. x12 = x11
* (1 – r) / [x31 + x11* (1 – r)]
x31 calculado antes do ciclo
31. x31 = 1 – x11* 02’. x12 = x11
* (1 – r) / [x31 + x11* (1 – r)]
32. x32 = 1 – x12
04. x13 = k x12 / [1 + (k – 1) x12]07. W3 = W1* x11* r / x13
01. W2 = W1* x31 / x32
Sequência com Ciclo Sequência sem Ciclo
Substituição Algébrica
![Page 165: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/165.jpg)
3.1 Equações Não-Lineares 3.1.1 Representação 3.1.2 Métodos Numéricos3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e representação 3.2.2 Resolução: partição, abertura, ordenação de equações
3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global 3.4.2 Estratégia Modular3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Dimensionamento 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.3 Dimensionamento e Simulação de Equipamentos
![Page 166: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/166.jpg)
3.3 Dimensionamento e Simulação dos Equipamentos
Adquirir familiaridade com os equipamentos antes de integrá-los no processo (livres de interações).
Motivação para estudar os equipamentos isolados
Montar as rotinas de dimensionamento e de simulaçãoque integram o programa de análise do processo.
Rever conhecimentos adquiridos em disciplinas anteriormente cursadas.
CIÊNCIAS BÁSICAS
FUNDAMENTOS
ENG. DE EQUIPAMENTOS
Analogia: estudar os instrumentos isoladamente antes de compor a melodia para a orquestra
![Page 167: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/167.jpg)
ENGENHARIA DE EQUIPAMENTOS
Projeto e Análise dos Equipamentosde Processo
ReatoresTrocadores de calorSeparadores
Torres de destilaçãoTorres de absorçãoExtratoresCristalizadoresFiltrosOutros...
Instrumentos de Controle Automático
CIÊNCIAS BÁSICAS
FUNDAMENTOS
ENG. DE EQUIPAMENTOS
Tratamento compartimentado!
![Page 168: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/168.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 169: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/169.jpg)
Dimensões dos principais Equipamentos.
Consumo de utilidadesmatérias primas e insumos
Especificaçõesde projeto
Modelo Matemático
previsão
Dimensões dos principais equipamentos
Consumo de utilidadesmatérias primas e insumos
Modelo Econômico
avaliaçãoLucro
No caso de processos químicos, a Análise consiste em prever e avaliar o desempenho de cada fluxograma gerado na Síntese, para fins de comparação
![Page 170: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/170.jpg)
1.6.4 Análise
Genericamente: análise significa
- decompor um todo em suas partes,
- depreender o comportamento do todo a partir do comportamento das partes.
PROJETO = SÍNTESE ANÁLISE
![Page 171: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/171.jpg)
W6
T6
W10 T10
W13 T13 W11
T11
W8
T8
W1
x11
T1
f11
f31
W7 T7
W5 T5
W3 x13
T3 f13 f23
W4 x14
T4 f14 f24
W12 T12
W12 T12
W14 T14
W2
x12
T2 f12 f32
EXTRATOR
Extrato
Rafinado
EVAPORADOR
CONDENSADORRESFRIADORMISTURADOR
BOMBA
1
2
3
4
5
67
8
9
10
11
12
13
14
15
Vd Ae
AcAr
Alimentação
Vapor
ÁguaÁgua
Benzeno
Benzeno
Produto
Condensado
W15 T15
Para analisar o Processo
![Page 172: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/172.jpg)
W10 T10
W13 T13
W12 T12
RESFRIADOR
10
11
12
13
Ar
Água
W13 T13
W8
T8
W5 T5
W12 T12
CONDENSADOR
58
9
Ac
Água
W10 T10
10
Benzeno
W6
T6
W7 T7
W3 x13
T3 f13 f23
W4 x14
T4 f14 f24
EVAPORADOR
4
67
Ae
Vapor
W5 T55
Benzeno
Produto
Condensado
3
Extrato
W1
x11
T1
f11
f31
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
Vd
W3 x13
T3 f13 f23
W15 T15
Fragmentando o Processo
W14 T14
MISTURADOR
14
15
Benzeno
W15 T15
13
![Page 173: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/173.jpg)
W10 T10
W13 T13
W12 T12
RESFRIADOR
10
11
12
13
Ar
Água
W13 T13
W1
x11
T1
f11
f31
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
Vd
W3 x13
T3 f13 f23
W15 T15
W8
T8
W5 T5
W12 T12
CONDENSADOR
58
9
Ac
Água
W10 T10
10
Benzeno
W6
T6
W7 T7
W3 x13
T3 f13 f23
W4 x14
T4 f14 f24
EVAPORADOR
4
67
Ae
Vapor
W5 T55
Benzeno
Produto
Condensado
3
Extrato
O Processo fragmentado...
W14 T14
MISTURADOR
14
15
Benzeno
W15 T15
13
![Page 174: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/174.jpg)
Segue a análise de cada equipamento através de problemas de dimensionamento e de simulação
![Page 175: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/175.jpg)
EXTRATOR
![Page 176: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/176.jpg)
Problema proposto: determinar o volume do decantador e a vazão de benzeno necessários para recuperar 60% do ácido benzóico presente a 0,2% nos 100.000 kg/h de alimentação, a 25 oC, com um tempo de residência de 5 min (0,0833 h). Determinar as concentrações das correntes de extrato e de rafinado. A temperatura do benzeno é 25 oC.
W3 x13
T3 f13 f23
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
f11
f31
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
T*15 = 25 oC
*= 0,0833 hr* = 0,60
Vd
W15 W15
Vd
DIMENSIONAMENTO DE UM EXTRATOR
Por enquanto, o extrator é apenas uma figura (ainda não existe)
Fluxograma
![Page 177: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/177.jpg)
01. Balanço Material do Ácido Benzóico: f11 - f12 - f13 = 002. Balanço Material do Benzeno: W15 - f23 = 003. Balanço Material da Água: f31 - f32 = 004. Relação de Equilíbrio Líquido-Líquido: x13 - k x12 = 005. Relação de Equilíbrio Líquido-Líquido: k – (3 + 0,04 Td) = 006. Balanço de Energia: (f11 Cp1 + f31 Cp3) (T1 - Td) + W15 Cp2l (T15 - Td) = 007. Equação de Dimensionamento: Vd - (f11 / 1 + W15 / 2 + f31 / 3) = 008. Fração Recuperada de Ácido Benzóico: r - f13 / f11 = 009. Fases em Equilíbrio T2 – Td = 010. Fases em Equilíbrio T3 – Td = 0
MODELO
W1
x11
T1
f11
f31
1
15
Alimentação
Extrato3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
Vd
W3 x13
T3 f13 f23
W15 T15
34. Vazão Total na Corrente 1: f11 + f31 - W1 = 035. Fração Mássica na Corrente 1: x11 - f11 / W1 = 036. Vazão Total na Corrente 2: f12 + f32 – W2 = 037. Fração Mássica na Corrente 2: x12 - f12 / W2 = 038. Vazão Total na Corrente 3: f13 + f23 – W3 = 039. Fração Mássica na Corrente 3: x13 - f13 / W3 = 0
![Page 178: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/178.jpg)
Balanço de InformaçãoV = 22N = 16C = 4G = 2 !
Metas de Projeto Máximo = 2
V = 22N = 16C = 4M = 2G = 0
W3 x13
T3 f13 f23
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
f11
f31
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
T*15 = 25 oC
Vd
W15 W15
Vd
*= 0,0833 hr* = 0,60
![Page 179: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/179.jpg)
Ordenação das Equações pelo Algoritmo A.O.E.
Quando ocorre um ciclo, o programa busca uma equação final de cima para baixo e toma a primeira que encontra (no caso, a eq 02)
![Page 180: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/180.jpg)
Montagem da rotina Dimensionar Extrator do programa BenzoDSO (ciclo eliminado)
f11 = x11 * W1 '35f13 = r * f11 '08f12 = f11 - f13 '01f31 = W1 - f11 '34f32 = f31 '03W2 = f12 + f32 '36x12 = f12 / W2 '37W = Cp1 * f11 + f31a = (25 / x12) * (W + f13 * Cp2l)b = W * (T15 + 75) - f13 * Cp2l * (T15 + 75 + 25 / x12)c = f13 * Cp2l * (T15 + 75)discr = Sqr(b ^ 2 + 4 * a * c)x13 = (discr + b) / (2 * a) '06' (Variável de abertura)W3 = f13 / x13 '39 (Início do Ciclo)k = x13 / x12 '04f23 = W3 - f13 '38Td = 25 * (k - 3) '05W15 = f23 '02 (Final do Ciclo)Vd = Tau * (f11 / Ro1 + W15 / Ro2 + f31 / Ro3) '07T2 = Td '09T3 = Td '10
Aqui aparece a sequência com a eq. 06 como final e x13 var. de abertura
![Page 181: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/181.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 182: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/182.jpg)
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
W3 x13
T3 f13 f23
T*15 = 25 oC
*= 0,0833 hr* = 0,60
Vd
W2 = 99.880 kg/hx12 = 0,0008
T2 = 25 oCf12 = 80 kg/hf32 = 99.800 kg/h
W3 = 37.490 kg/hx13 = 0,0032
T3 = 25 oCf13 = 120 kg/hf23 = 37.370 kg/h
f11 = 200 kg/hf31 = 99.800 kg/h
W15 = 37.370 kg/hW15
Vd = 11.855l
Extrator dimensionado
Problema proposto: determinar o volume do decantador e a vazão de benzeno necessários para recuperar 60% do ácido benzóico presente a 0,2% nos 100.000 kg/h de alimentação, a 25 oC, com um tempo de residência de 5 min (0,0833 h). Determinar as concentrações das correntes de extrato e de rafinado. A temperatura do benzeno é 25 oC.
Agora o extrator passa a existir, com um volume conhecido.
![Page 183: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/183.jpg)
Uma vez existindo e dimensionado para certas metas, pode-se prever o efeito de alterações nas correntes de entrada
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
1
15
Alimentação
Extrato
3
EXTRATOR
Rafinado
BOMBA
2
T*15 = 25 oC
*= 0,0833 hr* = 0,60
Vd
W2 = 99.880 kg/hx12 = 0,0008T2 = 25 oCf12 = 80 kg/hf32 = 99.800 kg/h
W3 = 37.490 kg/hx13 = 0,0032T3 = 25 oCf13 = 120 kg/hf23 = 37.370 kg/h
f11 = 200 kg/hf31 = 99.800 kg/h
W15 = 37.370 kg/hW15
Vd = 11.855l
Extrator dimensionado
sobre as metas originais e as variáveis de saída
EXEMPLO
SIMULAÇÃO !
![Page 184: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/184.jpg)
Problema proposto: determinar as vazões e as concentrações das correntes de extrato e de rafinado, a fração recuperada de ácido benzóico e o tempo de residência, caso o extrator de Vd = 11.855L fosse alimentado com 50.000 kg/h de benzeno, e não com os 37.370 kg/h de projeto (as demais condições de entrada permanecendo as mesmas de projeto).
SIMULAÇÃO DE UM EXTRATOR
r = =
W*15 = 50.000 kg/h
W3 x13
T3 f13 f23
W2
x12
T2 f12 f32
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
1
15
Alimentação
Extrato
3
EXTRATOR
Rafinado
BOMBA
2
V*d = 11.855 l
T*15 =25 oC
extrator para simulação
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
f11
f31
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
W3 x13
T3 f13 f23
T*15 = 25 oC
*= 0,0833 hr* = 0,60
Vd
W2 = 99.880 kg/hx12 = 0,0008
T2 = 25 oCf12 = 80 kg/hf32 = 99.800 kg/h
W3 = 37.490 kg/hx13 = 0,0032
T3 = 25 oCf13 = 120 kg/hf23 = 37.370 kg/h
f11 = 200 kg/hf31 = 99.800 kg/h
W15 = 37.370 kg/hW15
Vd = 11.855l
Extrator dimensionado
![Page 185: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/185.jpg)
01. Balanço Material do Ácido Benzóico: f11 - f12 - f13 = 002. Balanço Material do Benzeno: W15 - f23 = 003. Balanço Material da Água: f31 - f32 = 004. Relação de Equilíbrio Líquido-Líquido: x13 - k x12 = 005. Relação de Equilíbrio Líquido-Líquido: k – (3 + 0,04 Td) = 006. Balanço de Energia: (f11 Cp1 + f31 Cp3) (T1 - Td) + W15 Cp2l (T15 - Td) = 007. Equação de Dimensionamento: Vd - (f11 / 1 + W15 / 2 + f31 / 3) = 008. Fração Recuperada de Ácido Benzóico: r - f13 / f11 = 009. Fases em Equilíbrio T2 – Td = 010. Fases em Equilíbrio T3 – Td = 0
MODELO
W1
x11
T1
f11
f31
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
Vd
W3 x13
T3 f13 f23
W15 T15
34. Vazão Total na Corrente 1: f11 + f31 - W1 = 035. Fração Mássica na Corrente 1: x11 - f11 / W1 = 036. Vazão Total na Corrente 2: f12 + f32 – W2 = 037. Fração Mássica na Corrente 2: x12 - f12 / W2 = 038. Vazão Total na Corrente 3: f13 + f23 – W3 = 039. Fração Mássica na Corrente 3: x13 - f13 / W3 = 0
![Page 186: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/186.jpg)
Ordenação das Equações pelo Algoritmo AOE
Quando ocorre um ciclo, o programa examina as equações de cima para baixo e toma a primeira candidata como eq. Final (no caso, a eq 04)
![Page 187: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/187.jpg)
f23 = W15 '02f11 = W1 * x11 '35f31 = W1 - f11 '34f32 = f31 '03a = f11 * Cp1 + f31 * Cp3b = W15 * Cp2lTd = (a * T1 + b * T15) / (a + b) '06Tau = Vd / (f11 / Ro1 + W15 / Ro2 + f31 / Ro3) '07k = 3 + 0.04 * Td '05T2 = Td '09T3 = Td '10a = k - 1: 'Cells(24, 7) = ab = k * (f11 + f23) + f32 - f11: 'Cells(25, 7) = bc = f11 * f32: Cells(26, 7) = cdiscr = Sqr(b ^ 2 - 4 * a * c): 'Cells(27, 7) = discrf12 = (b - discr) / (2 * a) '04' Variável de Aberturaf13 = f11 - f12 '01 Início do CicloW2 = f12 + f32 '36x12 = f12 / W2 '37W3 = f13 + f23 '38x13 = f13 / W3 '39 Final de Ciclor = f13 / f11 '08
Montagem da rotina Simular Extrator do programa BenzoDSO (ciclo eliminado)
Aqui aparece a sequência com a eq. 04 como final e f12 var. de
abertura
![Page 188: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/188.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 189: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/189.jpg)
Problema proposto: determinar as vazões e as concentrações das correntes de extrato e de rafinado, a fração recuperada de ácido benzóico e o tempo de residência, caso o extrator de Vd = 11.855L fosse alimentado com 50.000 kg/h de benzeno, e não com os 37.370 kg/h de projeto (as demais condições de entrada permanecendo as mesmas de projeto).
SIMULAÇÃO DO EXTRATOR
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
f11
f31
1
15
Alimentação
Extrato
3
W2
x12
T2 f12 f32
EXTRATOR
Rafinado
BOMBA
2
W3 x13
T3 f13 f23
T*15 = 25 oC
*= 0,0833 hr* = 0,60
Vd
W2 = 99.880 kg/hx12 = 0,0008
T2 = 25 oCf12 = 80 kg/hf32 = 99.800 kg/h
W3 = 37.490 kg/hx13 = 0,0032
T3 = 25 oCf13 = 120 kg/hf23 = 37.370 kg/h
f11 = 200 kg/hf31 = 99.800 kg/h
W15 = 37.370 kg/hW15
Vd = 11.855l
Extrator dimensionado
W*1= 100.000 kg/h
x*11 = 0,002
T*1 = 25 oC
1
15
Alimentação
Extrato
3
EXTRATOR
Rafinado
BOMBA
2
V*d = 11.855 l
W*15 = 50.000 kg/h
T*15 =25 oC
r = 0,67
W3 = 50.133 kg/hx13 = 0,0026
T3 = 25 oCf13 = 133 kg/hf23 = 50.000 kg/h
W2 = 99.867 kg/hx12 = 0,0007
T2 = 25 oCf12 = 67 kg/hf32 = 99.800 kg/h
= 0,075 h
Extrator após Simulação
![Page 190: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/190.jpg)
RESFRIADOR
![Page 191: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/191.jpg)
Problema proposto: determinar a vazão de água de resfriamento e a área de troca térmica do resfriador necessárias para resfriar 36.345 kg/h de benzeno liquido saturado até 25 oC. A água se encontra a 15 oC e deve sair a 30 oC.
DIMENSIONAMENTO DE UM RESFRIADOR
Por enquanto, o resfriador é apenas uma figura (ainda não existe)
Fluxograma
W13 = T*
13 = 25 oC
W*10 = 36.345 kg/h
T*10 = 80 oC
W12 = T*
12 = 30 oC
10
11
12
13Ar
Água
W11 = T*
11 = 15 oC
![Page 192: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/192.jpg)
26. Balanço Material da Água: W11 - W12 = 027. Balanço Material do Benzeno: W10 - W13 = 028. Balanço de Energia na Corrente de Água: Qr - W11 Cp3 (T12 - T11) = 029. Balanço de Energia na Corrente de Benzeno: Qr - W10 Cp2l (T10 - T13) = 030. Equação de Dimensionamento: Qr - Ur Ar r = 031. Definição do T Médio Logarítmico (r ): r - [(T10 - T12) - (T13 - T11) ] / ln[(T10 - T12) / (T13 - T11)] = 0
MODELO
W10 T10
W13 T13
W12 T12
10
11
12
13
Ar
Água
W13 T13
![Page 193: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/193.jpg)
W13 = T*
13 = 25 oC
W*10 = 36.345 kg/h
T*10 = 80 oC
W12 = T*
12 = 30 oC
10
11
12
13Ar
Água
W11 = T*
11 = 15 oC
Balanço de Informação
V = 11N = 6C = 3M = 2G = 0 !
![Page 194: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/194.jpg)
Ordenação das Equações pelo Algoritmo AOE
![Page 195: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/195.jpg)
27. W13 = W10
29. Qr = W10 Cp2l (T10 - T13)28. W11 = Qr / (Cp3 (T12 - T11))26. W12 = W11
d1 = T10 - T12: d2 = T13 - T11
31. dr = (d1 - d2) / ln (d1 / d2)30. Ar = Qr / (Ur dr )
Montagem da rotina Dimensionar Resfriador do programa BenzoDSO
![Page 196: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/196.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 197: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/197.jpg)
W*10 = 36.345 kg/h
T*10 = 80 oC
W12 = 59.969 kg/hT*
12 = 30 oC
10
11
12
13Ar = 362 m2
Água
W11 = 59.969 kg/hT*
11 = 15 oC
W13 = 36.345 kg/hT*
13 = 25 oC
Agora o resfriador passa a existir, com uma área conhecida
Resfriador dimensionado
Problema proposto: determinar a vazão de água de resfriamento e a área de troca térmica do resfriador necessárias para resfriar 36.345 kg/h de benzeno liquido saturado até 25 oC. A água se encontra a 15 oC e deve sair a 30 oC.
![Page 198: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/198.jpg)
Uma vez existindo e dimensionado para certas metas, pode-se prever o efeito de alterações nas correntes de entrada
sobre as metas originais e as variáveis de saída
W*10 = 36.345 kg/h
T*10 = 80 oC
W12 = 59.969 kg/hT*
12 = 30 oC
10
11
12
13Ar = 362 m2
Água
W11 = 59.969 kg/hT*
11 = 15 oC
W13 = 36.345 kg/hT*
13 = 25 oC
EXEMPLO
SIMULAÇÃO !
![Page 199: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/199.jpg)
SIMULAÇÃO DO RESFRIADOR
T*10 = 80 oC
W12 T12
10
11
12
13A*
r = 362 m2
Água
W*11 = 59.969 kg/hT*11 = 15 oC
W13
T13
W*10 = 20.000 kg/h
Problema proposto: pretende-se determinar as temperaturas de saída do benzeno e da água, caso o resfriador projetado para 362 m2 veha a ser alimentado com 20.000 kg/h de benzeno ao invés de 36.345 kg/h, mantidas a vazão e a temperatura da água de resfriamento.
W*10 = 36.345 kg/h
T*10 = 80 oC
W12 = 59.969 kg/hT*
12 = 30 oC
10
11
12
13Ar = 362 m2
Água
W11 = 59.969 kg/hT*
11 = 15 oC
W13 = 36.345 kg/hT*
13 = 25 oC
resfriador dimensionado resfriador para simulação
![Page 200: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/200.jpg)
26. Balanço Material da Água: W11 - W12 = 027. Balanço Material do Benzeno: W10 - W13 = 028. Balanço de Energia na Corrente de Água: Qr - W11 Cp3 (T12 - T11) = 029. Balanço de Energia na Corrente de Benzeno: Qr - W10 Cp2l (T10 - T13) = 030. Equação de Dimensionamento: Qr - Ur Ar r = 031. Definição do T Médio Logarítmico (r ): r - [(T10 - T12) - (T13 - T11) ] / ln[(T10 - T12) / (T13 - T11)] = 0
MODELO
W10 T10
W13 T13
W12 T12
10
11
12
13
Ar
Água
W13 T13
![Page 201: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/201.jpg)
Ordenação das Equações pelo Algoritmo AOE
Logo de saída são ordenadas as duas EIU
As 4 demais equações formam um ciclo.Qualquer uma pode ser escolhida como final.
Por experiência, a única que permite e explicitação da variável de abertura é a 30.
![Page 202: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/202.jpg)
Ordenação das Equações pelo Algoritmo AOE
Este é o resultado apresentado pelo AOE com a equação 28 como final e r como variável abertura.
![Page 203: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/203.jpg)
Esta é a sequência com a equação 30 como final e Qr como variável de abertura.
![Page 204: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/204.jpg)
Resultou a rotina SimularResfriador
W12 = W11 '26W13 = W10 '27a =T10 – T11: b = 1/(W11*cp3) : c = 1/(W10*cp2l) : d = 1/(Ur * Ar) : e = exp ((c-d)/b)Qr = a*(e-1)/(c*e-b) ‘(30’) Var. Abert.T12 = T11 + Qr * a2 '28 T13 = T10 - Qr * a1 ‘29d1 = T10 - T12: d2 = T13 - T11If Abs(d1 - d2) < 0.00001 Then dr = d1 Else dr = (d1 - d2) / Log(d1 / d2) '31
![Page 205: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/205.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 206: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/206.jpg)
SIMULAÇÃO DO RESFRIADOR
Problema proposto: pretende-se determinar as temperaturas de saída do benzeno e da água, caso o resfriador projetado para 361 m2 fosse alimentado com 20.000 kg/h de benzeno ao invés de 36.345 kg/h, mantidas a vazão e a temperatura da água de resfriamento.
W*10 = 36.345 kg/h
T*10 = 80 oC
W12 = 59.969 kg/hT*
12 = 30 oC
10
11
12
13Ar = 362 m2
Água
W11 = 59.969 kg/hT*
11 = 15 oC
W13 = 36.345 kg/hT*
13 = 25 oC
W*10 = 20.000 kg/h
T*10 = 80 oC
W12 = 59.969 kg/hT12 = 24,5 oC
10
11
12
13A*
r = 362 m2
Água
W*11 = 59.969 kg/hT*11 = 15 oC
W13 = 20.000 kg/hT13 = 16,8 oC
Resultado do dimensionamento Resultado da simulação
![Page 207: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/207.jpg)
20. Balanço Material da Água: W8 - W9 = 021. Balanço Material do Benzeno: W5 - W10 = 022. Balanço de Energia na Corrente de Água: Qc - W8 Cp3 (T9 - T8) = 023. Balanço de Energia na Corrente de Benzeno: W5 [2 + Cp2g (T5 – T10)] - Qc = 024. Equação de Dimensionamento: Qc - Uc Ac c = 025. Definição do T Médio Logarítmico (c): c - [(T5 - T9) - (T10 - T8)]/ln[(T5 - T9)/(T10 - T8)] = 0
CONDENSADOR
W5 T5
W10 T10
W9 T9
5
8
9
10
Ar
Água
W8 T8
![Page 208: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/208.jpg)
DIMENSIONAMENTO DO CONDENSADOR
Problema proposto: determinar a vazão de água de resfriamento e a área de troca térmica necessárias para condensar 36.345 kg/h de benzeno de vapor saturado a líquido saturado. A água se encontra a 15 oC e deve sair a 30 oC .
W*5 = 36.345 kg/h
T*5 = 80 oC
W10 T*
10 = 80 oC
W9 T*
9 = 30 oC
5
8
9
10
Água
W8 T*
8 = 15 oC
Ac
W10 = 36.345 kg/hT*
10 = 80 oC
W9 = 228.101 kg/hT*
9 = 30 oC
5
8
9
10
Água
W8 = 228.101 kg/hT*
8 = 15 oC
Ac = 120 m2
W*5 = 36.345 kg/h
T*5 = 80 oC
V = 11N = 6C = 3M = 2G = 0 !
![Page 209: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/209.jpg)
20. Balanço Material da Água: W8 - W9 = 021. Balanço Material do Benzeno: W5 - W10 = 022. Balanço de Energia na Corrente de Água: Qc - W8 Cp3 (T9 - T8) = 023. Balanço de Energia na Corrente de Benzeno: W5 [2 + Cp2g (T5 – T10)] - Qc = 024. Equação de Dimensionamento: Qc - Uc Ac c = 025. Definição do T Médio Logarítmico (c): c - [(T5 - T9) - (T10 - T8)]/ln[(T5 - T9)/(T10 - T8)] = 0
CONDENSADOR
W5 T5
W10 T10
W9 T9
5
8
9
10
Ar
Água
W8 T8
![Page 210: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/210.jpg)
![Page 211: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/211.jpg)
21. W10 = W5
23. Qc = W5 2
d1 = T5 - T9: d2 = T10 - T8
22. W8 = Qc / (Cp3 * (T9 - T8))20. W9 = W8
25. dc = (d1 - d2) / ln (d1 / d2)24. Ac = Qc / (Uc * dc)
Resultando a rotina DimensionarCondensador
![Page 212: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/212.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 213: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/213.jpg)
SIMULAÇÃO DO CONDENSADOR
Problema proposto: determinar a vazão de água necessária para condensar 20.000 kg/h de benzeno, ao invés dos 36.345 kg/h para os quais foi calculada a área de 119 m2. Pretende-se que o benzeno deixe o condensador como líquido saturado. A água se encontra a 15 oC.
W*5 = 20.000 kg/h
T*5 = 80 oC
W10 = T*
10 = 80 oC
W9 = T9 =
5
8
9
10
Água
W8 = 228.101 kg/hT*
8 = 15 oC
A*c = 120 m2
W10 = 36.345 kg/hT*
10 = 80 oC
W9 = 228.101 kg/hT*
9 = 30 oC
5
8
9
10
Água
W8 = 228.101 kg/hT*
8 = 15 oC
Ac = 120 m2
W*5 = 36.345 kg/h
T*5 = 80 oC
resultado do dimensionamento
V = 11N = 6C = 5M = 1G = - 1 !!!
Pretendido na simulação
![Page 214: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/214.jpg)
V = 11N = 6C = 5M = 1G = -1 !!!
Este problema difere da simulação do extrator e do resfriador porque contem uma meta: o benzeno deve sair como líquido saturado, logo a 80 oC.
Daí: G = -1.
W*5 = 20.000 kg/h
T*5 = 80 oC
W10 = T*
10 = 80 oC
W9 = T9 =
5
8
9
10
Água
W8 = 228.101 kg/hT*
8 = 15 oC
A*c = 120 m2
![Page 215: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/215.jpg)
V = 11N = 6C = 5M = 1G = -1 !!!
W*5 = 20.000 kg/h
T*5 = 80 oC
W10 = T*
10 = 80 oC
W9 = T9 =
5
8
9
10
Água
W8 = 228.101 kg/hT*
8 = 15 oC
A*c = 120 m2
O problema só pode ser resolvido se alguma condição conhecida deixar de sê-lo.
Uma solução consiste em transformar W8 em incógnita. O seu valor calculado pode ser tomado como set-point de um sistema de controle que a utilize como variável de controle.
W*5 = 20.000 kg/h
T*5 = 80 oC
W10 = 20.000 kg/hT*
10 = 80 oC
W9 = 35.718 kg/hT9 = 67,7 oC
5
8
9
10
Água
W8 = 35.718 kg/hT*
8 = 15 oC
A*c = 120 m2
V = 11N = 5C = 5M = 1G = 0
![Page 216: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/216.jpg)
W10 = 36.345 kg/hT*
10 = 80 oC
W9 = 228.101 kg/hT*
9 = 30 oC
5
8
9
10
Água
W8 = 228.101 kg/hT*
8 = 15 oC
Ac = 120 m2
W*5 = 36.345 kg/h
T*5 = 80 oC
resultado do dimensionamento
W*5 = 20.000 kg/h
T*5 = 80 oC
W10 = 20.000 kg/hT*
10 = 80 oC
W9 = 35.727 kg/hT9 = 67,7 oC
5
8
9
10
Água
W8 = 35.727 kg/hT*
8 = 15 oC
A*c = 120 m2
resultado da simulação
![Page 217: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/217.jpg)
![Page 218: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/218.jpg)
21. W10 = W5
23. Qc = W5 2
24. dc = Qc / (Uc Ac)25. c - [(T5 - T9) - (T10 - T8)]/ln[(T5 - T9)/(T10 - T8)] = 0 (Resolução por Bisseção)22. W8 = Qc / (Cp3 (T9 - T8))20. W9 = W8
Resulta a rotina SimularCondensador
![Page 219: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/219.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 220: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/220.jpg)
11. Balanço Material do Ácido Benzóico: f13 - f14 = 012. Balanço Material do Benzeno: f23 - f24 - W5 = 013. Balanço Material do Vapor: W6 - W7 = 014. Balanço de Energia na Corrente de Vapor: W6 [3 + Cpv (T6 – T7)] - Qe = 015. Balanço de Energia na Corrente de Processo: Qe – [(f13Cp1 + f23Cp2l)(Te - T3) + W5 2] = 016. Equação de Dimensionamento: Qe - Ue Ae e = 017. Definição da Diferença de Temperatura (e): e - (T6 - Te) = 018. Fases em Equilíbrio T4 – Te = 019. Fases em Equilíbrio T5 – Te = 0
EVAPORADOR
W6
T6
W7 T7
W3 x13
T3 f13 f23
W4 x14
T4 f14 f24
4
67
Ae
Vapor
W5 T55
Benzeno
Produto
Condensado
3
Extrato
38. Vazão Total na Corrente 3: f13 + f23 – W3 = 039. Fração Mássica na Corrente 3: x13 - f13 /W3 = 040. Vazão Total na Corrente 4: f14 + f24 - W4 = 041. Fração Mássica na Corrente 4: x14 - f14/W4 = 0
![Page 221: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/221.jpg)
DIMENSIONAMENTO DO EVAPORADOR
Problema proposto: determinar a vazão de um vapor a 150 oC e a área de troca térmica necessárias para obter um concentrado com 10% de ácido benzóico, a partir de uma corrente com 37.545 kg/h de uma solução de 0,32% de ácido benzóico em
benzeno, a 25 oC. O condensado deve sair como líquido saturado a 150 oC . O evaporador opera a 1 atm.
W6
T*6 = 150 oC
W7 T*
7 = 150 oC
W*3 = 37.345 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 f23
W4 x*
14 = 0,10
T4 f14 f24
4
67
Ae
Vapor
W5 T55
Benzeno
Produto
Condensado
3 Te* = 80 oC
W6 = 8.569 kg/hT*
6 = 150 oCW7 = 8.569 kg/hT*
7 = 150 oC
W*3 = 37.545 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 = 120 kg/h f23 = 37.225 kg/h
W4 = 1.195 kg/h x*
14 = 0,10T4 = 80 oCf14 = 120 kg/h f24 = 1.076kg/h
4
67
Ae= 124 m2
Vapor
W5 = 36.150 kg/h T5 = 80 oC
5Benzeno
Produto
Condensado
3 Te* = 80 oC
V = 20N = 13C = 4M = 3G = 0 !
![Page 222: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/222.jpg)
![Page 223: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/223.jpg)
15. De = T6 - T35. f13 = W3 x13
09. f14 = f13
34. f23 = W3 - f13
37. W4 = f14 / x14
36. f24 = W4 - f14
10. W5 = f23 - f24
13. Qe = (f13 Cp1 + f23 Cp2l) (T - T3) + W5 L2
12. W6 = Qe / (L 3 + Cp3 (T6 - T7))11. W7 = W6
14. Ae = Qe / (Ue De)
Resulta a rotina DimensionarEvaporador
![Page 224: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/224.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 225: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/225.jpg)
SIMULAÇÃO DO EVAPORADOR
Problema proposto: determinar as vazões de vapor e de evaporado, a vazão e a concentração do concentrado, caso o evaporador, com os mesmos 124 m2 de área de projeto, fosse alimentado com 50.000 kg/h de solução e não mais com 37.545 kg/h. Pretende-se que o condensado saia como líquido saturado (150 oC).
V = 20N = 13C = 7M = 1G = -1 !!!
pretendido na simulaçãoresultado do dimensionamento
W6 = 8.615 kg/hT*
6 = 150 oCW7 = 8.615 kg/hT*
7 = 150 oC
W*3 = 37.545 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 = 120 kg/h f23 = 37.425 kg/h
W4 = 1.201 kg/h x*
14 = 0,10T4 = 80 oCf14 = 120 kg/h f24 = 1.081kg/h
4
67
Ae=124 m2
Vapor
W5 = 36.344 kg/h T5 = 80 oC
5Benzeno
Produto
Condensado
3 Te* = 80 oC
W6 = 8.615 kg/hT*
6 = 150 oCW7 = T*
7 = 150 oC
W*3 = 50.000 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 = 120 kg/h f23 = 37.425 kg/h
W4 = x14 = T4 =f14 = f24 =
4
67
Ae=124 m2
Vapor
W5 = T5 =5
Benzeno
Produto
Condensado
3 Te* = 80 oC
![Page 226: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/226.jpg)
V = 20N = 13C = 7M = 1G = -1 !!!
W5 = T5 =
W6 = 8.615 kg/hT*
6 = 150 oCW7 = T*
7 = 150 oC
W*3 = 50.000 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 = 120 kg/h f23 = 37.425 kg/h
W4 = x14 = T4 =f14 = f24 =
4
67
Ae=124 m2
Vapor
5Benzeno
Produto
Condensado
3 Te* = 80 oC
Uma solução consiste em transformar W6 em incógnita. O seu valor calculado pode ser tomado como set-point de um sistema de controle que a utilize como variável de controle.
Situação semelhante à da simulação do condensador
V = 20N = 13C = 6M = 1G = 0
W6 = 8.569 kg/hT*
6 = 150 oCW7 = 8.569 kg/hT*
7 = 150 oC
W*3 = 50.000 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 = 160 kg/h f23 = 49.840 kg/h
W4 = 17.177 kg/h x14 = 0,0093T4 = 80 oCf14 = 160 kg/h f24 = 17.017 kg/h
4
67 Vapor
W5 = 32.823 kg/h T5 = 80 oC
5Benzeno
Produto
Condensado
3 Te* = 80 oC
Ae=124 m2
![Page 227: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/227.jpg)
resultado do dimensionamento
W6 = 8.569 kg/hT*
6 = 150 oCW7 = 8.569 kg/hT*
7 = 150 oC
W*3 = 37.545 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 = 120 kg/h f23 = 37.425 kg/h
W4 = 1.201 kg/h x*
14 = 0,10T4 = 80 oCf14 = 120 kg/h f24 = 1.081kg/h
4
67
Ae=124 m2
Vapor
W5 = 36.344 kg/h T5 = 80 oC
5Benzeno
Produto
Condensado
3 Te* = 80 oC
W6 = 8.569 kg/hT*
6 = 150 oCW7 = 8.569 kg/hT*
7 = 150 oC
W*3 = 50.000 kg/h
x*13 = 0,0032
T*3 = 25 oC
f13 = 160 kg/h f23 = 49.840 kg/h
W4 = 17.177 kg/h x14 = 0,0093T4 = 80 oCf14 = 160 kg/h f24 = 17.017 kg/h
4
67 Vapor
W5 = 32.823 kg/h T5 = 80 oC
5Benzeno
Produto
Condensado
3 Te* = 80 oC
Ae=124 m2
resultado da simulação
![Page 228: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/228.jpg)
![Page 229: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/229.jpg)
15. De = T6 - T14. Qe = Ue Ae De 12. W6 = Qe / (l3 + Cpv * (T6 - T7))11.W7 = W6
35. f13 = W3 * x13
09. f14 = f13
34. f23 = W3 - f13
13. W5 = (Qe - (f13 * Cp1 + f23 * Cp2l) * (T - T3)) / l210. f24 = f23 - W536. W4 = f14 + f24
37. x14 = f14 / W4
Resulta a rotina SimularEvaporador
![Page 230: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/230.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 231: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/231.jpg)
3.1 Equações Não-Lineares 3.1.1 Representação 3.1.2 Métodos Numéricos3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e representação 3.2.2 Resolução: partição, abertura, ordenação de equações3.3 Dimensionamento e Simulação de Equipamentos
3.4.2 Estratégia Modular3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Dimensionamento 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global
![Page 232: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/232.jpg)
3.4 DIMENSIONAMENTO E SIMULAÇÃO DE PROCESSOS
Todas as equações são consideradas simultaneamente, independentemente dos equipamentos a que pertencem.
É a estratégia mais indicada para dimensionamento.
O Algoritmo de Ordenação de Equações é executado como se fosse para um equipamento isolado.
3.4.1 ESTRATÉGIA GLOBAL
Existem duas estratégias básicas:
- Estratégia Global- Estratégia Modular
![Page 233: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/233.jpg)
Dimensionamento do Processo – Estratégia Global
01. f11 - f12 - f13 = 002. W15 - f23 = 003. f31 - f32 = 004. k – (3 + 0,04 Td) = 005. k - x13 / x12= 006. (f11 Cp1 + f31 Cp3) (T1 - Td) + W15 Cp2l (T15 - Td) = 007. Vd - (f11 /1 + W15/2 + f31/3) = 008. r - f13/f11 = 009. T2 – Td = 010. T3 – Td = 0
11. f13 - f14 = 012. f23 - f24 - W5 = 013. W6 - W7 = 014. W6 [3 + Cpv (T6 – T7)] - Qe = 015. Qe – [(f13Cp1 + f23Cp2l)(Te - T3) + W5 2] = 016. Qe - Ue Ae e = 017. e - (T6- Te) = 018. T4 – Te = 019. T5 – Te = 0
20. W8 - W9 = 021. W5 - W10 = 022. Qc - W8 Cp3 (T9 - T8) = 023. W5 [2 + Cp2g (T5 – T10)] - Qc = 024. Qc - Uc Ac c = 025. c - [(T5 - T9) - (T10 - T8)]/ln[(T5 - T9)/(T10 - T8)] = 0
26. W11 - W12 = 027. W10 - W13 = 028. Qr - W11 Cp3 (T12 - T11) = 029. Qr - W10 Cp2l (T10 - T13) = 030. Qr - Ur Ar r = 031. r - [(T10 - T12) - (T13 - T11)]/ln[(T10 - T12)/(T13 - T11)] = 0
32. W13 + W14 - W15 = 033. W13 (T15 - T13) + W14 (T15 - T14) = 0
34. f11 + f31 - W1 = 035. x11 - f11 /W1 = 036. f12 + f22 – W2 = 037. x12 - f12/W2 = 038. f13 + f23 – W3 = 039. x13 - f13 /W3 = 040. f14 + f24 - W4 = 041. x14 - f14/W4 = 0
![Page 234: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/234.jpg)
f11 f12 f13 E V
1 * O * W15f23 5 k
2 O * f31 f32 9 T2
3 * O k 10 T3
4 * * O * * Td 17 De
5 O X T1T15 18 T4
6 * * X X O Vd 19 T5
7 * * * O X r 25 dc
8 * O X T2 35 f11
9 X O T3 8 f13
10 X O f14 1 f12
11 * O f24W5 11 f14
12 * * O W6W7 34 f31
13 * O T6 T7 Qe 3 f32
14 O X X * Te 4 f23
15 * * * * O X Ae
e 2 W15
16 * O * 6 T15
17 X X O T4 7 Vd
18 X O T5 36 W2
19 X O W8W9 37 x12
20 * O W10 38 W3
21 * O Qc T9 T8 39 x13
22 O * X X T10 41 W4
23 * * O X Ac
c 40 f24
24 * O * 12 W5
25 * X X X O W11W12 15 Qe
26 * O W13 14 W6
27 * O Qr T11T12 13 W7
28 O * X X T13 16 Ae
29 * X O * A r
r 21 W10
30 * O * 23 Qc
31 X X X * O W14T14 22 W8
32 * * O 20 W9
33 * * O * X W1 24 Ac
34 * O X x11 27 W13
35 O X X W2 32 W14
36 * * O x12 33 T14
37 * * O W3 31 dr
38 * * O x13 29 Qr
39 * * O W4 28 W11
40 * O * x14 26 W1241 * O X 30 Ar
Extrator
Evaporador
Correntes Multicomponentes
Condensador
Resfriador
Misturador
![Page 235: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/235.jpg)
Dimensionar Processo
(03) T3 = T2
(13) T4 = T5 (16) e = T6 - T5 (22) D1 = T5 - T9: D2 = T10 - T8 : c = (D1 - D2) / ln (D1 / D2)(32) f11 = W1 x11 (08) f13 = f11 r(31) f31 = W1 - f11 (01) f12 = f11 - f13 (09) f14 = f13 (03) f32 = f31 (04) f23 = f13 f32 / (k f12)(34) W4 = f14 / x14 (02) W15 = f23
(33) f24 = W4 - f14 (05) T15 = T2 - (f11 Cp1 + f31 Cp3) (T1 - T2) / (W15 Cp2l) (07) Vd = (f11 / 1 + W15 / 2 + f31 / 3) (10) W5 = f23 - f24 (14) Qe = (f13 Cp1 + f23 Cp2l) (T5 - T3) + W52
![Page 236: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/236.jpg)
(18) W10 = W5 (20) Qc = W5 (2 + Cp2l (T5 - T10)) (12) W6 = Qe / ( 3 + Cp3 (T6 - T7)) (15) Ae = Qe / (Ue e) (24) W13 = W10 (19) W8 = Qc / (Cp3 (T9 - T8)) (21) Ac = Qc / (Uc c) (11) W7 = W6 (29) W14 = W15 - W13 (17) W9 = W8 (30) T13 = T15 + W14 (T15 - T14) / W13 (26) Qr = W10 Cp2l (T10 - T13) (28) D1 = T10 - T12: D2 = T13 - T11: r = (D1 - D2) / ln (D1 / D2)(25) W11 = Qr / (Cp3 (T12 - T11)) (27) Ar = Qr / (Ur r) (23) W12 = W11
![Page 237: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/237.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 238: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/238.jpg)
3.1 Equações Não-Lineares 3.1.1 Representação 3.1.2 Métodos Numéricos3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e representação 3.2.2 Resolução: partição, abertura, ordenação de equações3.3 Dimensionamento e Simulação de Equipamentos3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global
3.5 Incerteza e Análise de Sensibilidade 3.5.1 Questionamento do Dimensionamento 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.4.2 Estratégia Modular
![Page 239: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/239.jpg)
3.4.2 Estratégia Modular
Para cada problema, os módulos são seqüenciados convenientemente segundo o fluxograma material do
processo.
Havendo a presença de reciclos no fluxograma, torna-se necessária a abertura de um certo número de correntes e a
inserção de um módulo promotor de convergência para cada uma.
É a estratégia mais indicada para simulação.
Utiliza módulos criados previamente para cada equipamento.Cada módulo contem as equações já ordenadas para
dimensionamento ou simulação (Seção 3.3).
![Page 240: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/240.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 241: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/241.jpg)
Simulação do Processo IlustrativoEstratégia Modular
![Page 242: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/242.jpg)
W6 =8.594 kg/hT*
6 = 150 oC
W10 =36.284 kg/hT*
10 = 80 oCW13 = 36.284 kg/hT13 = 25 oC
W11 = 59.969 kg/hT*
11 = 15 oCW8 = 232.603 kg/hT*
8 = 15 oC
W*1 = 150.000 kg/h
x*11 = 0,002
T*1 = 25 oC
f11 = 300 kg/hf31 = 149.700 kg/h
W7 = 8.594 kg/hT*
7 = 150 oC
W5 = 36.284 kg/hT*
5 = 80 oC
W3 = 37.477 kg/hx13 = 0,004
T3 = 25 oCf13 = 149 kg/hf23 = 37.328 kg/h
W4 = 1.130 kg/hx14 = 0,12
T4 = 80 oCf14 = 150 kg/hf24 = 1.080 kg/h
W12 = 59.969 kg/hT12 = 29 oC
W12 = 232.603 kg/hT12 = 29 oC
W*14 = 1.080 kg/h
T*14 = 25 oC
W2 = 149.850 kg/hx12 = 0,001
T2 = 25 oCf12 = 150 kg/hf32 = 149.700 kg/h
EXTRATOR
Extrato
Rafinado
EVAPORADOR
CONDENSADORRESFRIADORMISTURADOR
BOMBA
1
2
3
4
5
67
8
9
10
11
12
13
14
15
V*d = 11.859 l
= 0,0617 h
r = 0,50
A*e =
124 m2
A*c = 119 m2
A*r = 361 m2
W15 = 37.328 kg/hT13 = 25 oC
O fluxograma exibe um reciclo.
A cada iteração o módulo confere a convergência e atualiza o valor de W5
O valor inicial arbitrado para W5 pode ser aquele obtido noDimensionamento.
Implementa-se um módulo promotor de convergência: no caso, o deSubstituição Direta.
Seleciona-se uma corrente de abertura com o menor número possívelde variáveis (simplificar o gerenciamento da convergência): no caso, foi selecionada a corrente 5 (é preciso gerenciar apenas W5).
![Page 243: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/243.jpg)
Simulação do Processo Ilustrativo - Estratégia Modular
EXTRATOR
RESFRIADOR
MISTURADOR
CONDENSADOR
EVAPORADOR
SS
18. W10
20. Qc
19. c
22'. T9
21. W8
17. W9
24. W13
23. W12
25'. Qr
28. T13
27. T12
26. r
29. W15
30. T15
02. f23
32. f11
31. f31
03. f32
05. T2
07. 06. T3
01' f12
04. f13
08. r
W1
T1
x11
f11
f31
W15
T15
W45
T14
W13
T13
W10
T10
f13
f23
T3
W4
T4
x14
f14
f24
09. f14
13. T4
16. e
15. Qe
12. W6
14. W5
10. f24
11. W7
33. W4
34. x14
T5
T2
f12
f32
W5a
W5c
Repetição até convergir
|W5c – W5a| / W5a
erro relativo
![Page 244: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/244.jpg)
SUB SimularOProcesso'----------------------------------------------------------------------------INPUT "W5= "; W5cW5$ = "W5 = " + STR$(INT(W5c))NoDeIteracoes = 0DO W5a = W5c SimularOCondensador SimularOResfriador SimularOMisturador SimularOExtrator SimularOEvaporador MostrarOResultado NoDeIteracoes = NoDeIteracoes + 1 ErroRelativo = ABS(W5a - W5c) / W5a PausaSeQuizerLOOP UNTIL ConvergirEND SUB
PRINCIPAL
Simular
Simular
Simular
Simular
Simular
Simular
Dimensionar
Dimensionar
Dimensionar
Dimensionar
Dimensionar
DimensionarExtratorExtrator
EvaporadorEvaporador
Condensador Condensador
ResfriadorResfriador
MisturadorMisturador
ProcessoProcesso
Mostrar Lucro doCalcular Lucro doEmpreendimento
Processo
Resolver Problema
Otimizar
Empreendimento
![Page 245: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/245.jpg)
Resolver Problema
Otimizar Processo
Calcular Lucro
DimensionarExtrator
DimensionarEvaporador
DimensionarCondensador
DimensionarResfriador
DimensionarMisturador
SimularExtrator
SimularEvaporador
SimularCondensador
SimularResfriador
SimularMisturador
SimularProcesso
DimensionarProcesso
![Page 246: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/246.jpg)
Simulação de Processos com Estrutura Complexa
1 2 3 4 5 6 7 81* 2 3
4
5 6
7
8
9
10 11
12
13
14
Problema:Estabelecer uma estratégia de cálculo e implementá-la sob a forma de um algoritmo executável em computador.
Cada equipamento é representado por um módulo computacional em que as equações se encontram ordenadas para simulação.
A estratégia de cálculo é a ordem em que os equipamentos devem ser simulados.
![Page 247: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/247.jpg)
Simulação de Processos com Estrutura Complexa
1 2 3 4 5 6 7 81* 2 3
4
5 6
7
8
9
10 11
12
13
14
Procedimento:(a) identificação dos ciclos.(b) seleção das correntes de abertura(c ) construção do algoritmo de simulação
Dificuldade: os diversos reciclos
![Page 248: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/248.jpg)
(a) Identificação dos Ciclos
Pode-se utilizar o Método do Traçado de Percursos (labirinto)
1 2 3 4 5 6 7 81* 2 3
4
5 6
7
8
9
10 11
12
13
14
Percorre-se o fluxograma anotando, numa lista dupla, as correntes (LC) e os equipamentos visitados (LE).
Corrente: 1 2 3 4 Destino : 1 2 3 1
Um ciclo é identificado ao se chegar a um equipamento já visitado.
Equipamento 1 já visitado : ciclo 2 3 4
![Page 249: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/249.jpg)
(a) Identificação dos Ciclos
1 2 3 4 5 6 7 81* 2 3
4
5 6
7
8
9
10 11
12
13
14
ALGORITMO RESUMIDO
Colocar uma corrente conhecida na LC (toda corrente colocada fica “aberta”)Colocar o seu destino na LEREPETIR O destino já está na LE: registrar o ciclo formado por todas as correntes abertas após o primeiro registro do destino) e fechar as correntes até à última aberta; tomar a última corrente aberta e colocar seu destino na LE. O destino não está na LE: colocar as suas correntes de saída na LC; colocar o destino da primeira delas na LE.ATÉ NÃO MAIS HAVER CORRENTES ABERTAS
![Page 250: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/250.jpg)
MATRIZ CICLO - CORRENTE
2 3 4 5 6 7 8 9 10 11 12 131 1 1 12 1 1 13 1 1 1 1 1 14 1 1 1 15 1 1 1 1 1 1 1 6 1 1 1 1 1
Os Ciclos encontrados são registrados na
1 2 3 4 5 6 7 81* 2 3
4
5 6
7
8
9
10 11
12
13
14
![Page 251: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/251.jpg)
APLICAÇÃO AO PROBLEMA ILUSTRATIVO
![Page 252: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/252.jpg)
C: D:
3 4 5 6 711*
22 3
4
5 6
7
8
9
10 11
12
13
814
11
22
33
541
22
33
541
C: 1 2 3 5D: 1 2 3 4
765
86
1110 4
765
86
1110 4
7C: 1 2 3 5 6 8 11 D: 1 2 3 4 5 6 8
13 2
7C: 1 2 3 5 6 8 11 D: 1 2 3 4 5 6 8
13 2
C: 1 2 3 5 7 D: 1 2 3 4 7
1295
86
1110 4
C: 1 2 3 5 7 D: 1 2 3 4 7
1295
86
1110 4
12C: 1 2 3 5 7 9 8 11 D: 1 2 3 4 7 5 6 8
13 2
12C: 1 2 3 5 7 9 8 11 D: 1 2 3 4 7 5 6 8
13 2
C: 1 2 3 5 7 12 D: 1 2 3 4 7 8
13 2
C: 1 2 3 5 7 12 D: 1 2 3 4 7 8
13 2
Colocar uma corrente conhecida na LC (toda corrente colocada fica “aberta”)Colocar o seu destino na LEREPETIR O destino já está na LE: registrar o ciclo formado por todas as correntes abertas após o primeiro registro do destino) e fechar as correntes até à última aberta; tomar a última corrente aberta e colocar seu destino na LE. O destino não está na LE: colocar as suas correntes de saída na LC; colocar o destino da primeira delas na LE.ATÉ NÃO MAIS HAVER CORRENTES ABERTAS
![Page 253: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/253.jpg)
(b) Seleção das Correntes de Abertura
Matriz Ciclo - Corrente
ALGORITMOCalcular os elementos de CRepetir Identificar a corrente com o maior valor em C (pode ser a primeira encontrada) Inscrever a corrente em A Remover os ciclos abertos pela corrente (anular os elementos na linhas correspondentes) Atualizar C Até C = 0
2 3 4 5 6 7 8 9 10 11 12 131 1 1 12 1 1 13 1 1 1 1 1 14 1 1 1 15 1 1 1 1 1 1 1 6 1 1 1 1 1
1 4 1 3 2 3 4 1 2 2 1 3C
000000
A
![Page 254: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/254.jpg)
2 3 4 5 6 7 8 9 10 11 12 131 1 1 12 1 1 13 1 1 1 1 1 14 1 1 1 15 1 1 1 1 1 1 1 6 1 1 1 1 1
1 4 1 3 2 3 4 1 2 2 1 3C
000000
A
2 3 4 5 6 7 8 9 10 11 12 131 0 0 02 1 1 13 0 0 0 0 0 04 1 1 1 15 0 0 0 0 0 0 0 6 0 0 0 0 0
0 0 0 0 1 1 2 1 2 0 0 0C
300000
A
![Page 255: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/255.jpg)
2 3 4 5 6 7 8 9 10 11 12 131 0 0 02 0 0 03 0 0 0 0 0 04 0 0 0 05 0 0 0 0 0 0 0 6 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0C
380000
A
2 3 4 5 6 7 8 9 10 11 12 131 0 0 02 1 1 13 0 0 0 0 0 04 1 1 1 15 0 0 0 0 0 0 0 6 0 0 0 0 0
0 0 0 0 1 1 2 1 2 0 0 0C
300000
A
![Page 256: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/256.jpg)
(c) Construção do Algoritmo de Simulação
1 2 3 4 5 6 7 81* 2
4
5 6
7
8
9
10 11
12
13
143
Abrir C3
REPETIRSimular E3 (C4,C5)Simular E1 (C2)
REPETIRSimular E6 (C10,C11)Simular E4 (C6,C7 )Simular E7 (C9, C12)Simular E5 (C8)
ATÉ Convergir C8
Simular E8 (C13, C14)Simular E2 (C3)
ATÉ Convergir C3
Abrir C8
Corrente 1: única conhecida
![Page 257: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/257.jpg)
3.1 Equações Não-Lineares 3.1.1 Representação 3.1.2 Métodos Numéricos .3.2 Sistemas de Equações Não-Lineares 3.2.1 Estrutura e representação 3.2.2 Resolução: partição, abertura, ordenação de equações3.3 Dimensionamento e Simulação de Equipamentos3.4 Dimensionamento e Simulação de Processos 3.4.1 Estratégia Global 3.4.2 Estratégia Modular
3.5.1 Questionamento do Dimensionamento 3.5.2 Questionamento do Desempenho Futuro
3. ESTRATÉGIAS DE CÁLCULO
3.5 Incerteza e Análise de Sensibilidade
![Page 258: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/258.jpg)
ASSUNTO TRANSFERIDO PARA DEPOIS DO CAPÍTULO 4
![Page 259: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/259.jpg)
3.5 INCERTEZA E ANÁLISE DE SENSIBILIDADE
(a) modelos matemáticos: aproximações lineares, coeficientes constantes...
A Análise de Processos é executada em ambiente de muita incerteza.
A avaliação dos efeitos da incerteza é efetuada através da
(b) parâmetros físicos e econômicos: valores incertos (aproximados e variáveis).
Fontes de incerteza:
Análise de Sensibilidade
![Page 260: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/260.jpg)
(b) questionamento do desempenho futuro:
(a) questionamento do próprio dimensionamento:
A Análise de Sensibilidade avalia os efeitos da incerteza através de dois questionamentos ao final do dimensionamento,
Em que grau a incerteza nos parâmetros compromete o resultado do dimensionamento ?
Em que grau a incerteza nos parâmetros comprometerá as metas de projeto que serviram de base para o
dimensionamento?
![Page 261: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/261.jpg)
Fazem parte da Análise:
- as variáveis características do dimensionamento: dimensões.
- as variáveis características do desempenho do processo: variáveis de saída (metas de projeto).
- os parâmetros cujos valores são considerados incertos (variáveis conhecidas são aqui incorporadas ao conjunto dos parâmetros).
![Page 262: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/262.jpg)
F: variável do processo cujo valor é incerto devido à incerteza nos parâmetros . Exemplo: W3, A.
: vetor dos parâmetros (físicos e econômicos) e das variáveis especificadas cujos valores são incertos. Exemplo: Cp1, Cp3, U, W1, T1, T3.
Fundamento da Análise de Sensibilidade
Exemplo: Trocador de Calor
T1* = 80 oC
W1* = 30.000 kg/h
A = 265,6 m2
T 2* = 25 oC
W3 = 44.000 kg/h
T3* = 15 oC
T4* = 30 oC
0
TT
TTln
)TT()TT(.4
0UAQ.3
0)TT(CpWQ.2
0)TT(CpWQ.1
32
41
3241
3433
2111
![Page 263: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/263.jpg)
F: variável do processo cujo valor é incerto devido à incerteza nos parâmetros . Exemplo: W3, A.
S (F; i): Sensibilidade de F à incerteza no parâmetro i.
: vetor dos parâmetros (físicos e econômicos) e das variáveis especificadas cujos valores são incertos. Exemplo: Cp1, Cp3, U, W1, T1, T3.
Fundamento da Análise de Sensibilidade
i *
F
i
*ii
ii
)(F);F(S
Exemplo:
100
U
i
U
)(A)U;A(S
A Sensibilidade é função do parâmetro
![Page 264: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/264.jpg)
Conveniência: usar variáveis adimensionais F/F* e i / i*
Análise de Sensibilidade com Variáveis Adimensionais
)(F
)(F
)/(
)](F/)(F[)/;F/F(S
*i
*i
i
i*ii
*ii*
ii*
*i
1
Vantagens:(a) os valores independem das dimensões das variáveis e dos parâmetros.(b) as Sensibilidades podem ser comparadas, permitindo verificar a qual parâmetro a variável de interesse é mais sensível, e em que grau.
Nova definição de Sensibilidade:
6265
100
1001,U
)U(A
)U/U(
)]U(A/)U(A[)U/U;A/A(S
U*
***
Exemplo
![Page 265: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/265.jpg)
Sensibilidade de F/F* à incerteza em i / i*
1
F/F*
i / i *
)(F
)(F
)/(
)](F/)(F[)/;F/F(S
*
i
*
i
*ii
i
1
*
ii
*
ii*
ii
*
F
i i *
F*
![Page 266: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/266.jpg)
Utilizando um incremento de 1% para melhor aproximar a derivada
Em processos complexos é impossível obter a derivada aproximação linear
)(F
)(F
)/(
)](F/)(F[;
F
FS
*i
*i
i
i*ii
*ii
*i
i*
*i
1
i
*i
*i
*ii
*i
*i
*i
i
*ii
*i
*i
i* )(F
)(F)(F
)(F
)(F)(F,
F
FS
01,0/ *
ii
)F(
)F(ξ)ξF(1,01100
ξ,
F
FS
*i
*i
*i
*i
i*
![Page 267: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/267.jpg)
S(F/F*;i/ i*) estima a incerteza % em F diante de uma incerteza
de 1% em i
|S| > 1 : incerteza ampliada |S| < 1 : incerteza amortecida
)F(
)F(ξ)ξF(1,01100
ξ,
F
FS
*i
*i
*i
*i
i*
![Page 268: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/268.jpg)
S (T2;U) = 100 (24,828-25)/25 = - 0,686
S (T4;U) = 100(30,047-30)/30 = 0,156
S(A;U) = 100 (262,93-265,6)/265,6 = - 0,99
S(W3;U) = 0
QUESTIONAMENTO DO PROJETORe-dimensionamento com U = 101
QUESTIONAMENTO DO DESEMPENHO
Simulação com U = 101
DIMENSIONAMENTO ORIGINAL(BASE)
T1* = 80 oC
W1* = 30.000 kg/h
A = 265,6 m2
[U = 100]
T 2* = 25 oC
W3 = 44.000 kg/hT3
* = 15 oC
T4* = 30 oC
[U = 101]
A = 262,93 m2
T1* = 80 oC
W1* = 30.000 kg/h
T3* = 15 oC
W3 = 44.000 kg/h
T4* = 30 oC
T 2* = 25 oC
[U = 101]
T2 = 24,828 oC
T1* = 80 oC
W1* = 30.000 kg/h
T3* = 15 oC
W3* = 44.000 kg/h
T4 = 30,047 oC
A* = 265,6 m2
![Page 269: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/269.jpg)
Para um incremento de 1% em todos os parâmetros (aproximação da derivada):
n
1i*
i
*
ii*
ii
*
*
*
)/;F/F(S)(F
)(F)(F
);F(S01,0)/;F/F(S01,0)(F
)(F)01,1(F n
1i
*
ii
*
*
**
Ou seja, a Sensibilidade de F à incerteza em todos os parâmetros é a soma das Sensibilidades a cada parâmetro:
)F(ξ
)F(ξ)ξF(1,01100)ξ/ξ;S(F/F
*
****
A Sensibilidade de F à incerteza em todos os parâmetros, considerados simultaneamente, pode ser estimada:
![Page 270: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/270.jpg)
Questionamento do Projeto
Sensibilidades de W3, A e CT à incerteza em cada parâmetro e variável especificada e ao conjunto:
i S(W3; i) S(A; i) S(CT; i)
W1 1 1 0,93
T1 1,45 0,45 1,21
T3 1,01 0,56 0,88
Cp1 1 1 0,93
Cp3 - 1 0 - 0,78
U 0 - 1 - 0,13
Por serem adimensionais, as Sensibilidades podem ser comparadas.
S(F; ) 3,46 2,01 3,04
![Page 271: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/271.jpg)
Questionamento do Projeto
Sensibilidades de W3, A e CT à incerteza em cada parâmetro e ao conjunto de parâmetros:
i S(W3; i) S(A; i) S(CT; i)
S(F; ) 3,46 2,01 3,04
Caso os valores reais de todos os parâmetros estivessem positivamente afastados em 1% dos seus valores-base (usados no dimensionamento), os valores de corretos de W3, A e CT estariam afastados de seus valores-base (calculados no dimensionamento) nos percentuais acima.
![Page 272: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/272.jpg)
Questionamento do Desempenho
Sensibilidades de T2 e T4 à incerteza em cada parâmetro e variável especificada, e ao conjunto:
i S(T2; i) S(T4; i)
W1 0,80 0,32
T1 0,48 0,63
T3 0,48 0,37
W3 - 0,12 - 0,47
A - 0,68 0,17
Cp1 0,80 0,32
Cp3 - 0,12 - 0,47
U - 0,68 0,17
Por serem adimensionais, as Sensibilidades podem ser comparadas.
S(F; ) 0,96 1,04
![Page 273: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/273.jpg)
Questionamento do Desempenho
Sensibilidades de T2 e T4 à incerteza em cada parâmetro e variável especificada, e ao conjunto:
i S(T2; i) S(T4; i)
S(F; ) 0,96 1,04
Caso os valores reais de todos os parâmetros estivessem positivamente afastados em 1% dos seus valores-base (usados no dimensionamento), os valores reais esperados para T2 e T4 , durante a operação do trocador, estariam afastados de seus valores-base nos percentuais acima.
![Page 274: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/274.jpg)
FIM
![Page 275: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/275.jpg)
f11 f12 f13 W15 f23 f31 f32 x13 k x12 Td Vd T2 T3 W2 W3
1 1 1 1
2 1 1
3 1 1
4 1 1 1
5 1 1
6 1 1 1 1
7 1 1 1 1
8 1 1
9 1 1
10 1 1
11 1 1
12 1
13 1 1 1
14 1 1 1
15 1 1 1
16 1 1 1
ORDENAR EQUAÇÕESMONTAR MATRIZ
EXTRATOR: DIMENSIONAMENTO
![Page 276: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/276.jpg)
f11 f12 f13 f23 f31 f32 x13 k x12 Td tau r T2 T3 x11 W2 W3
1 1 1 1
2 1
3 1 1
4 1 1 1
5 1 1
6 1 1 1
7 1 1 1
8 1 1 1
9 1 1
10 1 1
11 1 1
12 1 1
13 1 1 1
14 1 1 1
15 1 1 1
16 1 1 1
ORDENAR EQUAÇÕESMONTAR MATRIZ
EXTRATOR: SIMULAÇÃO
![Page 277: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/277.jpg)
f13 f14 f23 f24 W5 W6 W7 Qe Ae De T4 T5 W4
11 1 1
12 1 1 1
13 1 1
14 1 1
15 1 1 1 1
16 1 1 1
17 1
18 1
19 1
38 1 1
39 1
40 1 1 1
41 1 1
ORDENAR EQUAÇÕESMONTAR MATRIZ
EVAPORADOR: DIMENSIONAMENTO
![Page 278: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/278.jpg)
f13 f14 f23 f24 W5 W6 W7 Qe De T4 T5 W4 x14
11 1 1
12 1 1 1
13 1 1
14 1 1
15 1 1 1 1
16 1 1
17 1
18 1
19 1
38 1 1
39 1
40 1 1 1
41 1 1 1
ORDENAR EQUAÇÕESMONTAR MATRIZ
EVAPORADOR: SIMULAÇÃO
![Page 279: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/279.jpg)
EXEMPLO: convergência pela Bisseção
31. x31 = 1 – x1132. x32 = 1 – x12a04. x13 = k x12a07. W3 = W11 x11 r / x1301. W2 = W1 x31 / x3202. W1 x11 – W2 x12a – W3 x13 = 0
x12a
32 04 07 0201x31
x32
x13
W3
W2BISS
x12a
f (x12)
Na ordenação das equações surgiu um ciclo. A eq. 02 foi escolhida como Final. Resultou x12 como variável de abertura
f (x12) = W1 x11 – W2 x12a – W3 x13 = 0
![Page 280: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/280.jpg)
f (x)
x
xi
xs
fs
x1
f1x2
f2
fi
Esquema de convergência pela Bisseção
32 04 07 0201x31
x32
x13
W3
W2BISS
x12a
f (x12)
f (x12) = W1 x11 – W2 x12a – W3 x13 = 0
Até convergir
![Page 281: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/281.jpg)
EXEMPLO: convergência pela Substituição Direta
31. x31 = 1 – x1132. x32 = 1 – x12a04. x13 = k x12a07. W3 = W11 x11 r / x1301. W2 = W1 x31 / x32 02. W1 x11 – W2 x12a – W3 x13 = 0
x12a
32 04 07 0201x31
x32
x13
W3
W2SD
x12a = x12c
x12c
Na ordenação das equações surgiu um ciclo. A eq. 02 foi escolhida como Final. Resultou x12 como variável de abertura
x12c = (W1 x11– W3 x13) / W2
![Page 282: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/282.jpg)
Um instrumento fundamental para a resolução de problemas
ALGORITMO
![Page 283: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/283.jpg)
ALGORITMO é uma seqüência inequívoca de ações bem definidas que conduzem sempre à solução de um problema
Assim, qualquer pessoa, ou mesmo um computador por ela programado, chegará sempre à solução do problema.
O exemplo mais trivial e prosaico de algoritmo é uma receita culinária. Um outro no campo da matemática é o da
extração da raiz quadrada de um número.
Algoritmos podem incluir etapas repetitivas (iterações) ou exigir decisões (lógica e comparações).
Algoritmos podem ser programas em computadores
Existem algoritmos complexos e poderosos capazes de gerar outros algoritmos (Inteligência Artificial)
![Page 284: CAPÍTULO 3 ESTRATÉGIAS DE CÁLCULO](https://reader035.fdocumentos.tips/reader035/viewer/2022062517/568139b8550346895da156e3/html5/thumbnails/284.jpg)
Origem dos Algoritmos
An algorithm is a procedure or formula for solving a problem. The word derives from the name of the mathematician, Mohammed ibn-Musa al-Khwarizmi, who was part of the royal court in Baghdad and who lived from about 780 to 850. Al-Khwarizmi's work is the likely source for the word algebra as well.
www.nist.gov/dads/html/algorithm.html