Asymptotic profiles of the constraint satisfying solution plotted with respect to heliocentric radius (s). Top panel: Profiles of T/T_*, (u_∞ − u)/(u_∞ − u_*), and an s^−2/7 power law for reference. Middle panel: Profiles of f/f_*, and σ/σ_*, with an s⁻² power law for reference. In both cases, the computed solutions are in good agreement with the representative power laws in the asymptotic region (r ≳ 10⁴ R_⊙). Bottom panel: Total energy flux associated with heat conduction (f), viscous stress (σu), bulk kinetic energy (ρu³/2), as well as saturation limits associated with free-streaming of electrons (p_ev_e) and ion thermal energy density (p_iu). The total energy flux, which is proportional to the sum of (ρu³/2 + f + σu), is shown in the black dashed-dotted curve.