Enunciado

Substâncias dissolvidas em um aquífero são transportadas por dois mecanismos separados. O processo pelo qual um soluto é transportado pelo grande movimento do fluxo da água subterrânea é chamado de advecção. Além disso, o soluto é espalhado pelas flutuações em escala pequena na velocidade da água subterrânea ao longo de caminhos tortuosos do fluxo dentro de poros individuais, um processo chamado de dispersão mecânica. A forma unidimensional da equação de advecção-dispersão para uma substância dissolvida não reativa em um meio poroso isotrópico homogêneo saturado sob um fluxo uniforme regular é

em que c(x, t) é a concentração do soluto, υ é a velocidade linear média da água subterrânea, D é o coeficiente da dispersão hidrodinâmica e L é o comprimento do aquífero. Suponha que as condições de contorno são

e que a condição inicial é

em que f(x) é a concentração inicial dada do soluto.

Usando a solução encontrada no item (e), faça o gráfico de c(x, t) em função de x para diversos valores de t, tais como t = 0,5, 1, 3, 6 e 10. Faça também o gráfico de c(x, t) em função de t para diversos valores de x. Note que o número de termos necessários para obter um gráfico preciso depende fortemente dos valores de t e de x..

Passo 1

Oi oi, meu xuxuzinho!! Tudo bem com você?? Vamos analisar o que a nossa questão nos pede e matar essa em 10 minutos?? Bora. Bom, nosso enunciado nos pede para analisarmos as deformações verticais de uma barra em suas condições de contorno. Parece muita coisa, mas é tranquilo. Partiu??

Passo 2

Ok, nosso primeiro passo vai ser reescrever a equação diferencial.

X n x = e v x 2 D sin ⁡ ( μ n x )  

X n ' x = e v x 2 D v 2 D sin ⁡ μ n x + μ n cos ⁡ μ n x

Com isso, teremos o seguinte sistema de equações:

1 T t d T d t = λ D 1 X x d 2 X d x 2 - ( v / D ) X x d X d x   = λ  

T ' + λ D T = 0 X ' ' - v D X ' + λ X = 0  

Lembrando que:

x = L

Passo 3

Pronto!! Agora é só fazermos umas matemágicas aí pra achar o resultado!!

X n x = e v x D sin ⁡ μ n x

c x , t = ∑ a n   e - λ o D x   e v x D sin ⁡ μ n x

Resolvendo essa equação, temos que:

a n = 2 ∫ d x sin ⁡ μ n x f x L - μ n - 1 sin ⁡ 2 μ n L

μ n = 2 n - 1 π 20

Resposta

μ n = 2 n - 1 π 20

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Vamos considerar agora um problema de autovalores de ordem mais alta. A análise das vibrações transversas de uma barra elástica uniforme se baseia na equação diferencial λ 4 - λ y = 0 , em que y é o d
Vamos considerar agora um problema de autovalores de ordem mais alta. A análise das vibrações transversas de uma barra elástica uniforme se baseia na equação diferencial λ 4 - λ y = 0 , em que y é o d
Vamos considerar agora um problema de autovalores de ordem mais alta. A análise das vibrações transversas de uma barra elástica uniforme se baseia na equação diferencial λ 4 - λ y = 0 , em que y é o d
Em cada um dos problemas de 14 a 18, determine se o problema de valores de contorno dado é autoadjunto. y ' ' + y ' + 2 y = 0 , y 0 = 0ey 1 = 0
Siga as instruções do Problema 2 usando f ( x )=x ( 1–x ) para 0x1 .