[2.0 pts] Modelo Dinâmico de Transferência de Oxigênio
Considere a concentração de oxigênio dissolvido,
, em um biorreator. Uma modelagem simplificada pode gerar:
onde:
-
é concentração de oxigênio dissolvido,
- K é o seu último número de matrícula.
d) Analise a Estabilidade (Estável, Instável ou Neutra). Justifique sua resposta.
Considere a concentração de oxigênio dissolvido,
Passo 1
O enunciado dá a equação diferencial linear com coeficientes constantes
onde
é a concentração e
é um parâmetro (o último número da matrícula). A pergunta é: analisar a estabilidade do sistema — isto é, dizer se o sistema é estável (todas as soluções tendem a zero quando
), instável (alguma solução cresce) ou neutro/marginal (não cresce nem decai para zero, por ex. oscilações sustentadas ou termo constante).
Informações relevantes extraídas do enunciado:
- É uma EDO linear homogênea de ordem 3 com coeficientes constantes.
- Coeficientes numéricos:
na derivada terceira, na segunda, na primeira, e na parte algébrica. - A estabilidade será determinada pelos autovalores (raízes do polinômio característico) — precisamos saber o sinal das partes reais dessas raízes em função de
.
Passo 2
Conceito e estratégia (explicação didática):
- Para uma EDO linear homogênea com coeficientes constantes, a solução é combinação de termos exponenciais
onde são as raízes do polinômio característico. A estabilidade assintótica exige que todas as partes reais de sejam negativas. - Aqui o polinômio característico é obtido trocando
por :
- Em vez de resolver explicitamente as raízes (o que dá conta, mas é trabalhoso e depende de
), existe um critério prático chamado critério de Routh–Hurwitz que permite verificar estabilidade (todas partes reais negativas) usando apenas os coeficientes do polinômio. Para um polinômio cúbico há uma condição simples e suficiente que vou apresentar e explicar.
Passo 3
Routh–Hurwitz para polinômio cúbico — explicação e fórmula:
- Considere o polinômio genérico cúbico
com . - O critério de Routh–Hurwitz diz que todas as raízes têm parte real negativa (estabilidade assintótica) se, e somente se, as seguintes condições se verificarem simultaneamente:
-
, -
, -
, -
.
-
- Intuição curta: essas condições garantem que, ao construir a tabela de Routh (que organiza combinações lineares dos coeficientes), a primeira coluna seja toda positiva, o que impede mudança de sinal e portanto raízes com parte real positiva.
- Vamos aplicar diretamente essa condição ao nosso polinômio, identificando
, , , .
Passo 4
Aplicação numérica das condições de Hurwitz ao polinômio dado:
- Checamos as condições fáceis:
-
(ok), -
(ok).
-
- Para
precisamos de . Portanto é necessário para estabilidade. - A condição adicional é
. Substituindo os valores:
Calcule o lado esquerdo:
Assim a condição torna-se
- Conclusão parcial: para que todas as raízes tenham parte real negativa (estabilidade assintótica) é necessário e suficiente que
Passo 5
Análise dos casos limítrofes e fora da faixa — interpretação detalhada (neutro / instável):
- Caso
: então . Pelo sinal alternado dos coeficientes ou pelo critério de Routh (a primeira coluna da tabela tem mudança de sinal) existe pelo menos uma raiz com parte real positiva. Logo o sistema é instável para . - Caso
: o polinômio fica
Uma raiz é
(modo constante) e as outras duas são raízes do quadrático
. O discriminante desse quadrático é
logo as duas outras raízes são complexas conjugadas com parte real negativa (porque os coeficientes do quadrático são positivos). Assim há uma raiz exatamente em zero e as demais têm parte real negativa → comportamento marginal (neutro): não há decaimento total (termo constante), mas também não há crescimento exponencial. Classificamos como neutra/marginal para
.
- Caso
: aqui a condição limite ocorre. Na tabela de Routh o elemento da linha de torna-se zero; isso corresponde a um par de raízes puramente imaginárias. Para ver isso formamos o polinômio auxiliar a partir da linha :
Com
obtemos
logo
(raízes puramente imaginárias) e a terceira raiz será real negativa (porque os coeficientes são positivos). Assim o sistema apresenta oscilações sustentadas (par de imaginários) e decaimento ausente para esse modo — isto é, comportamento neutro/marginal. Portanto
dá estabilidade marginal (neutra), não assintótica.
- Caso
: a quantidade (que aparece na primeira coluna da linha da tabela de Routh) torna-se negativa, o que implica mudança de sinal na primeira coluna e portanto existência de raiz com parte real positiva → sistema instável.
Passo 6
Resumo final das conclusões (didático e direto):
- A faixa de estabilidade assintótica é exatamente quando todos os coeficientes são positivos e
, o que deu . - Para
fora dessa faixa ocorrem situações de instabilidade (alguma raiz com parte real positiva) ou estabilidade marginal (raiz nula ou par imaginário).
Resposta
Para
: Instável.
Para : Neutra (marginal), há uma raiz
e os outros dois têm parte real negativa.
Para : Estável (assintoticamente estável), todas as raízes têm parte real negativa.
Para : Neutra (marginal), existe um par de raízes puramente imaginárias.
Para : Instável.
Para
Para
Para
Para