We perform 3D radiative hydrodynamic simulations to study convection in low-mass main-sequence stars with the aim of improving stellar models. Comparing models from a 0.90 M⊙ evolutionary track with 3D simulations reveals distinct differences between simulations and mixing length theory. The simulations show obvious structural differences throughout the superadiabatic layer where convection is inefficient at transporting energy. The discrepancy between MLT and simulation changes as the star evolves and the dynamical effects of turbulence increase. Further, the simulations reveal a T-tau relation that is sensitive to the strength of the turbulence, which is in contrast to 1D stellar models that use the same T-tau relation across the HR diagram.