Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers. (C) 2013 Elsevier B.V. All rights reserved.

Average Modelling of a Voltage Source Inverter with Dead-Time in a Synchronous Reference Frame

Voltage source inverters are an integral part of renewable power sources and smart grid systems. Computationally efficient and fairly accurate models for the voltage source inverter are required to carry out extensive simulation studies on complex power networks. Accuracy requires that the effect of dead-time be incorporated in the inverter model. The dead-time is essentially a short delay introduced between the gating pulses to the complementary switches in an inverter leg for the safety of power devices. As the modern voltage source inverters switch at fairly high frequencies, the dead-time significantly influences the output fundamental voltage. Dead-time also causes low-frequency harmonic distortion and is hence important from a power quality perspective. This paper studies the dead-time effect in a synchronous dq reference frame, since dynamic studies and controller design are typically carried out in this frame of reference. For the sake of computational efficiency, average models are derived, incorporating the dead-time effect, in both RYB and dq reference frames. The average models are shown to consume less computation time than their corresponding switching models, the accuracies of the models being comparable. The proposed average synchronous reference frame model, including effect of dead-time, is validated through experimental results.

Influence of soil grading on the characteristics of cement stabilised soil compacts

The paper deals with experimental investigations aiming at specifying optimum soil grading limits for the production of cement stabilised soil bricks (CSSB). Wide range of soil grading curves encompassing both fine and coarse grained soils were considered. Strength, durability and absorption characteristics of CSSB were examined considering 14 different types of soil grading curves and three cement contents. The investigations show that there is optimum clay content for the soil mix which yields maximum compressive strength for CSSB and the optimum clay content is about 10 and 14 % for fine grained and coarse grained soils respectively. Void ratio of the compacted specimens is the lowest at the optimum clay content and therefore possesses maximum strength at that point. CSSB using fine grained soils shows higher strength and better durability characteristics when compared to the bricks using coarse grained soils.

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.

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.

Why is the influence of soil macrofauna on soil structure only considered by soil ecologists?

These last twenty years have seen the development of an abundant literature on the influence of soil macrofauna on soil structure. Amongst these organisms, earthworms, termites and ants are considered to play a key role in regulating the physical, chemical and microbiological properties of soils. Due to these influential impacts, soil ecologists consider these soil macro-invertebrates as `soil engineers' and their diversity and abundance are nowadays considered as relevant bioindicators of soil quality by many scientists and policy makers. Despite this abundant literature, the soil engineering concept remains a `preach to the choir' and bioturbation only perceived as important for soil ecologists. We discussed in this article the main mechanisms by which soil engineers impact soil structure and proposed to classify soil engineers with respect to their capacity to produce biostructures and modify them. We underlined the lack of studies considering biostructure dynamics and presented recent techniques in this purpose. We discussed why soil engineering concept is mainly considered by soil ecologists and call for a better collaboration between soil ecologists and soil physicists. Finally, we summarized main challenges and questions that need to be answered to integrate soil engineers activities in soil structure studies. (C) 2014 Elsevier B.V. All rights reserved.

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.

Multi-scale modelling of Suzuki segregation in gamma' precipitates in Ni and Co-base superalloys

The high temperature strength of alloys with (gamma +gamma') microstructure is primarily due to the resistance of the ordered precipitate to cutting by matrix dislocations. Such shearing requires higher stresses since it involves the creation of a planar fault. Planar fault energy is known to be dependent on composition. This implies that the composition on the fault may be different from that in the bulk for energetic reasons. Such segregation (or desegregation) of specific alloying elements to the fault may result in Suzuki strengthening which has not been explored extensively in these systems. In this work, segregation (or desegregation) of alloying elements to planar faults was studied computationally in Ni-3(Al, Ti) and Co-3(W, Al) type gamma' precipitates. The composition dependence of APB energy and heat of mixing were evaluated from first principle electronic structure calculations. A phase field model incorporating the first principles results, was used to simulate the motion of an extended superdislocation under stress concurrently with composition evolution. Results reveal that in both systems, significant (de) segregation occurs on equilibration. On application of stress, solutes were dragged along with the APB in some cases. Additionally, it was also noted the velocity of the superdislocation under an applied stress is strongly dependent on atomic mobility (i. e. diffusivity).

A Generalist Predator Regulating Spread of a Wildlife Disease: Exploring Two Infection Transmission Scenarios

Ecoepidemiology is a well-developed branch of theoretical ecology, which explores interplay between the trophic interactions and the disease spread. In most ecoepidemiological models, however, the authors assume the predator to be a specialist, which consumes only a single prey species. In few existing papers, in which the predator was suggested to be a generalist, the alternative food supply was always considered to be constant. This is obviously a simplification of reality, since predators can often choose between a number of different prey. Consumption of these alternative prey can dramatically change their densities and strongly influence the model predictions. In this paper, we try to bridge the gap and explore a generic ecoepidemiological system with a generalist predator, where the densities of all prey are dynamical variables. The model consists of two prey species, one of which is subject to an infectious disease, and a predator, which consumes both prey species. We investigate two main scenarios of infection transmission mode: (i) the disease transmission rate is predator independent and (ii) the transmission rate is a function of predator density. For both scenarios we fulfil an extensive bifurcation analysis. We show that including a second dynamical prey in the system can drastically change the dynamics of the single prey case. In particular, the presence of a second prey impedes disease spread by decreasing the basic reproduction number and can result in a substantial drop of the disease prevalence. We demonstrate that with efficient consumption of the second prey species by the predator, the predator-dependent disease transmission can not destabilize interactions, as in the case with a specialist predator. Interestingly, even if the population of the second prey eventually vanishes and only one prey species finally remains, the system with two prey species may exhibit different properties to those of the single prey system.

Retrieval and Multi-scale Validation of Soil Moisture from Multi-temporal SAR Data in a Semi-Arid Tropical Region

The current study presents an algorithm to retrieve surface Soil Moisture (SM) from multi-temporal Synthetic Aperture Radar (SAR) data. The developed algorithm is based on the Cumulative Density Function (CDF) transformation of multi-temporal RADARSAT-2 backscatter coefficient (BC) to obtain relative SM values, and then converts relative SM values into absolute SM values using soil information. The algorithm is tested in a semi-arid tropical region in South India using 30 satellite images of RADARSAT-2, SMOS L2 SM products, and 1262 SM field measurements in 50 plots spanning over 4 years. The validation with the field data showed the ability of the developed algorithm to retrieve SM with RMSE ranging from 0.02 to 0.06 m(3)/m(3) for the majority of plots. Comparison with the SMOS SM showed a good temporal behaviour with RMSE of approximately 0.05 m(3)/m(3) and a correlation coefficient of approximately 0.9. The developed model is compared and found to be better than the change detection and delta index model. The approach does not require calibration of any parameter to obtain relative SM and hence can easily be extended to any region having time series of SAR data available.