-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsettings.ini
executable file
·77 lines (77 loc) · 4.51 KB
/
settings.ini
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
# PARAMETROS GLOBAIS
# Convenção: O primeiro grupo de partículas (part_0) deve ser o de referência (m = epsilon = sigma = 1)
# 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 = 1009 # número de partículas
Ntype = 2 # 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.00005 # passo de tempo
dimX = 100 # Dimensões da região de cálculo
dimY = 100 #
mesh = 10 10 # 300 60 # malha(elementos)
rcut = 3 # *sigma
wall = 'eeee' # condição Elastica ou Periodica nas paredes norte, sul, leste, oeste
termostato_nose_hoover = .false. #liga termostato Nose Hoover. Este funciona durante toda execução.
termostato_vel_scaling = .false. # liga termostato por velocity scaling
Mc = 1
Td = -1 #temperatura desejada global constante (-1 = não aplicar) kb = 1.38064852E-23, Td = (2/(3*kb))*E/epsilon(literatura)
temp_Td = 0 10000 2000 # Tempo de velocity scaling [inicial, final, iterações_por_scaling]
cold_cells = 1 12 1 30 # [cellx_ini, cellx_fin, celly_ini, celly_fin]
hot_cells = 190 200 1 30 # [cellx_ini, cellx_fin, celly_ini, celly_fin]
Td_hot = 2 # temperatura nas cold cells constante. Td_hot = -1 desliga termostato
Td_cold = 0.8 # v = sqrt(2*T) temperatura nas hot cells constante. Td_cold = -1 desliga termostato
force_lingradT = -1 # (NUMERO MULTIPLO DE 4 RECOMENDADO!) se > 0 irá dividir região no numero de subregiões entradas na direação X e será forçado um gradiente de temperatura linear. A região dos banhos será ignorada. A temperatura será controlado por um termostato velocity scaling
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
MField = 0 0 0 0 0 # Campo magnético vezes momento dipolo da particula magnética Hm. O primeiro numero número é a intensidade de m*H. O segundo e terceiro numeros são a direção do campo.. Os últimos dois números são a e b na função seno e heaviside de controlam o campo no tempo. Hm = H*m*sin(b*t)*heaviside(x-a)
print_TC = .false. # 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.
# se ismolecule = .true. então a partícula é levada em conta no termostato
# pr = constantes [A, B, alpha, beta] do lennard jones rugoso que simulará um cristal de atomos. Se o valor for dado, a particula poderá girar.
# o grupo com pr deve ser o último e deverá haver apenas um. Ver em rot_par.py como será o potencial
[par_0]
x = 'molp.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 = 1000
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 não funciona quando o termostato nose hoover está ligado
ismolecule = .true.
[par_1]
x = 'parp.csv'
v = 0 0 # velocidade global
v_file = '%%v_file_1.csv' # velocidade de cada partícula
m = 10
nome = 'g2'
sigma = 1 # Lennard Jones
epsilon = 1 # Lennard Jones
quantidade = 9
x_lockdelay = 0 # só vai poder mudar de posição a partir de t = x_lockdelay
rs = 5 # raio sólido. Posição de partículas na superfície de um sólido de raio rs
fric_term = 0 # fricção artifical
ismolecule = .false.
pr = 0 0.001 20 20 0.071847 # constantes [A, B, alpha, beta, ph] do lennard jones rugoso que simulará um cristal de atomos 43.72584 43.72584
#
#[par_2]
# x = 'p_g.csv'
# v = 0 0
# v_file = "%vp_g.csv"
# m = 20
# nome = 'g3'
# sigma = 1 # Lennard Jones
# epsilon = 1 # Lennard Jones
# quantidade = 20
# x_lockdelay = 300 # só vai poder mudar de posição a partir de t = x_lockdelay
# rs = 5
# fric_term = 0 # fricção artifical