mass_density * bulk_velocity**2.0)
bulk_reynolds_number = bulk_velocity * hydraulic_diameter / kinematic_viscosity
friction_velocity = (wall_shear_stress / mass_density)**0.5
viscous_length_scale = kinematic_viscosity / friction_velocity
friction_reynolds_number = half_height / viscous_length_scale
speed_of_sound = sd.liquid_water_speed_of_sound(temperature)
bulk_mach_number = bulk_velocity / speed_of_sound
friction_mach_number = friction_velocity / speed_of_sound
series_identifier = sd.add_series(
cursor,
flow_class=flow_class,
year=year,
study_number=study_number,
series_number=series_number,
number_of_dimensions=2,
coordinate_system=sd.RECTANGULAR_COORDINATE_SYSTEM,
)
sd.add_working_fluid_component(
cursor,
series_identifier,
sd.WATER_LIQUID,
)
sd.update_series_geometry(cursor, series_identifier,
sd.RECTANGULAR_GEOMETRY)
sd.set_series_value(
Re_bulk = bulk_velocity * diameter / kinematic_viscosity
volumetric_flow_rate = 0.25 * math.pi * diameter**2.0 * bulk_velocity
speed_of_sound = sd.sdfloat("inf")
if (working_fluid == "Air"):
speed_of_sound = sd.ideal_gas_speed_of_sound(temperature)
elif (working_fluid == "Water"):
speed_of_sound = sd.liquid_water_speed_of_sound(temperature)
Ma_bulk = bulk_velocity / speed_of_sound
series_identifier = sd.add_series(
cursor,
flow_class=flow_class,
year=year,
study_number=study_number,
series_number=series_number,
number_of_dimensions=2,
coordinate_system=sd.CYLINDRICAL_COORDINATE_SYSTEM,
)
if (working_fluid == "Air"):
sd.add_air_components(cursor, series_identifier)
elif (working_fluid == "Water"):
sd.add_working_fluid_component(
cursor,
series_identifier,
sd.WATER_LIQUID,
)
sd.update_series_geometry(cursor, series_identifier,