Length of near-wall plumes in turbulent convection

We present planforms of line plumes formed on horizontal surfaces in turbulent convection, along with the length of line plumes measured from these planforms, in a six decade range of Rayleigh numbers (10(5) < Ra < 10(11)) and at three Prandtl numbers (Pr = 0.7, 5.2, 602). Using geometric constraints on the relations for the mean plume spacings, we obtain expressions for the total length of near-wall plumes on horizontal surfaces in turbulent convection. The plume length per unit area (L(p)/A), made dimensionless by the near-wall length scale in turbulent convection (Z(w)), remains constant for a given fluid. The Nusselt number is shown to be directly proportional to L(p)H/A for a given fluid layer of height H. The increase in Pr has a weak influence in decreasing L(p)/A. These expressions match the measurements, thereby showing that the assumption of laminar natural convection boundary layers in turbulent convection is consistent with the observed total length of line plumes. We then show that similar relationships are obtained based on the assumption that the line plumes are the outcome of the instability of laminar natural convection boundary layers on the horizontal surfaces.

Planform structure of turbulent free convection on horizontal surfaces

The planform structure of turbulent free convection over a heated horizontal surface has been visualized and analyzed for different boundary conditions at the top and for different aspect ratios, for flux Rayleigh numbers ranging from 10 exp 8 - 10 exp 10. The different boundary conditions correspond to Rayleigh-Benard convection, open convection with evaporation at the top and with an imposed external flow on the heated boundary. Without the external flow the planform is one randomly oriented line plume. At large Ra, these line plumes seem to align along the diagonal, persumably due to a large-scale flow along as visualized in the side view. When the external flow is imposed, the line plumes clearly align in the direction of external flow. Flow visualization reveals that at these Ra, the shear tends to break the plumes which otherwise would reach the opposite boundary. (Author)

Cumulus clouds and convective boundary layers: a tropical perspective on two turbulent shear flows

This paper is a review prepared for the second Marseille Colloquium on the mechanics of turbulence, held in 2011, 50 years after the first. The review covers recent developments in our understanding of the large-scale dynamics of cumulus cloud flows and of the atmospheric boundary layer in the low-wind convective regime that is often encountered in the tropics. It has recently been shown that a variety of cumulus cloud forms and life cycles can be experimentally realized in the laboratory, with the transient diabatic plume taken as the flow model for a cumulus cloud. The plume is subjected to diabatic heating scaled to be dynamically similar to heat release from phase changes in clouds. The experiments are complemented by exact numerical solutions of the Navier-Stokes-Boussinesq equations for plumes with scaled off-source heating. The results show that the Taylor entrainment coefficient first increases with heating, reaches a positive maximum and then drops rapidly to zero or even negative values. This reduction in entrainment is a consequence of structural changes in the flow, smoothing out the convoluted boundaries in the non-diabatic plume, including the tongues engulfing the ambient flow. This is accompanied by a greater degree of mixedness in the core flow because of lower dilution by the ambient fluid. The cloud forms generated depend strongly on the history of the diabatic heating profile in the vertical direction. The striking effects of heating on the flow are attributable to the operation of the baroclinic torque due to the temperature field. The mean baroclinic torque is shown to peak around a quasi-cylindrical sheet situated midway between the axis of the flow and the edges. This torque is shear-enhancing and folds down the engulfment tongues. The increase in mixedness can be traced to an explosive growth in the enstrophy, triggered by a strong fluctuating baroclinic torque that acts as a source, especially at the higher wave numbers, thus enhancing the mixedness. In convective boundary layers field measurements show that, under conditions prevailing in the tropics, the eddy fluxes of momentum and energy do not follow the Monin-Obukhov similarity. Instead, the eddy momentum flux is found to be linear in the wind speed at low winds; and the eddy heat flux is, to a first approximation, governed by free convection laws, with wind acting as a small perturbation on a regime of free convection. A new boundary layer code, based on heat flux scaling rather than wall-stress scaling, shows promising improvements in predictive skills of a general circulation model.

Chromite-silicate chemistry of the Neoarchean Sittampundi Complex, southern India: Implications for subduction-related arc magmatism

