Compositional Nonlinear Solver Based On Trust Regions Of The Flux Function Along Key Tie-Lines