Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.

Retention characteristics of Cu2+, Pb2+, and Zn2+ from aqueous solutions by two types of low lime fly ashes

The technical feasibility of utilization of fly ash as a low-cost adsorbent for the removal of metals from water has been studied. For two types of fly ashes, the retention capacities of copper, lead, and zinc metal ions have been studied. Contact time, initial concentration, and pH have been varied and their effect on retention mechanism has been studied. The dominant mechanisms responsible for retention are found to be precipitation due to the presence of calcium hydroxide, and adsorption due to the presence of silica and alumina oxide surfaces in the fly ash. First-order kinetic plots have revealed that the rate constant increases with increase in the initial concentration and pH. Langmuir adsorption isotherms have been plotted to study the maximum adsorption capacities for metal ions considered under different conditions. X-ray diffraction studies revealed the formation of new peaks corresponding to respective metal ions precipitates under alkaline conditions.

Temporal Moments for Reactive Transport through Fractured Impermeable/Permeable Formations

The transport of reactive solutes through fractured porous formations has been analyzed. The transport through the porous block is represented by a general multiprocess nonequilibrium equation (MPNE), which, for the fracture, is represented by an advection-dispersion equation with linear equilibrium sorption and first-order transformation. An implicit finite-difference technique has been used to solve the two coupled equations. The transport characteristics have been analyzed in terms of zeroth, first, and second temporal moments of the solute in the fracture. The solute behavior for fractured impermeable and fractured permeable formations are first compared and the effects of various fracture and matrix transport parameters are analyzed. Subsequently, the transport through a fractured permeable formation is analyzed to ascertain the effect of equilibrium sorption, rate-limited sorption, and the multiprocess nonequilibrium transport process. It was found that the temporal moments were nearly identical for the fractured impermeable and permeable formations when both the diffusion coefficient and the first-order transformation coefficient were relatively large. The multiprocess nonequilibrium model resulted in a smaller mass recovery in the fracture and higher dispersion than the equilibrium and rate-limited sorption models. DOI: 10.1061/(ASCE)HE.19435584.0000586. (C) 2012 American Society of Civil Engineers.

Simulation of Influence of Bilayer Melting on Dynamics and Thermodynamics of Interfacial Water

We investigate the effect of bilayer melting transition on thermodynamics and dynamics of interfacial water using molecular dynamics simulation with the two-phase thermodynamic model. We show that the diffusivity of interface water depicts a dynamic crossover at the chain melting transition following an Arrhenius behavior until the transition temperature. The corresponding change in the diffusion coefficient from the bulk to the interface water is comparable with experimental observations found recently for water near 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles Phys. Chem. Chem. Phys. 13, 7732 (2011)]. The entropy and potential energy of interfacial water show distinct changes at the bilayer melting transition, indicating a strong correlation in the thermodynamic state of water and the accompanying first-order phase transition of the bilayer membrane. DOI: 10.1103/PhysRevLett.110.018303

Updating reliability models of statically loaded instrumented structures

The study extends the first order reliability method (FORM) and inverse FORM to update reliability models for existing, statically loaded structures based on measured responses. Solutions based on Bayes' theorem, Markov chain Monte Carlo simulations, and inverse reliability analysis are developed. The case of linear systems with Gaussian uncertainties and linear performance functions is shown to be exactly solvable. FORM and inverse reliability based methods are subsequently developed to deal with more general problems. The proposed procedures are implemented by combining Matlab based reliability modules with finite element models residing on the Abaqus software. Numerical illustrations on linear and nonlinear frames are presented. (c) 2012 Elsevier Ltd. All rights reserved.

Efficient Methods for Robust Classification Under Uncertainty in Kernel Matrices

In this paper we study the problem of designing SVM classifiers when the kernel matrix, K, is affected by uncertainty. Specifically K is modeled as a positive affine combination of given positive semi definite kernels, with the coefficients ranging in a norm-bounded uncertainty set. We treat the problem using the Robust Optimization methodology. This reduces the uncertain SVM problem into a deterministic conic quadratic problem which can be solved in principle by a polynomial time Interior Point (IP) algorithm. However, for large-scale classification problems, IP methods become intractable and one has to resort to first-order gradient type methods. The strategy we use here is to reformulate the robust counterpart of the uncertain SVM problem as a saddle point problem and employ a special gradient scheme which works directly on the convex-concave saddle function. The algorithm is a simplified version of a general scheme due to Juditski and Nemirovski (2011). It achieves an O(1/T-2) reduction of the initial error after T iterations. A comprehensive empirical study on both synthetic data and real-world protein structure data sets show that the proposed formulations achieve the desired robustness, and the saddle point based algorithm outperforms the IP method significantly.

Sensory-organ-like response determines the magnetism of zigzag-edged honeycomb nanoribbons

