How to build a coupling boundary conditions of a multi-region solver
 

Hi guys,
I tried to build a multi-region solver fo solving neutron diffusion equation. However, difficulties are encountered when the coupling boundary between regions is given. The default solver use the compressible::turbulentTemperatureCoupledBaffleMix ed boundary，how to modify it for satisfying the neutron flux volscalarfield?It is also necessary to meet the same continuity and heat flux.