Confining two-body systems provide a useful setting for examining the internal consistency of perturbative treatments of bound-state wavefunctions. In this work, we investigate the applicability of Dalgarno–Lewis perturbation theory beyond first order within a nonrelativistic Cornell-type potential, \( V(r)=br-\tfrac{4\alpha _s}{3r}+c \) , treating the linear confinement term as the parent Hamiltonian and incorporating the Coulombic interaction perturbatively so as to retain analytic control over the solutions. Going beyond standard first-order implementations, we carry out a systematic second-order Dalgarno–Lewis correction to the Airy-function ground state and analyze its impact on wavefunction moments. We show that while lower spatial moments are only mildly affected at first order, observables involving higher spatial moments are not fully stabilized at first order at the sub-percent level, indicating that first-order Dalgarno–Lewis theory does not by itself guarantee quantitative control of such quantities in confining two-body systems. Inclusion of second-order contributions is therefore necessary to ensure controlled perturbative accuracy and restore quantitative stability. As an explicit application, we evaluate effective overlap functions, their slopes and curvatures, as well as geometric and electromagnetic radii for representative heavy–light mesons (D, \(D_s\) , B, and \(B_s\) ). While the analysis is model dependent and is carried out within a nonrelativistic potential-model framework and does not rely on strict heavy-quark universality, the resulting magnitudes and systematic trends are compatible with those obtained in relativistic quark models, QCD sum-rule approaches, and continuum Dyson–Schwinger/Bethe–Salpeter studies. The present results provide semi-analytic benchmarks that clarify the role of higher-order wavefunction corrections in few-body bound-state problems.