Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

The behaviour of laminar stratified methane/air flames in counterflow

Flames propagating through a mixture with a gradient of equivalence ratio have been previously demonstrated to travel faster or slower than their equivalent premixed flames. The present study aims to numerically investigate the response of strained laminar methane-air flames to such gradients. The flames are simulated in a counterflow configuration where a premixed reactant stream at equivalence ratio φR opposes a hot equilibrium stream at equivalence ratio φP. Premixed and stratified flames are compared with respect to the equivalence ratio φ* and the corresponding gradient ∇φ* at the point of peak heat release rate, for three strain rates, a=50, 300 and 500s-1 and a range of φ*. The effect of different stratification levels is also investigated by varying the ratio of φP to φR, Θ. Results indicate that, as long as flames stabilize within the diffusion layer and Θ>1, increased heat release rate Q is seen throughout the progress variable space in comparison to the premixed state. In contrast, an attenuation of heat release rate is seen for Θ<1. The enhancement (or attenuation) of heat release varies monotonically with Θ. The effect of stratification on flame behavior becomes more pronounced as the strain rate increases. The present study reveals the mechanisms for the propagation of quasi-steady stratified flames under lean and rich conditions: stratified flames are primarily dominated by the diffusion of heat under lean conditions, and diffusion of H2 under rich conditions. Thanks to species and thermal support, stratified flames continue to burn beyond the premixed lean and rich flammability limits. Further investigation on the unsteady response of flames to the fluctuating equivalence ratio implies that the steady results represent the unsteady response well, as long as φ* and ∇φ* are similar in both steady and unsteady cases. © 2013 The Combustion Institute.

Microheterogeneous Thoria-Urania Fuels for Pressurized Water Reactors

A thorium-based fuel cycle for light water reactors will reduce the plutonium generation rate and enhance the proliferation resistance of the spent fuel. However, priming the thorium cycle with ²³⁵U is necessary, and the ²³⁵U fraction in the uranium must be limited to below 20% to minimize proliferation concerns. Thus, a once-through thorium-uranium dioxide (ThO2-UO2) fuel cycle of no less than 25% uranium becomes necessary for normal pressurized water reactor (PWR) operating cycle lengths. Spatial separation of the uranium and thorium parts of the fuel can improve the achievable burnup of the thorium-uranium fuel designs through more effective breeding of ²³³U from the ²³²Th. Focus is on microheterogeneous fuel designs for PWRs, where the spatial separation of the uranium and thorium is on the order of a few millimetres to a few centimetres, including duplex pellet, axially microheterogeneous fuel, and a checkerboard of uranium and thorium pins. A special effort was made to understand the underlying reactor physics mechanisms responsible for enhancing the achievable burnup at spatial separation of the two fuels. The neutron spectral shift was identified as the primary reason for the enhancement of burnup capabilities. Mutual resonance shielding of uranium and thorium is also a factor; however, it is small in magnitude. It is shown that the microheterogeneous fuel can achieve higher burnups, by up to 15%, than the reference all-uranium fuel. However, denaturing of the ²³³U in the thorium portion of the fuel with small amounts of uranium significantly impairs this enhancement. The denaturing is also necessary to meet conventional PWR thermal limits by improving the power share of the thorium region at the beginning of fuel irradiation. Meeting thermal-hydraulic design requirements by some of the microheterogeneous fuels while still meeting or exceeding the burnup of the all-uranium case is shown to be potentially feasible. However, the large power imbalance between the uranium and thorium regions creates several design challenges, such as higher fission gas release and cladding temperature gradients. A reduction of plutonium generation by a factor of 3 in comparison with all-uranium PWR fuel using the same initial ²³⁵U content was estimated. In contrast to homogeneously mixed U-Th fuel, microheterogeneous fuel has a potential for economic performance comparable to the all-UO2 fuel provided that the microheterogeneous fuel incremental manufacturing costs are negligibly small.

Vortical structures on three-dimensional shock control bumps