We present an analytical effective theory for the magnetic phase diagram for zigzag-edge terminated honeycomb nanoribbons described by a Hubbard model with an interaction parameter U. We show that the edge magnetic moment varies as ln U and uncover its dependence on the width W of the ribbon. The physics of this owes its origin to the sensory-organ-like response of the nanoribbons, demonstrating that considerations beyond the usual Stoner-Landau theory are necessary to understand the magnetism of these systems. A first-order magnetic transition from an antiparallel orientation of the moments on opposite edges to a parallel orientation occurs upon doping with holes or electrons. The critical doping for this transition is shown to depend inversely on the width of the ribbon. Using variational Monte Carlo calculations, we show that magnetism is robust to fluctuations. Additionally, we show that the magnetic phase diagram is generic to zigzag-edge terminated nanostructures such as nanodots. Furthermore, we perform first-principles modeling to show how such magnetic transitions can be realized in substituted graphene nanoribbons. DOI: 10.1103/PhysRevB.87.085412

Effect of Metal Ions (Ag, Co, Ni, and Pd) on the Visible Light Degradation of Rhodamine B by Carbon-Covered Alumina-Supported TiO2 in Aqueous Solutions

Metal-ion (Ag, Co, Ni, and Pd) doped TiO2 nanocatalysts were successfully embedded on carbon-covered alumina supports. The CCA-embedded catalysts were crystalline and had a high surface area compared to the free metal-ion doped titania nanocatalysts while they still retained the anatase phase of the core TiO2. These catalysts were photocatalytically active under solar light irradiation. Rhodamine B was used as a model pollutant and the reactivity followed a pseudo-first-order reaction kinetics. The reaction rate of the CCA-supported catalysts was Pd > Ag > Co > Ni. Among the ratios of the CCA:catalyst used, it was found that the 1:1 ratio had the fastest reaction rate, followed by the 1:2 ratio, while the 2:1 ratio exhibited the lowest reaction rate. The CCA/metal-ion doped titania were found to have photocatalytic activities comparable with those of CCA-supported titania.

Structure, tensile properties and cytotoxicity assessment of sebacic acid based biodegradable polyesters with ricinoleic acid

A new family of ricinoleic acid based polyesters was synthesized using catalyst free melt-condensation polymerization with sebacic acid, citric acid, mannitol and ricinoleic acid as precursors. The use of FT-IR and NMR characterisation techniques confirms the presence of ester linkages in the as-synthesized polymers. Depending on the precursor combination, their relative amount and the degree of curing, a broad range of elastic modulus (22-327 MPa) and tensile strength (0.7-12.7 MPa) can be obtained in the newly synthesized biopolymers. The polymers show rubbery behaviour at a physiological temperature (37 degrees C) and the contact angles of the synthesized polymers fall in the range of 42 degrees to 71 degrees, making them ideal substrates to study delivery of drugs through polymer scaffolds. The cytocompatibility assessment of the cured polymers confirmed good cell attachment and growth of smooth muscle cells (C2C12 myoblast cells). Importantly, oriented cell growth was observed after culturing myoblast cells for 3 days. The in vitro degradation in PBS indicates that the mild cured polymers follow a first order reaction kinetics and have degradation rate constants in the range of 0.009-0.038 h(-1), depending on the relative proportions of monomers. Overall, the results of our study indicate that the physical properties can be tailored by varying the composition of the monomers and curing conditions in the newly developed polyesters. Hence, they may be used as potential substrates for tissue engineering scaffolds and for localized drug delivery.

Organophotocatalysis in nanostructured soft gel materials as tunable reaction vessels: comparison with homogeneous and micellar solutions

Riboflavin tetraacetate-catalyzed aerobic photooxidation of 1-(4-methoxyphenyl)ethanol was investigated as a model reaction under blue visible light in different soft gel materials, aiming to establish their potential as reaction vessels for photochemical transformations. Three strategies involving different degrees of organization of the catalyst within the gel network were explored, and the results compared to those obtained in homogeneous and micellar solutions. In general, physical entrapment of both the catalyst and the substrate under optimized concentrations into several hydrogel matrices (including low-molecular-weight and biopolymer-based gels) allowed the photooxidation with conversions between 55 and 100% within 120 min (TOF similar to 0.045-0.08 min(-1); k(obs) similar to 0.011-0.028 min(-1)), albeit with first-order rates ca. 1-3-fold lower than in solution under comparable non-stirred conditions. Remarkably, the organogel made of a cyclohexane-based bisamide gelator in CH3CN not only prevented the photodegradation of the catalyst but also afforded full conversion in less than 60 min (TOF similar to 0.167 min(-1); k(obs) similar to 0.073 min(-1)) without the need of additional proton transfer mediators (e. g., thiourea) as it occurs in CH3CN solutions. In general, the gelators could be recycled without detriment to their gelation ability and reaction rates. Moreover, kinetics could be fine-tuned according to the characteristics of the gel media. For instance, entangled fibrillar networks with relatively high mechanical strength were usually associated with lower reaction rates, whereas wrinkled laminated morphologies seemed to favor the reaction. In addition, the kinetics results showed in most cases a good correlation with the aeration efficiency of the gel media.

