Modelling of Effect of Non-Uniform Current Density on the Performance of Soluble Lead Redox Flow Batteries

Soluble lead acid redox flow battery (SLRFB) offers a number of advantages. These advantages can be harnessed after problems associated with buildup of active material on. electrodes (residue) are resolved. A mathematical model is developed to understand residue formation in SLRFB. The model incorporates fluid flow, ion transport, electrode reactions, and non-uniform current distribution on electrode surfaces. A number of limiting cases are studied to conclude that ion transport and electrode reaction on anode simultaneously control battery performance. The model fits the reported cell voltage vs. time profiles very well. During the discharge cycle, the model predicts complete dissolution of deposited material from trailing edge side of the electrodes. With time, the active surface area of electrodes decreases rapidly. The corresponding increase in current density leads to precipitous decrease in cell potential before all the deposited material is dissolved. The successive charge-discharge cycles add to the residue. The model correctly captures the marginal effect of flow rate on cell voltage profiles, and identifies flow rate and flow direction as new variables for controlling residue buildup. Simulations carried out with alternating flow direction and a SLRFB with cylindrical electrodes show improved performance with respect to energy efficiency and residue buildup. (C) 2014 The Electrochemical Society. All rights reserved.

Modelling of optical transport behavior of organic photovoltaic devices with nano-pillar transparent conducting electrodes

Optical transport behavior of organic photo-voltaic devices with nano-pillar transparent electrodes is investigated in this paper in order to understand possible enhancement of their charge-collection efficiency. Modeling and simulations of optical transport due to this architecture show an interesting regime of length-scale dependent optical characteristics. An electromagnetic wave propagation model is employed with simulation objectives toward understanding the mechanism of optical scattering and waveguide effects due to the nano-pillars and effective transmission through the active layer. Partial filling of gaps between the nano-pillars due to the nano-fabrication process is taken into consideration. Observations made in this paper will facilitate appropriate design rules for nano-pillar electrodes. (C) 2014 AIP Publishing LLC.

Social insect colony as a biological regulatory system: modelling information flow in dominance networks

