Wall-mode instability in plane shear flow of viscoelastic fluid over a deformable solid

The linear stability analysis of a plane Couette flow of an Oldroyd-B viscoelastic fluid past a flexible solid medium is carried out to investigate the role of polymer addition in the stability behavior. The system consists of a viscoelastic fluid layer of thickness R, density rho, viscosity eta, relaxation time lambda, and retardation time beta lambda flowing past a linear elastic solid medium of thickness HR, density rho, and shear modulus G. The emphasis is on the high-Reynolds-number wall-mode instability, which has recently been shown in experiments to destabilize the laminar flow of Newtonian fluids in soft-walled tubes and channels at a significantly lower Reynolds number than that for flows in rigid conduits. For Newtonian fluids, the linear stability studies have shown that the wall modes become unstable when flow Reynolds number exceeds a certain critical value Re c which scales as Sigma(3/4), where Reynolds number Re = rho VR/eta, V is the top-plate velocity, and dimensionless parameter Sigma = rho GR(2)/eta(2) characterizes the fluid-solid system. For high-Reynolds-number flow, the addition of polymer tends to decrease the critical Reynolds number in comparison to that for the Newtonian fluid, indicating a destabilizing role for fluid viscoelasticity. Numerical calculations show that the critical Reynolds number could be decreased by up to a factor of 10 by the addition of small amount of polymer. The critical Reynolds number follows the same scaling Re-c similar to Sigma(3/4) as the wall modes for a Newtonian fluid for very high Reynolds number. However, for moderate Reynolds number, there exists a narrow region in beta-H parametric space, corresponding to very dilute polymer solution (0.9 less than or similar to beta < 1) and thin solids (H less than or similar to 1.1), in which the addition of polymer tends to increase the critical Reynolds number in comparison to the Newtonian fluid. Thus, Reynolds number and polymer properties can be tailored to either increase or decrease the critical Reynolds number for unstable modes, thus providing an additional degree of control over the laminar-turbulent transition.

Low frequency modulation of jets in quasigeostrophic turbulence

Quasigeostrophic turbulence on a beta-plane with a finite deformation radius is studied numerically, with particular emphasis on frequency and combined wavenumber-frequency domain analyses. Under suitable conditions, simulations with small-scale random forcing and large-scale drag exhibit a spontaneous formation of multiple zonal jets. The first hint of wave-like features is seen in the distribution of kinetic energy as a function of frequency; specifically, for progressively larger deformation scales, there are systematic departures in the form of isolated peaks (at progressively higher frequencies) from a power-law scaling. Concomitantly, there is an inverse flux of kinetic energy in frequency space which extends to lower frequencies for smaller deformation scales. The identification of these peaks as Rossby waves is made possible by examining the energy spectrum in frequency-zonal wavenumber and frequency-meridional wavenumber diagrams. In fact, the modified Rhines scale turns out to be a useful measure of the dominant meridional wavenumber of the modulating Rossby waves; once this is fixed, apart from a spectral peak at the origin (the steady jet), almost all the energy is contained in westward propagating disturbances that follow the theoretical Rossby dispersion relation. Quite consistently, noting that the zonal scale of the modulating waves is restricted to the first few wavenumbers, the energy spectrum is almost entirely contained within the corresponding Rossby dispersion curves on a frequency-meridional wavenumber diagram. Cases when jets do not form are also considered; once again, there is a hint of Rossby wave activity, though the spectral peaks are quite muted. Further, the kinetic energy scaling in frequency domain follows a -5/3 power-law and is distributed much more broadly in frequency-wavenumber diagrams. (C) 2015 AIP Publishing LLC.

Cryogenian magmatism and crustal reworking in the Southern Granulite Terrane, India

