R-parameter: A local truncation error based adaptive framework for finite volume compressible flow solvers

A residual-based strategy to estimate the local truncation error in a finite volume framework for steady compressible flows is proposed. This estimator, referred to as the -parameter, is derived from the imbalance arising from the use of an exact operator on the numerical solution for conservation laws. The behaviour of the residual estimator for linear and non-linear hyperbolic problems is systematically analysed. The relationship of the residual to the global error is also studied. The -parameter is used to derive a target length scale and consequently devise a suitable criterion for refinement/derefinement. This strategy, devoid of any user-defined parameters, is validated using two standard test cases involving smooth flows. A hybrid adaptive strategy based on both the error indicators and the -parameter, for flows involving shocks is also developed. Numerical studies on several compressible flow cases show that the adaptive algorithm performs excellently well in both two and three dimensions.

Modeling of Convection and Macrosegregation through Appropriate Consideration of Multiphase/Multiscale Phenomena during Alloy Solidification

Solidification processes are complex in nature, involving multiple phases and several length scales. The properties of solidified products are dictated by the microstructure, the mactostructure, and various defects present in the casting. These, in turn, are governed by the multiphase transport phenomena Occurring at different length scales. In order to control and improve the quality of cast products, it is important to have a thorough understanding of various physical and physicochemical phenomena Occurring at various length scales. preferably through predictive models and controlled experiments. In this context, the modeling of transport phenomena during alloy solidification has evolved over the last few decades due to the complex multiscale nature of the problem. Despite this, a model accounting for all the important length scales directly is computationally prohibitive. Thus, in the past, single-phase continuum models have often been employed with respect to a single length scale to model solidification processing. However, continuous development in understanding the physics of solidification at various length scales oil one hand and the phenomenal growth of computational power oil the other have allowed researchers to use increasingly complex multiphase/multiscale models in recent. times. These models have allowed greater understanding of the coupled micro/macro nature of the process and have made it possible to predict solute segregation and microstructure evolution at different length scales. In this paper, a brief overview of the current status of modeling of convection and macrosegregation in alloy solidification processing is presented.

Underwater Sustainability of the ``Cassie'' State of Wetting

A rough hydrophobic surface when immersed in water can result in a ``Cassie'' state of wetting in which the water is in contact with both the solid surface and the entrapped air. The sustainability of the entrapped air on such surfaces is important for underwater applications such as reduction of flow resistance in microchannels and drag reduction of submerged bodies such as hydrofoils. We utilize an optical technique based oil total internal reflection of light at the water-air interface to quantify the spatial distribution of trapped air oil such a surface and its variation with immersion time. With this technique, we evaluate the sustainability of the Cassie state on hydrophobic surfaces with four different kinds of textures. The textures studied are regular arrays of pillars, ridges, and holes that were created in silicon by a wet etching technique, and also a texture of random craters that was obtained through electrodischarge machining of aluminum. These surfaces were rendered hydrophobic with a self-assembled layer Of fluorooctyl trichlorosilane. Depending on the texture, the size and shape of the trapped air pockets were found to vary. However, irrespective of the texture, both the size and the number of air pockets were found to decrease with time gradually and eventually disappear, suggesting that the sustainability of the ``Cassie'' state is finite for all the microstructures Studied. This is possibly due to diffusion of air from the trapped air pockets into the water. The time scale for disappearance of air pockets was found to depend on the kind of microstructure and the hydrostatic pressure at the water-air interface. For the surface with a regular array of pillars, the air pockets were found to be in the form of a thin layer perched on top of the pillars with a large lateral extent compared to the spacing between pillars. For other surfaces studied, the air pockets are smaller and are of the same order as the characteristic length scale of the texture. Measurements for the surface with holes indicate that the time for air-pocket disappearance reduces as the hydrostatic pressure is increased.

Observation of the intense de-pairing effect in YBa2Cu3O7-delta due to the spin injection from La0.5Sr0.5CoO3

