Arquivo:Eq2.PNG

2021-05-29T02:13:55Z

Augustog:

Arquivo:Eq1.PNG

2021-05-29T02:13:14Z

Augustog:

Problema de Fermi-Pasta-Ulam

2021-05-29T02:03:40Z

Augustog: /* O Problema */

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas<ref name= FPU>https://www.respondeai.com.br/conteudo/fisica/mhs-e-mha/oscilacoes-acopladas/951</ref>.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F_{i} = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola (aqui considerado, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação ===

'''O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou'''<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref>

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Corda final.gif|center|thumb|500px|Simulação FPU para o problema proposto.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias fft.jpeg|thumb|700px|center|Porcentagem de cada modo vibracional ao longo do tempo. A escala do tempo está reduzida em 20 vezes.]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|500px|thumb|Energias por modo de oscilação. A escala do tempo está aumentada em 2 vezes.]]</li>
<li style="display: inline-block;"> [[Arquivo:Energias certo2.jpg|520px|thumb|Energias por modo de oscilação com fft. A escala do tempo está aumentada em 5 vezes.]]</li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados e o comportamento das energias através da simulação com dinâmica molecular, obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema <ref name=FPUT1></ref>.

A questão de simular a energia usando a porcentagens dos senos também faz sentido pensando ''quanta contribuição'' o modo normal de oscilação tem no movimento atual, como pudemos ver que o comportamento foi bem similar ao das energias calculadas.

== Implementação ==

A questão de variar os parâmetros do problema, como já mencionado previamente nesta wiki, e também apresentado nos gráficos e no gif, vimos que não interferia muito no caráter dos resultados (e o quanto isso poderia influenciar no paradoxo de equipartição de energia), porém o fizemos para gerar imagens melhores e o gif com menos de 2 MB.

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)
calcula_energias(pos)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/Fermi-Pasta-Ulam]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam

2021-05-29T01:58:05Z

Augustog:

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas<ref name= FPU>https://www.respondeai.com.br/conteudo/fisica/mhs-e-mha/oscilacoes-acopladas/951</ref>.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F_{i} = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola (aqui considerado, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação ===

'''O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou'''<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref>

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Corda final.gif|center|thumb|500px|Simulação FPU para o problema proposto.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias fft.jpeg|thumb|700px|center|Porcentagem de cada modo vibracional ao longo do tempo. A escala do tempo está reduzida em 20 vezes.]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|500px|thumb|Energias por modo de oscilação. A escala do tempo está aumentada em 2 vezes.]]</li>
<li style="display: inline-block;"> [[Arquivo:Energias certo2.jpg|520px|thumb|Energias por modo de oscilação com fft. A escala do tempo está aumentada em 5 vezes.]]</li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados e o comportamento das energias através da simulação com dinâmica molecular, obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema <ref name=FPUT1></ref>.

A questão de simular a energia usando a porcentagens dos senos também faz sentido pensando ''quanta contribuição'' o modo normal de oscilação tem no movimento atual, como pudemos ver que o comportamento foi bem similar ao das energias calculadas.

== Implementação ==

A questão de variar os parâmetros do problema, como já mencionado previamente nesta wiki, e também apresentado nos gráficos e no gif, vimos que não interferia muito no caráter dos resultados (e o quanto isso poderia influenciar no paradoxo de equipartição de energia), porém o fizemos para gerar imagens melhores e o gif com menos de 2 MB.

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)
calcula_energias(pos)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/Fermi-Pasta-Ulam]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam

2021-05-27T01:03:41Z

Augustog:

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação ===

'''O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou'''<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref>

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Corda (26).gif|center|thumb|500px|legenda.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias fft.jpeg|thumb|700px|center|Porcentagem de cada seno usando fft ]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|500px|thumb|Energias por modo de oscilação.]]</li>
<li style="display: inline-block;"> [[Arquivo:Energias certo2.jpg|520px|thumb|Energias por modo de oscilação com fft.]]</li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados e o comportamento das energias através da simulação com dinâmica molecular, obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema <ref name=FPUT1></ref>.

A questão de simular a energia usando a porcentagens dos senos também faz sentido pensando ''quanta contribuição'' o modo normal de oscilação tem no movimento atual, como pudemos ver que o comportamento foi bem similar ao das energias calculadas.

== Implementação ==