Understanding Neoproterozoic crustal evolution is fundamental to reconstructing the Gondwana supercontinent, which was assembled at this time. Here we report evidence of Cryogenian crustal reworking in the Madurai Block of the Southern Granulite Terrane of India. The study focuses on a garnet-bearing granite-charnockite suite, where the granite shows in situ dehydration into patches and veins of incipient charnockite along the contact with charnockite. The granite also carries dismembered layers of Mg-Al-rich granulite. Micro-textural evidence for dehydration of granite in the presence of CO2-rich fluids includes the formation of orthopyroxene by the breakdown of biotite, neoblastic zircon growth in the dehydration zone, at around 870 degrees C and 8kbar. The zircon U-Pb ages suggest formation of the granite, charnockite, and incipient charnockite at 836 +/- 73, 831 +/- 31, and 772 +/- 49Ma, respectively. Negative zircon epsilon Hf (t) (-5 to -20) values suggest that these rocks were derived from a reworked Palaeoproterozoic crustal source. Zircon grains in the Mg-Al-rich granulite record a spectrum of ages from ca. 2300 to ca. 500Ma, suggesting multiple provenances ranging from Palaeoproterozoic to mid-Neoproterozoic, with neoblastic zircon growth during high-temperature metamorphism in the Cambrian. We propose that the garnet-bearing granite and charnockite reflect the crustal reworking of aluminous crustal material indicated by the presence of biotite+quartz+aluminosilicate inclusions in the garnet within the granite. This crustal source can be the Mg-Al-rich layers carried by the granite itself, which later experienced high-temperature regional metamorphism at ca. 550Ma. Our model also envisages that the CO2 which dehydrated the garnet-bearing granite generating incipient charnockite was sourced from the proximal massive charnockite through advection. These Cryogenian crustal reworking events are related to prolonged tectonic activities prior to the final assembly of the Gondwana supercontinent.

A trade-off between maxima in efficiency and specific work output of super- and trans-critical CO2 Brayton cycles