We report experimental evidence for a huge pair breaking effect induced by spin polarized quasiparticles in a YBa2Cu3O7-delta/La0.5Sr0.5CoO3 bi-layer fabricated by pulsed laser deposition. The temperature dependent magnetization measurements show evidence for the presence of both ferromagnetic and diamagnetic phases in the bi-layer. The current dependent electrical transport studies in the bi-layer exhibit a significant reduction in the superconducting transition temperature with the increase in applied current as compared to a single YBa2Cu3O7-delta layer and it follows a I-2/3 dependence in accordance with the pair breaking effect. Here, we find that the current driven from a ferromagnetic electrode with low spin polarization, such as La0.5Sr0.5CoO3 (-11%), into the superconductor can act as a strong pair breaker. This indicates that the spin polarization of the injecting electrode is not the only criterion in determining the pair breaking effect, rather the transparency of the interface for the spin polarization may also be significant. More interestingly, the spin diffusion length for YBa2Cu3O7-delta has a much longer length scale than that reported earlier in the study of ferromagnetic/superconducting heterostructures.

Locating the seat of the solar dynamo

The hypothesis that the solar dynamo operates in a thin layer at the bottom of the convection zone is addressed. Recent work on the question whether the magnetic flux can be made to emerge at sunspot latitudes is reviewed. It is concluded that this hypothesis can fit the observational facts only if there is turbulence with a length scale of a few hundred kilometers in and around the dynamo region.

High-strength bulk Al-based bimodal ultrafine eutectic composite with enhanced plasticity

An in situ bulk ultrafine bimodal eutectic Al-Cu-Si composite was synthesized by solidification. This heterostructured composite with microstructural length scale hierarchy in the eutectic microstructure, which combines an ultrafine-scale binary cellular eutectic (alpha-Al + Al2Cu) and a nanometer-sized anomalous ternary eutectic (alpha-Al + Al2Cu + Si), exhibits high fracture strength (1.1 +/- 0.1 GPa) and large compressive plastic strain (11 +/- 2%) at room temperature. The improved compressive plasticity of the bimodal-nanoeutectic composite originates from homogeneous and uniform distribution of inhomogeneous plastic deformation (localized shear bands), together with strong interaction between shear bands in the spatially heterogeneous structure.

Eddy diffusivity: A single dispersion analysis of high resolution drifters in a tidal shallow estuary

In an estuary, mixing and dispersion resulting from turbulence and small scale fluctuation has strong spatio-temporal variability which cannot be resolved in conventional hydrodynamic models while some models employs parameterizations large water bodies. This paper presents small scale diffusivity estimates from high resolution drifters sampled at 10 Hz for periods of about 4 hours to resolve turbulence and shear diffusivity within a tidal shallow estuary (depth < 3 m). Taylor's diffusion theorem forms the basis of a first order estimate for the diffusivity scale. Diffusivity varied between 0.001 – 0.02 m2/s during the flood tide experiment. The diffusivity showed strong dependence (R2 > 0.9) on the horizontal mean velocity within the channel. Enhanced diffusivity caused by shear dispersion resulting from the interaction of large scale flow with the boundary geometries was observed. Turbulence within the shallow channel showed some similarities with the boundary layer flow which include consistency with slope of 5/3 predicted by Kolmogorov's similarity hypothesis within the inertial subrange. The diffusivities scale locally by 4/3 power law following Okubo's scaling and the length scale scales as 3/2 power law of the time scale. The diffusivity scaling herein suggests that the modelling of small scale mixing within tidal shallow estuaries can be approached from classical turbulence scaling upon identifying pertinent parameters.

Surface instability of confined elastic bilayers: Theory and simulations

The surface of a soft elastic film becomes unstable and forms a self-organized undulating pattern because of adhesive interactions when it comes in contact proximity with a rigid surface. For a single film, the pattern length scale lambda, which is governed by the minimization of the elastic stored energy, gives lambda similar to 3h, where h is the film thickness. Based on a linear stability analysis and simulations of adhesion and debonding, we consider the contact instability of an elastic bilayer, which provides greater flexibility in the morphological control of interfacial instability. Unlike the case of a single film, the morphology of the contact instability patterns, debonding distance, and debonding force in a bilayer can be controlled in a nonlinear way by varying the thicknesses and shear moduli of the films. Interestingly, the pattern wavelength in a bilayer can be greatly increased or decreased compared to a single film when the adhesive contact is formed by the stiffer or the softer of the two films, respectively. In particular, lambda as small as 0.5h can be obtained. This indicates a new strategy for pattern miniaturization in elastic contact lithography.