The vortical wake structure produced by a three-dimensional shock control bump (SCB) is thought to be useful for controlling transonic buffet on airfoils. However, at present the vorticity produced is relatively weak and the production mechanism is not well understood. Using a combined experimental and computational approach, a preliminary investigation on the wake vorticity for different bump geometries has been carried out. The structure of the wake for on and off-design conditions are considered, and the effects on the downstream boundary layer demonstrated. Three main vortical structures are observed: a primary vortex pair, weak inter-bump vortices and shear flow in the lambda-shock region. The effect of pressure gradients on vortex strength is examined and it is found that spanwise pressure gradients on the front section of the bump are the most significant parameter influencing vortex strength. © 2013 by S.P. Colliss et al.

Linear stability of miscible two-fluid flow down an incline

We show that miscible two-layer free-surface flows of varying viscosity down an inclined substrate are different in their stability characteristics from both immiscible two-layer flows, and flows with viscosity gradients spanning the entire flow. New instability modes arise when the critical layer of the viscosity transport equation overlaps the viscosity gradient. A lubricating configuration with a less viscous wall layer is identified to be the most stabilizing at moderate miscibility (moderate Peclet numbers). This also is in contrast with the immiscible case, where the lubrication configuration is always destabilizing. The co-existence that we find under certain circumstances, of several growing overlap modes, the usual surface mode, and a Tollmien-Schlichting mode, presents interesting new possibilities for nonlinear breakdown. © 2013 AIP Publishing LLC.

Computational analysis of liquid crystalline elastomer membranes: Changing Gaussian curvature without stretch energy

Liquid crystalline elastomers (LCEs) can undergo extremely large reversible shape changes when exposed to external stimuli, such as mechanical deformations, heating or illumination. The deformation of LCEs result from a combination of directional reorientation of the nematic director and entropic elasticity. In this paper, we study the energetics of initially flat, thin LCE membranes by stress driven reorientation of the nematic director. The energy functional used in the variational formulation includes contributions depending on the deformation gradient and the second gradient of the deformation. The deformation gradient models the in-plane stretching of the membrane. The second gradient regularises the non-convex membrane energy functional so that infinitely fine in-plane microstructures and infinitely fine out-of-plane membrane wrinkling are penalised. For a specific example, our computational results show that a non-developable surface can be generated from an initially flat sheet at cost of only energy terms resulting from the second gradients. That is, Gaussian curvature can be generated in LCE membranes without the cost of stretch energy in contrast to conventional materials. © 2013 Elsevier Ltd. All rights reserved.

Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables

We demonstrate how a prior assumption of smoothness can be used to enhance the reconstruction of free energy profiles from multiple umbrella sampling simulations using the Bayesian Gaussian process regression approach. The method we derive allows the concurrent use of histograms and free energy gradients and can easily be extended to include further data. In Part I we review the necessary theory and test the method for one collective variable. We demonstrate improved performance with respect to the weighted histogram analysis method and obtain meaningful error bars without any significant additional computation. In Part II we consider the case of multiple collective variables and compare to a reconstruction using least squares fitting of radial basis functions. We find substantial improvements in the regimes of spatially sparse data or short sampling trajectories. A software implementation is made available on www.libatoms.org.

Centrifuge modeling of seismic liquefaction effects on adjacent shallow foundations

Shallow foundations built on saturated deposits of granular soils in seismically active areas are, regardless of their static bearing capacity, critical structures during seismic events. A single centrifuge experiment involving shallow foundations situated atop a liquefiable soil deposit has been performed to identify the mechanisms involved in the interaction between liquefaction-induced effects on neighboring shallow foundations. Centrifuge test results indicate that liquefaction causes significant settlements of footings, which are affected by the presence of neighboring foundations and can be extremely damaging to the superstructure. The understanding of these interaction effects is very important, mainly in densely populated urban areas. The development of high excess pore-pressures, localized drainage in response to the high transient hydraulic gradients, and earthquake-induced vertical motions to the footings are also important effects that are discussed to assist in enhancing current understanding and ability to predict liquefaction effects on shallow foundations. © 2014 Taylor & Francis Group.

Laminar propagation of lean premixed flames ignited in stratified mixture

