Spatially resolved measurements of soot and gaseous precursors in ethylene counterflow diffusion flames up to 32 atm

We investigate the effect of pressure on both flame structure and soot formation in nitrogen diluted counterflow diffusion flames of ethylene in the 8–32atm pressure range. Capillary-probe gas sampling is performed to resolve spatially the profiles of gaseous species up to three-ring aromatics by GC/MS analysis and multi-color pyrometry is used to quantify the soot volume fraction and dispersion exponent. Self-similarity of flames is preserved by keeping constant mixture fraction and strain rate, so that profiles of concentrations and temperature, normalized with respect to their peak values, are unaffected by changes in pressure, once the axial coordinate is nondimensionalized with respect to the pressure-dependent diffusion length scale. When conditions are chosen so that the overall soot loading is approximately constant and compatible with the diagnostics, it is found that both the soot volume fraction and the profiles of key aromatics in the high-temperature nucleation region are virtually invariant. For it to happen, a twofold increase in pressure must be compensated by a ～100 K decrease in peak flame temperature and, therefore, in the temperature across the soot forming region. The implication is that from the perspective of the chemical kinetics of soot formation these two actions counterbalance each other. As pressure increases (and temperature decreases) the peak production rate of the high-temperature soot mechanism decreases and, further downstream, towards the particle stagnation plane, a low-temperature soot mechanism sets in, yielding an increase in soot H/C content. This mechanism is enhanced as the pressure is raised, causing a higher overall soot volume production rate in the 16atm flame and, especially, in the 32atm one. The role of C4/C2 species in the formation of C₆H₆ increases with increasing pressure and dominates over the recombination of propargyl radical at sufficiently high pressures. A comprehensive database is established for soot models at high pressures of relevance to applications.     

Study of MILD combustion using LES and advanced analysis tools

A cylindrical confined combustor operating under MILD condition is investigated using LES. The combustion and its interaction with turbulence are modeled using two reactor based models, PaSR and EDC. Results show that the Partially Stirred Reactor (PaSR) model yields improved estimation for mean temperature and species mole fractions compared to Eddy Dissipation Concept (EDC). LES data are analysed using advanced post-processing methods such as the chemical Tangential Stretching Rate (TSR), balance analysis and local Principle Component (PCA) analysis. TSR can identify chemical explosive (ignition-like) and contractive (burnt) regions. With the balance analysis of the convective, diffusive and reactive terms in temperature equation, regions with substantial heat release coming from ignition or flame are identified. The local PCA analysis classifies the whole domain into clusters (regions with specific features) and provides the leading species in each cluster. The three analyses correlate well with one another and it is observed that the most chemically active region locates upstream (in the near-field). Also, both autoignition and flame-like structures play equally important roles in MILD combustion.     

Parton distributions and new physics searches: the Drell–Yan forward–backward asymmetry as a case study

We present the first unquenched lattice-QCD calculation of the form factors for the decay \(B\rightarrow D^*\ell \nu \) at nonzero recoil. Our analysis includes 15 MILC ensembles with \(N_f=2+1\) flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from \(a\approx 0.15\) fm down to 0.045 fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valence b and c quarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element \(|V_{cb}|\). We obtain \(\left| V_{cb}\right| = (38.40 \pm 0.68_{\text {th}} \pm 0.34_{\text {exp}} \pm 0.18_{\text {EM}})\times 10^{-3}\). The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall \(\chi ^2\text {/dof} = 126/84\), which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict \(R(D^*) = 0.265 \pm 0.013\), which confirms the current tension between theory and experiment.

     

Torques on a nearly rigid body in a relativistic gravitational field

The effect of post-Newtonian potentials on the rotation of a perfect fluid, nearly rigid body is shown to consist of a precession and a torque. The frequency of the precession can be exactly represented by means of suitable differential operators. The relativistic torques, in the quadrupole approximation, depend on the instantaneous orientation of the principal axes of one body with respect to the position-like the classical torque-and velocity of the other. For a relatively low-mass body, such as a gyroscope, these velocity-dependent torques have no observable consequences.National Research Council Resident Research Associate.     

Master integrals for double real radiation emission in heavy-to-light quark decay

We evaluate analytically the master integrals for double real radiation emission in the \(b \rightarrow u W^*\) decay, where b and u are a massive and massless quark, respectively, while \(W^{*}\) is an off-shell charged weak boson. Since the W boson can subsequently decay in a lepton anti-neutrino pair, the results of the present paper constitute a further step toward a fully analytic computation of differential distributions for the semileptonic decay of a b quark at NNLO in QCD. The latter partonic process plays a crucial role in the study of inclusive semileptonic charmless decays of B mesons. Our results are expressed in terms of multiple polylogarithms of maximum weight four.     

Inductive Construction of the Loop Transform for Abelian Gauge Theories

We construct the loop transform in the case of Abelian gauge theories as a unitary operator given by the inductive limit of Fourier transforms on tori. We also show that its range, i.e. the space of kinematical states of the quantum loop representation, is the Hilbert space of square integrable complex valued functions on the group of hoops.     

Scaling, Optimality, and Landscape Evolution

A nonlinear model is studied which describes the evolution of a landscape under the effects of erosion and regeneration by geologic uplift by mean of a simple differential equation. The equation, already in wide use among geomorphologists and in that context obtained phenomenologically, is here derived by reparametrization invariance arguments and exactly solved in dimension d=1. Results of numerical simulations in d=2 show that the model is able to reproduce the critical scaling characterizing landscapes associated with natural river basins. We show that configurations minimizing the rate of energy dissipation (optimal channel networks) are stationary solutions of the equation describing the landscape evolution. Numerical simulations show that a careful annealing of the equation in the presence of additive noise leads to configurations very close to the global minimum of the dissipated energy, characterized by mean field exponents. We further show that if one considers generalized river network configurations in which splitting of the flow (i.e., braiding) and loops are allowed, the minimization of the dissipated energy results in spanning loopless configurations, under the constraints imposed by the continuity equations. This is stated in the form of a general theorem applicable to generic networks, suggesting that other branching structures occurring in nature may possibly arise as optimal structures minimizing a cost function.