cursor,
station_identifier,
sd.Q_BULK_TO_CENTER_LINE_VELOCITY_RATIO,
bulk_to_center_line_velocity_ratio,
averaging_system=sd.UNWEIGHTED_AVERAGING_SYSTEM,
)
for quantity in [
sd.Q_ROUGHNESS_HEIGHT,
sd.Q_INNER_LAYER_ROUGHNESS_HEIGHT,
sd.Q_OUTER_LAYER_ROUGHNESS_HEIGHT,
]:
sd.set_labeled_value(
cursor,
station_identifier,
quantity,
sd.WALL_POINT_LABEL,
0.0,
)
sd.set_labeled_value(
cursor,
station_identifier,
sd.Q_SHEAR_STRESS,
sd.WALL_POINT_LABEL,
wall_shear_stress,
averaging_system=sd.UNWEIGHTED_AVERAGING_SYSTEM,
)
sd.set_labeled_value(
cursor,
station_identifier,
station_number=station_number,
point_number=point_number,
point_label=sd.WALL_POINT_LABEL,
)
# In general, the surface is not well-described in this study at all.
# The data is consistent with a smooth surface, though.
for quantity in [
sd.Q_ROUGHNESS_HEIGHT,
sd.Q_INNER_LAYER_ROUGHNESS_HEIGHT,
sd.Q_OUTER_LAYER_ROUGHNESS_HEIGHT,
]:
sd.set_labeled_value(
cursor,
station_identifier,
quantity,
sd.WALL_POINT_LABEL,
0.0,
measurement_technique=sd.MT_ASSUMPTION,
)
sd.set_labeled_value(
cursor,
station_identifier,
sd.Q_MASS_DENSITY,
sd.WALL_POINT_LABEL,
mass_density,
averaging_system=sd.BOTH_AVERAGING_SYSTEMS,
measurement_technique=sd.MT_CALCULATION,
)
sd.set_labeled_value(
cursor,
averaging_system=sd.BOTH_AVERAGING_SYSTEMS,
measurement_technique=sd.MT_CALCULATION,
)
i += 1
if (is_rough_wall == False):
for quantity in [
sd.Q_ROUGHNESS_HEIGHT,
sd.Q_INNER_LAYER_ROUGHNESS_HEIGHT,
sd.Q_OUTER_LAYER_ROUGHNESS_HEIGHT,
]:
sd.set_labeled_value(
cursor,
station_identifier,
quantity,
sd.WALL_POINT_LABEL,
0.0,
measurement_technique=sd.MT_ASSUMPTION,
)
r_prof, u_prof = sd.get_twin_profiles(
cursor,
station_identifier,
sd.Q_TRANSVERSE_COORDINATE,
sd.Q_STREAMWISE_VELOCITY,
)
volumetric_flow_rate = -2.0 * math.pi * sd.integrate_using_trapezoid_rule(
r_prof, u_prof * r_prof)
mass_flow_rate = mass_density * volumetric_flow_rate
bulk_velocity = 4.0 * volumetric_flow_rate / (math.pi * diameter**2.0)