def find_zero_constrained_reactions(model):
    """Return list of reactions that are constrained to zero flux."""
    return [
        rxn for rxn in model.reactions if rxn.lower_bound == 0 and rxn.upper_bound == 0
    ]