Aerothermodynamic cycle analysis of a dual-spool, separate-exhaust turbofan engine with an interstage turbine burner

This study focuses on a specific engine, i.e., a dual-spool, separate-flow turbofan engine with an Interstage Turbine Burner (ITB). This conventional turbofan engine has been modified to include a secondary isobaric burner, i.e., ITB, in a transition duct between the high-pressure turbine and the low-pressure turbine. The preliminary design phase for this modified engine starts with the aerothermodynamics cycle analysis is consisting of parametric (i.e., on-design) and performance (i.e., off-design) cycle analyses. In parametric analysis, the modified engine performance parameters are evaluated and compared with baseline engine in terms of design limitation (maximum turbine inlet temperature), flight conditions (such as flight Mach condition, ambient temperature and pressure), and design choices (such as compressor pressure ratio, fan pressure ratio, fan bypass ratio etc.). A turbine cooling model is also included to account for the effect of cooling air on engine performance. The results from the on-design analysis confirmed the advantage of using ITB, i.e., higher specific thrust with small increases in thrust specific fuel consumption, less cooling air, and less NOx production, provided that the main burner exit temperature and ITB exit temperature are properly specified. It is also important to identify the critical ITB temperature, beyond which the ITB is turned off and has no advantage at all. With the encouraging results from parametric cycle analysis, a detailed performance cycle analysis of the identical engine is also conducted for steady-stateengine performance prediction. The results from off-design cycle analysis show that the ITB engine at full throttle setting has enhanced performance over baseline engine. Furthermore, ITB engine operating at partial throttle settings will exhibit higher thrust at lower specific fuel consumption and improved thermal efficiency over the baseline engine. A mission analysis is also presented to predict the fuel consumptions in certain mission phases. Excel macrocode, Visual Basic for Application, and Excel neuron cells are combined to facilitate Excel software to perform these cycle analyses. These user-friendly programs compute and plot the data sequentially without forcing users to open other types of post-processing programs.

Can the carbon isotopic composition of methane be reconstructed from multi-site firn air measurements?

Methane is a strong greenhouse gas and large uncertainties exist concerning the future evolution of its atmospheric abundance. Analyzing methane atmospheric mixing and stable isotope ratios in air trapped in polar ice sheets helps in reconstructing the evolution of its sources and sinks in the past. This is important to improve predictions of atmospheric CH4 mixing ratios in the future under the influence of a changing climate. The aim of this study is to assess whether past atmospheric δ13C(CH4) variations can be reliably reconstructed from firn air measurements. Isotope reconstructions obtained with a state of the art firn model from different individual sites show unexpectedly large discrepancies and are mutually inconsistent. We show that small changes in the diffusivity profiles at individual sites lead to strong differences in the firn fractionation, which can explain a large part of these discrepancies. Using slightly modified diffusivities for some sites, and neglecting samples for which the firn fractionation signals are strongest, a combined multi-site inversion can be performed, which returns an isotope reconstruction that is consistent with firn data. However, the isotope trends are lower than what has been concluded from Southern Hemisphere (SH) archived air samples and high-accumulation ice core data. We conclude that with the current datasets and understanding of firn air transport, a high precision reconstruction of δ13C of CH4 from firn air samples is not possible, because reconstructed atmospheric trends over the last 50 yr of 0.3–1.5 ‰ are of the same magnitude as inherent uncertainties in the method, which are the firn fractionation correction (up to ~2 ‰ at individual sites), the Kr isobaric interference (up to ~0.8 ‰, system dependent), inter-laboratory calibration offsets (~0.2 ‰) and uncertainties in past CH4 levels (~0.5 ‰).

A rapid and efficient ion-exchange chromatography for Lu–Hf, Sm–Nd, and Rb–Sr geochronology and the routine isotope analysis of sub-ng amounts of Hf by MC-ICP-MS