A questão de variar os parâmetros do problema, como já mencionado previamente nesta wiki, e também apresentado nos gráficos e no gif, vimos que não interferia muito no caráter dos resultados (e o quanto isso poderia influenciar no paradoxo de equipartição de energia), porém o fizemos para gerar imagens melhores e o gif com menos de 2 MB.

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)
calcula_energias(pos)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Arquivo:Corda (26).gif

2021-05-27T01:02:57Z

Augustog:

Arquivo:Energias fft.jpeg

2021-05-27T01:02:01Z

Augustog:

Problema de Fermi-Pasta-Ulam

2021-05-27T00:53:47Z

Augustog:

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação ===

'''O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou'''<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref>

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Trab2.gif|center|thumb|500px|legenda.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias vibracao3.png|thumb|700px|center|Porcentagem de cada seno usando fft ]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|500px|thumb|Energias por modo de oscilação.]]</li>
<li style="display: inline-block;"> [[Arquivo:Energias certo2.jpg|520px|thumb|Energias por modo de oscilação com fft.]]</li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados e o comportamento das energias através da simulação com dinâmica molecular, obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema <ref name=FPUT1></ref>.

A questão de simular a energia usando a porcentagens dos senos também faz sentido pensando ''quanta contribuição'' o modo normal de oscilação tem no movimento atual, como pudemos ver que o comportamento foi bem similar ao das energias calculadas.

== Implementação ==

A questão de variar os parâmetros do problema, como já mencionado previamente nesta wiki, e também apresentado nos gráficos e no gif, vimos que não interferia muito no caráter dos resultados (e o quanto isso poderia influenciar no paradoxo de equipartição de energia), porém o fizemos para gerar imagens melhores e o gif com menos de 2 MB.

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)
calcula_energias(pos)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Arquivo:Energias certo2.jpg

2021-05-27T00:31:52Z

Augustog:

Problema de Fermi-Pasta-Ulam

2021-05-27T00:14:05Z

Augustog:

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação ===

'''O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou'''<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref>

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Trab2.gif|center|thumb|500px|legenda.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias vibracao3.png|thumb|700px|center|Porcentagem de cada seno usando fft ]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> </li>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|thumb|500px|Energias por modo de oscilação.]] </li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados e o comportamento das energias através da simulação com dinâmica molecular, obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema <ref name=FPUT1></ref>.

A questão de simular a energia usando a porcentagens dos senos também faz sentido pensando ''quanta contribuição'' o modo normal de oscilação tem no movimento atual, como pudemos ver que o comportamento foi bem similar ao das energias calculadas.

== Implementação ==

A questão de variar os parâmetros do problema, como já mencionado previamente nesta wiki, e também apresentado nos gráficos e no gif, vimos que não interferia muito no caráter dos resultados (e o quanto isso poderia influenciar no paradoxo de equipartição de energia), porém o fizemos para gerar imagens melhores e o gif com menos de 2 MB.

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam

2021-05-27T00:13:17Z

Augustog:

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação: O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref> ===

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Trab2.gif|center|thumb|500px|legenda.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias vibracao3.png|thumb|700px|center|Porcentagem de cada seno usando fft ]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> </li>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|thumb|500px|Energias por modo de oscilação.]] </li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados e o comportamento das energias através da simulação com dinâmica molecular, obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema <ref name=FPUT1></ref>.

A questão de simular a energia usando a porcentagens dos senos também faz sentido pensando ''quanta contribuição'' o modo normal de oscilação tem no movimento atual, como pudemos ver que o comportamento foi bem similar ao das energias calculadas.

== Implementação ==

A questão de variar os parâmetros do problema, como já mencionado previamente nesta wiki, e também apresentado nos gráficos e no gif, vimos que não interferia muito no caráter dos resultados (e o quanto isso poderia influenciar no paradoxo de equipartição de energia), porém o fizemos para gerar imagens melhores e o gif com menos de 2 MB.

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam

2021-05-27T00:03:40Z

Augustog: /* Discussões */

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação: O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref> ===

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Trab2.gif|center|thumb|500px|legenda.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias vibracao3.png|thumb|700px|center|Porcentagem de cada seno usando fft ]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> </li>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|thumb|500px|Energias por modo de oscilação.]] </li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados e o comportamento das energias através da simulação com dinâmica molecular, obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema <ref=FPUT></ref>

