Is it possible to perform PyFluent operations within a PySystemCoupling solve loop?

[johnhavenar]

Hello All,

I'm seeking to couple Fluent with an external Python solver. It deals with a thermally-coupled system (Fluent providing heat flow at a surface; the Python solver providing wall temperature at the coupling interface).

As the domain of the external Python solver heats up, it produces gas which is ejected through the coupled interface. I'd like to take species mass fluxes of ejected gases from the Python solver and inject them into the Fluent domain. Is it a supported capability to within each System Coupling solution iteration perform a PyFluent operation on the domain which would modify a boundary condition?

Mockup code:

while needs_iteration:
    update_inputs()
    apply_species_flux_with_pyfluent()
    solve_coupled_system()
    update_outputs()