-
Notifications
You must be signed in to change notification settings - Fork 1
/
StableABL_precursor1.inp
128 lines (122 loc) · 8.44 KB
/
StableABL_precursor1.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
# --- 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 4005240.0 0.0
geometry.prob_hi = 651620.0 4040760.0 800.0
amr.n_cell = 4224 3552 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
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 351
sampling.p_hub.origin = 609380.0 4005240.0 0.0
sampling.p_hub.axis1 = 42000.0 0.0 0.0
sampling.p_hub.axis2 = 0.0 35000.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
#== END AMR-WIND INPUT ==