The Neoarchean layered anorthositic complex at Sittampundi in southern India is known for its chromitite layers that are mostly associated with anorthosite (An(90-100)). The chromitites contain FeAl-rich chromites concentrated in layers between amphibole-rich layers with a dominant mineralogy of amphibole-spinel-plagiocase+/-sapphirine. The chromite-rich layers contain only amphibole and plagioclase. Mineral compositions illustrated by X-ray composition maps and profiles show subtle chemical differences. The chrome spinels are of refractory grade with Cr2O3 and Al2O3 contents varying between 34-40 wt.% and 23-28 wt.%. The chromite compositions are noticeably different from those in layered igneous intrusions of the Bushveld-Stillwater type. The existence of original highly calcic plagioclase, FeAl-rich chromite, and magmatic amphibole is consistent with derivation from a parental magma of hydrous tholeiitic composition that was most likely generated in a supra-subduction zone arc setting. In terms of mineralogy and field relations, the Sittampundi chromitites are remarkably similar to anorthosite-hosted chromitites in the Neoarchean Fiskensset anorthositic complex, Greenland. We propose that the Sittampundi chromitites formed by partial melting of unusually aluminous harzburgite in a hydrated mantle wedge above a subduction zone. This melting process produced hydrous, aluminous basalt, which fractionated at depth to give rise to a variety of high-alumina basalt compositions from which the anorthositic complex with its cumulate chromite-rich and amphibole-rich layers formed within the magma chamber of a supra-subduction zone arc. (C) 2011 Elsevier B.V. All rights reserved.

Review on atomization and sprays of biofuels for IC engine applications

Ever increasing energy requirements, environmental concerns and energy security needs are strongly influencing engine researchers to consider renewable biofuels as alternatives to fossil fuels. Spray process being important in IC engine combustion, existing literature on various biofuel sprays is reviewed and summarized. Both experimental and computational research findings are reviewed in a detailed manner for compression ignition (CI) engine sprays and briefly for spark ignition (SI) engine sprays. The physics of basic atomization process of sprays from various injectors is included to highlight the most recent research findings followed by discussion highlighting the effect of physico-chemical properties on spray atomization for both biofuels and fossil fuels. Biodiesel sprays are found to penetrate faster and haw narrow spray plume angle and larger droplet sizes compared to diesel. Results of analytical and computational models are shown to be useful in shedding light on the actual process of atomization. However, further studies on understanding primary atomization and the effect of fuel properties on primary atomization are required. As far as secondary atomization is concerned, changes in regimes are observed to occur at higher air-jet velocities for biodiesel compared to those of diesel. Evaporating sprays revealed that the liquid length is longer for biodiesel. Pure plant oil sprays with potential use in CI engines may require alternative injector technology due to slower breakup as compared to diesel. Application of ethanol to gasoline engines may be feasible without any modifications to port fuel injection (PFI) engines. More studies are required on the application of alternative fuels to high pressure sprays used in Gasoline Direct Injection (GDI) engines.

Dynamics of the North American continent

The forces that cause deformation of western North America have been debated for decades. Recent studies, primarily based on analysis of crustal stresses in the western United States, have suggested that the deformation of the region is mainly controlled by gravitational potential energy (GPE) variations and boundary loads, with basal tractions due to mantle flow playing a relatively minor role. We address these issues by modelling the deviatoric stress field over western North America from a 3-D finite element mantle circulation model with lateral viscosity variations. Our approach takes into account the contribution from both topography and shallow lithosphere structure (GPE) as well as that from deeper mantle flow in one single model, as opposed to separate lithosphere and circulation models, as has been done so far. In addition to predicting the deviatoric stresses we also jointly fit the constraints of geoid, dynamic topography and plate motion both globally and over North America, in order to ensure that the forces that arise in our models are dynamically consistent. We examine the sensitivity of the dynamic models to different lateral viscosity variations. We find that circulation models that include upper mantle slabs yield a better fit to observed plate velocities. Our results indicate that a model of GPE variations coupled with mantle convection gives the best fit to the observational constraints. We argue that although GPE variations control a large part of the deformation of the western United States, deeper mantle tractions also play a significant role. The average deviatoric stress magnitudes in the western United States range 30-40 MPa. The cratonic region exhibits higher coupling to mantle flow than the rest of the continent. We find that a relatively strong San Andreas fault gives a better fit to the observational constraints, especially that of plate velocity in western North America.

Establishing a Relation between the Mass and the Spin of Stellar-Mass Black Holes

Stellar mass black holes (SMBHs), forming by the core collapse of very massive, rapidly rotating stars, are expected to exhibit a high density accretion disk around them developed from the spinning mantle of the collapsing star. A wide class of such disks, due to their high density and temperature, are effective emitters of neutrinos and hence called neutrino cooled disks. Tracking the physics relating the observed (neutrino) luminosity to the mass, spin of black holes (BHs) and the accretion rate ((M) over dot) of such disks, here we establish a correlation between the spin and mass of SMBHs at their formation stage. Our work shows that spinning BHs are more massive than nonspinning BHs for a given (M) over dot. However, slowly spinning BHs can turn out to be more massive than spinning BHs if (M) over dot at their formation stage was higher compared to faster spinning BHs.

Spatial variations of effective elastic thickness over the Ninetyeast Ridge and implications for its structure and tectonic evolution

