def test_regression_yaw(sample_inputs_fixture):
"""
Tandem turbines with the upstream turbine yawed
"""
sample_inputs_fixture.floris["wake"]["properties"][
"velocity_model"
] = VELOCITY_MODEL
sample_inputs_fixture.floris["wake"]["properties"][
"deflection_model"
] = DEFLECTION_MODEL
floris = Floris(input_dict=sample_inputs_fixture.floris)
# yaw the upstream turbine 5 degrees
floris.farm.turbines[0].yaw_angle = 5.0
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i in range(len(floris.farm.turbine_map.turbines)):
baseline = yawed_baseline
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
def test_regression_secondary_steering(sample_inputs_fixture):
"""
Tandem turbines with the upstream turbine yawed and yaw added recovery
correction enabled
"""
sample_inputs_fixture.floris["wake"]["properties"][
"velocity_model"
] = VELOCITY_MODEL
sample_inputs_fixture.floris["wake"]["properties"][
"deflection_model"
] = DEFLECTION_MODEL
floris = Floris(input_dict=sample_inputs_fixture.floris)
floris.farm.turbines[0].yaw_angle = 5.0
# Enable secondary steering
floris.farm.wake.deflection_model.use_secondary_steering = True
floris.farm.wake.velocity_model.calculate_VW_velocities = True
floris.farm.flow_field.reinitialize_flow_field()
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i in range(len(floris.farm.turbine_map.turbines)):
check = secondary_steering_baseline
assert test_results[i][0] == approx(check[i][0])
assert test_results[i][1] == approx(check[i][1])
assert test_results[i][2] == approx(check[i][2])
assert test_results[i][3] == approx(check[i][3])
def test_regression_tandem(sample_inputs_fixture):
"""
Tandem turbines
"""
sample_inputs_fixture.floris["wake"]["model_strings"]["velocity_model"] = VELOCITY_MODEL
sample_inputs_fixture.floris["wake"]["model_strings"]["deflection_model"] = DEFLECTION_MODEL
sample_inputs_fixture.floris["wake"]["model_strings"]["combination_model"] = COMBINATION_MODEL
floris = Floris.from_dict(sample_inputs_fixture.floris)
floris.initialize_domain()
floris.steady_state_atmospheric_condition()
n_turbines = floris.farm.n_turbines
n_wind_speeds = floris.flow_field.n_wind_speeds
n_wind_directions = floris.flow_field.n_wind_directions
velocities = floris.flow_field.u
yaw_angles = floris.farm.yaw_angles
test_results = np.zeros((n_wind_directions, n_wind_speeds, n_turbines, 4))
farm_avg_velocities = average_velocity(
velocities,
)
farm_cts = Ct(
velocities,
yaw_angles,
floris.farm.turbine_fCts,
floris.farm.turbine_type_map,
)
farm_powers = power(
floris.flow_field.air_density,
velocities,
yaw_angles,
floris.farm.pPs,
floris.farm.turbine_power_interps,
floris.farm.turbine_type_map,
)
farm_axial_inductions = axial_induction(
velocities,
yaw_angles,
floris.farm.turbine_fCts,
floris.farm.turbine_type_map,
)
for i in range(n_wind_directions):
for j in range(n_wind_speeds):
for k in range(n_turbines):
test_results[i, j, k, 0] = farm_avg_velocities[i, j, k]
test_results[i, j, k, 1] = farm_cts[i, j, k]
test_results[i, j, k, 2] = farm_powers[i, j, k]
test_results[i, j, k, 3] = farm_axial_inductions[i, j, k]
if DEBUG:
print_test_values(
farm_avg_velocities,
farm_cts,
farm_powers,
farm_axial_inductions,
)
assert_results_arrays(test_results[0], baseline)
def test_change_wind_direction(sample_inputs_fixture):
"""
Tandem turbines aligned to the wind direction, first calculated at 270 degrees wind
direction, and then calculated at 315 degrees wind direction.
"""
sample_inputs_fixture.floris["wake"]["properties"][
"velocity_model"
] = VELOCITY_MODEL
sample_inputs_fixture.floris["wake"]["properties"][
"deflection_model"
] = DEFLECTION_MODEL
floris = Floris(input_dict=sample_inputs_fixture.floris)
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i, turbine in enumerate(floris.farm.turbine_map.turbines):
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
floris.farm.wind_map.input_direction = [315.0]
floris.farm.wind_map.calculate_wind_direction()
floris.farm.turbine_map.reinitialize_turbines()
floris.farm.flow_field.reinitialize_flow_field()
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i, turbine in enumerate(floris.farm.turbine_map.turbines):
assert test_results[i][0] == approx(wd315_baseline[i][0])
assert test_results[i][1] == approx(wd315_baseline[i][1])
assert test_results[i][2] == approx(wd315_baseline[i][2])
assert test_results[i][3] == approx(wd315_baseline[i][3])
def test_gauss_to_curl_to_gauss(sample_inputs_fixture):
"""
Start with the Gauss wake model
Then, switch to Curl
Then, switch back to Gauss
"""
# Establish that the Gauss test passes
sample_inputs_fixture.floris["wake"]["properties"][
"velocity_model"] = "gauss_legacy"
sample_inputs_fixture.floris["wake"]["properties"][
"deflection_model"] = "gauss"
floris = Floris(input_dict=sample_inputs_fixture.floris)
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i in range(len(floris.farm.turbine_map.turbines)):
baseline = gauss_baseline
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
# Change the model to Curl, rerun calculate_wake, and compare to Curl
floris.farm.set_wake_model("curl")
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i in range(len(floris.farm.turbine_map.turbines)):
baseline = curl_baseline
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
# Change back to Gauss, rerun calculate_wake, and compare to gauss
floris.farm.set_wake_model("gauss_legacy")
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i in range(len(floris.farm.turbine_map.turbines)):
baseline = gauss_baseline
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
def test_regression_tandem(sample_inputs_fixture):
"""
Tandem turbines
"""
sample_inputs_fixture.floris["wake"]["properties"][
"velocity_model"
] = VELOCITY_MODEL
sample_inputs_fixture.floris["wake"]["properties"][
"deflection_model"
] = DEFLECTION_MODEL
floris = Floris(input_dict=sample_inputs_fixture.floris)
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i in range(len(floris.farm.turbine_map.turbines)):
check = baseline
assert test_results[i][0] == approx(check[i][0])
assert test_results[i][1] == approx(check[i][1])
assert test_results[i][2] == approx(check[i][2])
assert test_results[i][3] == approx(check[i][3])
def test_regression_multiple_calc_wake(sample_inputs_fixture):
"""
Verify turbine values stay the same with repeated (3x) calculate_wake calls.
"""
sample_inputs_fixture.floris["wake"]["properties"][
"velocity_model"
] = VELOCITY_MODEL
sample_inputs_fixture.floris["wake"]["properties"][
"deflection_model"
] = DEFLECTION_MODEL
floris = Floris(input_dict=sample_inputs_fixture.floris)
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i, turbine in enumerate(floris.farm.turbine_map.turbines):
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i, turbine in enumerate(floris.farm.turbine_map.turbines):
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
floris.farm.flow_field.calculate_wake()
test_results = turbines_to_array(floris.farm.turbine_map.turbines)
if DEBUG:
print_test_values(floris.farm.turbine_map.turbines)
for i, turbine in enumerate(floris.farm.turbine_map.turbines):
assert test_results[i][0] == approx(baseline[i][0])
assert test_results[i][1] == approx(baseline[i][1])
assert test_results[i][2] == approx(baseline[i][2])
assert test_results[i][3] == approx(baseline[i][3])
def test_regression_secondary_steering(sample_inputs_fixture):
"""
Tandem turbines with the upstream turbine yawed and secondary steering enabled
"""
sample_inputs_fixture.floris["wake"]["model_strings"]["velocity_model"] = VELOCITY_MODEL
sample_inputs_fixture.floris["wake"]["model_strings"]["deflection_model"] = DEFLECTION_MODEL
sample_inputs_fixture.floris["wake"]["enable_transverse_velocities"] = True
sample_inputs_fixture.floris["wake"]["enable_secondary_steering"] = True
sample_inputs_fixture.floris["wake"]["enable_yaw_added_recovery"] = False
floris = Floris.from_dict(sample_inputs_fixture.floris)
yaw_angles = np.zeros((N_WIND_DIRECTIONS, N_WIND_SPEEDS, N_TURBINES))
yaw_angles[:,:,0] = 5.0
floris.farm.yaw_angles = yaw_angles
floris.initialize_domain()
floris.steady_state_atmospheric_condition()
n_turbines = floris.farm.n_turbines
n_wind_speeds = floris.flow_field.n_wind_speeds
n_wind_directions = floris.flow_field.n_wind_directions
velocities = floris.flow_field.u
yaw_angles = floris.farm.yaw_angles
test_results = np.zeros((n_wind_directions, n_wind_speeds, n_turbines, 4))
farm_avg_velocities = average_velocity(
velocities,
)
farm_cts = Ct(
velocities,
yaw_angles,
floris.farm.turbine_fCts,
floris.farm.turbine_type_map,
)
farm_powers = power(
floris.flow_field.air_density,
velocities,
yaw_angles,
floris.farm.pPs,
floris.farm.turbine_power_interps,
floris.farm.turbine_type_map,
)
farm_axial_inductions = axial_induction(
velocities,
yaw_angles,
floris.farm.turbine_fCts,
floris.farm.turbine_type_map,
)
for i in range(n_wind_directions):
for j in range(n_wind_speeds):
for k in range(n_turbines):
test_results[i, j, k, 0] = farm_avg_velocities[i, j, k]
test_results[i, j, k, 1] = farm_cts[i, j, k]
test_results[i, j, k, 2] = farm_powers[i, j, k]
test_results[i, j, k, 3] = farm_axial_inductions[i, j, k]
if DEBUG:
print_test_values(
farm_avg_velocities,
farm_cts,
farm_powers,
farm_axial_inductions,
)
assert_results_arrays(test_results[0], secondary_steering_baseline)
