Questão 1 - Na solução numérica da equação de Blasius
utilizando o método de Runge-Kutta de quarta ordem quais são as condições de contorno utilizadas na solução numérica? Explique cada uma delas. (1,0 ponto).
Passo 1
O enunciado pergunta: quais são as condições de contorno usadas na solução numérica da equação de Blasius
quando se emprega o método de Runge–Kutta de quarta ordem, e pede que cada condição seja explicada. A única informação extra do enunciado é que a solução é numérica (ou seja, transforma-se o problema numa integração a partir de um ponto inicial).
Passo 2
Apresento as condições de contorno clássicas usadas no problema de Blasius (vou somente listá‑las aqui; a explicação vem a seguir):
Passo 3
Explicação física e intuitiva de cada condição:
- Para
: em problemas de camada limite a função é a função de corrente adimensionalizada (normalizada por e escalas de comprimento). O valor corresponde à condição de impermeabilidade da parede: a componente normal da velocidade na superfície é nula (não há penetração através da parede). - Para
: a derivada é proporcional à componente tangencial da velocidade do fluxo adimensionalizada. Assim representa a condição de não deslizamento (no‑slip) na parede: a velocidade do fluido junto à parede iguala a velocidade da parede (aqui a parede está fixa, então velocidade tangencial é zero). - Para
: longe da parede (no infinito em coordenada adimensional ) a velocidade tangencial deve igualar a velocidade livre . Como a adimensionalização foi feita por , isso dá . Essa é a condição de correspondência com o escoamento exterior invíscido.
Passo 4
Por que isso é um desafio numérico e o que se faz na prática: a equação é de terceira ordem, portanto formalmente precisa de três condições. Duas condições são dadas em
(ponto inicial), mas a terceira está em infinito, não no ponto inicial — isso impede aplicar diretamente um método de passo único que exige todos os dados no mesmo ponto. A estratégia numérica típica é transformar o problema de contorno numa família de problemas de valor inicial (método de shooting): assume‑se um valor inicial desconhecido
(e
é um parâmetro a determinar), e integra‑se a partir de
até um
grande usando RK4. Ajusta‑se
até que a condição no infinito seja satisfeita numericamente, isto é até que
.
Passo 5
Como converter para um sistema de primeira ordem (passo necessário para aplicar RK4) e como as condições aparecem: defina
O sistema equivalente é
As condições iniciais em
ficam
onde
é a incógnita que será ajustada. O RK4 integra esse sistema para um dado
. Após a integração até
calcula‑se o resíduo
e usa‑se um método de root‑finding (por exemplo bissecção, secante ou Newton) sobre
para levar
a zero dentro de uma tolerância.
Passo 6
Algumas dicas práticas e números úteis (isto ajuda quando for implementar na prova ou no laboratório):
- Em vez de
usa‑se um finito (tipicamente entre e , dependendo de tolerância), pois para Blasius a solução já chegou muito próximo ao regime assintótico nesse intervalo. - Escolha uma tolerância razoável para
como ou . - Um chute inicial para
costuma ser perto de ; o valor conhecido (referência) é aproximadamente . - O método RK4 exige que em cada passo se calcule quatros incrementos (k1…k4) por componente; para um sistema
a iteração básica é (em notação compacta):
Aplicando isso ao sistema acima e repetindo até
obtém‑se
para comparar com
.
Passo 7
Significado adicional de
: além de ser o parâmetro de ajuste,
tem interpretação física — está ligado à tensão de cisalhamento na parede (skin friction). Por isso encontrar
é também encontrar a intensidade do atrito na parede no problema de Blasius.