Lucaso em 22h44min de 5 de fevereiro de 2024

2024-02-05T22:44:36Z

← Edição anterior		Edição das 19h44min de 5 de fevereiro de 2024
Linha 9:		Linha 9:
	<center><math>		<center><math>
	\frac{u_i^{n+1} - u_i^n}{\Delta t} + v\frac{\Delta u_i^{n} - u_{i-1}^{n}}{\Delta x} = 0		\frac{u_i^{n+1} - u_i^n}{\Delta t} + v\frac{\Delta u_i^{n} - u_{i-1}^{n}}{\Delta x} = 0
	</math></center>(15)		</math></center>(1)

	Para <math>v < 0</math>:		Para <math>v < 0</math>:
Linha 15:		Linha 15:
	<center><math>		<center><math>
	\frac{u_i^{n+1} - u_i^n}{\Delta t} + v\frac{u_{i+1}^{n} - u_{i}^{n}}{\Delta x} = 0		\frac{u_i^{n+1} - u_i^n}{\Delta t} + v\frac{u_{i+1}^{n} - u_{i}^{n}}{\Delta x} = 0
	</math></center>(16)		</math></center>(2)

	Combinando (15) e (16):		Combinando (1) e (2):

	<center><math>		<center><math>
	u_i^{n+1} = u_i^{n} - \Delta t (v^{+}u_x^{-} + v^{-}u_x^{+})		u_i^{n+1} = u_i^{n} - \Delta t (v^{+}u_x^{-} + v^{-}u_x^{+})
	</math></center>(17)		</math></center>(3)

	Onde		Onde

Lucaso em 22h14min de 4 de fevereiro de 2024

2024-02-04T22:14:52Z

← Edição anterior		Edição das 19h14min de 4 de fevereiro de 2024
Linha 155:		Linha 155:
	<center>		<center>
	{\|		{\|
	\|[[Arquivo: amplitude up-w.png\|thumb\|upright=0.0\|center\|Solução pelo método ~~Lax-Wendroff~~\|600px]]		\|[[Arquivo: amplitude up-w.png\|thumb\|upright=0.0\|center\|Solução pelo método upwind differencing\|600px]]
	\|}		\|}
	</center>		</center>

Lucaso: Criou página com 'Foram apresentados métodos numéricos eficientes com propósitos gerais para resolver a equação de advecção. Entretanto, é necessário restringir nossa abordagem a formas de onda suficientemente suaves para evitar oscilações espúrias tanto nas bordas de ataque quanto nas bordas de queda da forma de onda. Uma abordagem reconhecida para suprimir desvios nas extremidades inicial e final de uma forma de onda afiada é a adoção de um esquema conhecido como "esquema...'

2024-02-04T21:41:15Z

Criou página com 'Foram apresentados métodos numéricos eficientes com propósitos gerais para resolver a equação de advecção. Entretanto, é necessário restringir nossa abordagem a formas de onda suficientemente suaves para evitar oscilações espúrias tanto nas bordas de ataque quanto nas bordas de queda da forma de onda. Uma abordagem reconhecida para suprimir desvios nas extremidades inicial e final de uma forma de onda afiada é a adoção de um esquema conhecido como "esquema...'

Página nova

Foram apresentados métodos numéricos eficientes com propósitos gerais para resolver a equação de advecção. Entretanto, é necessário restringir nossa abordagem a formas de onda suficientemente suaves para evitar oscilações espúrias tanto nas bordas de ataque quanto nas bordas de queda da forma de onda. Uma abordagem reconhecida para suprimir desvios nas extremidades inicial e final de uma forma de onda afiada é a adoção de um esquema conhecido como "esquema de diferenças upwind".

Nesse esquema, as diferenças espaciais são inclinadas na direção "upwind", ou seja, na direção de onde o fluxo advectivo se origina.

Portanto, a versão upwind do simples esquema de diferenças explícitas é escrita da seguinte forma:

Para <math>v > 0</math>:

<center><math>
\frac{u_i^{n+1} - u_i^n}{\Delta t} + v\frac{\Delta u_i^{n} - u_{i-1}^{n}}{\Delta x} = 0
</math></center>(15)

Para <math>v < 0</math>:

<center><math>
\frac{u_i^{n+1} - u_i^n}{\Delta t} + v\frac{u_{i+1}^{n} - u_{i}^{n}}{\Delta x} = 0
</math></center>(16)

Combinando (15) e (16):

<center><math>
u_i^{n+1} = u_i^{n} - \Delta t (v^{+}u_x^{-} + v^{-}u_x^{+})
</math></center>(17)

Onde

<center><math>
\begin{cases}
v^{+} &= max(v,0) \\
v^{-} &= min(v,0) \\
\end{cases}
</math></center>

e

<center><math>
\begin{cases}
u_x^{-} &= \frac{u_i^{n} - u_{i-1}^{n}}{\Delta x} \\
u_x^{+} &= \frac{u_{i+1}^{n} - u_{i}^{n}}{\Delta x}
\end{cases}
</math></center>

O esquema upwind é estável quando a seguinte condição Courant-Friedrichs-Lewy (CFL) é atendida:

<center><math>
|v\frac{\Delta t}{\Delta x}| \leq 1
</math></center>

= Implementação do método =

* Condição inicial: <math>u(x,0) = 1-\cos(x)</math>;

* Condições de contorno para bordas cíclicas.

<source lang = "python">

# Solução pelo método Upwind para equação de advecção

def Upwindad(L, tf, v, Nx, Nt):
"""
Parâmetros:
- L: comprimento
- tf: tempo final
- v: velocidade de propagação
- Nx: número de pontos na direção espacial
- Nt: número de pontos na direção temporal

Retorna:
- Matriz com a solução da equação da onda
"""

dx = L / (Nx - 1)
dt = tf / (Nt - 1)
r = v * dt / dx

u = np.zeros((Nt, Nx+1))

# Condição inicial: u(x,0) = f(x)
x = np.linspace(0, L, Nx+1)
u[0,:] = 1 - np.cos(x) # Função que descreve a perturbação da onda

# Condições de contorno borda infinita:
xpos = np.zeros(Nx+1)
xneg = np.zeros(Nx+1)

for i in range(0, Nx+1):
xpos[i] = i+1
xneg[i] = i-1
xpos[Nx] = 0
xneg[0] = Nx

# Iteração no tempo
for n in range(0, Nt - 1):
for i in range(0, Nx+1):
v_mais = np.maximum(v, 0)
v_menos = np.minimum(v, 0)
u_x_mais = 1/dx * (u[n, int(xpos[i])] - u[n, i])
u_x_menos = 1/dx * (u[n, i] - u[n, int(xneg[i])])

u[n+1, i] = u[n, i] - dt * (v_mais * u_x_menos + v_menos * u_x_mais)

return u

</source>

<source lang = "python">

# Parâmetros
L = 2*np.pi
tf =1
v = -1 # -1. muda direção de propagação
Nx = 100
Nt = 500

solv4 = Upwindad(L, tf, v, Nx, Nt)

listX = np.linspace(0, L, Nx+1)
listT = np.linspace(0, tf, Nt)

X, T = np.meshgrid(listX, listT)

plt.figure(figsize=(10, 6))
plt.pcolormesh(X, T, solv4, cmap='viridis', shading='auto')
plt.colorbar(label='Amplitude(u)')
plt.xlabel('Posição (x)')
plt.ylabel('Tempo (t)')
plt.title('Solução Upwind da Equação da advecção (1D)', fontsize=16)
plt.show()

</source>

<center>
{|
|[[Arquivo: up-w.png|thumb|upright=0.0|center|Solução pelo método upwind differencing|600px]]
|}
</center>

<source lang = "python">

# Teste: Plota todas as cuvas amplitude por posição de todos os tempos:

for tt in range(len(listT-1)):
amplitudes_tt = solv4[ tt,:]
plt.plot(listX, amplitudes_tt)

plt.title('Amplitude em Função da Posição')
plt.xlabel('Posição (x)')
plt.ylabel('Amplitude (u)')
plt.legend()
plt.grid(True)
plt.show()

</source>

<center>
{|
|[[Arquivo: amplitude up-w.png|thumb|upright=0.0|center|Solução pelo método Lax-Wendroff|600px]]
|}
</center>

Upwind Differencing - Histórico de revisão

Lucaso em 22h44min de 5 de fevereiro de 2024

Lucaso em 22h14min de 4 de fevereiro de 2024