We present new data on the strength of oceanic lithosphere along the Ninetyeast Ridge (NER) from two independent methods: spectral analysis (Bouguer coherence) using the fan wavelet transform technique, and spatial analysis (flexure inversion) with the convolution method. The two methods provide effective elastic thickness (T-e) patterns that broadly complement each other, and correlate well with known surface structures and regional-scale features. Furthermore, our study presents a new high resolution database on the Moho configuration, which obeys flexural isostasy, and exhibit regional correlations with the T-e variations. A continuous ridge structure with a much lower T-e value than that of normal oceanic lithosphere provides strong support for the hotspot theory. The derived T-e values vary over the northern (higher T-e similar to 10-20 km), central (anomalously low T-e similar to 0-5 km), and southern (low T-e similar to 5 km) segments of the NER. The lack of correlation of the T-e value with the progressive aging of the lithosphere implies differences in thermo-mechanical setting of the crust and underlying mantle in different parts of the NER, again indicating diversity in their evolution. The anomalously low T-e and deeper Moho (similar to 22 km) estimates of the central NER (between 0.5 degrees N and 17 degrees S) are attributed to the interaction of a hotspot with the Wharton spreading ridge that caused significant thermal rejuvenation and hence weakening of the lithosphere. The higher mechanical strength values in the northern NER (north of 0.5 degrees N) may support the idea of off-ridge emplacement and a relatively large plate motion at the time of volcanism. The low T-e and deeper Moho (similar to 22 km) estimates in the southern part (south of 17 degrees S) suggest that the lithosphere was weak and therefore younger at the time of volcanism, and this supports the idea that the southern NER was emplaced on the edge of the Indian plate. (C) 2013 Elsevier B.V. All rights reserved.

Imprints of Archean to Neoproterozoic crustal processes in the Madurai Block, Southern India

Madurai Block, the largest crustal block in the Southern Granulite Terrane (SGT) of Peninsular India, preserves the imprints of multistage tectonic evolution. Here, we present U-Pb and Hf isotope data on zircons from a charnockite-granite suite in the north-western part of this block. The oscillatory zoning, and the LREE to HREE enriched patterns of the zircons with positive Ce and negative Eu anomalies suggest that the zircon cores are of magmatic origin, with ages in the range of 2634-2435 Ma implying Neoarchean-Paleoproterozoic magmatism followed by subsequent metamorphism and protocontinent formation in the north-western part of the Madurai Block. A regional 550-500 Ma metamorphic overprint is also preserved in the zircons coinciding with the final amalgamation of the Gondwana supercontinent. The Hf isotopic data suggest that the granite and charnockite were derived from isotopically heterogeneous juvenile crustal domains and the charnockites show a significant contribution of mantle-derived components. Therefore, the Hf isotopic data reflect mixing of crustal and mantle-derived sources for the generation of Neoarchean crust in the north-western Madurai Block, possibly in a suprasubduction zone setting during continent building processes. (c) 2014 Elsevier Ltd. All rights reserved.

Current kinematics and dynamics of Africa and the East African Rift System

Although the East African Rift System (EARS) is an archetype continental rift, the forces driving its evolution remain debated. Some contend buoyancy forces arising from gravitational potential energy (GPE) gradients within the lithosphere drive rifting. Others argue for a major role of the diverging mantle flow associated with the African Superplume. Here we quantify the forces driving present-day continental rifting in East Africa by (1) solving the depth averaged 3-D force balance equations for 3-D deviatoric stress associated with GPE, (2) inverting for a stress field boundary condition that we interpret as originating from large-scale mantle tractions, (3) calculating dynamic velocities due to lithospheric buoyancy forces, lateral viscosity variations, and velocity boundary conditions, and (4) calculating dynamic velocities that result from the stress response of horizontal mantle tractions acting on a viscous lithosphere in Africa and surroundings. We find deviatoric stress associated with lithospheric GPE gradients are similar to 8-20 MPa in EARS, and the minimum deviatoric stress resulting from basal shear is similar to 1.6 MPa along the EARS. Our dynamic velocity calculations confirm that a force contribution from GPE gradients alone is sufficient to drive Nubia-Somalia divergence and that additional forcing from horizontal mantle tractions overestimates surface kinematics. Stresses from GPE gradients appear sufficient to sustain present-day rifting in East Africa; however, they are lower than the vertically integrated strength of the lithosphere along most of the EARS. This indicates additional processes are required to initiate rupture of continental lithosphere, but once it is initiated, lithospheric buoyancy forces are enough to maintain rifting.

Silicon isotopes: from cosmos to benthos

