arrumando_visc

anotações.txt

@@ -35,3 +35,16 @@ se aglutinarem mais ou menos.
 
 Ek =  T'* epsilon 
 T' = Ek/epsilon (numero grande) = Ek'
+
+
+Raios atômicos:
+https://en.wikipedia.org/wiki/Atomic_radii_of_the_elements_(data_page)
+para o carbono usar uma média entre sigma e raio = 0.335, 0.170 (van der Waals)
+https://www.sciencedirect.com/topics/chemistry/rhombohedral-graphite
+
+Usar o fator de empacotamento (em geral de 0.64 a 0.74 BCC, FCC, HCF)
+
+Ideias: 
+
+Variar densidade (número)
+Variar Temperatura

lennard.f90

@@ -3099,6 +3099,7 @@ program main
                     end do
                     ! print*, "AAA"
                     call vec2csv(rFUp,N,5,'rF_u_P',j,t,nimpre,start)
+
 
                 end if 
                 ! print*, "B"

part_gen_pequeno.py

@@ -24,13 +24,13 @@
 import matplotlib.pyplot as plt
 from random import randint
 
-ff = 4.0/1 # formato da região ff = x/y
-Lx = 2000
-arquivo2 = 'bin_peq2_5x1.csv'
-arquivo1 = 'bin_peq1_5x1.csv'
+ff = 1 # formato da região ff = x/y
+Lx = 400
+arquivo2 = 'poucas1.csv'
+arquivo1 = 'poucas2.csv'
 
-ratio = 0.8
-N1 = 5000
+ratio = 0
+N1 = 500
 N2 = 100
 
 L1 = Lx
@@ -104,6 +104,7 @@
 print('Ntot = {}'.format(Ntot))
 plt.scatter(p1[:,0],p1[:,1], label='pequeno2')
 
+plt.legend()
 plt.show()
 
 

settings.ini

@@ -3,18 +3,18 @@
 # Os outros parâmetros são adimensionais. Checar anotações.txt T' = T*k_b/epsilon
 # Argonio: epsilon/kb = 120K, sigma = 0.341nm 
 [global]
-    nimpre  = 200 # quntidade de saidas 
-    N = 5035 # número de partículas 
-    Ntype = 2 # número de tipos de partícula
-    t_fim = 400  # tempo de execução 
+    nimpre  = 5000 # quntidade de saidas 
+    N = 2000 # número de partículas 
+    Ntype = 1 # número de tipos de partícula
+    t_fim = 50  # tempo de execução 
     nimpre_init = 0 # começa a simular a partir no tempo referente ao snapshot (nimpre). 
     dt = 0.0001 # passo de tempo
     dimX = 2000 # Dimensões da região de cálculo
     dimY = 500 #
     mesh = 150 30 # malha(elementos)
     rcut = 3 # *sigma
     wall = 'ppee' # condição Elastica ou Periodica nas paredes norte, sul, leste, oeste
-    Td = -1 # 2.897189213660259e+24 #temperatura desejada global constante (-1 = não aplicar) kb = 1.38064852E-23, Td = (2/(3*kb))*E/epsilon(literatura)
+    Td = -1 #temperatura desejada global constante (-1 = não aplicar) kb = 1.38064852E-23, Td = (2/(3*kb))*E/epsilon(literatura)
     temp_Td = 0 100 100 # Tempo de velocity scaling [inicial, final, iterações_por_scaling] 
     cold_cells = 1 5 1 20  # [cellx_ini, cellx_fin, celly_ini, celly_fin]
     hot_cells = 95 100 1 20 # [cellx_ini, cellx_fin, celly_ini, celly_fin]
@@ -23,34 +23,34 @@
     vd = 2 2 # velocidade distribuida com Maxwell Boltzmann = sqrt(Td') 
     NMPT = 1000 # Número máximo de partículas que poderá mudar de process. Se não souber estimar, usar -1 ou o valor de N.
     GField = 0 0 # Campo de força gravitacional que afeta todas as partículas
-    print_TC = 0 # imprimir dados para calcular coeficiente de transporte
+    print_TC = 1 # imprimir dados para calcular coeficiente de transporte
 # PARTICLE PROPERTIES
 # v é velocidade global, entrar dois número correspondentes ao vetor vx vy que serão aplicadas em todas as partículas
 # se v_file é dada, então ela será usada para dar a velocidade às partículas. Deve ser arquivo .csv. Deixar "%" se não usar. 
 [par_0]
-    x = 'bin_peq1_5x1.csv' 
+    x = 'gas.csv' 
     v = 0 0 # velocidade global
     v_file = '%v_file_0.csv' # velocidade de cada partícula
     m = 1
     nome = 'g1'
     sigma = 1 # Lennard Jones
     epsilon = 1 # Lennard Jones
-    quantidade = 4935
+    quantidade = 2000
     x_lockdelay = 0 # só vai poder mudar de posição a partir de t = x_lockdelay
     rs = 0 # raio sólido. Posição de partículas na superfície de um sólido de raio rs
     fric_term = 0 # fricção artifical				
-[par_1]
-    x = 'bin_peq2_5x1.csv' 
-    v = 0 0 # velocidade global
-    v_file = '%v_file_0.csv' # velocidade de cada partícula
-    m = 1
-    nome = 'g1'
-    sigma = 1 # Lennard Jones
-    epsilon = 1 # Lennard Jones
-    quantidade = 100
-    x_lockdelay = 0 # só vai poder mudar de posição a partir de t = x_lockdelay
-    rs = 1 # raio sólido. Posição de partículas na superfície de um sólido de raio rs
-    fric_term = 0 # fricção artifical				
+#[par_1]
+#    x = 'bin_peq2_5x1.csv' 
+#    v = 0 0 # velocidade global
+#    v_file = '%v_file_0.csv' # velocidade de cada partícula
+#    m = 1
+#    nome = 'g1'
+#    sigma = 1 # Lennard Jones
+#    epsilon = 1 # Lennard Jones
+#    quantidade = 100
+#    x_lockdelay = 0 # só vai poder mudar de posição a partir de t = x_lockdelay
+#    rs = 1 # raio sólido. Posição de partículas na superfície de um sólido de raio rs
+#    fric_term = 0 # fricção artifical				
 #[par_2]
 #    x = 'p_g.csv'
 #    v = 0 0