-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathin.CS_3021_150km
103 lines (76 loc) · 3.76 KB
/
in.CS_3021_150km
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
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
# SPARTA input file for satellite CS_3021, for an altitude of 150.0km
print ""
print "***** Running SPARTA simulation for CS_3021, at h=150km *****"
print ""
seed 12345
dimension 3
global gridcut 0.1 comm/sort yes surfmax 10000 splitmax 1000
boundary o o o
create_box -0.5500 0.5500 -0.3500 0.3500 -0.3500 0.3500
create_grid 63 40 40
balance_grid rcb part
global nrho 3.0712e+15 fnum 5.3478e+10 vstream -3475.5100 0.0 0.0 temp 172.4040
species ../atmo.species CO2 N2 Ar CO O O2
mixture atmo CO2 frac 0.6583
mixture atmo N2 frac 0.1192
mixture atmo Ar frac 0.0491
mixture atmo CO frac 0.0780
mixture atmo O frac 0.0853
mixture atmo O2 frac 0.0101
collide vss atmo ../atmo.vss
read_surf ../data/data.CS_3021 trans 0.1500 0 0
surf_collide 1 diffuse 293.15 0.7721
surf_modify all collide 1
region sat block -0.5000 0.2500 -0.3500 0.3500 -0.3500 0.3500
fix in emit/face atmo xhi zhi zlo yhi ylo
timestep 3.7679e-06
compute forces surf all all fx fy fz
fix avg ave/surf all 1 33 50 c_forces[*] ave running
compute sum_force reduce sum f_avg[*]
compute n grid all all n
fix n_avg ave/grid all 1 33 50 c_n[*]
compute nrho grid all all nrho
fix nrho_avg ave/grid all 1 33 50 c_nrho[*]
compute massrho grid all all massrho
fix massrho_avg ave/grid all 1 33 50 c_massrho[*]
compute u grid all all u
fix u_avg ave/grid all 1 33 50 c_u[*]
compute avg_ppc reduce ave f_n_avg
compute T thermal/grid all all temp
fix T_avg ave/grid all 1 33 50 c_T[*]
compute knudsen lambda/grid f_nrho_avg f_T_avg CO2 kall
stats 50
stats_style step cpu wall np nscoll nexit c_sum_force[*] c_avg_ppc
dump 0 grid all 100 ../results_sparta/CS_3021/vals_150km_0.*.dat id f_n_avg f_nrho_avg f_massrho_avg f_u_avg f_T_avg c_knudsen[*]
write_grid ../results_sparta/CS_3021/grid_150km_0.dat
run 750
print "Refinement level 1"
timestep 3.8801e-06
adapt_grid all refine coarsen value c_knudsen[2] 5 50 combine min thresh less more cells 2 2 2
scale_particles all 5
global fnum 1.0696e+10
undump 0
dump 1 grid all 50 ../results_sparta/CS_3021/vals_150km_1.*.dat id f_n_avg f_nrho_avg f_massrho_avg f_u_avg f_T_avg c_knudsen[*]
balance_grid rcb part
write_grid ../results_sparta/CS_3021/grid_150km_1.dat
run 250
print "Refinement level 2"
timestep 1.9400e-06
adapt_grid all refine coarsen value c_knudsen[2] 5 50 combine min thresh less more cells 2 2 2 region sat one
scale_particles all 10
global fnum 1.0696e+09
undump 1
dump 2 grid all 50 ../results_sparta/CS_3021/vals_150km_2.*.dat id f_n_avg f_nrho_avg f_massrho_avg f_u_avg f_T_avg c_knudsen[*]
balance_grid rcb part
write_grid ../results_sparta/CS_3021/grid_150km_2.dat
run 200
print "Refinement level 3"
timestep 9.7002e-07
adapt_grid all refine coarsen value c_knudsen[2] 5 50 combine min thresh less more cells 2 2 2 region sat one
scale_particles all 5
global fnum 2.1391e+08
undump 2
dump 3 grid all 50 ../results_sparta/CS_3021/vals_150km_3.*.dat id f_n_avg f_nrho_avg f_massrho_avg f_u_avg f_T_avg c_knudsen[*]
balance_grid rcb part
write_grid ../results_sparta/CS_3021/grid_150km_3.dat
run 300