Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Local Structural and Environmental Factors Define the Efficiency of an RNA Pseudoknot Involved in Programmed Ribosomal Frameshift Process

In programmed -1 ribosomal frameshift, an RNA pseudoknot stalls the ribosome at specific sequence and restarts translation in a new reading frame. A precise understanding of structural characteristics of these pseudoknots and their PRF inducing ability has not been clear to date. To investigate this phenomenon, we have studied various structural aspects of a -1 PRF inducing RNA pseudoknot from BWYV using extensive molecular dynamics simulations. A set of functional and poorly functional forms, for which previous mutational data were available, were chosen for analysis. These structures differ from each other by either single base substitutions or base-pair replacements from the native structure. We have rationalized how certain mutations in RNA pseudoknot affect its function; e.g., a specific base substitution in loop 2 stabilizes the junction geometry by forming multiple noncanonical hydrogen bonds, leading to a highly rigid structure that could effectively resist ribosome-induced unfolding, thereby increasing efficiency. While, a CG to AU pair substitution in stem 1 leads to loss of noncanonical hydrogen bonds between stems and loop, resulting in a less stable structure and reduced PRF inducing ability, inversion of a pair in stem 2 alters specific base-pair geometry that might be required in ribosomal recognition of nucleobase groups, negatively affecting pseudoknot functioning. These observations illustrate that the ability of an RNA pseudoknot to induce -1 PRF with an optimal rate depends on several independent factors that contribute to either the local conformational variability or geometry

Zinc oxide based photocatalysis: tailoring surface-bulk structure and related interfacial charge carrier dynamics for better environmental applications

As an alternative to the gold standard TiO2 photocatalyst, the use of zinc oxide (ZnO) as a robust candidate for wastewater treatment is widespread due to its similarity in charge carrier dynamics upon bandgap excitation and the generation of reactive oxygen species in aqueous suspensions with TiO2. However, the large bandgap of ZnO, the massive charge carrier recombination, and the photoinduced corrosion-dissolution at extreme pH conditions, together with the formation of inert Zn(OH)(2) during photocatalytic reactions act as barriers for its extensive applicability. To this end, research has been intensified to improve the performance of ZnO by tailoring its surface-bulk structure and by altering its photogenerated charge transfer pathways with an intention to inhibit the surface-bulk charge carrier recombination. For the first time, the several strategies, such as tailoring the intrinsic defects, surface modification with organic compounds, doping with foreign ions, noble metal deposition, heterostructuring with other semiconductors and modification with carbon nanostructures, which have been successfully employed to improve the photoactivity and stability of ZnO are critically reviewed. Such modifications enhance the charge separation and facilitate the generation of reactive oxygenated free radicals, and also the interaction with the pollutant molecules. The synthetic route to obtain hierarchical nanostructured morphologies and study their impact on the photocatalytic performance is explained by considering the morphological influence and the defect-rich chemistry of ZnO. Finally, the crystal facet engineering of polar and non-polar facets and their relevance in photocatalysis is outlined. It is with this intention that the present review directs the further design, tailoring and tuning of the physico-chemical and optoelectronic properties of ZnO for better applications, ranging from photocatalysis to photovoltaics.

Investigation of the separated region ahead of three-dimensional protuberances on plates and cones in hypersonic flows with laminar boundary layers

Heat transfer rate and pressure measurements were made upstream of surface pro-tuberances on a flat plate and a sharp cone subjected to hypersonic flow in a conventional shock tunnel. Heat flux was measured using platinum thin-film sensors deposited on macor substrate and the pressure measurements were made using fast acting piezoelectric sensors. A distinctive hot spot with highest heat flux was obtained near the foot of the protuberance due to heavy vortex activity in the recirculating region. Schlieren flow visualization was used to capture the shock structures and the separation distance ahead of the protrusions was quantitatively measured for varying protuberance heights. A computational analysis was conducted on the flat plate model using commercial computational fluid dynamics software and the obtained trends of heat flux and pressure were compared with the experimental observation. Experiments were also conducted by physically disturbing the laminar boundary layer to check its effect on the magnitude of the hot spot heat flux. In addition to air, argon was also used as test gas so that the Reynolds number can be varied. (C) 2014 AIP Publishing LLC.