The development and improvement of MC-ICP-MS instruments have fueled the growth of Lu–Hf geochronology over the last two decades, but some limitations remain. Here, we present improvements in chemical separation and mass spectrometry that allow accurate and precise measurements of 176Hf/177Hf and 176Lu/177Hf in high-Lu/Hf samples (e.g., garnet and apatite), as well as for samples containing sub-nanogram quantities of Hf. When such samples are spiked, correcting for the isobaric interference of 176Lu on 176Hf is not always possible if the separation of Lu and Hf is insufficient. To improve the purification of Hf, the high field strength elements (HFSE, including Hf) are first separated from the rare earth elements (REE, including Lu) on a first-stage cation column modified after Patchett and Tatsumoto (Contrib. Mineral. Petrol., 1980, 75, 263–267). Hafnium is further purified on an Ln-Spec column adapted from the procedures of Münker et al. (Geochem., Geophys., Geosyst., 2001, DOI: 10.1029/2001gc000183) and Wimpenny et al. (Anal. Chem., 2013, 85, 11258–11264) typically resulting in Lu/Hf < 0.0001, Zr/Hf < 1, and Ti/Hf < 0.1. In addition, Sm–Nd and Rb–Sr separations can easily be added to the described two-stage ion-exchange procedure for Lu–Hf. The isotopic compositions are measured on a Thermo Scientific Neptune Plus MC-ICP-MS equipped with three 1012 Ω resistors. Multiple 176Hf/177Hf measurements of international reference rocks yield a precision of 5–20 ppm for solutions containing 40 ppb of Hf, and 50–180 ppm for 1 ppb solutions (=0.5 ng sample Hf 0.5 in ml). The routine analysis of sub-ng amounts of Hf will facilitate Lu–Hf dating of low-concentration samples.

Contrasting PT paths within the Barchi-Kol UHP terrain (Kokchetav Complex): Implications for subduction and exhumation of continental crust

The Barchi-Kol terrain is a classic locality of ultrahigh-pressure (UHP) metamorphism within the Kokchetav metamorphic belt. We provide a detailed and systematic characterization of four metasedimentary samples using dominant mineral assemblages, mineral inclusions in zircon and monazite, garnet zonation with respect to major and trace elements, and Zr-in-rutile and Ti-in-zircon temperatures. A typical diamond-bearing gneiss records peak conditions of 49 ± 4 kbar and 950–1000 °C. Near isothermal decompression of this rock resulted in the breakdown of phengite associated with a pervasive recrystallization of the rock. The same terrain also contains mica schists that experienced peak conditions close to those of the diamond-bearing rocks, but they were exhumed along a cooler path where phengite remained stable. In these rocks, major and trace element zoning in garnet has been completely equilibrated. A layered gneiss was metamorphosed at UHP conditions in the coesite field, but did not reach diamond-facies conditions (peak conditions: 30 kbar and 800–900 °C). In this sample, garnet records retrograde zonation in major elements and also retains prograde zoning in trace elements. A garnet-kyanite-micaschist that reached significantly lower pressures (24 ± 2 kbar, 710 ± 20 °C) contains garnet with major and trace element zoning. The diverse garnet zoning in samples that experienced different metamorphic conditions allows to establish that diffusional equilibration of rare earth element in garnet likely occurs at ~900–950 °C. Different metamorphic conditions in the four investigated samples are also documented in zircon trace element zonation and mineral inclusions in zircon and monazite. U-Pb geochronology of metamorphic zircon and monazite domains demonstrates that prograde (528–521 Ma), peak (528–522 Ma), and peak to retrograde metamorphism (503–532 Ma) occurred over a relatively short time interval that is indistinguishable from metamorphism of other UHP rocks within the Kokchetav metamorphic belt. Therefore, the assembly of rocks with contrasting P-T trajectories must have occurred in a single subduction-exhumation cycle, providing a snapshot of the thermal structure of a subducted continental margin prior to collision. The rocks were initially buried along a low geothermal gradient. At 20–25 kbar they underwent near isobaric heating of 200 °C, which was followed by continued burial along a low geothermal gradient. Such a step-wise geotherm is in good agreement with predictions from subduction zone thermal models.

Chemical and germanium isotopic compositions of sediments and sedimentary rocks from ODP Sites 129-801 and 185-1149

