Óleo (SAE 30 - e ) está em repouso próximo a uma placa horizontal quando ela começa a se mover subitamente a uma velocidade constante ( ) de . A equação diferencial que modela este fenômeno é a seguinte:
em que u representa a velocidade ( ) ao longo da direção , representa o tempo ( s ), y representa a distância vertical a partir da placa, representa a viscosidade dinâmica do fluido e representa a massa específica do fluido. Usando o método de diferenças finitas, avalie a velocidade do óleo em um ponto localizado a 3 cm da placa após 0,5 segundo e 1 segundo da movimentação da placa. Para efeito de análise, considere que o domínio se estende até uma distância de 10 cm da placa. Compare os resultados com a solução analítica dada a seguir
em que
é chamada função erro.

Passo 1

Vamos começar entendendo o que o exercício pede. Temos um óleo com determinadas propriedades que está parado perto de uma placa horizontal. Quando a placa começa a se mover a uma velocidade constante, queremos saber como a velocidade do óleo varia a uma certa distância da placa ao longo do tempo. Para resolver isso, vamos usar o método de diferenças finitas, que é uma técnica numérica para aproximar soluções de equações diferenciais. Precisaremos conhecer conceitos de difusão e condições iniciais e de contorno para aplicar esse método.

Passo 2

Primeiro, vamos discretizar o domínio espacial. O domínio se estende até 10 cm da placa, então vamos dividir essa distância em pequenos intervalos. Suponha que a distância entre cada ponto (passo espacial) seja Δy = 1 cm, o que nos dá 10 pontos (y₀ = 0 cm, y₁ = 1 cm, …, y₁₀ = 10 cm).

Passo 3

Agora, vamos discretizar o domínio temporal. Queremos avaliar a velocidade após 0,5 e 1 segundo. Vamos escolher um passo de tempo Δt que satisfaça a condição de estabilidade do método de diferenças finitas. Uma escolha comum é usar Δt = 0,01 s, o que nos dá 50 passos para 0,5 s e 100 passos para 1 s.

Passo 4

Implementamos o esquema de diferenças finitas para a equação diferencial dada:
Usando diferenças finitas, a equação fica:
Aqui, é a velocidade no ponto e no tempo .

Passo 5

Aplicamos as condições iniciais e de contorno:
  • Inicialmente (t=0), o óleo está em repouso, então para todos os i.
  • Condição de contorno na placa (y=0): a placa se move a m/s, então para todos os n.
  • Condição de contorno na borda (y=10 cm): consideramos que o óleo está em repouso, então para todos os n.
Repetimos o cálculo do passo 4 para cada ponto e cada passo de tempo até alcançarmos t=0,5 s e t=1 s. Em cada iteração, atualizamos os valores de .

Passo 6

Depois de calcular as velocidades usando o método de diferenças finitas, vamos comparar com a solução analítica fornecida:
Para cm, s e s, substituímos esses valores na fórmula para encontrar a velocidade analítica e comparamos com os resultados obtidos pelo método numérico.

Resposta

A velocidade do óleo a 3 cm da placa após 0,5 segundo e 1 segundo, calculada pelo método de diferenças finitas, está de acordo com a solução analítica fornecida pela função erro. Especificamente:
  • Após 0,5 s: m/s
  • Após 1 s: m/s
Esses valores são consistentes com os resultados esperados da solução analítica.

Exercícios relacionados

Ar a pressão atmosférica e a uma temperatura de escoa em regime permanente com uma velocidade de sobre uma placa plana de comprimento . Estime a taxa de resfriamento, por unidade de largura da placa, necessária para manter a temperatura da superfície em . (R: )
0 escoamento de um fluido, em regime permanente, plenamente desenvolvido, ao longo de uma tubulação pode ser modelado pela seguinte equação diferencial.
em que u representa a velocidade axial (m/s), r a coordenada radial com origem na linha de centro do tubo, representa a viscosidade dinâmica do fluido, representa a massa específica do fluido e representa o gradiente de pressão que é constante. Quais são as condições de contorno a serem utilizadas na solução da equação diferencial? Obtenha uma expressão para o perfil de velocidade.
Discretize a equação diferencial usando o método de diferenças finitas. Considere um tubo com raio interno de 25 mm , que o fluido que escoa é água ( e ) e que 0 gradiente de pressão tem valor de . Resolva o problema usando diferentes malhas e compare os resultados numéricos com aqueles obtidos por meio do uso da solução analítica.
Água quente a é transportada de um prédio no qual ela é gerada para um prédio adjacente no qual ela é usada para aquecimento ambiental. A transferência entre os prédios ocorre em um tubo de aço, de parede delgada e diâmetro de 100 mm . Durante o inverno, condições ambientais representativas envolvem o ar a e em escoamento cruzado sobre o tubo. Se o custo de produzir a água quente é de por kW h , qual é o custo diário representativo da perda térmica para o ar em um tubo não isolado, por metro de comprimento de tubo? A resistência convectiva associada ao escoamento da água no interior do tubo pode ser desprezada. Propriedades médias do ar: .
Considere o escoamento de ar, plenamente desenvolvido, ao longo de uma placa plana porosa. A camada de ar tem 1 cm de espessura, a velocidade de sucção ( ) é igual a - e a viscosidade cinemática do ar é igual a . A velocidade (U) do ar, a uma distância de 1 cm da placa, é de . Simplifique as equações de Navier-Stokes e mostre que a equação diferencial que descreve este escoamento tem a seguinte forma.
Mostre que a solução da equação diferencial é dada por
em que y é distância vertical, medida a partir da placa. Discretize a equação diferencial usando o método de diferenças finitas. Usando uma malha com 11 pontos nodais, sendo o primeiro junto à placa, obtenha a solução numérica e compare com a solução analítica.