Teste conv: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Heitor (discussão | contribs)
Sem resumo de edição
Heitor (discussão | contribs)
Sem resumo de edição
 
(Uma revisão intermediária pelo mesmo usuário não está sendo mostrada)
Linha 4: Linha 4:
|}
|}


<math>f'(x)=\lim_{\Delta x\to 0}{f(x+\Delta x)-f(x)\over \Delta x} \frac{\partial}{\partial t} \int_{V}f(x,t) = -\int_{S}\mathbf{J} \cdot d\mathbf{s} + \int_{V}S(x,t)dV</math>
:<math>S = \int_a^b f(x)\, dx = \lim_{N\to\infty} \sum_{i=0}^{N} f(x_i) \Delta x </math>


<math> \frac{\partial}{\partial t} \int_{V}f(x,t) = -\int_{S}\mathbf{J} \cdot d\mathbf{s} + \int_{V}S(x,t)dV</math>
<math>f'(x)=\lim_{\Delta x\to 0}{f(x+\Delta x)-f(x)\over \Delta x^3}</math>
 
<math>f'(x)=\lim_{\Delta x\to 0}{f(x+\Delta x)-f(x)\over \Delta x^2} </math>
 
 
<math> x^3 </math>
 
<math> x^2 </math>
 
<math>f'(x)=\lim_{\Delta x\to 0}{f(x+\Delta x)-f(x)\over \Delta x} </math>
 
<math> S = \frac{\partial}{\partial t} \int_{V}f(x,t) = -\int_{S}\mathbf{J} \cdot d\mathbf{s} + \int_{V}S(x,t)dV </math>


<math>\frac{\partial}{\partial t} \int_{V}f(x,t) = -\int_{V} (\nabla \cdot \mathbf{J}) dV + \int_{V}S(x,t)dV</math>
<math>\frac{\partial}{\partial t} \int_{V}f(x,t) = -\int_{V} (\nabla \cdot \mathbf{J}) dV + \int_{V}S(x,t)dV</math>
Linha 73: Linha 84:


<math>0 < k \leq \frac{1}{2}</math>
<math>0 < k \leq \frac{1}{2}</math>
TESTE!!

Edição atual tal como às 12h53min de 19 de setembro de 2017

f′(x)=limΔx→0f(x+Δx)−f(x)Δx (Eq. 1)
S=∫abf(x)dx=limN→∞∑i=0Nf(xi)Δx

f′(x)=limΔx→0f(x+Δx)−f(x)Δx3

f′(x)=limΔx→0f(x+Δx)−f(x)Δx2


x3

x2

f′(x)=limΔx→0f(x+Δx)−f(x)Δx

S=∂∂t∫Vf(x,t)=−∫S𝐉⋅d𝐬+∫VS(x,t)dV

∂∂t∫Vf(x,t)=−∫V(∇⋅𝐉)dV+∫VS(x,t)dV

∫V(∂∂tf(x,t)+∇⋅𝐉−S(x,t))dV=0

∂f∂t=−∇⋅𝐉+S(x,t)

Lei de Fick:

𝐉=−D∇f

Onde D é a constante de difusão.

∂f∂t=−∇⋅(D∇f)+S(x,t)

∂f∂t=D∇2f+S(x,t)+v∂f∂x

Equação da difusão:

∂f∂t=D∇2f

Em uma dimensão:

∂f∂t=D∂2f∂t2

FTCS (Foward Time Central Space):

∂f∂t→f(x,t+Δt)−f(x,t)Δt

∂2f∂x2→f(x,t+Δt)+f(x−Δx,t)−2f(x,t)Δx2

(Escrever a equação em termos numéricos...)

Teste de establilidade do método FTCS:

Um dos modos de Fourier da solução:

fjn=Aneiqjh

h=Δx

fjn+1=An+1eiqjh

fj+1n=Aneiq(j+1)h

fj−1n=Aneiq(j−1)h

An+1eiqjh=k(Aneiq(j+1)h+Aneiq(j−1)h−2Aneiqjh)+Aneiqjh

k=DΔtΔx2

|An+1An|≤1

An+1An=1+k(eiqh+e−iqh−2)

An+1An=1+2k[cos(qh)−1]

An+1An=1+4ksen2(qh2)

ξ=|An+1An|=|1+4ksen2(qh2)|

Na pior hipótese, o seno quadrado é 1.

ξ=|1+4k|≤1

0<k≤12

TESTE!!