Quantitative estimation of density variation in high-speed flows through inversion of the measured wavefront distortion

A simple method employing an optical probe is presented to measure density variations in a hypersonic flow obstructed by a test model in a typical shock tunnel. The probe has a plane light wave trans-illuminating the flow and casting a shadow of a random dot pattern. Local slopes of the distorted wavefront are obtained from shifts of the dots in the pattern. Local shifts in the dots are accurately measured by cross-correlating local shifted shadows with the corresponding unshifted originals. The measured slopes are suitably unwrapped by using a discrete cosine transform based phase unwrapping procedure and also through iterative procedures. The unwrapped phase information is used in an iterative scheme for a full quantitative recovery of density distribution in the shock around the model through refraction tomographic inversion. Hypersonic flow field parameters around a missile shaped body at a free-stream Mach number of 5.8 measured using this technique are compared with the numerically estimated values. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

Cross-correlation-aided transport in stochastically driven accretion flows

The origin of linear instability resulting in rotating sheared accretion flows has remained a controversial subject for a long time. While some explanations of such non-normal transient growth of disturbances in the Rayleigh stable limit were available for magnetized accretion flows, similar instabilities in the absence of magnetic perturbations remained unexplained. This dichotomy was resolved in two recent publications by Chattopadhyay and co-workers Mukhopadhyay and Chattopadhyay, J. Phys. A 46, 035501 (2013); Nath et al., Phys. Rev. E 88, 013010 (2013)] where it was shown that such instabilities, especially for nonmagnetized accretion flows, were introduced through interaction of the inherent stochastic noise in the system (even a ``cold'' accretion flow at 3000Kis too ``hot'' in the statistical parlance and is capable of inducing strong thermal modes) with the underlying Taylor-Couette flow profiles. Both studies, however, excluded the additional energy influx (or efflux) that could result from nonzero cross correlation of a noise perturbing the velocity flow, say, with the noise that is driving the vorticity flow (or equivalently the magnetic field and magnetic vorticity flow dynamics). Through the introduction of such a time symmetry violating effect, in this article we show that nonzero noise cross correlations essentially renormalize the strength of temporal correlations. Apart from an overall boost in the energy rate (both for spatial and temporal correlations, and hence in the ensemble averaged energy spectra), this results in mutual competition in growth rates of affected variables often resulting in suppression of oscillating Alfven waves at small times while leading to faster saturations at relatively longer time scales. The effects are seen to be more pronounced with magnetic field fluxes where the noise cross correlation magnifies the strength of the field concerned. Another remarkable feature noted specifically for the autocorrelation functions is the removal of energy degeneracy in the temporal profiles of fast growing non-normal modes leading to faster saturation with minimum oscillations. These results, including those presented in the previous two publications, now convincingly explain subcritical transition to turbulence in the linear limit for all possible situations that could now serve as the benchmark for nonlinear stability studies in Keplerian accretion disks.

Tectono-geomorphic and environmental set-up deduced during deposition of Mio-Pleistocene sediments in NW Himalaya, India

