Linear stability of planar premixed flames: reactive Navier-Stokes equations with finite activation energy and arbitrary Lewis number

A numerical shooting method for performing linear stability analyses of travelling waves is described and applied to the problem of freely propagating planar premixed flames. Previous linear stability analyses of premixed flames either employ high activation temperature asymptotics or have been performed numerically with finite activation temperature, but either for unit Lewis numbers (which ignores thermal-diffusive effects) or in the limit of small heat release (which ignores hydrodynamic effects). In this paper the full reactive Navier-Stokes equations are used with arbitrary values of the parameters (activation temperature, Lewis number, heat of reaction, Prandtl number), for which both thermal-diffusive and hydrodynamic effects on the instability, and their interactions, are taken into account. Comparisons are made with previous asymptotic and numerical results. For Lewis numbers very close to or above unity, for which hydrodynamic effects caused by thermal expansion are the dominant destablizing mechanism, it is shown that slowly varying flame analyses give qualitatively good but quantitatively poor predictions, and also that the stability is insensitive to the activation temperature. However, for Lewis numbers sufficiently below unity for which thermal-diffusive effects play a major role, the stability of the flame becomes very sensitive to the activation temperature. Indeed, unphysically high activation temperatures are required for the high activation temperature analysis to give quantitatively good predictions at such low Lewis numbers. It is also shown that state-insensitive viscosity has a small destabilizing effect on the cellular instability at low Lewis numbers

Effect of thermal expansion on the linear stability of planar premixed flames for a simple chain-branching model: The high activation energy asymptotic limit

The linear stability of freely propagating, adiabatic, planar premixed ames is investigated in the context of a simple chain-branching chemistry model consisting of a chain-branching reaction step and a completion reaction step. The role of chain-branching is governed by a crossover temperature. Hydrodynamic effects, induced by thermal expansion, are taken into account and the results compared and contrasted with those from a previous purely thermal-di�usive constant density linear stability study. It is shown that when thermal expansion is properly accounted for, a region of stable ames predicted by the constant density model disappears, and instead the ame is unstable to a long-wavelength cellular instability. For a pulsating mode, however, thermal expansion is shown to have only a weak e�ect on the critical fuel Lewis number required for instability. These e�ects of thermal expansion on the two-step chain-branching ame are shown to be qualitatively similar to those on the standard one-step reaction model. Indeed, as found by constant density studies, in the limit that the chain-branching crossover temperature tends to the adiabatic ame temperature, the two-step model can be described to leading order by the one-step model with a suitably de�ned e�ective activation energy

Partially quenched chiral perturbation theory without Φ0\Phi_0

This paper completes the argument that lattice simulations of partially quenched QCD can provide quantitative information about QCD itself, with the aid of partially quenched chiral perturbation theory. A barrier to doing this has been the inclusion of $\Phi_0$, the partially quenched generalization of the $\eta'$, in previous calculations in the partially quenched effective theory. This invalidates the low energy perturbative expansion, gives rise to many new unknown parameters, and makes it impossible to reliably calculate the relation between the partially quenched theory and low energy QCD. We show that it is straightforward and natural to formulate partially quenched chiral perturbation theory without $\Phi_0$, and that the resulting theory contains the effective theory for QCD without the $\eta'$. We also show that previous results, obtained including $\Phi_0$, can be reinterpreted as applying to the theory without $\Phi_0$. We contrast the situation with that in the quenched effective theory, where we explain why it is necessary to include $\Phi_0$. We also compare the derivation of chiral perturbation theory in partially quenched QCD with the standard derivation in unquenched QCD. We find that the former cannot be justified as rigorously as the latter, because of the absence of a physical Hilbert space. Finally, we present an encouraging result: unphysical double poles in certain correlation functions in partially quenched chiral perturbation theory can be shown to be a property of the underlying theory, given only the symmetries and some plausible assumptions.Comment: 45 pages, no figure

Physical Results from Unphysical Simulations

