In normal tees, SIFs are static numbers because of the simplicity of the geometry.
In the proposed geometry, the SIFs would be variable based on the forces and moments coming in from all sides. This would be computationally expensive, and is not available and may not be for quite some time.
The best I could imagine you could do would be to model a single branch in NozzlePro to calculate flexibilities for it, plug them into CAESAR, and then model the cross in a 3d FEA software, applying loads as calculated from CAESAR.
Historically, I would suspect several have taken the SIFs for normal tees, doubled them, and hope for the best.