Detailed pedofacies characterization along-with lithofacies investigations of the Mio-Pleistocene Siwalik sediments exposed in the Ramnagar sub-basin have been studied so as to elucidate variability in time and space of fluvial processes and the role of intra- and extra-basinal controls on fluvial sedimentation during the evolution of the Himalayan foreland basin (HFB). Dominance of multiple, moderately to strongly developed palaeosol assemblages during deposition of Lower Siwalik (similar to 12-10.8 Ma) sediments suggest that the HFB was marked by Upland set-up of Thomas et al. (2002). Activity of intra-basinal faults on the uplands and deposition of terminal fans at different times caused the development of multiple soils. Further, detailed pedofacies along-with lithofacies studies indicate prevalence of stable tectonic conditions and development of meandering streams with broad floodplains. However, the Middle Siwalik (similar to 10.8-4.92 Ma) sub-group is marked by multistoried sandstones and minor mudstone and mainly weakly developed palaeosols, indicating deposition by large braided rivers in the form of megafans in a Lowland set-up of Thomas et al. (2002). Significant change in nature and size of rivers from the Lower to Middle Siwalik at similar to 10 Ma is found almost throughout of the basin from Kohat Plateau (Pakistan) to Nepal because the Himalayan orogeny witnessed its greatest tectonic upheaval at this time leading to attainment of great heights by the Himalaya, intensification of the monsoon, development of large rivers systems and a high rate of sedimentation, hereby a major change from the Upland set-up to the Lowland set-up over major parts of the HFB. An interesting geomorphic environmental set-up prevailed in the Ramnagar sub-basin during deposition of the studied Upper Siwalik (similar to 4.92 to <1.68 Ma) sediments as observed from the degree of pedogenesis and the type of palaeosols. In general, the Upper Siwalik sub-group in the Ramnagar sub-basin is subdivided from bottom to top into the Purmandal sandstone (4.92-4.49 Ma), Nagrota (4.49-1.68 Ma) and Boulder Conglomerate (<1.68 Ma) formations on the basis of sedimentological characters and change in dominant lithology. Presence of mudstone, a few thin gravel beds and dominant sandstone lithology with weakly to moderately developed palaeosols in the Purmandal sandstone Fm. indicates deposition by shallow braided fluvial streams. The deposition of mudstone dominant Nagrota Fm. with moderately to some well developed palaeosols and a zone of gleyed palaeosols with laminated mudstones and thin sandstones took place in an environment marked by numerous small lakes, water-logged regions and small streams in an environment just south of the Piedmont zone, perhaps similar to what is happening presently in the Upland region/the Upper Gangetic plain. This area is locally called the `Trai region' (Pascoe, 1964). Deposition of Boulder Conglomerate Fm. took place by gravelly braided river system close to the Himalayan Ranges. Activity along the Main Boundary Fault led to progradation of these environments distal-ward and led to development of on the whole a coarsening upward sequence. (C) 2014 Elsevier B.V. All rights reserved.

Kinetic theory for sheared granular flows

Rapid granular flows are far-from-equilibrium-driven dissipative systems where the interaction between the particles dissipates energy, and so a continuous supply of energy is required to agitate the particles and facilitate the rearrangement required for the flow. This is in contrast to flows of molecular fluids, which are usually close to equilibrium, where the molecules are agitated by thermal fluctuations. Sheared granular flows form a class of flows where the energy required for agitating the particles in the flowing state is provided by the mean shear. These flows have been studied using the methods of kinetic theory of gases, where the particles are treated in a manner similar to molecules in a molecular gas, and the interactions between particles are treated as instantaneous energy-dissipating binary collisions. The validity of the assumptions underlying kinetic theory, and their applicability to the idealistic case of dilute sheared granular flows are first discussed. The successes and challenges for applying kinetic theory for realistic dense sheared granular flows are then summarised. (C) 2014 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Environmental monitoring of bacterial contamination and antibiotic resistance patterns of the fecal coliforms isolated from Cauvery River, a major drinking water source in Karnataka, India

The present study focuses prudent elucidation of microbial pollution and antibiotic sensitivity profiling of the fecal coliforms isolated from River Cauvery, a major drinking water source in Karnataka, India. Water samples were collected from ten hotspots during the year 2011-2012. The physiochemical characteristics and microbial count of water samples collected from most of the hotspots exhibited greater biological oxygen demand and bacterial count especially coliforms in comparison with control samples (p <= 0.01). The antibiotic sensitivity testing was performed using 48 antibiotics against the bacterial isolates by disk-diffusion assay. The current study showed that out of 848 bacterial isolates, 93.51 % (n=793) of the isolates were found to be multidrug-resistant to most of the current generation antibiotics. Among the major isolates, 96.46 % (n=273) of the isolates were found to be multidrug-resistant to 30 antibiotics and they were identified to be Escherichia coli by 16S rDNA gene sequencing. Similarly, 93.85 % (n=107), 94.49 % (n=103), and 90.22 % (n=157) of the isolates exhibited multiple drug resistance to 32, 40, and 37 antibiotics, and they were identified to be Enterobacter cloacae, Pseudomonas trivialis, and Shigella sonnei, respectively. The molecular studies suggested the prevalence of blaTEM genes in all the four isolates and dhfr gene in Escherichia coli and Sh. sonnei. Analogously, most of the other Gram-negative bacteria were found to be multidrug-resistant and the Gram-positive bacteria, Staphylococcus spp. isolated from the water samples were found to be methicillin and vancomycin-resistant Staphylococcus aureus. This is probably the first study elucidating the bacterial pollution and antibiotic sensitivity profiling of fecal coliforms isolated from River Cauvery, Karnataka, India.

HYDROMAGNETICS OF ADVECTIVE ACCRETION FLOWS AROUND BLACK HOLES: REMOVAL OF ANGULAR MOMENTUM BY LARGE-SCALE MAGNETIC STRESSES

We show that the removal of angular momentum is possible in the presence of large-scale magnetic stresses in geometrically thick, advective, sub-Keplerian accretion flows around black holes in steady state, in the complete absence of alpha-viscosity. The efficiency of such an angular momentum transfer could be equivalent to that of alpha-viscosity with alpha = 0.01-0.08. Nevertheless, the required field is well below its equipartition value, leading to a magnetically stable disk flow. This is essentially important in order to describe the hard spectral state of the sources when the flow is non/sub-Keplerian. We show in our simpler 1.5 dimensional, vertically averaged disk model that the larger the vertical-gradient of the azimuthal component of the magnetic field is, the stronger the rate of angular momentum transfer becomes, which in turn may lead to a faster rate of outflowing matter. Finding efficient angular momentum transfer in black hole disks via magnetic stresses alone, is very interesting when the generic origin of alpha-viscosity is still being explored.

SUPPORTING ENVIRONMENTALLY-BENIGN DESIGN ELUCIDATING ENVIRONMENTAL IMPACT PROPAGATION IN CONCEPTUAL DESIGN PHASE BY SAPPHIRE MODEL OF CAUSALITY

Conceptual Design Phase is the most critical for design decisions and their impact on the Environment. It is also a phase of many `unknowns' making it flexible and allowing exploration of many solutions. Thus, it is a challenge to determine the most Environmentally-benign Solution or Concept to be translated in to a `good' product. The SAPPhIRE Model captures the various levels of abstractions present in Conceptual Design by Outcomes and defines a Solution-variant as a set of verifiable and quantifiable Outcomes. The Causality explains the propagation of Environmental Impact across Outcomes at varying levels of abstraction, suggesting that the Environmental Impact of an Outcome at a certain level can be represented as a collation of Environmental Impact information of all the Outcomes at each of its subsequent lower levels of abstraction. Thus a ball-park impact value can be associated with the higher-levels of abstraction, thereby supporting design decisions taken earlier on in Conceptual Design directing towards Environmentally-benign Design.

