From nonlinear models and direct numerical simulations we report on several new findings of relevance to the single-mode Rayleigh-Taylor (RT) instability driven by time-varying acceleration histories. The incompressible, Direct Numerical Simulations (DNS) were performed in two-and three-dimensions, and at a range of density ratios of the fluid combinations (characterized by the Atwood number). We investigated several acceleration histories, including acceleration profiles of and with corresponding 3D drag buoyancy model solutions derived in this article. Our generalization of the RT problem to study variable g(t) affords us the opportunity to investigate the appropriate scaling for bubble and spike amplitudes under these conditions. We consider two candidates, the displacement Z and width s 2 , but find the appropriate scaling is dependent on the density ratios between the fluids -at low density ratios, bubble and spike amplitudes are explained by both s 2 and Z, while at large density differences the displacement collapses the spike data. Finally, for all the acceleration profiles studied here, spikes enter a free-fall regime at lower Atwood numbers than predicted by all the models.2