Influence of turbulent scalar mixing physics on premixed flame propagation

The influence of reactive scalar mixing physics on turbulent premixed flame propagation is studied, within the framework of turbulent flame speed modelling, by comparing predictive ability of two algebraic flame speed models: one that includes all relevant physics and the other ignoring dilatation effects on reactive scalar mixing. This study is an extension of a previous work analysing and validating the former model. The latter is obtained by neglecting modelling terms that include dilatation effects: a direct effect because of density change across the flame front and an indirect effect due to dilatation on turbulence-scalar interaction. An analysis of the limiting behaviour shows that neglecting the indirect effect alters the flame speed scaling considerably when is small and the scaling remains unaffected when is large. This is evident from comparisons of the two models with experimental data which show that the quantitative difference between the two models is as high as 66% at but only 4% at . Furthermore, neglecting the direct effect results in a poor prediction of turbulent flame speed for all values of , and both effects are important for practically relevant values of this velocity ratio.Peer Reviewe

On velocity and reactive scalar spectra in turbulent premixed flames

AbstractKinetic energy and reactive scalar spectra in turbulent premixed flames are studied from compressible three-dimensional direct numerical simulations (DNS) of a temporally evolving rectangular slot-jet premixed flame, a statistically one-dimensional configuration. The flames correspond to a lean premixed hydrogen–air mixture at an equivalence ratio of 0.7, preheated to 700 K and at 1 atm, and three DNS are considered with a fixed jet Reynolds number of 10 000 and a jet Damköhler number varying between 0.13 and 0.54. For the study of spectra, motivated by the need to account for density change, which can be locally strong in premixed flames, a new density-weighted definition for two-point velocity/scalar correlations is proposed. The density-weighted two-point correlation tensor retains the essential properties of its constant-density (incompressible) counterpart and recovers the density-weighted Reynolds stress tensor in the limit of zero separation. The density weighting also allows the derivation of balance equations for velocity and scalar spectrum functions in the wavenumber space that illuminate physics unique to combusting flows. Pressure–dilatation correlation is a source of kinetic energy at high wavenumbers and, analogously, reaction rate–scalar fluctuation correlation is a high-wavenumber source of scalar energy. These results are verified by the spectra constructed from the DNS data. The kinetic energy spectra show a distinct inertial range with a $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}-5/3$ scaling followed by a ‘diffusive–reactive’ range at higher wavenumbers. The exponential drop-off in this range shows a distinct inflection in the vicinity of the wavenumber corresponding to a laminar flame thickness, $\delta _L$, and this is attributed to the contribution from the pressure–dilatation term in the energy balance in wavenumber space. Likewise, a clear spike in spectra of major reactant species (hydrogen) arising from the reaction-rate term is observed at wavenumbers close to $\delta _L$. It appears that in the inertial range classical scaling laws for the spectra involving the Kolmogorov scale are applicable, but in the high-wavenumber range where chemical reactions have a strong signature the laminar flame thickness produces a better collapse. It is suggested that a full scaling should perhaps involve the Kolmogorov scale, laminar flame thickness, Damköhler number and Karlovitz number.This is the accepted manuscript. The final version is available from CUP at http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=9317552&fileId=S0022112014003929

Topology of pocket formation in turbulent premixed flames

© 2018 The Combustion Institute. In turbulent premixed flame propagation, the formation of isolated pockets of reactants or products is associated with flame pinch-off events which cause rapid changes in the flame surface area. Previous topological analysis of these phenomena has been carried out based on Morse theory and Direct Numerical Simulation (DNS) in two spatial dimensions. The present work extends the topological analysis to three dimensions with emphasis on the formation and subsequent burnout of reactant pockets. Singular behaviour observed previously for terms of the Surface Density Function (SDF) transport equation in the two-dimensional case is shown to occur also in three dimensions. Further singular behaviour is observed in the displacement speed close to reactant pocket burnout. The theory is compared against DNS data from hydrogen-air flames and good agreement is obtained

Large Eddy Simulation of a dual swirl gas turbine combustor: Flame/flow structures and stabilisation under thermoacoustically stable and unstable conditions

