The present paper presents a numerical solution method for the two-dimensional elliptic partial differential equations modeling a premixed laminar flame freely propagating at constant pressure between parallel plates having a constant, cool wall temperature. The technique developed is an efficient, semidirect (semi-iterative) method termed "multistep damped Picard iteration." Picard linearization of the nonlinear terms together with variable damping as the solution progressed and a sparse matrix elliptic solver allowed rapid solution to be obtained as a limit of a convergent sequence of solutions. The flame-quenching limit was determined by iterating the solution with respect to the distance separating the plates and approaching the limit state to within a specific criterion. Global kinetics and average transport properties were used for a simulated propane-air mixture; radiation loss was neglected throughout the analysis. The method was used to study the effect of the "ignition temperature" on the flame parameters as well as to reveal the two-dimensional structure of the flame at the quenching limit.
Nomenclature= cross-sectional area, Sec. IV B, cm 2 = specific heat at constant pressure of the mixture, cal/g K = diffusion coefficient for the fuel species cm 2 /s = hydraulic (equivalent) diameter of the parallel plates, cm = activation energy, cal/mole = mass flow rate (mass burning rate), g/cm 2 s; eigenvalue of the flame equations = heat of combustion of the fuel, -cal/g = frequency factor (reaction rate constant), units depend on reaction order = distance separating the parallel plates, cm = reaction order = perimeter, Sec. IV B, cm = pressure, atm = quenching Peclet number = quenching distance, cm = universal gas constant, 1.986 cal/mole K = surface of the upstream boundary, cm 2 , Sec. Ill B = surface of the downstream boundary cm 2 , Sec. Ill B = temperature, K = ignition (cutoff) temperature, K = ambient temperature, K = x component of the mixture velocity, cm/s = volume in Sec. Ill B cm 3 = flame speed in Eq. (22), cm/s = velocity of the mixture at the downstream boundary (_L toS out ),cm/s = rate of consumption of fuel per unit volume and time, g/cm 3 s = mass of reactant (fuel) per unit mass of mixture -mass of reactant (fuel) in unburned mixture per unit mass of that mixture A = thermal conductivity of the mixture cal/cm s K p = local mixture density, g/cm P 0 = density of unburned mixture, g/cm 3 Pom = density of the mixture at the downstream boundary, Sec. Ill B Subscript 0 = upstream boundary (reactants or cold)