-
Notifications
You must be signed in to change notification settings - Fork 6
/
sample.inp
167 lines (141 loc) · 5.98 KB
/
sample.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
159
160
161
162
163
164
165
166
#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# SIMULATION STOP #
#.......................................#
time.stop_time = 20000.0 # Max (simulated) time to evolve
time.max_step = -1 # Max number of time steps
#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# TIME STEP COMPUTATION #
#.......................................#
time.fixed_dt = 0.25 # Use this constant dt if > 0
time.cfl = 0.95 # CFL factor
#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# INPUT AND OUTPUT #
#.......................................#
time.plot_interval = 1000 # Steps between plot files
time.checkpoint_interval = 10000 # Steps between checkpoint files
#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# PHYSICS #
#.......................................#
incflo.gravity = 0.0 0.0 -9.81 # Gravitational force (3D)
incflo.density = 1.0 # Reference density
incflo.use_godunov = 1
incflo.godunov_type = "ppm_nolim"
transport.viscosity = 1.872e-5
transport.laminar_prandtl = 0.7
transport.turbulent_prandtl = 0.3333
turbulence.model = OneEqKsgsM84
incflo.physics = ABL Actuator
ICNS.source_terms = BoussinesqBuoyancy CoriolisForcing ABLForcing ActuatorForcing
TKE.source_terms = KsgsM84Src
BoussinesqBuoyancy.reference_temperature = 300.0
CoriolisForcing.east_vector = 1.0 0.0 0.0
CoriolisForcing.north_vector = 0.0 1.0 0.0
CoriolisForcing.latitude = 41.5
CoriolisForcing.rotational_time_period = 86164.0900027328
ABLForcing.abl_forcing_height = 20.0
incflo.velocity = 10 0 0 #3.5355339059327373 3.5355339059327373 0.0
ABL.reference_temperature = 300.0
ABL.temperature_heights = 0.0 650.0 750.0 1000.0
ABL.temperature_values = 300.0 300.0 308.0 308.75
ABL.perturb_temperature = false
ABL.cutoff_height = 50.0
ABL.perturb_velocity = true
ABL.perturb_ref_height = 50.0
ABL.kappa = 0.40
ABL.surface_roughness_z0 = 0.0005
ABL.surface_temp_rate = -0.32
ABL.normal_direction = 2
ABL.stats_output_frequency = 2
ABL.stats_output_format = netcdf
#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# ADAPTIVE MESH REFINEMENT #
#.......................................#
amr.n_cell = 576 576 384 # Grid cells at coarsest AMRlevel
amr.max_level = 0 # Max AMR level in hierarchy
#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# GEOMETRY #
#.......................................#
geometry.prob_lo = 0. 0. 0. # Lo corner coordinates
geometry.prob_hi = 1500. 1500. 1000. # Hi corner coordinates
geometry.is_periodic = 1 1 0 # Periodicity x y z (0/1)
# Boundary conditions
zlo.type = "wall_model"
zlo.temperature_type = "wall_model"
zhi.type = "slip_wall"
zhi.temperature_type = "fixed_gradient"
zhi.temperature = 0.003
zlo.tke_type = "fixed_gradient"
incflo.verbose = 0
# MLMG settings
nodal_proj.mg_rtol = 1.0e-6
nodal_proj.mg_atol = 1.0e-12
mac_proj.mg_rtol = 1.0e-6
mac_proj.mg_atol = 1.0e-12
diffusion.mg_rtol = 1.0e-6
diffusion.mg_atol = 1.0e-12
temperature_diffusion.mg_rtol = 1.0e-10
temperature_diffusion.mg_atol = 1.0e-13
# Activate data probe sampling
incflo.post_processing = sampling
# Frequency of output for the data
sampling.output_frequency = 100
# Type of probes to output
sampling.labels = p_h p_f
# Fields to output
sampling.fields = velocity temperature
# Definitions for each probe
sampling.p_h.type = PlaneSampler
sampling.p_h.axis1 = 1500.0 0.0 0.0
sampling.p_h.axis2 = 0.0 1500.0 0.0
sampling.p_h.origin = 0.0 0.0 20.0
sampling.p_h.num_points = 76 76
sampling.p_h.normal = 0.0 0.0 1.0
sampling.p_h.offsets = 0.0 20.0 40.0
sampling.p_f.type = PlaneSampler
sampling.p_f.num_points = 21 21
sampling.p_f.origin = 500.0 500.0 90.0
sampling.p_f.axis1 = 500.0 0.0 0.0
sampling.p_f.axis2 = 0.0 500.0 0.0
sampling.p_f.normal = 0.0 0.0 1.0
sampling.p_f.offsets = 0.0 20.0 40.0
#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# MESH TAGGING #
#.......................................#
tagging.labels = s1 f1 g1
tagging.s1.type = CartBoxRefinement
tagging.s1.static_refinement_def = refinebox.txt
tagging.f1.type = FieldRefinement
tagging.f1.field_name = density
tagging.f1.grad_error = 0.1 0.1 0.1
tagging.g1.type = GeometryRefinement
tagging.g1.shapes = b1
tagging.g1.b1.type = box
tagging.g1.b1.origin = 50.0 500.0 50.0
tagging.g1.b1.xaxis = 300.0 0.0 0.0
tagging.g1.b1.yaxis = 0.0 500.0 0.0
tagging.g1.b1.zaxis = 0.0 0.0 100.0
#tagging.g1.c1.type = cylinder
#tagging.g1.c1.start = 500.0 500.0 250.0
#tagging.g1.c1.end = 500.0 500.0 750.0
#tagging.g1.c1.outer_radius = 300.0
#tagging.g1.c1.inner_radius = 275.0
Actuator.labels = WTG01 WTG02 WTG03
Actuator.type = TurbineFastLine
Actuator.TurbineFastLine.rotor_diameter = 126.0
Actuator.TurbineFastLine.hub_height = 90.0
Actuator.TurbineFastLine.num_points_blade = 64
Actuator.TurbineFastLine.num_points_tower = 12
Actuator.TurbineFastLine.epsilon = 10.0 10.0 10.0
Actuator.TurbineFastLine.epsilon_tower = 5.0 5.0 5.0
Actuator.TurbineFastLine.openfast_start_time = 0.0
Actuator.TurbineFastLine.openfast_stop_time = 1.0
Actuator.TurbineFastLine.nacelle_drag_coeff = 0.0
Actuator.TurbineFastLine.nacelle_area = 0.0
Actuator.TurbineFastLine.output_frequency = 10
Actuator.TurbineFastLine.density = 1.225
Actuator.WTG01.base_position = 750 750 0
Actuator.WTG01.openfast_input_file = "fast_inp/nrel5mw.fst"
Actuator.WTG02.base_position = 850 850 0
Actuator.WTG02.openfast_input_file = "fast_inp/nrel5mw.fst"
Actuator.WTG03.base_position = 650 650 0
Actuator.WTG03.openfast_input_file = "fast_inp/nrel5mw.fst"