A laboratory gas turbine model combustor with dual-swirler configuration is in- vestigated using Large Eddy Simulation (LES) with a flamelet subgrid combus- tion model. Two partially premixed methane/air flames with different equivalence ratio and thermal power are simulated: one stably burning with an elongated V- shape and another undergoing pronounced thermoacoustic oscillations exhibiting a flat shape. Additionally, both flames feature a hydrodynamic instability in the form of a precessing vortex core (PVC). Detailed comparisons between experi- mental and LES results show that the different flow and reaction zone structures in these two flames are reproduced well. The various flow dynamics resulting from the PVC and thermoacoustic oscillations are also captured accurately in the simulation. Further analyses on the lifted swirl flame stabilisation using phase averaged statistics at the PVC frequencies reveal that the PVC-induced stagnation points provide an anchoring mechanism for both the stable and unstable flames, although in the latter case large self-excited pressure oscillations are present. It is found that the PVC is significantly influenced by these oscillations, being axially stretched and compressed at high and low pressures, respectively. However, the formation of flame leading edge due to the PVC is robust during these unstable processes and the azimuthal movement of the leading point is found to be strongly correlated with the rotation of the PVC in both flames, further confirming the vital role of the PVC in the stabilisation process of these lifted swirl flames

Anteroposterior chest radiograph vs. chest CT scan in early detection of pneumothorax in trauma patients

Pneumothorax is a common complication following blunt chest wall trauma. In these patients, because of the restrictions regarding immobilization of the cervical spine, Anteroposterior (AP) chest radiograph is usually the most feasible initial study which is not as sensitive as the erect chest X-ray or CT chest for detection of a pneumothorax. We will present 3 case reports which serve for better understanding of the entity of occult pneumothorax. The first case is an example of a true occult pneumothorax where an initial AP chest X-ray revealed no evidence of pneumothorax and a CT chest immediately performed revealed evidence of pneumothorax. The second case represents an example of a missed rather than a truly occult pneumothorax where the initial chest radiograph revealed clues suggesting the presence of pneumothorax which were missed by the reading radiologist. The third case emphasizes the fact that "occult pneumothorax is predictable". The presence of subcutaneous emphesema and pulmonary contusion should call for further imaging with CT chest to rule out pneumothorax. Thoracic CT scan is therefore the "gold standard" for early detection of a pneumothorax in trauma patients. This report aims to sensitize readers to the entity of occult pneumothorax and create awareness among intensivists and ER physicians regarding the proper diagnosis and management

Anatomy and origin of authochthonous late Pleistocene forced regression deposits, east Coromandel inner shelf, New Zealand: implications for the development and definition of the regressive systems tract