We calculate various properties of pseudoscalar mesons in partially quenched QCD using chiral perturbation theory through next-to-leading order. Our results can be used to extrapolate to QCD from partially quenched simulations, as long as the latter use three light dynamical quarks. In other words, one can use unphysical simulations to extract physical quantities - in this case the quark masses, meson decay constants, and the Gasser-Leutwyler parameters L_4-L_8. Our proposal for determining L_7 makes explicit use of an unphysical (yet measurable) effect of partially quenched theories, namely the double-pole that appears in certain two-point correlation functions. Most of our calculations are done for sea quarks having up to three different masses, except for our result for L_7, which is derived for degenerate sea quarks.Comment: 26 pages, 12 figures (discussion on discretization errors at end of sec. IV clarified; minor improvements in presentation; results unchanged

Steady non-ideal detonations in cylindrical sticks of expolsives

Numerical simulations of detonations in cylindrical rate-sticks of highly non-ideal explosives are performed, using a simple model with a weakly pressure dependent rate law and a pseudo-polytropic equation of state. Some numerical issues with such simulations are investigated, and it is shown that very high resolution (hundreds of points in the reaction zone) are required for highly accurate (converged) solutions. High resolution simulations are then used to investigate the qualitative dependences of the detonation driving zone structure on the diameter and degree of confinement of the explosive charge. The simulation results are used to show that, given the radius of curvature of the shock at the charge axis, the steady detonation speed and the axial solution are accurately predicted by a quasi-one-dimensional theory, even for cases where the detonation propagates at speeds significantly below the Chapman-Jouguet speed. Given reaction rate and equation of state models, this quasi-one-dimensional theory offers a significant improvement to Wood-Kirkwood theories currently used in industry

Heavy-Meson Observables at One-Loop in Partially Quenched Chiral Perturbation Theory

I present one-loop level calculations of the Isgur-Wise functions for B -> D^{(*)} + e + nu, of the matrix elements of isovector twist-2 operators in B and D mesons, and the matrix elements for the radiative decays D^* -> D + gamma in partially quenched heavy quark chiral perturbation theory. Such expressions are required in order to extrapolate from the light quark masses used in lattice simulations of the foreseeable future to those of nature.Comment: 13 pages, 3 fig

Chiral Perturbation Theory for the Quenched Approximation of QCD

[This version is a minor revision of a previously submitted preprint. Only references have been changed.] We describe a technique for constructing the effective chiral theory for quenched QCD. The effective theory which results is a lagrangian one, with a graded symmetry group which mixes Goldstone bosons and fermions, and with a definite (though slightly peculiar) set of Feynman rules. The straightforward application of these rules gives automatic cancellation of diagrams which would arise from virtual quark loops. The techniques are used to calculate chiral logarithms in $f_K/f_\pi$, $m_\pi$, $m_K$, and the ratio of $\langle{\bar s}s\rangle$ to $\langle{\bar u}u\rangle$. The leading finite-volume corrections to these quantities are also computed. Problems for future study are described.Comment: 14 page

Applications of Partially Quenched Chiral Perturbation Theory

Partially quenched theories are theories in which the valence- and sea-quark masses are different. In this paper we calculate the nonanalytic one-loop corrections of some physical quantities: the chiral condensate, weak decay constants, Goldstone boson masses, B_K and the K+ to pi+ pi0 decay amplitude, using partially quenched chiral perturbation theory. Our results for weak decay constants and masses agree with, and generalize, results of previous work by Sharpe. We compare B_K and the K+ decay amplitude with their real-world values in some examples. For the latter quantity, two other systematic effects that plague lattice computations, namely, finite-volume effects and unphysical values of the quark masses and pion external momenta are also considered. We find that typical one-loop corrections can be substantial.Comment: 22 pages, TeX, refs. added, minor other changes, version to appear in Phys. Rev.

Kaon B parameter from quenched Lattice QCD

We present results of a large-scale simulation for the Kaon B parameter $B_K$ in quenched lattice QCD with the Kogut-Susskind quark action. Calculating $B_K$ at 1% statistical accuracy for seven values of lattice spacing in the range $a\approx 0.24-0.04$ fm on lattices up to $56^3\times 96$, we verify a quadratic $a$ dependence of $B_K$ theoretically predicted. Strong indications are found that, with our level of accuracy, $\alpha_{\bar{MS}}(1/a)^2$ terms arising from our one-loop matching procedure have to be included in the continuum extrapolation. We present $B_K$(NDR, 2 GeV)=0.628(42) as our final value, as obtained by a fit including the $\alpha_{\bar{MS}}(1/a)^2$ term.Comment: 8 pages, Latex(revtex, epsf), 2 epsf figure

High Precision determination of the pi, K, D and D_s decay constants from lattice QCD

We determine $D$ and $D_s$ decay constants from lattice QCD with 2% errors, 4 times better than experiment and previous theory: $f_{D_s}$ = 241(3) MeV, $f_D$ = 207(4) MeV and $f_{D_s}/f_D$ = 1.164(11). We also obtain $f_K/f_{\pi}$ = 1.189(7) and $(f_{D_s}/f_D)/(f_K/f_{\pi})$ = 0.979(11). Combining with experiment gives $V_{us}$=0.2262(14) and $V_{cs}/V_{cd}$ of 4.43(41). We use a highly improved quark discretization on MILC gluon fields that include realistic sea quarks fixing the $u/d, s$ and $c$ masses from the $\pi$, $K$, and $\eta_c$ meson masses. This allows a stringent test against experiment for $D$ and $D_s$ masses for the first time (to within 7 MeV).Comment: 4 pages, 2 figures. Published version - changes from original include a more extensive discussion of errors and an error budget table covering more quantities. There are very small changes in some of the values reporte