Modelling the impact of extensive irrigation on the groundwater resources

Drastic groundwater resource depletion due to excessive extraction for irrigation is a major concern in many parts of India. In this study, an attempt was made to simulate the groundwater scenario of the catchment using ArcSWAT. Due to the restriction on the maximum initial storage, the deep aquifer component in ArcSWAT was found to be insufficient to represent the excessive groundwater depletion scenario. Hence, a separate water balance model was used for simulating the deep aquifer water table. This approach is demonstrated through a case study for the Malaprabha catchment in India. Multi-site rainfall data was used to represent the spatial variation in the catchment climatology. Model parameters were calibrated using observed monthly stream flow data. Groundwater table simulation was validated using the qualitative information available from the field. The stream flow was found to be well simulated in the model. The simulated groundwater table fluctuation is also matching reasonably well with the field observations. From the model simulations, deep aquifer water table fluctuation was found very severe in the semi-arid lower parts of the catchment, with some areas showing around 60m depletion over a period of eight years. Copyright (c) 2012 John Wiley & Sons, Ltd.

Note on modelling directionality in piezoresistivity

When computing the change in electrical resistivity of a piezoresistive cubic material embedded in a deforming structure, the piezoresistive and the stress tensors should be in the same coordinate system. While the stress tensor is usually calculated in a coordinate system aligned with the principal axes of a regular structure, the specified piezoresistive coefficients may not be in that coordinate system. For instance, piezoresistive coefficients are usually given in an orthogonal cartesian coordinate system aligned with the <100> crystallographic directions and designers sometimes deliberately orient a crystallographic direction other than <100> along the principal directions of the structure to increase the gauge factor. In such structures, it is advantageous to calculate the piezoresistivity tensor in the coordinate system along which the stress tensors are known rather than the other way around. This is because the transformation of stress will have to be done at every point in the structure but piezoresistivity tensor needs to be transformed only once. Here, using tensor transformation relations, we show how to calculate the piezoresistive tensor along any arbitrary Cartesian coordinate system from the piezoresistive coefficients for the <100> coordinate system. Some of the software packages that simulate the piezoresistive effect do not have interfaces for calculation of the entire piezoresistive tensor for arbitrary directions. This warrants additional work for the user because not considering the complete piezoresisitive tensor can lead to large errors. This is illustrated with an example where the error is as high as 33%. Additionally, for elastic analysis, we used hybrid finite element formulation that estimates stresses more accurately than displacement-based formulation. Therefore, as shown in an example where the change in resistance can be calculated analytically, the percentage error of our piezoresistive program is an order of magnitude lower relative to displacement-based finite element method.

A Novel Design Method for SOGI-PLL for Minimum Settling Time and Low Unit Vector Distortion

Phase-locked loops (PLLs) are necessary in grid connected systems to obtain information about the frequency, amplitude and phase of the grid voltage. In stationary reference frame control, the unit vectors of PLLs are used for reference generation. It is important that the PLL performance is not affected significantly when grid voltage undergoes amplitude and frequency variations. In this paper, a novel design for the popular single-phase PLL topology, namely the second-order generalized integrator (SOGI) based PLL is proposed which achieves minimum settling time during grid voltage amplitude and frequency variations. The proposed design achieves a settling time of less than 27.7 ms. This design also ensures that the unit vectors generated by this PLL have a steady state THD of less than 1% during frequency variations of the grid voltage. The design of the SOGI-PLL based on the theoretical analysis is validated by experimental results.

Analytical Closed-Form Expressions for Harmonic Distortion Corresponding to Novel Switching Sequences for Neutral-Point-Clamped Inverters

Analytical closed-form expressions for harmonic distortion factors corresponding to various pulsewidth modulation (PWM) techniques for a two-level inverter have been reported in the literature. This paper derives such analytical closed-form expressions, pertaining to centered space-vector PWM (CSVPWM) and eight different advanced bus-clamping PWM (ABCPWM) schemes, for a three-level neutral-point-clamped (NPC) inverter. These ABCPWM schemes switch each phase at twice the nominal switching frequency in certain intervals of the line cycle while clamping each phase to one of the dc terminals over certain other intervals. The harmonic spectra of the output voltages, corresponding to the eight ABCPWM schemes, are studied and compared experimentally with that of CSVPWM over the entire modulation range. The measured values of weighted total harmonic distortion (WTHD) of the line voltage V-WTHD are used to validate the analytical closed-form expressions derived. The analytical expressions, pertaining to two of the ABCPWM methods, are also validated by measuring the total harmonic distortion (THD) in the line current I-THD on a 2.2-kW constant volts-per-hertz induction motor drive.

Modelling the effects of chorus species composition and caller density on acoustic masking interference in multispecies choruses of crickets and katydids

Natural multispecies acoustic choruses such as the dusk chorus of a tropical rain forest consist of simultaneously signalling individuals of different species whose calls travel through a common shared medium before reaching their `intended' receivers. This causes masking interference between signals and impedes signal detection, recognition and localization. The levels of acoustic overlap depend on a number of factors, including call structure, intensity, habitat-dependent signal attenuation and receiver tuning. In addition, acoustic overlaps should also depend on caller density and the species composition of choruses, including relative and absolute abundance of the different calling species. In this study, we used simulations to examine the effects of chorus species relative abundance and caller density on the levels of effective heterospecific acoustic overlap in multispecies choruses composed of the calls of five species of crickets and katydids that share the understorey of a rain forest in southern India. We found that on average species-even choruses resulted in higher levels of effective heterospecific acoustic overlap than choruses with strong dominance structures. This effect was found consistently across dominance levels ranging from 0.4 to 0.8 for larger choruses of forty individuals. For smaller choruses of twenty individuals, the effect was seen consistently for dominance levels of 0.6 and 0.8 but not 0.4. Effective acoustic overlap (EAO) increased with caller density but the manner and extent of increase depended both on the species' call structure and the acoustic context provided by the composition scenario. The Phaloria sp. experienced very low levels of EAO and was highly buffered to changes in acoustic context whereas other species experienced high FAO across contexts or were poorly buffered. These differences were not simply predictable from call structures. These simulation-based findings may have important implications for acoustic biodiversity monitoring and for the study of acoustic masking interference in natural environments. (C) 2013 Elsevier B.V. All rights reserved.

Octahedral distortion induced magnetic anomalies in LaMn0.5Co0.5O3 single crystals

Single crystals of LaMn0.5Co0.5O3 belonging to the ferromagnetic-insulator and distorted perovskite class were grown using a four-mirror optical float zone furnace. The as-grown crystal crystallizes into an orthorhombic Pbnm structure. The spatially resolved 2D Raman scan reveals a strain-induced distribution of transition metal (TM)-oxygen (O) octahedral deformation in the as-grown crystal. A rigorous annealing process releases the strain, thereby generating homogeneous octahedral distortion. The octahedra tilt by reducing the bond angle TM-O-TM, resulting in a decline of the exchange energy in the annealed crystal. The critical behavior is investigated from the bulk magnetization. It is found that the ground state magnetic behavior assigned to the strain-free LaMn0.5Co0.5O3 crystal is of the 3D Heisenberg kind. Strain induces mean field-like interaction in some sites, and consequently, the critical exponents deviate from the 3D Heisenberg class in the as-grown crystal. The temperature-dependent Raman scattering study reveals strong spin-phonon coupling and the existence of two magnetic ground states in the same crystal. (C) 2014 AIP Publishing LLC.

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.