A Comprehensive Study of Effects of Mixing and Chemical Kinetic Models on Predictions of n-heptane Jet Ignitions with the PDF Method

The composition Probability Density Function (PDF) model is coupled with a Reynolds-averaged k???ε turbulence model and three computationally efficient, yet widely used chemical mechanisms to simulate transient n-heptane spray injection and ignition in a high temperature and high density ambient fluid. Molecular diffusion is modelled by three mixing models, namely the interaction by exchange with the mean (IEM), modified Curl (MC) and Euclidean minimum spanning trees (EMST) models. The liquid phase is modelled by a discrete phase model (DPM). This represents among the first applications of the PDF method in practical diesel engine conditions. A non-reacting case is first considered, with the focus on the ability of the model to capture the spray structure, e.g., vapour penetration and liquid length, fuel mixture fraction and its variance. Reacting cases are then investigated to compare and evaluate the three different chemical mechanisms and the three mixing models. It is concluded that the EMST mixing model in conjunction with a reduced chemical kinetic model (Lu et al., Combust Flame 156(8):1542–1551, 2009) performs the best among the options considered. The sensitivity of the results to the choice of the mixing constant is also studied to understand its effect on the flame ignition and stabilisation. Finally, the PDF model is compared to a well-mixed model that assumes turbulent fluctuations are negligible, which has been widely used in the diesel spray combustion community. Significant structural differences in the modelled flame are revealed comparing the PDF method with the well-mixed model. Quantitatively, the PDF model exhibits excellent agreement with the measurements and shows much better results than the well-mixed model in all ambient O₂ and temperature conditions.     

Sample introduction in gas chromatography using a coiled wire filament

A simple device for field sampling and concentration of analytes for subsequent introduction into an injection port for gas chromatographic (GC) analysis has been developed. It consists of a tiny, coiled platinum wire filament (CWF) that is attached to a retractable plunger wire, which fits inside a syringe needle housing. Sampling is accomplished by dipping the end of the CWF in a liquid sample, which is drawn into the wire coil by capillary action, and introducing it into the injection port either before or after allowing the solvent to evaporate. The CWF can be used with or without a nonvolatile chemical coating. A major advantage of this sampling device is that nonvolatile sample matrix components remain on the wire coil, reducing the required injection port and liner cleaning frequency and contamination of the head of the chromatographic column. The coil itself can be easily cleaned between analyses by rinsing and/or burning off residual material in a small flame. The sampling coil facilitates specifically designed chemical reactions in the injection port, such as thermochemolysis and methylation. Applications demonstrated in this work include: (1) direct introduction of samples with little or no pre-treatment, (2) simultaneous thermochemolysis and methylation of lipid-containing samples such as bacteria and bacterial endospores for analysis of biomarkers, and (3) solid phase micro-extraction (SPME) using temporary wire coatings. The CWF allowed for significant reduction in sample preparation time, in most cases to less than a few minutes. The peak shapes examined for polycyclic aromatic hydrocarbon analytes (PAHs) were significantly better (asymmetry factors <1.3) when using the CWF sampling technique compared to splitless and on-column injection techniques (asymmetry factors >1.3). Extraction efficiencies for SPME (especially for high boiling point components such as PAHs) improved by an average of 2.5 times when using the CWF compared to the performance of commercially available SPME fibers. Coiled wire filaments and GC injection port liners were used for more than 100 Bacillus endospore thermochemolysis methylation analyses without the need for cleaning or replacement.     

Estimation of three-dimensional flame surface densities from planar images in turbulent premixed combustion

Turbulence motions are, by nature, three-dimensional while planar imaging techniques, widely used in turbulent combustion, give only access to two-dimensional information. For example, to extract flame surface densities, a key ingredient of some turbulent combustion models, from planar images implicitly assumes an instantaneously two-dimensional flow, neglecting the unresolved flame front wrinkling. The objective here is to estimate flame surface densities from two-dimensional measurements assuming that (1) the flow is statistically two dimensional; (2) the measuring plane is a plane of symmetry of the mean flow, either by translation (homogeneous third direction as in slot burners for example) or by rotation (axi-symmetrical flows such as jets) and (3) flame movements in transverse directions are similar. The unknown flame front wrinkling is then modelled from known quantities. An excellent agreement is achieved against direct numerical simulation (DNS) data where all three-dimensional quantities are known, but validations in other conditions (larger DNS, experiments) are required.     

Observation of e(+)e(-) annihilation into the C = +1 hadronic final states rho(0)rho(0) and phirho(0)

We report the first observation of e(+)e(-) annihilation into states of positive C parity, namely, rho(0)rho(0) and phirho(0). The two states are observed in the pi(+)pi(-)pi(+)pi(-) and K(+)K(-)pi(+)pi(-) final states, respectively, in a data sample of 225 fb(-1) collected by the BABAR experiment at the Positron-Electron Project II e(+)e(-) storage rings at energies near sqrt[s]=10.58 GeV. The distributions of costheta(*), where theta(*) is the center-of-mass polar angle of the phi meson or the forward rho(0) meson, suggest production by two-virtual-photon annihilation. We measure cross sections within the range |costheta(*)|<0.8 of sigma(e(+)e(-)-->rho(0)rho(0))=20.7+/-0.7(stat)+/-2.7(syst) fb and sigma(e(+)e(-)-->phirho(0))=5.7+/-0.5(stat)+/-0.8(syst) fb.     

A DNS evaluation of mixing and evaporation models for TPDF modelling of nonpremixed spray flames