Combustion in stratified mixtures is envisaged in practical energy systems such as direct-injection spark-ignited (DISI) car engines, gas turbines, for reducing CO2 and pollutant emissions while protecting their efficiency. The mixture gradients change the fundamental properties of the flame, especially by a difference in temperature and composition between the burnt gases and those of a flame consuming a homogeneous mixture. This paper presents an investigation of the properties of the flame propagating in a lean homogeneous mixture after ignition in a richer mixture according to the magnitude of the stratification. Three magnitudes of stratification are investigated. The local flame burning velocity is determined by an original PIV algorithm developed previously. The local equivalence ratio in the fresh gases is measured from anisole PLIF. From the simultaneous PIV-PLIF measurements, the flame burning velocities conditioned on the local stretch rate and equivalence ratio in fresh gases are measured. The flame propagating through the homogeneous lean mixture has properties depending on the ignition conditions in the stratified layer. The flame propagating in the lean mixture is back-supported longer for ignition under the richer condition. The change of stretch sensitivity and burning velocity of the flame in the lean mixture is measured over time for the three magnitudes of mixture stratification investigated. The ignition in richer mixtures compensates for the nonequidiffusion effect of lean propane flame and sustains its robustness to stretch. The flame propagation in the lean homogeneous mixture is enhanced by ignition in a richer stratified layer, as much by their robustness to stretch as by an increase in the flame speed or the burning velocity. The decay time of this influence of the stratification, called memory effect, is determined. © 2013 The Combustion Institute.

Factors affecting leaf morphology: a case study of Ranunculus natans C. A. Mey. (Ranunculaceae) in the arid zone of northwest China

In order to examine the role of environmental factors affecting foliar morphology, we performed a case study of leaf morphological variation of Ranunculus natans found in the arid zone of northwest China. We found that foliar phenotypic variation differed significantly between populations. We described substantial positive correlations between altitude and leaf area (LA) as well as leaf perimeter (LP), and also between longitude and number of teeth, along with dissection index (DI). The pH, conductivity, and salinity of the environment caused a significant decrease in both LA and LP. Ranked in terms of their impacts on leaf morphology, the six selected factors were: altitude > pH > conductivity > salinity > longitude > latitude. We found that foliar morphological variations are functional responses to water-quantity factors (e.g., altitude and longitude at regional scales) and water-availability relation factors (e.g., pH, conductivity, and salinity at local scales), rather than to temperature-relation factors (latitude). Therefore, altitude and longitude, along with pH, conductivity, and salinity, are the main factors that significantly influence foliar morphology in the arid zone of China. We found that main factors played major roles in plant phenotypic plasticity in a complex ecosystem, although different combinations and interactions of environmental and geographical factors in each local environment may obscure the general trends in trait changes along environmental gradients.

Electrostatic force microscopy and potentiometry of realistic nanostructured systems

We investigate the dependency of electrostatic interaction forces on applied potentials in electrostatic force microscopy (EFM) as well as in related local potentiometry techniques such as Kelvin probe microscopy (KPM). The approximated expression of electrostatic interaction between two conductors, usually employed in EFM and KPM, may loose its validity when probe-sample distance is not very small, as often realized when realistic nanostructured systems with complex topography are investigated. In such conditions, electrostatic interaction does not depend solely on the potential difference between probe and sample, but instead it may depend on the bias applied to each conductor. For instance, electrostatic force can change from repulsive to attractive for certain ranges of applied potentials and probe-sample distances, and this fact cannot be accounted for by approximated models. We propose a general capacitance model, even applicable to more than two conductors, considering values of potentials applied to each of the conductors to determine the resulting forces and force gradients, being able to account for the above phenomenon as well as to describe interactions at larger distances. Results from numerical simulations and experiments on metal stripe electrodes and semiconductor nanowires supporting such scenario in typical regimes of EFM investigations are presented, evidencing the importance of a more rigorous modeling for EFM data interpretation. Furthermore, physical meaning of Kelvin potential as used in KPM applications can also be clarified by means of the reported formalism. © 2009 American Institute of Physics.

Spatial distribution of benthic algae in the Gangqu River, Shangrila, China