== Implementação ==

Usamos XX partículas, com modo de oscilação YY

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

'''Falar um pouco mais da implementação?'''

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam

2021-05-27T00:02:42Z

Augustog: /* Resultados */

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação: O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref> ===

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

Segue abaixo a evolução do sistema ao longo do tempo, apenas para gerar o .gif considerando:

<math>
\begin{cases}

\alpha = 1,2 \\
k = 0,95 \\
m = 1,05 \\
N = 30 \\
t_{max} = 4000 \\

\end{cases}
</math>

[[Arquivo:Trab2.gif|center|thumb|500px|legenda.]]

Já os resultados das energias, consideramos todas constantes iguais à 1, mais de 100 partículas e o número de iterações está em cada imagem. procedemos de duas formas. Calculamos a porcentagem de cada modo normal de oscilação pela fft, o que representa um comportamento muito similar às energias por modo:

[[Arquivo:Energias vibracao3.png|thumb|700px|center|Porcentagem de cada seno usando fft ]]

E também calculamos pela equação de energia citada acima, obtivemos os seguintes comportamentos:

<center><div><ul>
<li style="display: inline-block;"> </li>
<li style="display: inline-block;"> [[Arquivo:Energias certo1.jpeg|thumb|500px|Energias por modo de oscilação.]] </li>
</ul></div></center>

== Discussões ==

Como o intuito era replicar os resultados através da simulação com dinâmica molecular obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema ('''Citar estes artigos''')

== Implementação ==

Usamos XX partículas, com modo de oscilação YY

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

'''Falar um pouco mais da implementação?'''

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Arquivo:Energias certo1.jpeg

2021-05-26T23:58:46Z

Augustog:

Problema de Fermi-Pasta-Ulam

2021-05-26T23:33:14Z

Augustog: /* Resultados */

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref> sobre sistemas dinâmicos não lineares. As análises serão sobre a solução dos modos de vibração comparados à solução analítica para poucas massas e também sobre a energia do sistema para os modos de oscilação, enquanto o sistema evolui no tempo.

== O Problema ==

O Problema proposto constitui-se de simulações em uma rede de partículas ligadas entre si através de molas que obedecem a Lei de Hooke com uma correção não-linear quadrática ou cúbica <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Springs.png|thumb|437px|Figura 1. Problema de Fermi-Pasta-Ulam, molas acopladas.]] </li>
</li>
</ul></div></center>

A lei de forças que rege o comportamento deste sistema é:

<math> F = -k \Delta x - \alpha k \Delta x^2 - \beta k \Delta x^3 </math>.

Onde <math> \Delta x </math> e a deformação a cada 2 massas acopladas (<math> x_{i+1} - x_i </math>), <math> k </math> é a constante elástica da mola, <math> \alpha </math> é um parâmetro de deformação arbitrário que controla a correção não linear quadrática e <math> \beta </math> é o parâmetro que controla a correção cúbica. Importante ressaltar que se <math> \alpha </math> é possuir assumir um valor não nulo, real, <math> \beta </math> é igual a zero no nosso sistema, ou vice-versa. Não estamos analisando correções quadráticas somadas com correções cúbicas neste trabalho.

=== Motivação: O paradoxo Enrico Fermi, John R. Pasta, Stanislaw M. Ulam e Mary Tsingou<ref name= FPUT1> https://www.scielo.br/j/rbef/a/SkRCy5fdnGbhfxNjpx5BkRD/?format=pdf&lang=pt - O problema Fermi-Pasta-Ulam-Tsingou: Equiparticão de
energia vista através de simulações computacionais </ref> ===

A premissa inicial do paradoxo de FPUT consiste no Teorema da Equipartição de Energia. O sistema consistia em uma corrente de partículas, com as extremidades fixas, que interagiam entre seus vizinhos somente com forças elásticas (as forças teriam um termo linear como a Força de Hooke e mais um termo não-linear, podendo ser quadrático ou cúbico). Era esperado que a energia total fosse distribuída igualmente entre as partículas. No caso em questão, a distribuição de energia entre as partículas pode ser descrita através dos seus modos normais de vibração.