Silicon is the second most abundant element on the Earth and one of the more abundant elements in our Solar System. Variations in the relative abundance of the stable isotopes of Si (Si isotope fractionation) in different natural reservoirs, both terrestrial (surface and deep Earth) as well as extra-terrestrial (e.g. meteorites, lunar samples), are a powerful tracer of present and past processes involving abiotic as well as biotic systems. The versatility of the Si isotope tracer is reflected in its wide-ranging applications from understanding the origin of early Solar System objects, planetary differentiation, Moon formation, mantle melting and magma differentiation on the Earth, ancient sea-water composition, to modern-day weathering, clay formation and biological fractionation on land as well as in the oceans. The application of Si isotopes as tracers of natural processes started over six decades ago and its usage has seen a sudden increase over the last decade due to improvements in mass spectrometry, particularly the advent of multi-collector inductively coupled plasma mass spectrometers, which has made Si isotope measurements safe and relatively easy while simultaneously improving the accuracy and precision of measurements.

India-Madagascar paleo-fit based on flexural isostasy of their rifted margins

The present study contributes new constraints on, and definitions of, the reconstructed plate margins of India and Madagascar based on flexural isostasy along the Western Continental Margin of India (WCMI) and the Eastern Continental Margin of Madagascar (ECMM). We have estimated the nature of isostasy and crustal geometry along the two margins, and have examined their possible conjugate structure. Here we utilize elastic thickness (Te) and Moho depth data as the primary basis for the correlation of these passive margins. We employ the flexure inversion technique that operates in spatial domain in order to estimate the spatial variation of effective elastic thickness. Gravity inversion and flexure inversion techniques are used to estimate the configuration of the Moho/Crust-Mantle Interface that reveals regional correlations with the elastic thickness variations. These results correlate well with the continental and oceanic segments of the Indian and African plates. The present study has found a linear zone of anomalously low-Te (1-5 km) along the WCMI (similar to 1680 km), which correlates well with the low-Te patterns obtained all along the ECMM. We suggest that the low-Te zones along the WCMI and ECMM represent paleo-rift inception points of lithosphere thermally and mechanically weakened by the combined effects of the Marion hotspot and lithospheric extension due to rifting. We have produced an India-Madagascar paleo-fit representing the initial phase of separation based on the Te estimates of the rifted conjugate margins, which is confirmed by a close-fit correlation of Moho geometry and bathymetry of the shelf margins. The matching of tectonic lineaments, lithologies and geochronological belts between India and Madagascar provide an additional support for the present plate reconstruction. (C) 2014 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

A Short Review on Lu-Hf Isotope System in Zircon: Implications for Crustal Evolution

Zircon has been recognized as the unaltered part of the Earth's history which preserves nearly 4 billion year record of earth's evolution. Zircon preserves igneous and metamorphic processes during its formation and remains unaffected by sedimentary processes and crustal recycling. U-Pb and Lu-Hf in zircon work as geochronometer and geochemical tracer respectively. Zircon provide valuable information about the source composition of the rocks and the intrinsic details of an unseen crust-mantle processes. The world wide data of U-Pb and Lu-Hf isotope systems in zircon reveal crustal evolution through geological history. Moreover, the U-Pb age pattern of zircons show distinct peaks attributed to preservation of crustal rocks or mountain building during supercontinent assembly. The histogram of continental crust preservation shows that nearly one-third of continental crust was formed during the Archean, almost 20% was formed during Paleoproterozoic and 14% in last 400 Ma.

Dynamos driven by weak thermal convection and heterogeneous outer boundary heat flux

We use numerical dynamo models with heterogeneous core-mantle boundary (CMB) heat flux to show that lower mantle lateral thermal variability may help support a dynamo under weak thermal convection. In our reference models with homogeneous CMB heat flux, convection is either marginally supercritical or absent, always below the threshold for dynamo onset. We find that lateral CMB heat flux variations organize the flow in the core into patterns that favour the growth of an early magnetic field. Heat flux patterns symmetric about the equator produce non-reversing magnetic fields, whereas anti-symmetric patterns produce polarity reversals. Our results may explain the existence of the geodynamo prior to inner core nucleation under a tight energy budget. Furthermore, in order to sustain a strong geomagnetic field, the lower mantle thermal distribution was likely dominantly symmetric about the equator. (C) 2015 Elsevier B.V. All rights reserved.

稳态等离子体推进器羽流的粒子模拟

本文利用结合了particle-in-cell方法的直接模拟Monte Carlo方法模拟了稳态等离子体推进器的羽流场.为了更真实地再现实验设备中的流动,模拟中包含了背压气体,并在计算中直接模拟,改变了以往的中性粒子平衡态的处理方式.为了匹配试验结果,修正了推进器出口处气流偏转角分布,并建议最大偏转半角取为20°.计算结果与现有的实验结果符合较好.