A new technique for the precise and accurate determination of Ge stable isotope compositions has been developed and applied to silicate rocks and biogenic opal. The analyses were performed using a continuous flow hydride generation system coupled to a MC-ICP-MS. Samples have been purified through anion- and cation-exchange resins to separate Ge from matrix elements and eliminate potential isobaric interferences. Variations of 74Ge/70Ge ratios are expressed as d74Ge values relative to our internal standard and the long-term external reproducibility of the data is better than 0.2? for sample size as low as 15 ng of Ge. Data are presented for igneous and sedimentary rocks, and the overall variation is 2.4? in d74Ge, representing 12 times the uncertainty of the measurements and demonstrating that the terrestrial isotopic composition of Ge is not unique. Co-variations of 74Ge/70Ge, 73Ge/70Ge and 72Ge/70Ge ratios follow a mass-dependent behaviour and imply natural isotopic fractionation of Ge by physicochemical processes. The range of d74Ge in igneous rocks is only 0.25? without systematic differences among continental crust, oceanic crust or mantle material. On this basis, a Bulk Silicate Earth reservoir with a d74Ge of 1.3+/-0.2? can be defined. In contrast, modern biogenic opal such as marine sponges and authigenic glauconite displayed higher d74Ge values between 2.0? and 3.0?. This suggests that biogenic opal may be significantly enriched in light isotopes with respect to seawater and places a lower bound on the d74Ge of the seawater to +3.0?.This suggests that seawater is isotopically heavy relative to Bulk Silicate Earth and that biogenic opal may be significantly fractionated with respect to seawater. Deep-sea sediments are within the range of the Bulk Silicate Earth while Mesozoic deep-sea cherts (opal and quartz) have d74Ge values ranging from 0.7? to 2.0?. The variable values of the cherts cannot be explained by binary mixing between a biogenic component and a detrital component and are suggestive of enrichment in the light isotope of diagenetic quartz. Further work is now required to determine Ge isotope fractionation by siliceous organisms and to investigate the effect of diagenetic processes during chert lithification.

Pollen distribution in marine sediment samples along the south-western African coast

The distribution of pollen in marine sediments is used to reconstruct pathways of terrigenous input to the oceans and provides a record of vegetation change on adjacent continents. The wind transport routes of aeolian pollen is comprehensively illustrated by clusters of trajectories. Isobaric, 4-day backward trajectories are calculated using the modelled wind-field of ECHAM3, and are clustered on a seasonal basis to estimate the main pathways of aeolian particles to sites of marine cores in the south-eastern Atlantic. Trajectories and clusters based on the modelled wind-field of the Last Glacial Maximum hardly differ from those of the present-day. Trajectory clusters show three regional, and two seasonal patterns, determining the pathways of aeolian pollen transport into the south-eastern Atlantic ocean. Mainly, transport out of the continent occurs during austral fall and winter, when easterly and south-easterly winds prevail. South of 25°S, winds blow mostly from the west and southwest, and aeolian terrestrial input is very low. Generally, a good latitudinal correspondence exists between the distribution patterns of pollen in marine surface sediments and the occurrence of the source plants on the adjacent continent. The northern Angola Basin receives pollen and spores from the Congolian and Zambezian forests mainly through river discharge. The Zambezian vegetation zone is the main source area for wind-blown pollen in sediments of the Angola Basin, while the semi-desert and desert areas are the main sources for pollen in sediments of the Walvis Basin and on the Walvis Ridge. A transect of six marine pollen records along the south-western African coast indicates considerable changes in the vegetation of southern Africa between glacial and interglacial periods. Important changes in the vegetation are the decline of forests in equatorial Africa and the north of southern Africa and a northward shift of winter rain vegetation along the western escarpment.

Geochemical composition and crystallization conditions of melts from the Mid Atlantic Ridge near Ascension Island (7-11°S)

