-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathdefinitions_small_unstable.h
109 lines (95 loc) · 4.22 KB
/
definitions_small_unstable.h
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
#define PHYSICS MHD
#define DIMENSIONS 2
#define COMPONENTS 3
#define GEOMETRY CYLINDRICAL
#define BODY_FORCE NO
#define COOLING NO
#define RECONSTRUCTION LINEAR
#define TIME_STEPPING RK2
#define DIMENSIONAL_SPLITTING YES
#define NTRACER 0
#define USER_DEF_PARAMETERS 9
/* -- physics dependent declarations -- */
#define EOS PVTE_LAW
#define ENTROPY_SWITCH NO
#define DIVB_CONTROL NO
#define BACKGROUND_FIELD NO
#define RESISTIVITY ALTERNATING_DIRECTION_IMPLICIT
#define THERMAL_CONDUCTION NO
#define VISCOSITY NO
#define ROTATING_FRAME NO
/* -- user-defined parameters (labels) -- */
#define ETAX_GAU 0
#define ETAY_GAU 1
#define ETAZ_GAU 2
#define KAPPA_GAUBOB 3
#define TWALL 4
#define T0 5
#define DENS0 6
#define VZ0 7
#define ALPHA_J 8
/* [Beg] user-defined constants (do not change this line) */
#define UNIT_DENSITY 1.67382e-6
#define UNIT_LENGTH 1.0e-3
#define UNIT_VELOCITY 1.0e6
#define TC_SATURATED_FLUX NO
/* [End] user-defined constants (do not change this line) */
/* -- supplementary constants (user editable) -- */
#define INITIAL_SMOOTHING NO
#define WARNING_MESSAGES YES
#define PRINT_TO_FILE YES
#define INTERNAL_BOUNDARY YES
#define SHOCK_FLATTENING NO
#define CHAR_LIMITING NO
#define LIMITER DEFAULT
#define ASSIGN_VECTOR_POTENTIAL NO
#define UPDATE_VECTOR_POTENTIAL NO
#define SHOW_TIME_STEPS YES
/* ---------------------------------------------------- */
/* Additional constants */
#define UNIT_ETA (4*CONST_PI/(CONST_c*CONST_c)*UNIT_VELOCITY*UNIT_LENGTH)
#define UNIT_KAPPA (UNIT_DENSITY*UNIT_VELOCITY*UNIT_LENGTH*CONST_kB/CONST_mp)
/* Ema's additional macros */
#define SPLIT_DIFF_ADV_ADV_DIFF
// #define PROFILE_GPROF_STOPSTEP 8
// #define FREEZE_FLUID
/* Either give value in ]0,0.5[
or comment if you prefer to use Peaceman-Rachford scheme */
// #define FRACTIONAL_THETA 0.3
// #define SPLIT_IMPLICIT
/* For a pseudo P-R algoritm: if FRACT==0.5 you have the usual P-R
otherwise you unbalance the scheme towards the implcit or explicit part
(keep it in ]0,1[). If you do not define it it will be set to 0.5*/
#define FRACT 0.5
// to set the order of the ADI scheme, allowed values: YES, NO, RANDOM, PERMUTE.
#define FIRST_JDIR_THEN_IDIR NO
// #define TEST_ADI
#define JOULE_EFFECT_AND_MAG_ENG (YES && RESISTIVITY==ALTERNATING_DIRECTION_IMPLICIT)
//Keep it YES for now. If YES: power flux is computed inside adi schemes (if NO, outside)
#define POW_INSIDE_ADI YES
/* Macros to impose T (B) on walls also for advection (unphisical!)
(if NO, conduction and B diffusion can be modeled only via ADI scheme) */
#define IMPOSE_TWALL NO
#define IMPOSE_BWALL NO
/* Number of subcycles performed by ADI scheme*/
#define NSUBS_ADI 1
// #define COMMON_RATIO_NSUBS_ADI 1.5
/* Decide whether the electrode must be set as a hom-Neumann boundary*/
#define ELECTR_B_NEUM
/* To set to 0 the mag field in a region outside capillary*/
// #define FLATTEN_B_OUTCAP
// #define PRESUBS_RES
// #define ELECTR_T_NEUM
/* ---------------------------------------------------- */
/* ---------------------------------------------------- */
/* CAPILLARY GEOMETRY SETTINGS */
#define RCAP 0.07
#define DZCAP 0.2 /*the electrodes are wide DZCAP cm*/
#define ZCAP 0.8 /*the capillary is long 2*ZCAP cm and wide 2*RCAP cm*/
/* ---------------------------------------------------- */
/* ---------------------------------------------------- */
/* ADDITTIONAL OUTPUT */
#define WRITE_J YES
/*------------------------------------------------------*/
/* FOR DEBUG */
#define DEBUG_EMA