Investigation of conditional source-term estimation applied to a non-premixed turbulent flame

A turbulent combustion model, Conditional Source-term Estimation (CSE) is applied to a non-premixed turbulent jet methane flame. The conditional chemical source terms are determined on the basis of first order closure and the conditional averaged species concentrations are obtained by inverting an integral equation. The Tikhonov method is implemented for regularisation. Detailed chemistry is tabulated using the trajectory generated low-dimensional manifold method. Radiation due to the gaseous species is included. Reynolds Averaged Navier–Stokes calculations are performed using two different turbulence models. The objectives of the paper are (i) assessment of the impact of the main numerical parameters in CSE and (ii) comparison of the CSE numerical predictions with available experimental data and results from previous simulations for the selected flame. The number of CSE domains and the number of points in each CSE domain are shown to have a significant impact on the results if not selected appropriately. The present CSE calculations always converge to unique and stable predictions. The corrected k–ε model yields mixture fraction profiles in good agreement with the experimental data values for axial locations in the first half of the flame. Farther downstream, the RNG k–ε model performs better. Overall, the current predictions for the mixture fraction are in good agreement with the experimental data. The predicted temperatures using CSE and the k–ε turbulence model with a modified value of C_ε1 = 1.47 are found to be in very good agreement with the experimental data. Further, the current CSE results are of comparable quality with previous simulations using the flamelet model and conditional moment closure. Future work may include further investigation on optimal determination of the regularisation parameter and alternative regularisation techniques, soot modelling within the CSE formulation, and improved formulation of radiation.     

Inverse analysis and regularisation in conditional source-term estimation modelling

Conditional Source-term Estimation (CSE) obtains the conditional species mass fractions by inverting a Fredholm integral equation of the first kind. In the present work, a Bayesian framework is used to compare two different regularisation methods: zeroth-order temporal Tikhonov regulatisation and first-order spatial Tikhonov regularisation. The objectives of the current study are: (i) to elucidate the ill-posedness of the inverse problem; (ii) to understand the origin of the perturbations in the data and quantify their magnitude; (iii) to quantify the uncertainty in the solution using different priors; and (iv) to determine the regularisation method best suited to this problem. A singular value decomposition shows that the current inverse problem is ill-posed. Perturbations to the data may be caused by the use of a discrete mixture fraction grid for calculating the mixture fraction PDF. The magnitude of the perturbations is estimated using a box filter and the uncertainty in the solution is determined based on the width of the credible intervals. The width of the credible intervals is significantly reduced with the inclusion of a smoothing prior and the recovered solution is in better agreement with the exact solution. The credible intervals for temporal and spatial smoothing are shown to be similar. Credible intervals for temporal smoothing depend on the solution from the previous time step and a smooth solution is not guaranteed. For spatial smoothing, the credible intervals are not dependent upon a previous solution and better predict characteristics for higher mixture fraction values. These characteristics make spatial smoothing a promising alternative method for recovering a solution from the CSE inversion process.     

NO2 sensing properties of YSZ-based sensor using NiO and Cr-doped NiO sensing electrodes at high temperature

Nickel oxide and chromium-doped nickel oxide (Ni_0.95Cr_0.03O_1−δ ) were prepared by thermal decomposition of nitrates. The obtained NiO and Ni_0.95Cr_0.03O_1−δ samples were utilized as sensing electrodes (SEs) in yttria-stabilized zirconia (YSZ)-based sensors for detection of NO₂ at 800 °C under wet condition (5 vol.% H₂O). While the mixed-potential-type planar sensor attached with NiO-SE gave rather large NO₂ sensitivity, the sensor attached with Ni_0.95Cr_0.03O_1−δ -SE exhibited fast recovery rate with an acceptable sensitivity. The Δemf (electromotive force) of the sensors varied linearly with NO₂ concentration in the examined range of 50–400 ppm on a logarithmic scale. Based on the results of measurements for polarization, complex impedance and gas phase catalysis, the fast recovery was attributable to the high rate for the anodic reaction of O₂ at the Ni_0.95Cr_0.03O_1−δ /YSZ interface, and the lower NO₂ sensitivity was caused by both the high rate for the anodic reaction of O₂ and the high degree for the gas phase conversion of NO₂ to NO.     

