# Replay de la simulación del fluido monoatomico en equilibrio # # Christopher Paredes 2015 # Esta simulación permite calcular la función de distribución radial para longitudes mas largas. ### Definición de las variables en SI de unidades ### variable kb equal 1.3806504*10^(-23) # J/K variable T_lj equal 120 # K , variable N equal ${N_particulas} variable masa equal 6.69042650*10^(-26) # Kg variable rho equal 1374 # Kg/m^3 variable T equal 94.4 # K , Temperatura de la simulación. variable sigma equal 3.4*10^(-10) # A variable epsilon equal ${kb}*${T_lj} # J 1.65678048*10^(-21) variable cutoff equal 2.25*${sigma} # m variable skin equal 0.5*${sigma} # m variable L equal 6 # variable escala equal (10.229/${L})*${sigma}# m variable dt equal 10^(-14) # s variable tDump equal 20*dt variable intDump equal ${N_pasos}/${N_configuraciones} ### Geometría del problema y condiciones iniciales ### units si atom_style atomic dimension 3 boundary p p p lattice fcc ${escala} region caja_simulacion block 0 ${L} 0 ${L} 0 ${L}# Red fcc compuesta de 216 celdas unitarias con vacancia para 4 átomos cada una. create_box 1 caja_simulacion create_atoms 1 box mass 1 ${masa} timestep ${dt} #> Computo del MSD compute MSD2 all msd fix result all print ${intDump} "$(step) $(c_MSD2[4])" file "./datos/dump/msd_rerun.dump" screen no timestep ${dt} thermo ${intDump} thermo_style custom step temp etotal ke pe epair rerun ${f_configuraciones} dump x y z vx vy vz