A análise do problema gerou um paradoxo que começaria a ser respondido somente 10 anos depois, o que ajudou no desenvolvimento das teorias de sólitons e do caos. Pretendia-se observar a distribuição uniforme de energia entre os diversos modos normais de vibração com o passar do tempo (ao longo das iterações da simulação computacional). Isso significaria que o sistema alcançou um equilíbrio térmico e seria uma exemplificação computacional do Teorema de Equipartição de Energia. Caso as forças entre as partículas fossem estritamente lineares, isso não ocorreria, pois a energia alocada em cada modo não conseguiria acessar outros modos. Imaginava-se que uma componente não-linear na força tornaria acessível qualquer modo de vibração, porém não foi o observado.

A princípio, foi observada a tendência do sistema de distribuir a energia. O primeiro modo de vibração antes estimulado, perdeu energia ao longo do tempo, a qual começou a se alocar nos modos de energia mais baixos. Entretanto, por um descuido, deixaram a simulação decorrer por um tempo maior do que era planejado. Ao retornar ao laboratório para corrigir tal erro, se depararam com um resultado inesperado. A energia, que supostamente deveria estar igualmente partilhada entre os modos de vibração, estava quase completamente alocada no primeiro modo de vibração. De fato, somente 3% da energia não estava presente no primeiro modo. Devido a esta observação, deixaram a simulação correr por ainda mais tempo. Notaram, então, que existia
um ciclo, no qual a energia saía do primeiro modo de vibração, começava a se distribuir nos modos mais baixos, para, por fim, voltar quase que inteiramente para o primeiro modo de vibração. Contudo, em 2015, relatou-se que o sistema FPUT poderia atingir equipartição de energia pelo menos entre modos normais livres (interação entre três fônons).

== Discretização ==

A discretização deste problema gira em torno de abrir a equação das forças, e com o termo de aceleração, iterar o movimento das partículas a partir disso <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>. Partimos do problema com correção quadrática, ou seja, <math> \beta = 0 </math>. Partindo de:

<math> F = -k \Delta x - \alpha k \Delta x^2 </math>,

substituímos pelas variáveis discretas:

<math> m \ddot{x_j} = -k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) - \alpha k \left( (x_{j+1} - x_j) - (x_j - x_{j-1}) \right) </math>,

Chegamos em:

<math> m \ddot{x_j} = k \left( x_{j+1} - 2x_{j} + x_{j-1} \right) \left[ 1 + \alpha \left( x_{j+1} - x_{j-1} \right) \right] </math>,

Em que <math> \ddot{x_j} </math> é a aceleração da j-ésima partícula, com ela conseguimos integrar o movimento das partículas.

A Energia do sistema é calculada para cada ciclo de oscilação, porém para obtermos os resultados dos modos de oscilação e compararmos com os estudos atuais e o original de de FPUT, calculamos a energia dos primeiros modos de vibração da corda para demonstrar o comportamento visivelmente periódico destas energias. É possível calcular a energia dos modos de vibração através da equação:

<math> E_{total} = \sum_{i=2}^{N} \frac{\omega_{i}^{2}A_{i}^{2} + \dot{A_{i}^{2}}}{2} </math>,

onde <math> A_i </math> corresponde à:

<math> A_{i} = C\sin{\left( i\theta \right )} </math>

Conseguimos resolver <math>C</math> com a equação geral de <math>x(t)</math>:

<math> x_{i}(t) = \sum_{m=1}^{N} C_m \sin{\left (\frac{im\pi}{N+1} \right )}\cos{\left (wmt \right )} </math>,

e as freqûencias <math> \omega_i </math> são dadas por: <math> \omega_m = 2\omega_0 \sin\left (\frac{m\pi}{2(N + 1)}\right ) </math>.

<math>N</math> é o número de partículas e <math>\omega_0 = \frac{k}{massa} </math>

== Resultados ==

O sistema foi iniciado com o modo normal de oscilação 1 (nenhum nodo), com velocidades iniciais igual a zero de cada partícula. O movimento começa apenas pelas forças já presentes entre cada partículas, pelas equações apresentadas acima.

[[Arquivo:Estado inicial das particulas.png|center|thumb|500px|legenda.]]

NOSSO GIF DEU MAIS DE 7 MEGA -----------------> VER COMO REDUZIR O GIF JÁ NO CÓDIGO!!!! ELE TA BEM RUIM ASSIM

[[Arquivo:Trab2.gif|center|thumb|500px|legenda.]]