This study consisted of sampling benthic algae at 32 sites in the Gangqu River, an important upstream tributary of the Yangtze River. Our aims were to characterize the benthic algae communities and relationships with environmental variables. Among the 162 taxa observed, Achnanthes linearis and Achnanthes lanceolata var. elliptica were the dominant species (17.10% and 14.30% of the total relative abundance, respectively). Major gradients and principal patterns of variation within the environmental variables were detected by principal component analysis (PCA). Then non-metric multidimensional scaling (NMS) divided all the sites into three groups, which were validated by multi-response permutation procedures (MRPP). Canonical correspondence analysis (CCA) indicated that three environmental variables (TN, TDS, and TP) significantly affected the distribution of benthic algae. Weighted averaging regression and cross-calibration produced strong models for predicting TN and TDS concentration, which enabled selection of algae taxa as potentially sensitive indicators of certain TN and TDS levels: for TN, Achnanthes lanceolata, Achnanthes lanceolata var. elliptica, and Cymbella ventricosa var. semicircularis; for TDS, Cocconeis placentula, Cymbella alpina var. minuta, and Fragilaria virescens. The present study represents an early step in establishing baseline conditions. Further monitoring is suggested to gain a better understanding of this region.

Influence of Potamogeton crispus growth on nutrients in the sediment and water of Lake Tangxunhu

An incubation experiment was performed on Potamogeton crispus (P. crispus) using sediment collected from Lake Tangxunhu in the center of China, in order to determine the effects of plant growth on Fe, Si, Cu, Zn, Mn, Mg, P, and Ca concentrations in the sediments and overlying waters. After 3 months of incubation, Ca, Mg, and Si concentrations in the water column were significantly lower, and P and Cu concentrations were significantly higher than in unplanted controls. The effect of P. crispus growth on sediment pore waters and water-extractable elements varied. Concentrations of Ca, Mg, Si, Fe, Cu, and Zn were significantly higher, and P was significantly lower, than in pore waters of the control. Water-extracted concentrations of Fe, Mg, and Si in the sediments were lower, and P was higher, than in the control. Presence of P. crispus generally enhanced concentration gradients of elements between pore waters and overlying waters but not for P. The growth of P. crispus was associated with an increase in water pH and formation of root plaques, resulting in complex effects on the sediment nutritional status.

Spatiotemporal heterogeneity of plankton communities in Lake Donghu, China, as revealed by PCR-denaturing gradient gel electrophoresis and its relation to biotic and abiotic factors

The 16S and 18S rRNA genes of planktonic organisms derived from five stations with nutrient gradients in Lake Donghu, China, were studied by PCR-denaturing gradient gel electrophoresis (DGGE) fingerprinting, and the relationships between the genetic diversity of the plankton community and biotic/abiotic factors are discussed. The concentrations of total nitrogen (TN), total phosphorus (TP), NH4-N and As were found to be significantly related (P < 0.05) to morphological composition of the plankton community. Both chemical and morphological analyses suggested that temporal heterogeneity was comparatively higher than spatial heterogeneity in Lake Donghu. Although the morphological composition was not identical to the DGGE fingerprints in characterizing habitat similarity, the two strongest eutrophic stations (I and II) were always initially grouped into one cluster. Canonical correspondence analysis suggested that the factors strongly correlated with the first two ordination axes were seasonally different. The concentrations of TN and TP and the densities of rotifers and crustaceans were generally the main factors related to the DGGE patterns of the plankton communities. The study suggested that genetic diversity as depicted by metagenomic techniques (such as PCR-DGGE fingerprinting) is a promising tool for ecological study of plankton communities and that such techniques are likely to play an increasingly important role in assessing the environmental conditions of aquatic habitats.

Sparse recovery of complex phase-encoded velocity images using iterative thresholding

In this paper we propose a new algorithm for reconstructing phase-encoded velocity images of catalytic reactors from undersampled NMR acquisitions. Previous work on this application has employed total variation and nonlinear conjugate gradients which, although promising, yields unsatisfactory, unphysical visual results. Our approach leverages prior knowledge about the piecewise-smoothness of the phase map and physical constraints imposed by the system under study. We show how iteratively regularizing the real and imaginary parts of the acquired complex image separately in a shift-invariant wavelet domain works to produce a piecewise-smooth velocity map, in general. Using appropriately defined metrics we demonstrate higher fidelity to the ground truth and physical system constraints than previous methods for this specific application. © 2013 IEEE.