The objective of this study, which is divided into two parts, is twofold: to address long-standing challenges in the sensing of atmospheric turbulence-induced wavefront aberrations under strong scintillation conditions via a comparative analysis of several basic scintillation-resistant wavefront sensing (SR-WFS) architectures and iterative phase retrieval (IPR) techniques (Part I, this paper), and to develop a framework for the potential integration of SR-WFS techniques into practical closed-loop non-astronomical atmospheric adaptive optics (AO) systems (Part II). In this paper, we consider basic SR-WFS mathematical models and phase retrieval algorithms, tradeoffs in sensor design and phase retrieval technique implementation, and methodologies for WFS parameter optimization and performance assessment. The analysis is based on wave-optics numerical simulations imitating realistic turbulence-induced phase aberrations and intensity scintillations, as well as optical field propagation inside the SR-WFSs. Several potential issues important for the practical implementation of SR-WFS and IPR techniques, such as the requirements for phase retrieval computational grid resolution, tolerance with respect to optical element misalignments, and the impact of camera noise and input light non-monochromaticity, are also considered. The results demonstrate that major wavefront sensing requirements desirable for AO operation under strong intensity scintillations can potentially be achieved by transitioning to novel SR-WFS architectures, based on iterative phase retrieval techniques.