A compressible numerical model is applied for three-dimensional (3-D) gravity wave (GW) packets undergoing momentum deposition, self-acceleration (SA), breaking, and secondary GW (SGW) generation in the presence of highly-structured environments enabling thermal and/or Doppler ducts, such as a mesospheric inversion layer (MIL), tidal wind (TW), or combination of MIL and TW. Simulations reveal that ducts can strongly modulate GW dynamics. Responses modeled here include reflection, trapping, suppressed transmission, strong local instabilities, reduced SGW generations, higher altitude SGW responses, and induced large-scale flows. Instabilities that arise in ducts experience strong dissipation after they emerge, while trapped smaller-amplitude and smaller-scale GWs can survive in ducts to much later times. Additionally, GW breaking and its associated dynamics enhance the local wind along the GW propagation direction in the ducts, and yield layering in the wind field. However, these dynamics do not yield significant heat transport in the ducts. The failure of GW breaking to induce stratified layers in the temperature field suggests that such heat transport might not be as strong as previously assumed or inferred from observations and theoretical assessments. The present numerical simulations confirm previous finding that MIL generation may not be caused by the breaking of a transient high-frequency GW packet alone.