'''TEM QUE MUDAR AS FIGURAS DE BAIXO:'''

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Energias vibracao3.png|thumb|500px|legenda.]] </li>
<li style="display: inline-block;"> [[Arquivo:Energias vibracao2.png|thumb|500px|legenda.]] </li>
</ul></div></center>

CASO PRECISE:
<math>
\begin{cases}
X(x = 0,y) = X{0} = ..... \\
Y(x = L,y) = Y(x,y = 0) = Y(x,y = L) = .....\\
\end{cases}
</math>

== Discussões ==

Como o intuito era replicar os resultados através da simulação com dinâmica molecular obtivemos resultados muito parecidos comparando estudos já realizados sobre este problema ('''Citar estes artigos''')

== Implementação ==

Usamos XX partículas, com modo de oscilação YY

<source lang="haskell" line='line'>
### Exemplo da iteração do movimento utilizando forcacom correção quadrática ###
### código em python

def aceleracao(pos,alpha,k,massa):
size = len(pos)
acel = [0.0 for i in range(size)]
for i in range(1,size-1):
acel[i] = (k/massa) * ((pos[i+1] + pos[i-1] - 2*pos[i]) * ( 1.0 + alpha*(pos[i+1]-pos[i-1]) ) )
return acel

def velocidade(velo, acel, dt):
size = len(velo)
new_velo = [0.0 for i in range(size)]
for i in range(size):
new_velo[i] = velo[i] + 0.5*acel[i]*dt
return new_velo

def posicao(pos, velo, dt):
size = len(posY)
new_posY = [0.0 for i in range(size)]
for i in range(size):
new_posY[i] = posY[i] + new_veloY[i]*dt
return new_posY

N = número de partículas
dt = 0.2
x = np.linspace(0, x_final, dt)
pos = np.sin( 2*x*pi / (N*dt))

while t < tmax: # Loop temporal

plt.scatter(x,pos) # plotagem dos gráficos
acel = aceleracao(pos_old,alpha,k,massa)
velo = velocidade(veloY, new_aceY,dt)
pos = posicao(posY_old, veloY, dt)

pos_old = pos.copy() #sem ".copy()" o python usa o mesmo endereço de memória para 2 variáveis
t = t + td

gera_gif()
</source>

'''Falar um pouco mais da implementação?'''

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Arquivo:Trab2.gif

2021-05-26T23:32:54Z

Augustog: Augustog enviou uma nova versão de "Arquivo:Trab2.gif"

Arquivo:Trab2.gif

2021-05-26T23:32:19Z

Augustog: Augustog enviou uma nova versão de "Arquivo:Trab2.gif"

Arquivo:Trab2.gif

2021-05-26T23:32:11Z

Augustog: Augustog enviou uma nova versão de "Arquivo:Trab2.gif"

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref> sobre sistemas dinâmicos não lineares, abordando como a energia do sistema se comporta de acordo com os modos de vibração de uma corda.

== O Problema ==

drele <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. acessos em 23 maio 2021. Epub 10-Out-2016. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref>

dele

<math> teste = teste </math>

dele doli <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>

== TITULO 1 ==

<center><math> \nabla^2\Phi = 0 </math>.</center> equações

== TITULO 2 ==

=== SUBTITULOS ===

'''negrito''', ''Simultaneous OverRelaxation''

<math>
\begin{cases}
\Phi(x = 0,y) = \Phi_{0} = 1 \\
\Phi(x = L,y) = \Phi(x,y = 0) = \Phi(x,y = L) = 0\\
\end{cases}
</math>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Problema1-image.PNG|thumb|left|437px|Problema da borda carregada eletricamente.]] </li>
<li style="display: inline-block;"> [[Arquivo:Escorregador analitico.png|thumb|right|400px|Gráfico da solução analítica somando até o termo n=199.]] </li>
</ul></div></center>

<source lang="haskell" line='line'>
### Exemplo da evolução temporal no método de relaxação ###
### Exemplo para o algoritmo de jacobi, Equação de Laplace ###
# P é a matriz do potencial no tempo n
# Q é a matriz do potencial no tempo n+1

while t < tmax: # Loop temporal

for i in range(1,L+1): # Loop em x
for j in range(1,L+1): # Loop em y
Q[i][j] = (P[i+1][j] + P[i-1][j] + P[i][j+1] + P[i][j-1])/4

P = Q.copy()
t = t + td

plt.plot(x,y,P) # plotagem dos gráficos
</source>

