Asymptotic and numerical methods are used to study several classes of singularly perturbed boundary value problems for which the underlying homogeneous operators have exponentially small eigenvalues. Examples considered include the familiar boundary layer resonance problems and some extensions and certain linearized equations associated with metastable internal layer motion. For the boundary layer resonance problems, a systematic projection method, motivated by the work of De Groen , is used to analytically calculate high-order asymptotic solutions. This method justifies and extends some previous results obtained from the variational method of Grasman and Matkowsky . A numerical approach, based on an integral equation formulation, is used to accurately compute boundary layer resonance solutions and their associated exponentially small eigenvalues. For various examples, the numerical results are shown to compare very favorably with two-term asymptotic results. Finally, some Sturm-Liouville operators with exponentially small spectral gap widths are studied. One such problem is applied to analyzing metastable internal layer motion for a certain forced Burgers equation.