Several operational aspects for thermal power plants in general are non-intuitive and involve simultaneous optimization of a number of operational parameters. In the case of solar operated power plants, it is even more difficult due to varying heat source temperatures induced by variability in insolation levels. This paper introduces a quantitative methodology for load regulation of a CO2 based Brayton cycle power plant using the `thermal efficiency and specific work output' coordinate system. The analysis shows that a transcritical CO2 cycle offers more flexibility under part load performance than the supercritical cycle in case of non-solar power plants. However, for concentrated solar power, where efficiency is important, supercritical CO2 cycle fares better than transcritical CO2 cycle. A number of empirical equations relating heat source temperature, high side pressure with efficiency and specific work output are proposed which could assist in generating control algorithms. (C) 2015 Elsevier B.V. All rights reserved.

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.

Acoustic response of vortex breakdown modes in a coaxial isothermal unconfined swirling jet

The present experimental work is concerned with the study of amplitude dependent acoustic response of an isothermal coaxial swirling jet. The excitation amplitude is increased in five distinct steps at the burner's Helmholtz resonator mode (i.e., 100 Hz). Two flow states are compared, namely, sub-critical and super-critical vortex breakdown (VB) that occur before and after the critical conical sheet breakdown, respectively. The geometric swirl number is varied in the range 2.14-4.03. Under the influence of external pulsing, global response characteristics are studied based on the topological changes observed in time-averaged 2D flow field. These are obtained from high resolution 2D PIV (particle image velocimetry) in the longitudinal-mid plane. PIV results also illustrate the changes in the normalized vortex core coordinates (r(vcc)/(r(vcc))(0) (Hz), y(vcc)/(y(vcc))(0) (Hz)) of internal recirculation zone (IRZ). A strong forced response is observed at 100 Hz (excitation frequency) in the convectively unstable region which get amplified based on the magnitude of external forcing. The radial extent of this forced response region at a given excitation amplitude is represented by the acoustic response region (b). The topological placement of the responsive convectively unstable region is a function of both the intensity of imparted swirl (characterized by swirl number) and forcing amplitude. It is observed that for sub-critical VB mode, an increase in the excitation amplitude till a critical value shifts the vortex core centre (particularly, the vortex core moves downstream and radially outwards) leading to drastic fanning-out/widening of the IRZ. This is accompanied by similar to 30% reduction in the recirculation velocity of the IRZ. It is also observed that b < R (R: radial distance from central axis to outer shear layer-OSL). At super-critical amplitudes, the sub-critical IRZ topology transits back (the vortex core retracts upstream and radially inwards) and finally undergoes a transverse shrinkage ((r(vcc))/(r(vcc))(0 Hz) decreases by similar to 20%) when b >= R. In contrast, the vortex core of super-critical breakdown mode consistently spreads radially outwards and is displaced further downstream. Finally, the IRZ fans-out at the threshold excitation amplitude. However, the acoustic response region b is still less than R. This is explained based on the characteristic geometric swirl number (S-G) of the flow regimes. The super-critical flow mode with higher S-G (hence, higher radial pressure drop due to rotational effect which scales as Delta P similar to rho u theta(2) and acts inwards towards the center line) compared to sub-critical state imposes a greater resistance to the radial outward spread of b. As a result, the acoustic energy supplied to the super-critical flow mode increases the degree of acoustic response at the pulsing frequency and energizes its harmonics (evident from power spectra). As a disturbance amplifier, the stronger convective instability mode within the flow structure of super-critical VB causes the topology to widen/fan-out severely at threshold excitation amplitude. (C) 2015 AIP Publishing LLC.

Experimental determination and theoretical correlation for the solubilities of dicarboxylic acid esters in supercritical carbon dioxide

The synthesis of high molecular weight esters such as bis (2-ethylhexyl) sebacate is of significance for its use as a lubricant. This ester is synthesized by the transesterification of dimethyl sebacate with 2-ethylhexanol. Therefore, the solubilities of bis (2-ethylhexyl) sebacate and dimethyl sebacate were determined at 308-328 K at pressures of 10-18 MPa in supercritical carbon dioxide. The solubility of dimethyl sebacate was always higher than bis (2-ethylhexyl) sebacate at a given temperature and pressure. The Mendez-Teja model was used to verify the self-consistency of data. Further, a new semi-empirical model with three parameters was developed using the solution theory coupled with Wilson activity coefficient. This model was used to correlate the experimental data of this work and solubilities of many high molecular weight esters reported in the literature. (C) 2015 Elsevier B.V. All rights reserved.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Quantitative metabolic profiling of NMR spectral signatures of branched chain amino acids in blood serum

Branched Chain Amino Acids (BCAAs) are related to different aspects of diseases like pathogenesis, diagnosis and even prognosis. While in some diseases, levels of all the BCAAs are perturbed; in some cases, perturbation occurs in one or two while the rest remain unaltered. In case of ischemic heart disease, there is an enhanced level of plasma leucine and isoleucine but valine level remains unaltered. In `Hypervalinemia', valine is elevated in serum and urine, but not leucine and isoleucine. Therefore, identification of these metabolites and profiling of individual BCAA in a quantitative manner in body-fluid like blood plasma/serum have long been in demand. H-1 NMR resonances of the BCAAs overlap with each other which complicates quantification of individual BCAAs. Further, the situation is limited by the overlap of broad resonances of lipoprotein with the resonances of BCAAs. The widely used commercially available kits cannot differentially estimate the BCAAs. Here, we have achieved proper identification and characterization of these BCAAs in serum in a quantitative manner employing a Nuclear Magnetic Resonance-based technique namely T-2-edited Correlation Spectroscopy (COSY). This approach can easily be extended to other body fluids like bile, follicular fluids, saliva, etc.

PREDICTION OF TRANSITIONAL AND SEPARATED BOUNDARY LAYERS IN A COMPRESSOR CASCADE

Numerical simulations were performed of experiments from a cascade of stator blades at three low Reynolds numbers representative of flight conditions. Solutions were assessed by comparing blade surface pressures, velocity and turbulence intensity along blade normals at several stations along the suction surface and in the wake. At Re = 210,000 and 380,000 the laminar boundary layer over the suction surface separates and reattaches with significant turbulence fluctuations. A new 3-equation transition model, the k-k(L)-omega model, was used to simulate this flow. Predicted locations of the separation bubble, and profiles of velocity and turbulence fluctuations on blade-normal lines at various stations along the blade were found to be quite close to measurements. Suction surface pressure distributions were not as close at the lower Re. The solution with the standard k-omega SST model showed significant differences in all quantities. At Re = 640,000 transition occurs earlier and it is a turbulent boundary layer that separates near the trailing edge. The solution with the Reynolds stress model was found to be quite close to the experiment in the separated region also, unlike the k-omega SST solution. Three-dimensional computations were performed at Re = 380,000 and 640,000. In both cases there were no significant differences between the midspan solution from 3D computations and the 2D solutions. However, the 3D solutions exhibited flow features observed in the experiments the nearly 2D structure of the flow over most of the span at 380,000 and the spanwise growth of corner vortices from the endwall at 640,000.