Social insects provide an excellent platform to investigate flow of information in regulatory systems since their successful social organization is essentially achieved by effective information transfer through complex connectivity patterns among the colony members. Network representation of such behavioural interactions offers a powerful tool for structural as well as dynamical analysis of the underlying regulatory systems. In this paper, we focus on the dominance interaction networks in the tropical social wasp Ropalidia marginata-a species where behavioural observations indicate that such interactions are principally responsible for the transfer of information between individuals about their colony needs, resulting in a regulation of their own activities. Our research reveals that the dominance networks of R. marginata are structurally similar to a class of naturally evolved information processing networks, a fact confirmed also by the predominance of a specific substructure-the `feed-forward loop'-a key functional component in many other information transfer networks. The dynamical analysis through Boolean modelling confirms that the networks are sufficiently stable under small fluctuations and yet capable of more efficient information transfer compared to their randomized counterparts. Our results suggest the involvement of a common structural design principle in different biological regulatory systems and a possible similarity with respect to the effect of selection on the organization levels of such systems. The findings are also consistent with the hypothesis that dominance behaviour has been shaped by natural selection to co-opt the information transfer process in such social insect species, in addition to its primal function of mediation of reproductive competition in the colony.

The irregularities of the sunspot cycle and their theoretical modelling

The 11-year sunspot cycle has many irregularities, the most prominent amongst them being the grand minima when sunspots may not be seen for several cycles. After summarizing the relevant observational data about the irregularities, we introduce the flux transport dynamo model, the currently most successful theoretical model for explaining the 11-year sunspot cycle. Then we analyze the respective roles of nonlinearities and random fluctuations in creating the irregularities. We also discuss how it has recently been realized that the fluctuations in meridional circulation also can be a source of irregularities. We end by pointing out that fluctuations in the poloidal field generation and fluctuations in meridional circulation together can explain the occurrences of grand minima.

Beyond classical dynamic structural plasticity using mesh-free modelling techniques

The problem of modelling the transient response of an elastic-perfectly-plastic cantilever beam, carrying an impulsively loaded tip mass, is,often referred to as the Parkes cantilever problem 25]; The permanent deformation of a cantilever struck transversely at its tip, Proc. R. Soc. A., 288, pp. 462). This paradigm for classical modelling of projectile impact on structures is re-visited and updated using the mesh-free method, smoothed particle hydrodynamics (SPH). The purpose of this study is to investigate further the behaviour of cantilever beams subjected to projectile impact at its tip, by considering especially physically real effects such as plastic shearing close to the projectile, shear deformation, and the variation of the shear strain along the length and across the thickness of the beam. Finally, going beyond macroscopic structural plasticity, a strategy to incorporate physical discontinuity (due to crack formation) in SPH discretization is discussed and explored in the context of tip-severance of the cantilever beam. Consequently, the proposed scheme illustrates the potency for a more refined treatment of penetration mechanics, paramount in the exploration of structural response under ballistic loading. The objective is to contribute to formulating a computational modelling framework within which transient dynamic plasticity and even penetration/failure phenomena for a range of materials, structures and impact conditions can be explored ab initio, this being essential for arriving at suitable tools for the design of armour systems. (C) 2014 Elsevier Ltd. All rights reserved.

An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

In this paper we present a massively parallel open source solver for Richards equation, named the RichardsFOAM solver. This solver has been developed in the framework of the open source generalist computational fluid dynamics tool box OpenFOAM (R) and is capable to deal with large scale problems in both space and time. The source code for RichardsFOAM may be downloaded from the CPC program library website. It exhibits good parallel performances (up to similar to 90% parallel efficiency with 1024 processors both in strong and weak scaling), and the conditions required for obtaining such performances are analysed and discussed. These performances enable the mechanistic modelling of water fluxes at the scale of experimental watersheds (up to few square kilometres of surface area), and on time scales of decades to a century. Such a solver can be useful in various applications, such as environmental engineering for long term transport of pollutants in soils, water engineering for assessing the impact of land settlement on water resources, or in the study of weathering processes on the watersheds. (C) 2014 Elsevier B.V. All rights reserved.

Modelling of transient heat conduction with diffuse interface methods

We present a survey on different numerical interpolation schemes used for two-phase transient heat conduction problems in the context of interface capturing phase-field methods. Examples are general transport problems in the context of diffuse interface methods with a non-equal heat conductivity in normal and tangential directions to the interface. We extend the tonsorial approach recently published by Nicoli M et al (2011 Phys. Rev. E 84 1-6) to the general three-dimensional (3D) transient evolution equations. Validations for one-dimensional, two-dimensional and 3D transient test cases are provided, and the results are in good agreement with analytical and numerical reference solutions.

Thermally Induced Demixing in an LCST Mixture in the Presence of Densely Grafted Nanoparticles: Tuning the Graft Chain Length To Induce Thermodynamic Miscibility

The demixing in an LCST mixture of PS/PVME (polystyrene/poly(vinyl methyl ether)) was probed here by melt rheology in the presence of gold nanoparticles which were densely coated with varying graft lengths of PS. The graft density for the gold nanoparticles coated with 3 kDa PS was ca. Sigma = 1.7 chains/nm(2), and that for 53 kDa PS was ca. Sigma = 1.2 chains/nm(2). The evolution of morphology, as the blends transit through the metastable and the unstable envelopes of the phase diagram, and the localization of the gold nanoparticles upon demixing were monitored using in situ hot-stage AFM and confocal Raman imaging. Interestingly, gold nanoparticles coated with 3 kDa polystyrene (PS(3 kDa)-g-nAu) were localized in the PVME phase, whereas gold nanoparticles coated with 53 kDa polystyrene (PS(53 kDa)-g-nAu) were localized in the PS phase of the blend. While the localization of PS(3 kDa)-g-nAu in the PVME phase can be expected to be of entropic origin due to expulsion from the PS phase as R-g,R-matrix chains > R-g,R-grafted chains (where R-g is the radius of gyration of the polymer chain), the localization of PS(53 kDa)-g-nAu in the PS phase is believed to be facilitated by favorable melt/graft interactions. The latter nanoparticles also delayed the demixing by 12 degrees C with respect to the neat mixture. The observed changes were addressed in context to enthalpic interactions between the grafted PS and the free PS, the entropic losses (deformational entropic losses on blending, translational entropic loss of the free PS, and the conformational entropic loss of the grafted PS), and the interface of the grafted and the free chains.

Electromagnetic shielding materials and coatings derived from gelation of multiwall carbon nanotubes in an LCST mixture

Thermally induced demixing in an LCST mixture, polystyrene (PS)/polyvinyl methyl ether] (PVME), was used as a template to design materials with high electrical conductivity. This was facilitated by gelation of multiwall carbon nanotubes (MWNTs) in a given phase of the blends. The MWNTs were mixed in the miscible blends and the thermodynamic driven demixing further resulted in selective localization in the PVME phase of the blends. This was further confirmed by atomic force microscopy (AFM). The time dependent gelation of MWNTs at shallow quench depth, evaluated using isochronal temperature sweep by rheology, was studied by monitoring the melt electrical conductivity of the samples in situ by an LCR meter coupled to a rheometer. By varying the composition in the mixture, several intricate shapes like gaskets and also coatings capable of attenuating the EM radiation in the microwave frequency can be derived. For instance, the PVME rich mixtures can be molded in the form of a gasket, O-ring and other intricate shapes while the PS rich mixtures can be coated onto an insulating polymer to enhance the shielding effectiveness (SE) for EM radiation. The SE of the various materials was analyzed using a vector network analyzer in both the X-band (8.2 to 12 GHz) and the K-u-band (12 to 18 GHz) frequency. The improved SE upon gelation of MWNTs in the demixed blends is well evident by comparing the SE before and after demixing. A reflection loss of -35 dB was observed in the blends with 2 wt% MWNTs. Further, by coating a layer of ca. 0.15 mm of PS/PVME/MWNT, a SE of -15 dB at 18 GHz could be obtained.

Mixture of Fuels Approach for the Synthesis of SrFeO3-delta Nanocatalyst and Its Impact on the Catalytic Reduction of Nitrobenzene

A modified solution combustion approach was applied in the synthesis of nanosize SrFeO3-delta (SFO) using single as well as mixture of citric acid, oxalic acid, and glycine as fuels with corresponding metal nitrates as precursors. The synthesized and calcined powders were characterized by Fourier transform infrared spectroscopy (FT-IR), X-ray diffraction (XRD), thermogravimetric analysis and derivative thermogravimetric analysis (TG-DTG), scanning electron microscopy, transmission electron microscopy, N-2 physisorption methods, and acidic strength by n-butyl amine titration methods. The FT-IR spectra show the lower-frequency band at 599 cm(-1) corresponds to metal-oxygen bond (possible Fe-O stretching frequencies) vibrations for the perovskite-structure compound. TG-DTG confirms the formation temperature of SFO ranging between 850-900 degrees C. XRD results reveal that the use of mixture of fuels in the preparation has effect on the crystallite size of the resultant compound. The average particle size of the samples prepared from single fuels as determined from XRD was similar to 50-35 nm, whereas for samples obtained from mixture of fuels, particles with a size of 30-25 nm were obtained. Specifically, the combination of mixture of fuels for the synthesis of SFO catalysts prevents agglomeration of the particles, which in turn leads to decrease in crystallite size and increase in the surface area of the catalysts. It was also observed that the present approach also impacted the catalytic activity of the SFO in the catalytic reduction of nitrobenzene to azoxybenzene.

MULTI-PITCH TRACKING USING GAUSSIAN MIXTURE MODEL WITH TIME VARYING PARAMETERS AND GRATING COMPRESSION TRANSFORM

Grating Compression Transform (GCT) is a two-dimensional analysis of speech signal which has been shown to be effective in multi-pitch tracking in speech mixtures. Multi-pitch tracking methods using GCT apply Kalman filter framework to obtain pitch tracks which requires training of the filter parameters using true pitch tracks. We propose an unsupervised method for obtaining multiple pitch tracks. In the proposed method, multiple pitch tracks are modeled using time-varying means of a Gaussian mixture model (GMM), referred to as TVGMM. The TVGMM parameters are estimated using multiple pitch values at each frame in a given utterance obtained from different patches of the spectrogram using GCT. We evaluate the performance of the proposed method on all voiced speech mixtures as well as random speech mixtures having well separated and close pitch tracks. TVGMM achieves multi-pitch tracking with 51% and 53% multi-pitch estimates having error <= 20% for random mixtures and all-voiced mixtures respectively. TVGMM also results in lower root mean squared error in pitch track estimation compared to that by Kalman filtering.

The key role of polymer grafted nanoparticles in the phase miscibility of an LCST mixture

Blends of bromo-terminated polystyrene (PS-Br) and poly(vinyl methylether) (PVME) exhibit lower critical solution temperatures. In this study, PS-Br was designed by atom transfer radical polymerization and was converted to thiol-capped polystyrene (PS-SH) by reacting with thiourea. The silver nanoparticles (nAg) were then decorated with covalently bound PS-SH macromolecules to improve the phase miscibility in the PS-Br-PVME blends. Thermally induced demixing in this model blend was followed in the presence of polystyrene immobilized silver nanoparticles (PS-g-nAg). The graft density of the PS macromolecules was estimated to be ca. 0.78 chains per nm(2). Although the matrix and the grafted molecular weights were similar, PS-g-nAg particles were expelled from the PS phase and were localized in the PVME phase of the blends. This was addressed with respect to intermediate graft density and favourable PS-PVME contacts from microscopic interactions point of view. Interestingly, blends with 0.5 wt% PS-g-nAg delayed the spinodal decomposition temperature in the blends by ca. 18 degrees C with respect to the control blends. The scale of cooperativity, as determined by differential scanning calorimetry, increased only marginally in the case of PS-g-nAg; however, it increased significantly in the presence of bare nAg particles.

Identifying regions for conservation of sloth bears through occupancy modelling in north-eastern Karnataka, India

In the absence of information on species in decline with contracting ranges, management should emphasize remaining populations and protection of their habitats. Threatened by anthropogenic pressure including habitat degradation and loss, sloth bears (Melursus ursinus) in India have become limited in range, habitat, and population size. We identified ecological and anthropogenic determinants of occurrence within an occupancy framework to evaluate habitat suitability of non-protected regions (with sloth bears) in northeastern Karnataka, India. We employed a systematic sampling methodology to yield presence absence data to examine a priori hypotheses of determinants that affected occupancy. These covariates were broadly classified as habitat or anthropogenic factors. Mean number of termite mounds and trees positively influenced sloth bear occupancy, and grazing pressure expounded by mean number of livestock dung affected it negatively. Also, mean percentage of shrub coverage had no impact on bear inhabitance. The best fitting model further predicted habitats in Bukkasagara, Agoli, and Benakal reserved forests to have 38%, 75%, and 88%, respectively, of their sampled grid cells with high occupancies (>0.70) albeit little or no legal protection. We recommend a conservation strategy that includes protection of vegetation stand-structure, maintenance of soil moisture, and enrichment of habitat for the long-term welfare of this species.

Modelling of soldier fly halteres for gyroscopic oscillations

Nature has evolved a beautiful design for small-scale vibratory rategyro in the form of dipteran halteres that detect body rotations via Coriolis acceleration. In most Diptera, including soldier fly, Hermetia illucens, halteres are a pair of special organs, located in the space between the thorax and the abdomen. The halteres along with their connecting joint with the fly's body constitute a mechanism that is used for muscle-actuated oscillations of the halteres along the actuation direction. These oscillations lead to bending vibrations in the sensing direction (out of the haltere's actuation plane) upon any impressed rotation due to the resulting Coriolis force. This induced vibration is sensed by the sensory organs at the base of the haltere in order to determine the rate of rotation. In this study, we evaluate the boundary conditions and the stiffness of the anesthetized halteres along the actuation and the sensing direction. We take several cross-sectional SEM (scanning electron microscope) images of the soldier fly haltere and construct its three dimensional model to get the mass properties. Based on these measurements, we estimate the natural frequency along both actuation and sensing directions, propose a finite element model of the haltere's joint mechanism, and discuss the significance of the haltere's asymmetric cross-section. The estimated natural frequency along the actuation direction is within the range of the haltere's flapping frequency. However, the natural frequency along the sensing direction is roughly double the haltere's flapping frequency that provides a large bandwidth for sensing the rate of rotation to the soldier flies.

REDUCED ORDER MODELLING OF COMBUSTION INSTABILITY IN A BACKWARD FACING STEP COMBUSTOR

This paper develops a fully coupled time domain Reduced Order Modelling (ROM) approach to model unsteady combustion dynamics in a backward facing step combustor The acoustic field equations are projected onto the canonical acoustic eigenmodes of the systems to obtain a coupled system of modal evolution equations. The heat release response of the flame is modelled using the G-equation approach. Vortical velocity fluctuations that arise due to shear layer rollup downstream of the step are modelled using a simplified 1D-advection equation whose phase speed is determined from a linear, local, temporal stability analysis of the shear layer just downstream of the step. The hydrodynamic stability analysis reveals a abrupt change in the value of disturbance phase speed from unity for Re < Re-crit to 0.5 for Re > Re-crit, where Remit for the present geometry was found to be approximate to 10425. The results for self-excited flame response show highly wrinkled flame shapes that are qualitatively similar to those seen in prior experiments of acoustically forced flames. The effect of constructive and destructive interference between the two contributions to flame surface wrinkling results in high amplitude wrinkles for the case when K-c -> 1.