Validation of flow simulation and gas combustion sub-models for the CFD-based prediction of NOx formation in biomass grate furnaces

While reasonably accurate in simulating gas phase combustion in biomass grate furnaces, CFD tools based on simple turbulence–chemistry interaction models and global reaction mechanisms have been shown to lack in reliability regarding the prediction of NO_x formation. Coupling detailed NO_x reaction kinetics with advanced turbulence–chemistry interaction models is a promising alternative, yet computationally inefficient for engineering purposes. In the present work, a model is proposed to overcome these difficulties. The model is based on the Realizable k–? model for turbulence, Eddy Dissipation Concept for turbulence–chemistry interaction and the HK97 reaction mechanism. The assessment of the sub-models in terms of accuracy and computational effort was carried out on three laboratory-scale turbulent jet flames in comparison with the experimental data. Without taking NO_x formation into account, the accuracy of turbulence modelling and turbulence–chemistry interaction modelling was systematically examined on Sandia Flame D and Sandia CO/H₂/N₂ Flame B to support the choice of the associated models. As revealed by the Large Eddy Simulations of the former flame, the shortcomings of turbulence modelling by the Reynolds averaged Navier–Stokes (RANS) approach considerably influence the prediction of the mixing-dominated combustion process. This reduced the sensitivity of the RANS results to the variations of turbulence–chemistry interaction models and combustion kinetics. Issues related to the NO_x formation with a focus on fuel bound nitrogen sources were investigated on a NH₃-doped syngas flame. The experimentally observed trend in NO_x yield from NH₃ was correctly reproduced by HK97, whereas the replacement of its combustion subset by that of a detailed reaction scheme led to a more accurate agreement, but at increased computational costs. Moreover, based on results of simulations with HK97, the main features of the local course of the NO_x formation processes were identified by a detailed analysis of the interactions between the nitrogen chemistry and the underlying flow field.     

Raman scattering and room temperature ferromagnetism in Co-doped SrTiO3 particles

Raman scattering has been used to study the influence of cobalt, an effective dopant to obtain SrTiO₃ magnetic oxide, on the lattice dynamics of SrTiO₃. It is found that Co doping increases the lattice defects and induces a Raman vibration mode of 690 cm⁻¹. On the other hand, the ferromagnetism dependence on the x and annealing temperature was clearly and coherently observed in SrTi_1−xCo_xO₃ (x = 0, 0.01, 0.03 and 0.05) nanoparticles. It is found that the ferromagnetism of SrTi_1−xCo_xO₃ nanoparticles is weakly related to crystal deformation and oxygen vacancies in SrTiO₃. So, F-center model can explain the origin of the ferromagnetism in the prepared Co-doped SrTiO₃ samples. At the same time, the finding of large room-temperature ferromagnetism (1.6 emu/g) in this system would stimulate further interest in the area of more complicated ternary oxides.     

The effect of hydrogen addition on flammability limit and NOx emission in ultra-lean counterflow CH4/air premixed flames

The effect of hydrogen addition to ultra lean counterflow CH₄/air premixed flames on the extinction limits and the characteristics of NO_x emission was investigated by numerical simulation. Detailed chemistry and complex thermal and transport properties were employed. The results show that the addition of hydrogen can significantly enlarge the flammable region and extend the flammability limit to lower equivalence ratios. If the equivalence ratio is kept constant, the addition of hydrogen increases the emission of NO in a flame due to the enhancement in the rate of the NNH or N₂O intermediate NO formation routes. The addition of hydrogen causes a monotonic decrease in the formation of NO₂ and N₂O, except flames near the extinction limits, where the emission of NO₂ and N₂O first increases, and then decreases with the increase in the fraction of hydrogen. Overall, hydrogen enrichment technology allows stable combustion under ultra lean conditions, resulting in significant CO₂ and NO emission reduction.     