ELEVATION OF HEAT REJECTION TEMPERATURE IN TRANSCRITICAL CONDENSING CYCLES USING CO2+PROPANE MIXTURES

In this paper, a detailed thermodynamic performance analysis of a transcritical condensing (TC) cycle is performed with pure CO2 and a blend of 48.5 % propane with 51.5 % CO2 as working fluids. A realistic thermodynamic model is used incorporating irreversibilities in turbo-machineries and heat exchangers. The Key finding is that the addition of propane elevates the heat rejection temperature, but does not impair any of the performance indicators. Such a fluid may be useful for power generation in concentrated solar power applications by using which a hike of up to 2 % can be realized in the thermal efficiency of a power plant.

Local structure of boundary layer transition in experiments with a single streamwise vortex

Transition induced by an isolated streamwise vortex embedded in a flat plate boundary layer was studied experimentally. The vortex was created by a gentle hill with a Gaussian profile that spanned on half of the width of a flat plate mounted in a low turbulence wind tunnel. PIV and hot-wire anemometry data were taken. Transition occurs as a non-inclined shear layer breaks up into a sequence of vortices, close to the boundary layer edge. The passing frequency of these vortices scales with square of the freestream velocity, similar to that in single-roughness induced transition. Quadrant analysis of streamwise and wall-normal velocity fluctuations show large ejection events in the outer layer. (C) 2015 Elsevier Inc. All rights reserved.

Intertwined vorticity and elastodynamics in flapping wing propulsion

We performed numerical experiments on a one-dimensional elastic solid oscillating in a two-dimensional viscous incompressible fluid with the intent of discerning the interplay of vorticity and elastodynamics in flapping wing propulsion. Perhaps for the first time, we have established the role of foil deflection topology and its influence on vorticity generation, through spatially and temporally evolving foil slope and curvature. Though the frequency of oscillation of the foil has a definite role, it is the phase relation between foil slope and pressure that determines thrust or drag. Similarly, the phase difference between flapping velocity, and pressure and inertial forces, determine the power input to the foil, and in turn drives propulsive efficiency. At low frequencies of oscillation, the sympathetic slope and curvature of deformation of the foil allow generation of leading-edge vortices that do not separate; they cause substantial rise in pressure between the leading edge and mid-chord. The circulatory component of pressure is determined primarily by the leading-edge vortex and therefore thrust too is predominantly circulatory in origin at low frequencies. In the intermediate and high-frequency range, thrust and drag on the foil spatially alternate and non-circulatory forces dominate over circulatory and viscous forces. For the mass ratios we simulated, thrust due to flapping varies quadratically as a function of Strouhal number or trailing-edge flapping velocity; further, the trailing edge flapping velocities peak at the same set of frequencies where the thrust is also a maximum. Propulsive efficiency, on the other hand, is roughly a mirror image of the thrust variation with respect to Strouhal number. Given that most instances of flapping propulsion in nature are primarily through distributed muscular actuation that enables precise control of deformation shape, leading to high thrust and efficiency, the results presented here are pointers towards understanding some of the mechanisms that drive thrust and propulsive efficiency.

Solubilities of 10-undecenoic acid and geraniol in supercritical carbon dioxide

The solubilities of 10-undecenoic acid and geraniol in supercritical carbon dioxide were measured at 308, 313,323 and 333K, and at pressures of 10-18 MPa. Solubilities (in mole fraction),ranged from 0.4 x 10(-3) to 17.4 x 10(-3) for 10-undecenoic acid and 2.7 x 10(-3) to 25 x 10(-3) for geraniol, respectively. The AARD was around 11% and 5% for these models for 10-undecenoic acid and geraniol, respectively. The solubilities of both compounds showed retrograde behavior wherein the solubilities decrease with temperature at isobaric conditions. The solubility of geraniol was higher than 10-undecenoic acid at all investigated temperatures and pressures. The data were found to be self consistent based on the Mendez-Santiago model. New models based on association theory using van Laar and Margules activity coefficient models for solute in liquid phase were derived, and used to correlate the solubilities. (C) 2015 Elsevier B.V. All rights reserved.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.