-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrun_AWBS_analysis.sh
executable file
·123 lines (110 loc) · 5 KB
/
run_AWBS_analysis.sh
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
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
#! /bin/bash
NPROC=2
## Working Efield general diffusive regime.
#RS=7 # Increasing spatial resolution improves the stability the best!
#MINGSAFE=2000
#F1ORDER=3
#F0ORDER=2
#SIGMA=1e15 # Kn 1e-5, safe for Z=2
##SIGMA=1e16 # Kn 1e-6, safe for Z=4
#SIGMASAFE=1e18 # Kn 1e-8, safe for any Z
RS=7
MINGSAFE=2000
F1ORDER=3
F0ORDER=2
NI=1e20
TMAX=1000
TMIN=100
TGRAD=180
XPOINT=0.046775 # in cm
SIGMA=1.2e12 # Zbar = 1 -> Kn 1e-2
SIGMASAFE=1.2e18 # Zbar = 1 -> Kn 1e-8
ZBAR=1
#ZBAR=2
#ZBAR=5
#ZBAR=10
#ZBAR=20
#ZBAR=50
#ZBAR=100
L=0.1
PROBLEM=5
#if [false] ; then
### Pascal's setting for nonlocal test
## First, diffusive, case sets sigma 1e5x higher, which assures SH solution.
## C7E
mpirun -np $NPROC C7 -p $PROBLEM -m data/segment01.mesh -rs $RS -tf 0.0 -ok $F1ORDER -ot $F0ORDER -no-vis -fa -print -Tmax $TMAX -Tmin $TMIN -sigma $SIGMASAFE -Tgrad $TGRAD -Z $ZBAR -ni $NI -L $L -xp $XPOINT -minG $MINGSAFE -S0 1.0 -E0 1.0
cp results/tmp/C7_1_profiles.* results/fe_analysis/Ecorrect_data/
cp results/tmp/C7_1_fe_point.txt results/fe_analysis/Ecorrect_data/fe_point_Ecorrect.txt
## C7*
mpirun -np $NPROC C7 -p $PROBLEM -m data/segment01.mesh -rs $RS -tf 0.0 -ok $F1ORDER -ot $F0ORDER -no-vis -fa -print -Tmax $TMAX -Tmin $TMIN -sigma $SIGMASAFE -Tgrad $TGRAD -Z $ZBAR -ni $NI -L $L -xp $XPOINT -minG $MINGSAFE
cp results/tmp/C7_1_profiles.* results/fe_analysis/Emimic_data/
cp results/tmp/C7_1_fe_point.txt results/fe_analysis/Emimic_data/fe_point_Emimic.txt
cd results/fe_analysis
python C7_AWBS_SH_analysis.py -N $NPROC -Z $ZBAR -s $SIGMASAFE -n $NI -xp $XPOINT --Emimic --Ecorrect --AWBSoriginal
cd ../..
#fi
#mpirun -np $NPROC C7 -p $PROBLEM -m data/segment01.mesh -rs $RS -tf 0.0 -ok $F1ORDER -ot $F0ORDER -no-vis -fa -print -Tmax $TMAX -Tmin $TMIN -sigma $SIGMA -Tgrad $TGRAD -Z $ZBAR -ni $NI -L $L -xp $XPOINT -minG $MINGSAFE -S0 1.0 -E0 1.0
#cd results/fe_analysis
#cp ../tmp/C7_1_fe_point.txt fe_point_Emimic.txt
#python C7_AWBS_SH_analysis.py -N $NPROC -Z $ZBAR -s $SIGMA -n $NI -xp $XPOINT --Ecorrect
#cd ../..
#if [false] ; then
## Nonlocal solution very well corresponding to Pascal's solution with Aladin.
## P1 closure.
## C7*
mpirun -np $NPROC C7 -p $PROBLEM -m data/segment01.mesh -rs $RS -tf 0.0 -ok $F1ORDER -ot $F0ORDER -no-vis -fa -print -Tmax $TMAX -Tmin $TMIN -sigma $SIGMA -Tgrad $TGRAD -Z $ZBAR -ni $NI -L $L -xp $XPOINT -minG $MINGSAFE
cp results/tmp/C7_1_profiles.* results/fe_analysis/Emimic_data/
cp results/tmp/C7_1_fe_point.txt results/fe_analysis/Emimic_data/fe_point_Emimic.txt
## M1 closure.
## C7*
mpirun -np $NPROC C7 -p $PROBLEM -m data/segment01.mesh -rs $RS -tf 0.0 -ok $F1ORDER -ot $F0ORDER -no-vis -fa -print -Tmax $TMAX -Tmin $TMIN -sigma $SIGMA -Tgrad $TGRAD -Z $ZBAR -ni $NI -L $L -xp $XPOINT -minG $MINGSAFE -M1
cp results/tmp/C7_1_profiles.* results/fe_analysis/Ecorrect_data/
cp results/tmp/C7_1_fe_point.txt results/fe_analysis/Ecorrect_data/fe_point_Ecorrect.txt
#cp results/tmp/C7_1_fe_point.txt results/fe_analysis/Emimic_data/fe_point_Emimic.txt
cd results/fe_analysis
python C7_AWBS_SH_analysis.py -N $NPROC -Z $ZBAR -s $SIGMA -n $NI -xp $XPOINT --Emimic --Ecorrect --labelEmimic C7P1 --labelEcorrect C7M1
cd ../..
#fi
## Diffusive asymptotic test. The SH, AWBS (original and corrected),
## and C7 (proper and mimic Efield) calculations are compared.
## In reality, the resulting Knudsen number is just on the diffusive limit
## in the following case.
#mpirun -np 8 C7 -p 8 -m data/segment01.mesh -rs 6 -tf 0.0 -ok 2 -ot 1 -no-vis -fa -print -Tmax 800.5 -Tmin 799.5 -Tgrad 2.3 -S0 1.0 -E0 1.0 -sigma 1e11 -Z 4 -ne 5e20 -M1 -minG 2000
## A pure diffusion case.
## Converged numerical flux from -minG 200. Err 1e-5 -minG 50.
RS=7
F1ORDER=4
F0ORDER=3
NI=5e20
ZBAR=100
MING=25
if [ $ZBAR -gt 10 ] ; then
MING=2000
fi
TMAX=800.5
TMIN=799.5
#TGRAD=2.3
#TGRAD=80
TGRAD=300
SIGMA=8.77e9 ## Kn=1.0e-4 nonlocality limit.
#SIGMA=8.77e8 ## Kn=1.0e-3 q=0.8*qSH E0=0.91
#SIGMA=1.75e8 ## Kn=5.0e-3 q=0.4*qSH E0=0.66
#SIGMA=8.45e15 ## Kn 1.0e-10
#SIGMA=8.45e11 ## Kn 1.0e-6
#SIGMA=1.69e11 ## Kn 5.0-6
#SIGMA=0.845e11 ## Kn 1.0-5
#SIGMA=0.535e11 ## Kn 1.6-5 qC7E zero
#SIGMA=0.845e10 ## Kn 1.0-4 qC7E negative
#ZBAR=100
#MING=$MINGSAFE
L=0.1
XPOINT=0.05
PROBLEM=5
mpirun -np $NPROC C7 -p $PROBLEM -m data/segment01.mesh -rs $RS -tf 0.0 -ok $F1ORDER -ot $F0ORDER -no-vis -fa -print -Tmax $TMAX -Tmin $TMIN -Tgrad $TGRAD -sigma $SIGMA -Z $ZBAR -ni $NI -L $L -xp $XPOINT -minG $MING
cp results/tmp/C7_1_profiles.* results/fe_analysis/Emimic_data/
cp results/tmp/C7_1_fe_point.txt results/fe_analysis/Emimic_data/fe_point_Emimic.txt
mpirun -np $NPROC C7 -p $PROBLEM -m data/segment01.mesh -rs $RS -tf 0.0 -ok $F1ORDER -ot $F0ORDER -no-vis -fa -print -Tmax $TMAX -Tmin $TMIN -Tgrad $TGRAD -S0 1.0 -E0 1.0 -sigma $SIGMA -Z $ZBAR -ni $NI -L $L -M1 -xp $XPOINT -minG $MING
cp results/tmp/C7_1_profiles.* results/fe_analysis/Ecorrect_data/
cp results/tmp/C7_1_fe_point.txt results/fe_analysis/Ecorrect_data/fe_point_Ecorrect.txt
cd results/fe_analysis
python C7_AWBS_SH_analysis.py -N $NPROC -Z $ZBAR -s $SIGMA -n $NI -xp $XPOINT --Ecorrect --Emimic --AWBSoriginal