Characterization of a continuous CO2 laser-welded Fe-Cu dissimilar couple

Continuous CO2 laser welding of an Fe-Cu dissimilar couple in a butt-weld geometry at different process conditions is studied. The process conditions are varied to identify and characterize the microstructural features that are independent of the welding mode. The study presents a characterization of the microstructure and mechanical properties of the welds. Detailed microstructural analysis of the weld/base-metal interface shows features that are different on the two sides of the weld. The iron side can grow into the weld with a local change in length scale, whereas the interface on the copper side indicates a barrier to growth. The interface is jagged, and a banded microstructure consisting of iron-rich layers could be observed next to the weld/Cu interface. The observations suggest that solidification initiates inside the melt, where iron and copper are mixed due to convective flow. The transmission electron microscopy (TEM) of the weld region also indicates the occasional presence of droplets of iron and copper. The microstructural observations are rationalized using arguments drawn from a thermodynamic analysis of the Fe-Cu system.

Electronic phase separation in correlated oxides: The phenomenon, its present status and future prospects

Many transition metal oxide materials of high chemical purity are not necessarily monophasic. Thus, single crystals of chemically pure rare earth manganites and cobaltates of the general formula Ln1-xAxMO3 (Ln=rare earth metal, A=alkaline earth metal, M=Mn, Co) exhibit the phenomenon of electronic phase separation wherein phases of different electronic and magnetic properties coexist. Such phase separation, the length scale of which can vary anywhere between a few nanometers to microns, gives distinct signatures in X-ray and neutron diffraction patterns, electrical and magnetic properties, as well as in NMR and other spectroscopies. While the probe one employs to investigate electronic phase separation depends on the length scale, it is noteworthy that direct imaging of the inhomogeneities has been accomplished. Some understanding of this phenomenon has been possible on the basis of some of the theoretical models, but we are far from unraveling the varied aspects of this new phenomenon. Herein, we present the highlights of experimental techniques and theoretical approaches, and comment on the future outlook for this fascinating phenomenon

Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age-old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucual role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large-scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean-field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some of the current issues of the solar simulations are put forward.

Heat transfer correlations for PCM-based heat sinks with plate fins

Characterization of melting process in a Phase Change Material (PCM)-based heat sink with plate fin type thermal conductivity enhancers (TCEs) is numerically studied in this paper. Detailed parametric investigations are performed to find the effect of aspect ratio of enclosure and the applied heat flux on the thermal performance of the heat sinks. Various non-dimensional numbers, such as Nusselt number (Nu), Rayleigh number (Ra), Stefan number (Ste) and Fourier number (Fo) based on a characteristic length scale, are identified as important parameters. The half fin thickness and the fin height are varied to obtain a wide range of aspect ratios of an enclosure. It is found that a single correlation of Nu with Ra is not applicable for all aspect ratios of enclosure with melt convection taken into account. To find appropriate length scales, enclosures with different aspect ratios are divided into three categories, viz. (a) shallow enclosure, (b) rectangular enclosure and (c) tall enclosure. Accordingly, an appropriate characteristic length scale is identified for each type of enclosure and correlation of Nu with Ra based on that characteristic length scale is developed. (C) 2010 Elsevier Ltd. All rights reserved.

Locating the seat of the solar dynamo

The hypothesis that the solar dynamo operates in a thin layer at the bottom of the convection zone is addressed. Recent work on the question whether the magnetic flux can be made to emerge at sunspot latitudes is reviewed. It is concluded that this hypothesis can fit the observational facts only if there is turbulence with a length scale of a few hundred kilometers in and around the dynamo region.