Simulation of Deformation Texture Evolution During Multi Axial Forging of Interstitial Free Steel

Bulk texture measurement of multi-axial forged body center cubic interstitial free steel performed in this study using x-ray and neutron diffraction indicated the presence of a strong {101}aOE (c) 111 > single texture component. Viscoplastic self-consistent simulations could successfully predict the formation of this texture component by incorporating the complicated strain path followed during this process and assuming the activity of {101}aOE (c) 111 > slip system. In addition, a first-order estimate of mechanical properties in terms of highly anisotropic yield locus and Lankford parameter was also obtained from the simulations.

Longitudinal localization of tropical intraseasonal variability

Intraseasonal time-scales play an important role in tropical variability. Two modes that contribute significantly to tropical intraseasonal variability (ISV) are the eastward-propagating MaddenJulian Oscillation (MJO), and westward-moving moist equatorial Rossby waves. This note reports on a correspondence between the longitudinal gradient of mean tropical precipitable water (PW), and the geographical regions of genesis, and convective activity, of both these large-scale tropical systems. Our finding is based on an analysis of PW from the MERRA reanalysis product. The data indicate that the mean tropical PW has a dominant wavenumber two (three) structure in longitude in the Northern (Southern) Hemisphere. Departures from a longitudinally homogeneous state are attributed to the influence of subtropical anticyclones, and are accentuated by monsoonal regions of both hemispheres. This mean structure results in a sharply localized longitudinal gradient of PW. Remarkably, regions with positive gradients (such as the Northern and Southern Hemisphere western Indian Ocean), i.e. they have larger PW to the east, are the very zones that are implicated in the formation, and show high levels of convective activity, of the eastward-moving MJO. On the other hand, regions with negative gradients (such as the Southern Hemisphere central Pacific) are the very regions where genesis, and maxima in variance, of westward-moving moist equatorial Rossby waves are known to occur. Apart from providing a first-order longitudinal footprint of the convective phase of these systems, this correspondence reinforces the role of the mean climatic state in tropical ISV. Copyright (c) 2012 Royal Meteorological Society

A risk-estimation-based formulation for speech enhancement and its relation to Wiener filtering

The goal of speech enhancement algorithms is to provide an estimate of clean speech starting from noisy observations. The often-employed cost function is the mean square error (MSE). However, the MSE can never be computed in practice. Therefore, it becomes necessary to find practical alternatives to the MSE. In image denoising problems, the cost function (also referred to as risk) is often replaced by an unbiased estimator. Motivated by this approach, we reformulate the problem of speech enhancement from the perspective of risk minimization. Some recent contributions in risk estimation have employed Stein's unbiased risk estimator (SURE) together with a parametric denoising function, which is a linear expansion of threshold/bases (LET). We show that the first-order case of SURE-LET results in a Wiener-filter type solution if the denoising function is made frequency-dependent. We also provide enhancement results obtained with both techniques and characterize the improvement by means of local as well as global SNR calculations.

Phase synchronization in brain networks derived from correlation between probabilities of recurrences in functional MRI data

It is increasingly being recognized that resting state brain connectivity derived from functional magnetic resonance imaging (fMRI) data is an important marker of brain function both in healthy and clinical populations. Though linear correlation has been extensively used to characterize brain connectivity, it is limited to detecting first order dependencies. In this study, we propose a framework where in phase synchronization (PS) between brain regions is characterized using a new metric ``correlation between probabilities of recurrence'' (CPR) and subsequent graph-theoretic analysis of the ensuing networks. We applied this method to resting state fMRI data obtained from human subjects with and without administration of propofol anesthetic. Our results showed decreased PS during anesthesia and a biologically more plausible community structure using CPR rather than linear correlation. We conclude that CPR provides an attractive nonparametric method for modeling interactions in brain networks as compared to standard correlation for obtaining physiologically meaningful insights about brain function.

Quantum phases and phase transitions of frustrated hard-core bosons on a triangular ladder

Kinetically frustrated bosons at half filling in the presence of a competing nearest-neighbor repulsion support a wide supersolid regime on the two-dimensional triangular lattice. We study this model on a two-leg ladder using the finite-size density-matrix renormalization-group method, obtaining a phase diagram which contains three phases: a uniform superfluid (SF), an insulating charge density wave (CDW) crystal, and a bond ordered insulator (BO). We show that the transitions from SF to CDW and SF to BO are continuous in nature, with critical exponents varying continuously along the phase boundaries, while the transition from CDW to BO is found to be first order. The phase diagram is also found to contain an exactly solvable Majumdar Ghosh point, and reentrant SF to CDW phase transitions.