Phase equilibria simulations were performed on naturally quenched basaltic glasses to determine crystallization conditions prior to eruption of magmas at the Mid-Atlantic Ridge (MAR) east of Ascension Island (7°11°S).The results indicate that midocean ridge basalt (MORB) magmas beneath different segments of the MAR have crystallized over a wide range of pressures (100-900MPa). However, each segment seems to have a specific crystallization history. Nearly isobaric crystallization conditions (100-300MPa) were obtained for the geochemically enriched MORB magmas of the central segments, whereas normal (N)-MORB magmas of the bounding segments are characterized by polybaric crystallization conditions (200-900MPa). In addition, our results demonstrate close to anhydrous crystallization conditions of N-MORBs, whereas geochemically enriched MORBs were successfully modeled in the presence of 0.4-1wt% H2O in the parental melts.These estimates are in agreement with direct (Fourier transform IR) measurements of H2O abundances in basaltic glasses and melt inclusions for selected samples. Water contents determined in the parental melts are in the range 0.04-0.09 and 0.30-0.55 wt% H2O for depleted and enriched MORBs, respectively. Our results are in general agreement (within ±200MPa) with previous approaches used to evaluate pressure estimates in MORB. However, the determination of pre-eruptive conditions of MORBs, including temperature and water content in addition to pressure, requires the improvement of magma crystallization models to simulate liquid lines of descent in the presence of small amounts of water. KEY WORDS: MORB; Mid-Atlantic Ridge; depth of crystallization; water abundances; phase equilibria calculations; cotectic crystallization; pressure estimates; polybaric fractionation

Application of capillary electrophoresis-electrospray ionization mass spectometry in the determination of molecular diversity.

By means of capillary electrophoresis coupled online to electrospray ionization MS, a library of theoretically 171 disubstituted xanthene derivatives was analyzed. The method allowed the purity and makeup of the library to be determined: 160 of the expected compounds were found to be present, and 12 side-products were also detected in the mixture. Due to the ability of capillary electrophoresis to separate analytes on the basis of charge, most of the xanthene derivatives could be resolved by simple capillary electrophoresis-MS procedures even though 124 of the 171 theoretical compounds were isobaric with at least one other molecule in the mixture. Any remaining unresolved peaks were resolved by MS/MS experiments. The method shows promise for the analysis of small combinatorial libraries with fewer than 1000 components.

Pitfalls in the Evaluation of the Thermodynamic Consistency of Experimental VLE Data Sets

The thermodynamic consistency of almost 90 VLE data series, including isothermal and isobaric conditions for systems of both total and partial miscibility in the liquid phase, has been examined by means of the area and point-to-point tests. In addition, the Gibbs energy of mixing function calculated from these experimental data has been inspected, with some rather surprising results: certain data sets exhibiting high dispersion or leading to Gibbs energy of mixing curves inconsistent with the total or partial miscibility of the liquid phase, surprisingly, pass the tests. Several possible inconsistencies in the tests themselves or in their application are discussed. Related to this is a very interesting and ambitious initiative that arose within the NIST organization: the development of an algorithm to assess the quality of experimental VLE data. The present paper questions the applicability of two of the five tests that are combined in the algorithm. It further shows that the deviation of the experimental VLE data from the correlation obtained by a given model, the basis of some point-to-point tests, should not be used to evaluate the quality of these data.

Phase equilibria of the water + 1-butanol + toluene ternary system at 101.3 kPa

Isobaric vapour–liquid and vapour–liquid–liquid equilibrium data for the water + 1-butanol + toluene ternary system were measured at 101.3 kPa with a modified VLE 602 Fischer apparatus. In addition, the liquid–liquid equilibrium data at 313.15 K were measured and compared with data from other authors at different temperatures. The system exhibits a ternary heterogeneous azeotrope whose temperature and composition have been determined by interpolation. The thermodynamic consistency of the experimental vapour–liquid and vapour–liquid–liquid data was checked by means of the Wisniak’s Li/Wi consistency test. Moreover, the vapour–liquid and the liquid–liquid equilibrium correlation for the ternary system with NRTL and UNIQUAC models, together with the prediction made with the UNIFAC model, were studied and discussed.

Improving some thermodynamic consistency tests for VLE data evaluation

