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
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>

Edição das 21h54min de 12 de agosto de 2015

f(x)=limΔx0f(x+Δx)f(x)Δx (Eq. 1)
S=abf(x)dx=limNi=0Nf(xi)Δx

f(x)=limΔx0f(x+Δx)f(x)Δx3

f(x)=limΔx0f(x+Δx)f(x)Δx2


x3

x2

f(x)=limΔx0f(x+Δx)f(x)Δx

S=tVf(x,t)=S𝐉d𝐬+VS(x,t)dV

tVf(x,t)=V(𝐉)dV+VS(x,t)dV

V(tf(x,t)+𝐉S(x,t))dV=0

ft=𝐉+S(x,t)

Lei de Fick:

𝐉=Df

Onde D é a constante de difusão.

ft=(Df)+S(x,t)

ft=D2f+S(x,t)+vfx

Equação da difusão:

ft=D2f

Em uma dimensão:

ft=D2ft2

FTCS (Foward Time Central Space):

ftf(x,t+Δt)f(x,t)Δt

2fx2f(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

fj1n=Aneiq(j1)h

An+1eiqjh=k(Aneiq(j+1)h+Aneiq(j1)h2Aneiqjh)+Aneiqjh

k=DΔtΔx2

|An+1An|1

An+1An=1+k(eiqh+eiqh2)

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<k12