-
Notifications
You must be signed in to change notification settings - Fork 1
/
StableABL_precursor2_hypre.inp
158 lines (149 loc) · 9.21 KB
/
StableABL_precursor2_hypre.inp
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
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
# --- Simulation time control parameters ---
time.stop_time = 30000.0 # Max (simulated) time to evolve [s]
time.max_step = -1
time.fixed_dt = 0.125 # Fixed timestep size (in seconds). If negative, then time.cfl is used
time.checkpoint_interval = 8000
incflo.physics = ABL # List of physics models to include in simulation.
incflo.verbose = 3
io.check_file = chk
incflo.use_godunov = true
incflo.godunov_type = weno_z
turbulence.model = OneEqKsgsM84
TKE.source_terms = KsgsM84Src
nodal_proj.mg_rtol = 1e-06
nodal_proj.mg_atol = 1e-12
#mac_proj.mg_rtol = 1e-06
#mac_proj.mg_atol = 1e-12
diffusion.mg_rtol = 1e-06
diffusion.mg_atol = 1e-12
temperature_diffusion.mg_rtol = 1e-10
temperature_diffusion.mg_atol = 1e-13
incflo.gravity = 0.0 0.0 -9.81 # Gravitational acceleration vector (x,y,z) [m/s^2]
incflo.density = 1.17 # Fluid density [kg/m^3]
transport.viscosity = 1.853e-05 # Fluid dynamic viscosity [kg/m-s]
transport.laminar_prandtl = 0.7 # Laminar prandtl number
transport.turbulent_prandtl = 0.3333 # Turbulent prandtl number
# --- Geometry and Mesh ---
geometry.prob_lo = 609380.0 3999000.0 0.0
geometry.prob_hi = 651620.0 4047000.0 800.0
amr.n_cell = 4224 4800 80 # Number of cells in x, y, and z directions
amr.max_level = 2
geometry.is_periodic = 1 1 0
zlo.type = wall_model
zlo.temperature_type = wall_model
zlo.tke_type = zero_gradient
zhi.type = slip_wall
zhi.temperature_type = fixed_gradient
zhi.temperature = 0.00125
# --- ABL parameters ---
ICNS.source_terms = ABLForcing BoussinesqBuoyancy CoriolisForcing ABLMeanBoussinesq
ABL.stats_output_frequency = 1
ABL.stats_output_format = netcdf
incflo.velocity = -0.8759152146139654 10.011756715822044 0.0
ABLForcing.abl_forcing_height = 91.0
ABL.kappa = 0.4
ABL.normal_direction = 2
ABL.surface_roughness_z0 = 0.1
ABL.reference_temperature = 303.803
ABL.surface_temp_rate = -0.07
ABL.mo_beta_m = 16.0 # Monin-Obukhov Beta m parameter
ABL.mo_gamma_m = 5.0 # Monin-Obukhov Gamma m parameter
ABL.mo_gamma_h = 5.0 # Monin-Obukhov Gamma h parameter
ABL.random_gauss_mean = 0.0
ABL.random_gauss_var = 1.0
CoriolisForcing.latitude = 36.57236455833334
CoriolisForcing.rotational_time_period = 86164.0900027328
CoriolisForcing.north_vector = 0.0 1.0 0.0
CoriolisForcing.east_vector = 1.0 0.0 0.0
BoussinesqBuoyancy.reference_temperature = 303.803
ABL.temperature_heights = 0 500.0 600.0 800.0
ABL.temperature_values = 305.803 305.803 313.803 314.053
ABL.perturb_velocity = true
ABL.perturb_ref_height = 50.0
ABL.Uperiods = 4.0
ABL.Vperiods = 4.0
ABL.deltaU = 1.0
ABL.deltaV = 1.0
ABL.perturb_temperature = false
time.plot_interval = 8000
io.plot_file = plt
io.KE_int = -1
#---- tagging defs ----
tagging.labels = Farm_level_0_center_630500_4023000_100 KP AF BR
tagging.Farm_level_0_center_630500_4023000_100.type = GeometryRefinement
tagging.Farm_level_0_center_630500_4023000_100.shapes = Farm_level_0_center_630500_4023000_100
tagging.Farm_level_0_center_630500_4023000_100.level = 0
tagging.Farm_level_0_center_630500_4023000_100.Farm_level_0_center_630500_4023000_100.type = box
tagging.Farm_level_0_center_630500_4023000_100.Farm_level_0_center_630500_4023000_100.origin = 650132.351164615 4010864.8944071657 0.0
tagging.Farm_level_0_center_630500_4023000_100.Farm_level_0_center_630500_4023000_100.xaxis = -2510.085391132557 28690.40730504227 0.0
tagging.Farm_level_0_center_630500_4023000_100.Farm_level_0_center_630500_4023000_100.yaxis = -36859.20382939459 -3224.7624816633547 0.0
tagging.Farm_level_0_center_630500_4023000_100.Farm_level_0_center_630500_4023000_100.zaxis = 0.0 0.0 450.0
tagging.KP.type = GeometryRefinement
tagging.KP.shapes = KP
tagging.KP.level = 1
tagging.KP.KP.type = box
tagging.KP.KP.origin = 648532.2646106855 4024309.1361471326 0.0
tagging.KP.KP.xaxis = -958.7131702242406 10958.1416790092 0.0
tagging.KP.KP.yaxis = -17931.50456565142 -1568.8033694578482 0.0
tagging.KP.KP.zaxis = 0.0 0.0 300.0
tagging.AF.type = GeometryRefinement
tagging.AF.shapes = AF
tagging.AF.level = 1
tagging.AF.AF.type = box
tagging.AF.AF.origin = 636732.7228190051 4010017.8883415214 0.0
tagging.AF.AF.xaxis = -1551.3722209083166 17732.26562603307 0.0
tagging.AF.AF.yaxis = -17931.50456565142 -1568.8033694578482 0.0
tagging.AF.AF.zaxis = 0.0 0.0 300.0
tagging.BR.type = GeometryRefinement
tagging.BR.shapes = BR
tagging.BR.level = 1
tagging.BR.BR.type = box
tagging.BR.BR.origin = 625297.7417857626 4031880.5919449204 0.0
tagging.BR.BR.xaxis = -488.07215938688614 5578.690309313775 0.0
tagging.BR.BR.yaxis = -12950.531075192692 -1133.0246557195571 0.0
tagging.BR.BR.zaxis = 0.0 0.0 300.0
#---- postprocessing defs ----
incflo.post_processing = sampling
sampling.type = Sampling
sampling.output_frequency = 2000
sampling.fields = velocity
#---- sample defs ----
sampling.labels = p_hub
sampling.p_hub.type = PlaneSampler
sampling.p_hub.num_points = 421 481
sampling.p_hub.origin = 609380.0 3999000.0 0.0
sampling.p_hub.axis1 = 42240.0 0.0 0.0
sampling.p_hub.axis2 = 0.0 48000.0 0.0
sampling.p_hub.normal = 0.0 0.0 1.0
sampling.p_hub.offsets = 60 91
#---- extra params ----
amr.blocking_factor = 16
amr.max_grid_size = 64
mac_proj.mg_rtol = 1.0e-11
mac_proj.mg_atol = 1.0e-9
mac_proj.do_semicoarsening = true
mac_proj.bottom_solver = hypre
mac_proj.bottom_verbose = 0
mac_proj.max_coarsening_level = 0
mac_proj.bottom_rtol = 1.0e-12
mac_proj.bottom_atol = 1.0e-12
hypre.hypre_solver = GMRES
hypre.hypre_preconditioner = BoomerAMG
hypre.verbose = 0
hypre.bamg_verbose = 0
hypre.num_krylov = 20
hypre.bamg_max_levels = 5
hypre.bamg_num_sweeps = 1
#hypre.bamg_relax_type = 11
hypre.bamg_coarsen_type = 8
hypre.bamg_interp_type = 6
hypre.bamg_cycle_type = 1
hypre.bamg_relax_order = 0
hypre.bamg_trunc_factor = 0.25
hypre.bamg_agg_num_levels = 1
hypre.bamg_agg_interp_type = 7
hypre.bamg_agg_pmax_elmts = 3
hypre.bamg_pmax_elmts = 3
hypre.bamg_strong_threshold = 0.57
hypre.recompute_preconditioner = 0
#== END AMR-WIND INPUT ==