[[Arquivo:Escorregador jacobi.png|center|thumb|500px|Solução numérica do problema da borda carregada.]]

<center><div><ul>
<li style="display: inline-block; vertical-align: top;"> [[Arquivo:Erro gauss 01.png|thumb|center|500px|Erro relativo médio para a solução de Gauss-Seidel para várias iterações.]] </li>
</ul></div></center>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam

2021-05-23T16:18:32Z

Augustog:

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é replicar os resultados do problema proposto por Fermi-Pasta-Ulam em 1953 <ref name= FPU> http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref> sobre sistemas dinâmicos não lineares, abordando como a energia do sistema se comporta de acordo com os modos de vibração de uma corda.

== O Problema ==

drele <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. acessos em 23 maio 2021. Epub 10-Out-2016. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref>

dele

dele doli <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>

== TITULO 1 ==

<center><math> \nabla^2\Phi = 0 </math>.</center> equações

== TITULO 2 ==

=== SUBTITULOS ===

'''negrito''', ''Simultaneous OverRelaxation''

<math>
\begin{cases}
\Phi(x = 0,y) = \Phi_{0} = 1 \\
\Phi(x = L,y) = \Phi(x,y = 0) = \Phi(x,y = L) = 0\\
\end{cases}
</math>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Problema1-image.PNG|thumb|left|437px|Problema da borda carregada eletricamente.]] </li>
<li style="display: inline-block;"> [[Arquivo:Escorregador analitico.png|thumb|right|400px|Gráfico da solução analítica somando até o termo n=199.]] </li>
</ul></div></center>

<source lang="haskell" line='line'>
### Exemplo da evolução temporal no método de relaxação ###
### Exemplo para o algoritmo de jacobi, Equação de Laplace ###
# P é a matriz do potencial no tempo n
# Q é a matriz do potencial no tempo n+1

while t < tmax: # Loop temporal

for i in range(1,L+1): # Loop em x
for j in range(1,L+1): # Loop em y
Q[i][j] = (P[i+1][j] + P[i-1][j] + P[i][j+1] + P[i][j-1])/4

P = Q.copy()
t = t + td

plt.plot(x,y,P) # plotagem dos gráficos
</source>

[[Arquivo:Escorregador jacobi.png|center|thumb|500px|Solução numérica do problema da borda carregada.]]

<center><div><ul>
<li style="display: inline-block; vertical-align: top;"> [[Arquivo:Erro gauss 01.png|thumb|center|500px|Erro relativo médio para a solução de Gauss-Seidel para várias iterações.]] </li>
</ul></div></center>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam

2021-05-23T16:12:26Z

Augustog: Criou página com ''''Grupo: Augusto M Giani e Henrique Padovani''' O objetivo deste trabalho é ... == Referências == xesque <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E....'

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é ...

== Referências ==

xesque <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. acessos em 23 maio 2021. Epub 10-Out-2016. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref>

dele <ref name=FPU>http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

dele doli <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>

== TITULO 1 ==

<center><math> \nabla^2\Phi = 0 </math>.</center> equações

== TITULO 2 ==

=== SUBTITULOS ===

'''negrito''', ''Simultaneous OverRelaxation''

<math>
\begin{cases}
\Phi(x = 0,y) = \Phi_{0} = 1 \\
\Phi(x = L,y) = \Phi(x,y = 0) = \Phi(x,y = L) = 0\\
\end{cases}
</math>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Problema1-image.PNG|thumb|left|437px|Problema da borda carregada eletricamente.]] </li>
<li style="display: inline-block;"> [[Arquivo:Escorregador analitico.png|thumb|right|400px|Gráfico da solução analítica somando até o termo n=199.]] </li>
</ul></div></center>

<source lang="haskell" line='line'>
### Exemplo da evolução temporal no método de relaxação ###
### Exemplo para o algoritmo de jacobi, Equação de Laplace ###
# P é a matriz do potencial no tempo n
# Q é a matriz do potencial no tempo n+1

while t < tmax: # Loop temporal

for i in range(1,L+1): # Loop em x
for j in range(1,L+1): # Loop em y
Q[i][j] = (P[i+1][j] + P[i-1][j] + P[i][j+1] + P[i][j-1])/4

P = Q.copy()
t = t + td

plt.plot(x,y,P) # plotagem dos gráficos
</source>

