Modeling the Subsurface Evolution of Active Region Flux Tubes