High-fidelity modeling of nuclear reactors requires the solution of nonlinear coupled multi-physics stiff problems with widely varying time and length scales that need to be resolved correctly. A numerical method that converges the implicit nonlinear terms to a small tolerance is often referred to as nonlinearly consistent (or tightly coupled). This nonlinear consistency is still lacking in the...