[[Arquivo:Escorregador jacobi.png|center|thumb|500px|Solução numérica do problema da borda carregada.]]

<center><div><ul>
<li style="display: inline-block; vertical-align: top;"> [[Arquivo:Erro gauss 01.png|thumb|center|500px|Erro relativo médio para a solução de Gauss-Seidel para várias iterações.]] </li>
</ul></div></center>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>

Problema de Fermi-Pasta-Ulam-Tsingou

2021-05-23T16:10:07Z

Augustog:

Trabalhos 2020/2

2021-05-23T16:07:14Z

Augustog: /* Problema de Fermi-Pasta-Ulam-Tsingou */

===[[Modelo de Keller-Segel para relação população-economia]] ===

===[[Equações de Laplace e Poisson]] ===

===[[Difusão ambipolar em plasmas]] ===

===[[ Equação de Cahn-Hilliard ]]===

===[[ Modelo de Fitzhugh-Nagumo para o potencial de ação em neurônios ]]===

===[[ Belousov-Zhabotinsky ]]===

===[[Modelo de Bornholdt para simulação de mercados financeiros artificiais]] ===

===[[Simulação de Micélio de Fungo]] ===

===[[Problema de Fermi-Pasta-Ulam]] ===

Problema de Fermi-Pasta-Ulam-Tsingou

2021-05-23T16:06:22Z

Augustog:

'''Grupo: Augusto M Giani e Henrique Padovani'''

O objetivo deste trabalho é ...

== Referências ==

xesque <ref name=andrade>ANDRADE, D. X.; ANJOS, P. H. R.; ASSIS, P. E. G.. Sobre a conexão entre alguns modelos físicos não-lineares. Rev. Bras. Ensino Fís., São Paulo , v. 39, n. 1, e1307, 2017 . Disponível em <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1806-11172017000100407&lng=pt&nrm=iso>. http://dx.doi.org/10.1590/1806-9126-rbef-2016-0083.</ref>

dele <ref name=FPU>http://www.physics.utah.edu/~detar/phys6720/handouts/fpu/FermiCollectedPapers1965.pdf - Fermi, Pasta, Ulam, Studies of non linear problems</ref>

dele doli <ref name=wiki>https://en.wikipedia.org/wiki/Fermi%E2%80%93Pasta%E2%80%93Ulam%E2%80%93Tsingou_problem</ref>

== TITULO 1 ==

<center><math> \nabla^2\Phi = 0 </math>.</center> equações

== TITULO 2 ==

=== SUBTITULOS ===

'''negrito''', ''Simultaneous OverRelaxation''

<math>
\begin{cases}
\Phi(x = 0,y) = \Phi_{0} = 1 \\
\Phi(x = L,y) = \Phi(x,y = 0) = \Phi(x,y = L) = 0\\
\end{cases}
</math>

<center><div><ul>
<li style="display: inline-block;"> [[Arquivo:Problema1-image.PNG|thumb|left|437px|Problema da borda carregada eletricamente.]] </li>
<li style="display: inline-block;"> [[Arquivo:Escorregador analitico.png|thumb|right|400px|Gráfico da solução analítica somando até o termo n=199.]] </li>
</ul></div></center>

<source lang="haskell" line='line'>
### Exemplo da evolução temporal no método de relaxação ###
### Exemplo para o algoritmo de jacobi, Equação de Laplace ###
# P é a matriz do potencial no tempo n
# Q é a matriz do potencial no tempo n+1

while t < tmax: # Loop temporal

for i in range(1,L+1): # Loop em x
for j in range(1,L+1): # Loop em y
Q[i][j] = (P[i+1][j] + P[i-1][j] + P[i][j+1] + P[i][j-1])/4

P = Q.copy()
t = t + td

plt.plot(x,y,P) # plotagem dos gráficos
</source>

[[Arquivo:Escorregador jacobi.png|center|thumb|500px|Solução numérica do problema da borda carregada.]]

<center><div><ul>
<li style="display: inline-block; vertical-align: top;"> [[Arquivo:Erro gauss 01.png|thumb|center|500px|Erro relativo médio para a solução de Gauss-Seidel para várias iterações.]] </li>
</ul></div></center>

== Link para Códigos ==

Fizemos no ambiente Colab em ''.ipynb'', segue link do github:[https://github.com/padovanih/equacao-de-laplace]

== Referências ==

<references/>