Thermodynamics Conference 2013 (Statistical Mechanics and Thermodynamics Group of the Royal Society of Chemistry), The University of Manchester, 3-6 September 2013.

Molecular dynamics simulation of organic-inorganic nanocomposites: Layering behavior and interlayer structure of organoclays

Isothermal-isobaric (NPT) molecular dynamics simulation has been performed to investigate the layering behavior and structure of nanoconfined quaternary alkylammoniums in organoclays. This work is focused on systems consisting of two clay layers and a number of alkylammoniums, and involves the use of modified Dreiding force field. The simulated basal spacings of organoclays agree satisfactorily with the experimental results in the literature. The atomic density profiles in the direction normal to the clay surface indicate that the alkyl chains within the interlayer space of montmorillonite exhibit an obvious layering behavior. The headgroups of long alkyl chains are distributed within two layers close to the clay surface, whereas the distributions of methyl and methylene groups are strongly dependent on the alkyl chain length and clay layer charge. Monolayer, bilayer, and pseudo-trilayer structures are found in organoclays modified with single long alkyl chains, which are identical to the structural models based on the measured basal spacings. A pseudo-quadrilayer structure, for the first time to our knowledge, is also identified in organoclays with double long alkyl chains. In the mixture structure of paraffin-type and multilayer, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer in pseudo-trilayer as well as next nearest layer in pseudo-quadrilayer.

Molecular dynamics simulation of the structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays

The structural and dynamic properties of dioctadecyldimethylammoniums (DODDMA) intercalated into 2:1 layered clays are investigated using isothermal-isobaric (NPT) molecular dynamics (MD) simulation. The simulated results are in reasonably good agreement with the available experimental measurements, such as X-ray diffraction (XRD), atom force microscopy (AFM), Fourier transform infrared (FTIR), and nuclear magnetic resonance (NMR) spectroscopies. The nitrogen atoms are found to be located mainly within two layers close to the clay surface whereas methylene groups form a pseudoquadrilayer structure. The results of tilt angle and order parameter show that interior two-bond segments of alkyl chains prefer an arrangement parallel to the clay surface, whereas the segments toward end groups adopt a random orientation. In addition, the alkyl chains within the layer structure lie almost parallel to the clay surface whereas those out of the layer structure are essentially perpendicular to the surface. The trans conformations are predominant in all cases although extensive gauche conformations are observed, which is in agreement with previous simulations on n-butane. Moreover, an odd-even effect in conformation distributions is observed mainly along the chains close to the head and tail groups. The diffusion constants of both nitrogen atoms and methylene groups in these nanoconfined alkyl chains increase with the temperature and methelene position toward the tail groups.

Phase equilibrium studies at moderate pressures

The theory of vapour-liquid equilibria is reviewed, as is the present status or prediction methods in this field. After discussion of the experimental methods available, development of a recirculating equilibrium still based on a previously successful design (the modified Raal, Code and Best still of O'Donnell and Jenkins) is described. This novel still is designed to work at pressures up to 35 bar and for the measurement of both isothermal and isobaric vapour-liquid equilibrium data. The equilibrium still was first commissioned by measuring the saturated vapour pressures of pure ethanol and cyclohexane in the temperature range 77-124°C and 80-142°C respectively. The data obtained were compared with available literature experimental values and with values derived from an extended form of the Antoine equation for which parameters were given in the literature. Commissioning continued with the study of the phase behaviour of mixtures of the two pure components as such mixtures are strongly non-ideal, showing azeotopic behaviour. Existing data did not exist above one atmosphere pressure. Isothermal measurements were made at 83.29°C and 106.54°C, whilst isobaric measurements were made at pressures of 1 bar, 3 bar and 5 bar respectively. The experimental vapour-liquid equilibrium data obtained are assessed by a standard literature method incorporating a themodynamic consistency test that minimises the errors in all the measured variables. This assessment showed that reasonable x-P-T data-sets had been measured, from which y-values could be deduced, but that the experimental y-values indicated the need for improvements in the design of the still. The final discussion sets out the improvements required and outlines how they might be attained.