High-resolution seismic reflection data from the east Coromandel coast, New Zealand, provide details of the sequence stratigraphy beneath an autochthonous, wave dominated inner shelf margin during the late Quaternary (0-140 ka). Since c. 1 Ma, the shelf has experienced limited subsidence and fluvial sediment input, producing a depositional regime characterised by extensive reworking of coastal and shelf sediments during glacio-eustatic sea-level fluctuations. It appears that only one complete fifth-order (c. 100 000 yr) depositional sequence is preserved beneath the inner shelf, the late Pleistocene Waihi Sequence, suggesting any earlier Quaternary sequences were mainly cannibalised into successively younger sequences. The predominantly Holocene-age Whangamata Sequence is also evident in seismic data and modern coastal deposits, and represents an incomplete depositional sequence in its early stages of formation. A prominent aspect of the sequence stratigraphy off parts of the east Coromandel coast is the presence of forced regressive deposits (FRDs) within the regressive systems tract (RST) of the late Pleistocene Waihi Sequence. The FRDs are interpreted to represent regressive barrier-shoreface sands that were sourced from erosion and onshore reworking of underlying Pleistocene sediments during the period of slow falling sea level from isotope stages 5 to 2 (c. 112-18 ka). The RST is volumetrically the most significant depositional component of the Waihi Sequence; the regressive deposits form a 15-20 m thick, sharp-based, tabular seismic unit that downsteps and progrades continuously across the inner shelf. The sequence boundary for the Waihi Sequence is placed at the most prominent, regionally correlative, and chronostratigraphically significant surface, namely an erosional unconformity characterised in many areas by large incised valleys that was generated above the RST. This unconformity is interpreted as a surface of maximum subaerial erosion generated during the last glacial lowstand (c. 18 ka). Although the base of the RST is associated with a prominent regressive surface of erosion, this is not used as the sequence boundary as it is highly diachronous and difficult to identify and correlate where FRDs are not developed. The previous highstand deposits are limited to subaerial barrier deposits preserved behind several modern Holocene barriers along the coast, while the transgressive systems tract is preserved locally as incised-valley fill deposits beneath the regressive surface of erosion at the base of the RST. Many documented late Pleistocene RSTs have been actively sourced from fluvial systems feeding the shelf and building basinward-thickening, often stacked wedges of FRDs, for which the name allochthonous FRDs is suggested. The Waihi Sequence RST is unusual in that it appears to have been sourced predominantly from reworking of underlying shelf sediments, and thus represents an autochthonous FRD. Autochthonous FRDs are also present on the Forster-Tuncurry shelf in southeast Australia, and may be a common feature in other shelf settings with low subsidence and low sediment supply rates, provided shelf gradients are not too steep, and an underlying source of unconsolidated shelf sediments is available to source FRDs. The preservation potential of such autochthonous FRDs in ancient deposits is probably low given that they are likely to be cannibalised during subsequent sea-level falls

Neodymium isotope constraints on provenance, dispersal, and climate-driven supply of Zambezi sediments along the Mozambique Margin during the past ∼45,000 years

Marine sediments deposited off the Zambezi River that drains a considerable part of the southeast African continent provide continuous records of the continental climatic and environmental conditions. Here we present time series of neodymium (Nd) isotope signatures of the detrital sediment fraction during the past ~45,000 years, to reconstruct climate-driven changes in the provenance of clays deposited along the Mozambique Margin. Coherent with the surface current regime, the Nd isotope distribution in surface sediments reveals mixing of the alongshore flowing Zambezi suspension load with sediments supplied by smaller rivers located further north. To reconstruct past changes in sediment provenances, Nd isotope signatures of clays that are not significantly fractionated during weathering processes have been obtained from core 64PE304-80, which was recovered just north of the Zambezi mouth at 1329 m water depth. Distinctly unradiogenic clay signatures (ENd values <214.2) are found during the Last Glacial Maximum, Heinrich Stadial 1, and Younger Dryas. In contrast, the Nd isotope record shows higher, more radiogenic isotope signatures during Marine Isotope Stage 3 and between ~15 and ~5 ka BP, the latter coinciding with the timing of the northern hemisphere African Humid Period. The clay-sized sediment fraction with the least radiogenic Nd isotope signatures was deposited during the Holocene, when the adjacent Mozambique Shelf became completely flooded. In general, the contribution of the distinctly unradiogenic Zambezi suspension load has followed the intensity of precession-forced monsoonal precipitation and enhanced during periods of increased southern hemisphere insolation and high-latitude northern hemispheric climate variability

Theoretical study on the electronic, structural, properties and reactivity of a series of mono-, di-, tri- and tetrachlorothiophenes as well as corresponding radical cation forms as monomers for conducting polymers

In this paper, electrical and structural properties of mono-, di-, tri- and tetrachlorothiophenes and their radical cations have been studied using the density functional theory and B3LYP method with 6-311++G** basis set. The effects of the number and position of the substituent of chlorine atoms on the properties of the thiophene ring for all chlorothiophenes and their radical cations have been studied. Vibrational frequencies, nuclear chemical shielding constants, spin-density distribution, size and direction of dipole moment vector, ionization potential, electric polarizabilities and NICS values of these compounds have been calculated as well. The analysis of these data showed that double bonds in 3-chlorothiophene are more delocalized and it is the best possible candidate monomer among all chlorothiophenes for the synthesis of corresponding conducting polymers with modified characteristics