In the present work, nonpremixed temporally evolving planar spray jet flames are simulated using both direct numerical simulation (DNS) and the composition transported probability density function (TPDF) method. The objective is to assess the performance of various mixing and evaporation source term distribution models which are required to close the PDF transport equation in spray flames. Quantities which would normally be provided to the TPDF solver by spray models and turbulence models are provided from the DNS: the mean flow velocity, turbulent diffusivity, mixing frequency, and cell-mean evaporation source term. Two cases with different Damköhler numbers (Da) are considered. The low Da case (Da-) features extinction followed by reignition while extinction in the high Da case (Da+) is insignificant. The TPDF modelling considers two mixing models: interaction by exchange with the mean (IEM) and Euclidean minimum spanning trees (EMST). Three models for distribution of the evaporation source terms are considered: EQUAL which distributes them in proportion to notional particles’ mass weight, NEW which creates new particles of pure fuel, and SAT which distributes the sources preferentially to notional particles close to saturation. It is found that the IEM model overpredicts the extinction when used with any evaporation model for both Da- and Da+ cases. The EMST model captures well the trend for extinction and reignition for the Da- case when it is coupled with the EQUAL evaporation model, but it overpredicts the extinction when coupled with the NEW or SAT evaporation model. For the Da+ case, all evaporation models reasonably capture the flame dynamics when coupled with EMST. The flame temperature in the mixture fraction space was examined to further assess the model performance. In general the EMST model results in narrow PDFs with little conditional fluctuation, while the IEM model produces bimodal PDFs with burning and partial extinction branches.     

Observation of the bottomonium ground state in the decay Upsilon(3S)-->gammaetab

We report the results of a search for the bottomonium ground state etab(1S) in the photon energy spectrum with a sample of (109+/-1) million of Upsilon(3S) recorded at the Upsilon(3S) energy with the BABAR detector at the PEP-II B factory at SLAC. We observe a peak in the photon energy spectrum at Egamma=921.2(-2.8)+2.1(stat)+/-2.4(syst) MeV with a significance of 10 standard deviations. We interpret the observed peak as being due to monochromatic photons from the radiative transition Upsilon(3S)-->gammaetab(1S). This photon energy corresponds to an etab(1S) mass of 9388.9(-2.3)+3.1(stat)+/-2.7(syst) MeV/c2. The hyperfine Upsilon(1S)-etab(1S) mass splitting is 71.4(-3.1)+2.3(stat)+/-2.7(syst) MeV/c2. The branching fraction for this radiative Upsilon(3S) decay is estimated to be [4.8+/-0.5(stat)+/-1.2(syst)]x10(-4).     

Exclusive branching-fraction measurements of semileptonic tau decays into three charged hadrons, into phipi(-)nu tau, and into phi K(-)nu tau

Using a data sample corresponding to an integrated luminosity of 342 fb(-1) collected with the BABAR detector at the SLAC PEP-II electron-positron storage ring operating at a center-of-mass energy near 10.58 GeV, we measure B(tau(-)--> pi(-)pi(-)pi+nu(tau)(ex.K(S0))=(8.83+/-0.01+/-0.13)%, B(tau(-) -->K(-)pi(-)pi+nu tau(ex.K(S0))=(0.273+/-0.002+/-0.009)%, B(tau(-) -->K(-)pi(-)K+nu tau)=(0.1346+/-0.0010+/-0.0036)%, and B(tau(-) -->K(-)K(-)K+nu tau)=(1.58+/-0.13+/-0.12)x10;{-5}, where the uncertainties are statistical and systematic, respectively. These include significant improvements over previous measurements and a first measurement of B(tau(-) -->K(-)K(-)K+nu tau) in which no resonance structure is assumed. We also report a first measurement of B(tau(-) -->var phi(-)nu tau)=(3.42+/-0.55+/-0.25)x10(-5), a new measurement of B(tau(-) -->var phi K(-)nu tau)=(3.39+/-0.20+/-0.28)x10(-5) and a first upper limit on B(tau(-) -->K(-)K(-)K+nu tau(ex.var phi)).     

Observation of Y(3940) --> J/psiomega in B --> J/psiomegaK at BABAR

We present a study of the decays B;{0,+}-->J/psiomegaK;{0,+} using 383x10;{6} BB[over ] events obtained with the BABAR detector at PEP-II. We observe Y(3940)-->J/psiomega, with mass 3914.6_{-3.4};{+3.8}(stat)+/-2.0(syst) MeV/c;{2}, and width 34_{-8};{+12}(stat)+/-5(syst) MeV. The ratio of B0 and B+ decay to YK is 0.27_{-0.23};{+0.28}(stat)-0.01+0.04(syst), and the relevant B0 and B+ branching fractions are reported.     

Search for CPT and lorentz violation in B0-B[over ]0 oscillations with dilepton events

We report results of a search for CPT and Lorentz violation in B(0)-B[over ](0) oscillations using inclusive dilepton events from 232 x 10(6) Upsilon(4S)-->BB[over ] decays recorded by the BABAR detector at the PEP-II B Factory at SLAC. We find 2.8sigma significance, compatible with no signal, for variations in the complex CPT violation parameter z at the Earth's sidereal frequency and extract values for the quantities Deltaa(micro) in the general Lorentz-violating standard-model extension. The spectral powers for variations in z over the frequency range 0.26 yr(-1) to 2.1 solar day(-1) are also compatible with no signal.     

Measurement of the B0-B-0 oscillation frequency with inclusive dilepton events

The B0-B-0 oscillation frequency has been measured with a sample of 23 x 10(6) BB- pairs collected with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. In this sample, we select events in which both B mesons decay semileptonically and use the charge of the leptons to identify the flavor of each B meson. A simultaneous fit to the decay time difference distributions for opposite- and same-sign dilepton events gives deltamd = 0.493+/-0.012(stat)+/-0.009(syst) ps-1.