Fingerprinting environmental conditions and related stress using stable isotopic composition of rice (Oryza sativa L.) grain organic matter

Sea level rise (SLR) is a primary factor responsible for inundation of low-lying coastal regions across the world, which in turn governs the agricultural productivity. In this study, rice (Oryza sativa L.) cultivated seasonally in the Kuttanad Wetland, a SLR prone region on the southwest coast of India, were analysed for oxygen, hydrogen and carbon isotopic ratios (delta O-18, delta H-2 and delta C-13) to distinguish the seasonal environmental conditions prevalent during rice cultivation. The region receives high rainfall during the wet season which promotes large supply of fresh water to the local water bodies via the rivers. In contrast, during the dry season reduced river discharge favours sea water incursion which adversely affects the rice cultivation. The water for rice cultivation is derived from regional water bodies that are characterised by seasonal salinity variation which co-varies with the delta O-18 and delta H-2 values. Rice cultivated during the wet and the dry season bears the isotopic imprints of this water. We explored the utility of a mechanistic model to quantify the contribution of two prominent factors, namely relative humidity and source water composition in governing the seasonal variation in oxygen isotopic composition of rice grain OM. delta C-13 values of rice grain OM were used to deduce the stress level by estimating the intrinsic water use efficiency (WUEi) of the crop during the two seasons. 1.3 times higher WUE, was exhibited by the same genotype during the dry season. The approach can be extended to other low lying coastal agro-ecosystems to infer the growth conditions of cultivated crops and can further be utilised for retrieving paleo-environmental information from well preserved archaeological plant remains. (c) 2015 Elsevier Ltd. All rights reserved.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.