Tópico 17 · Script LIGGGHTS
4/4
Dissecação do script: tempo, saída, execução e rotação
8. Gerenciamento do tempo (integração de Verlet)
timestep 0.00001
thermo_style custom step atoms ke
thermo 1000
thermo_modify lost ignore norm no
- timestep define o passo de tempo da integração (associado ao Tempo de Rayleigh). 0.00001 passo de 10 microssegundos (Δt), usado para prever o deslocamento a cada iteração.
- thermo_style custom formato próprio para os dados exibidos no terminal/log durante a simulação. step imprime o número do passo atual. atoms imprime a quantidade de partículas presentes. ke imprime a energia cinética global.
- thermo 1000 envia os dados definidos em thermo_style a cada 1000 passos, monitorando a simulação sem sobrecarregar a saída.
- thermo_modify customiza os dados de saída globais. lost ignore impede erro fatal quando partículas deixam o domínio (fronteiras fixas); o cálculo continua ignorando as perdidas. norm no desativa a normalização dos valores; saídas ficam em termos brutos (não divididos pela quantidade de partículas).
dump dmp all custom 2000 post/dump*.liggghts id type x y z vx vy vz radius
dump dmpstl all stl 2000 post/tambor*.stl
- dump dmp registra o estado das partículas para visualização (ex. ParaView); "dmp" é o nome do usuário. all inclui todas as partículas. custom formato personalizado, colunas listadas ao final.
- 2000 grava a cada 2000 passos. post/dump*.liggghts diretório e nome de saída; o curinga * é substituído pelo número do passo (ex. dump100000.liggghts).
- id type x y z vx vy vz radius colunas exportadas — identificador, tipo do material, coordenadas, velocidades e raio.
- dump dmpstl grava a geometria (malha) em formato STL. stl formato de saída poligonal, registrando a rotação da malha ao longo do tempo.
- post/tambor*.stl diretório e nome de saída, curinga substituído pelo número do passo.
10. Execução da simulação
run 100000
unfix ins
fix mover_tambor all move/mesh mesh cad1 rotate origin 0. 0. 0. axis 0. 1. 0. period 0.4188
run 500000
- run 100000 aciona o solucionador de integração, executando o ciclo de Newton e Verlet por 100.000 passos; com timestep de 0,00001s, corresponde a 1 segundo simulado — tempo para a chuva de partículas assentar antes de iniciar a rotação.
- unfix ins remove a operação de inserção ("ins") do ciclo, interrompendo a criação de novas partículas.
- mover_tambor nome do usuário para o movimento, adicionado após o assentamento do leito. move/mesh mesh cad1 impõe movimento prescrito à malha "cad1" (o tambor). rotate o movimento é rotacional.
- origin 0. 0. 0. ponto por onde passa o eixo de rotação (origem do domínio). axis 0. 1. 0. vetor do eixo de rotação (eixo Y). period 0.4188 tempo (s) para uma volta completa (2π rad), correspondendo a ~15 rad/s.
- run 500000 executa mais 500.000 passos com o tambor girando; com timestep 0,00001s, equivale a 5 segundos simulados de mistura/avalanche. Ao final, os arquivos do dump ficam disponíveis para leitura no ParaView.