An integral method is presented for determining the evolution of molecular mixing, finite rate chemical reactions, and local extinction in diffusion layers under the effect of an unsteady strain rate. The partial differential equations governing the reactant, product, and temperature profiles are used to derive ordinary differential equations governing the evolution of moments for the product and temperature profiles and for the reactant gradient profiles. The actual profiles enter these equations only through integral moments resulting from the reaction rate terms (referred to as "reaction integrals"). As a consequence, it is possible to accurately track the evolution of the profile moments, and thereby determine global properties of the layer such as burning rates and extinction conditions, using remarkably simple representations for the actual profiles to evaluate the reaction integrals. Here these profile shapes are specified as self-similar families of curves parameterized by just a few degrees of freedom, which then evolve from the moment equations. Results for combustion in isolated strained diffusion layers, as well as for consumption of a burning fuel strip, are generally within a few percent of the results from finite difference solutions of the full equations.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/29480/1/0000566.pd