Field and temperature induced colossal strain in Gd5(SixGe1−x)4

Gd₅(Si_xGe_1−x)₄, known for its giant magnetocaloric effect, also exhibits a colossal strain of the order of 10,000 ppm for a single crystal near its coupled first-order magnetic-structural phase transition, which occurs near room temperature for the compositions 0.41≤x≤0.575. Such colossal strain can be utilised for both magnetic sensor and actuator applications. In this study, various measurements have been carried out on strain as a function of magnetic field strength and as a function of temperature on single crystal Gd₅Si₂Ge₂ (x=0.5), and polycrystalline Gd₅Si_1.95Ge_2.05 (x=0.487) and Gd₅Si_2.09Ge_1.91 (x=0.52). Additionally a giant magnetostriction/thermally induced strain of the order of 1800 ppm in polycrystalline Gd₅Si_2.09Ge_1.91 was observed at its first order phase transition on varying temperature using a Peltier cell without the use of bulky equipment such as cryostat or superconducting magnet.     

Air-broadening parameters in the ν3 band of NO2 using a multispectrum fitting technique

Air-broadened linewidths, pressure-induced shift coefficients and their temperature dependences were retrieved for over 1000 transitions in the ν₃ band of ¹⁴N¹⁶O₂ at 6 μm. In addition, precise line center positions and relative intensities were also determined. The results were obtained by fitting simultaneously 27 spectra recorded at high resolution (0.002-0.006 cm⁻¹) with two Fourier transform spectrometers and gas sample temperatures ranging from 206 to 298 K. It was necessary to modify the multispectrum fitting software to accommodate constraints on the retrieved parameters of closely spaced spin-split doublets in order to successfully determine their broadening and shift parameters. The variations of the widths, shifts and their temperature dependences with the quantum numbers were investigated. Subsets of the observed linewidths were reproduced to within 3% using an empirical smoothing function.     

Temperature mapping in bread dough using SE and GE two-point MRI methods: experimental and theoretical estimation of uncertainty

Two-dimensional (2D)-SE, 2D-GE and tri-dimensional (3D)-GE two-point T(1)-weighted MRI methods were evaluated in this study in order to maximize the accuracy of temperature mapping of bread dough during thermal processing. Uncertainties were propagated throughout each protocol of measurement, and comparisons demonstrated that all the methods with comparable acquisition times minimized the temperature uncertainty to similar extent. The experimental uncertainties obtained with low-field MRI were also compared to the theoretical estimations. Some discrepancies were reported between experimental and theoretical values of uncertainties of temperature; however, experimental and theoretical trends with varying parameters agreed to a large extent for both SE and GE methods. The 2D-SE method was chosen for further applications on prefermented dough because of its lower sensitivity to susceptibility differences in porous media. It was applied for temperature mapping in prefermented dough during chilling prior to freezing and compared locally to optical fiber measurements.     

The methane Raman spectrum from 1200 to 5500 cm: A first step toward temperature diagnostic using methane as a probe molecule in combustion systems

We present a study of the spontaneous Raman spectra of ¹²CH₄ from 1200 to 5500 cm⁻¹ at various temperatures. This study is of interest from a fundamental as well as from a practical point of view with regards to the temperature diagnostic in hydrocarbon combustion. The present investigation shows that the spontaneous ¹²CH₄ Raman spectra are very sensitive to temperature and that the complexity of methane spectra is not an obstacle to use methane as a probe molecule in laser-diagnostic techniques. Our study consists in determining the polarisability parameters of methane ¹²CH₄, unknown at the present time, from spontaneous Raman spectra recorded at pressure less than one atmosphere to minimize the problem of line-mixing. This diagnostic method technique based on spontaneous Raman spectroscopy may be applied to temperature measurements in a uni-element rocket combustion chamber using methane-oxygen propellants.     

Formation pathways of HO2 and OH changing as a function of temperature in photolytically initiated oxidation of dimethyl ether

