The main result of this paper is a general Holder estimate in a class of singularly perturbed elliptic systems. This estimate is applied to the study of the well-known Burke–Schuman approximation in flame theory. After reviewing some classical cases (equidiffusional case; high activation energy approximation) we turn to the non-equidiffusional case, and to nonlinear diffusion models. The limiting problems are nonlinear elliptic equations; they have Holder or Lipschitz maximal global regularity. A natural question is then whether this regularity is kept uniformly throughout the approximation process, and we show that this is true in general.