The time resolved product formation in oxidation of dimethyl ether (DME) has been studied between 298-625 K and 20-90 torr total pressure. Near-infrared frequency modulation spectroscopy (FMS) with Herriott type multi pass optics and UV absorption spectroscopy (UV) were conducted in the same cell. The reaction was initiated by pulsed photolysis in a mixture of Cl₂, O₂, and DME via CH₃OCH₂ radical formation. The reaction process was investigated through FMS measurement of HO₂ and OH, and UV measurement of CH₃OCH₂O₂. The yields of HO₂ and OH are obtained by comparison with reference mixtures, Cl₂, O₂, and CH₃OH for HO₂, and Cl₂, O₂, CH₃OH, and NO for OH, which convert 100% of initial Cl to HO₂ and OH. The CH₃OCH₂O₂ yield is also obtained. It was found that the HO₂ yield increases sharply over 500 K mainly with a longer time constant than that of R + O₂ reaction, while a prompt component exists throughout the temperature range at a few percent yield. OH was found to be produced promptly at a yield considerably larger than that known for the simplest alkanes. The CH₃OCH₂O₂ profile has a prompt rise followed by a gradual decay whose rate is consistent with the slow HO₂ formation. The species profiles were successfully predicted with a model constructed by modifying the existing one to suit the reduced pressure condition. After modification, it was inferred that the HO₂ formation over 500 K is secondary from HCHO + OH and HCO + O₂ and a part of HCO is formed directly from the O₂ adduct, whereas the HO₂ formation below 500 K is governed by CH₃OCH₂O₂ chemistry. The HCO forming pathway via isomerization-decomposition of the O₂ adduct, which was not included in the former models, was supported by our quantum-chemical calculations.     

Ground states of the hydrogen molecule and its molecular ion in the presence of a magnetic field using the variational Monte Carlo method

ABSTRACT

By using the variational Monte Carlo (VMC) method, we calculated the 1sσ_g-state energies, the dissociation energies, and the binding energies of the hydrogen molecule and its molecular ion in the presence of an aligned magnetic field regime between 0 and 10?a.u. The present calculations are based on using two types of compact and accurate trial wave functions, which are put forward for consideration in calculating energies in the absence of a magnetic field. The obtained results are compared with the most recent accurate values. We conclude that the applications of the VMC method can be successfully extended to cover the case of molecules under the effect of a magnetic field.     

Lorentz half-width, pressure-induced shift and speed-dependent coefficients in oxygen-broadened CO2 bands at 6227 and 6348 cm using a constrained multispectrum analysis

The McMath-Pierce Fourier transform spectrometer located at the National Solar Observatory (NSO) on Kitt Peak, Arizona, was used to record infrared high resolution absorption spectra of CO₂ spectra broadened by O₂. These spectra were analyzed to measure O₂-broadened half-width coefficients, O₂-induced pressure-shift coefficients and speed dependent parameters for transitions in the 30013←00001 and 30012←00001 bands of ¹⁶O¹²C¹⁶O located near 6227 and 6348 cm⁻¹, respectively. All spectra were obtained at room temperature using the long path, 6 m base path White cell available at NSO. A multispectrum nonlinear least-squares fitting algorithm employing Voigt line shapes modified to include line mixing and speed dependence was used to fit simultaneously a total of 19 spectra in the 6120-6280 cm⁻¹ (30013←00001) and 6280-6395 cm⁻¹ (30012←00001) spectral regions. 16 of the 19 spectra analyzed in this work were self broadened and three spectra were lean mixtures of CO₂ in O₂. The volume mixing ratios of CO₂ in the three spectra varied between 0.06 and 0.1. Lorentz half-width and pressure-induced shift coefficients were measured for all transitions in the P(50)-R(50) range in both vibrational bands. The results obtained from present analysis have been compared with measurements available in the literature for self-, air-, oxygen- and argon-broadening. No significant differences were observed between the broadening and shift coefficients of the two bands. The N₂-broadened half-width and pressure-shift coefficients were computed from measured air- and O₂-broadened width and shift coefficients.